Impedance Imaging of the Thorax: Why it s difficult, and what we are doing about it?

Similar documents
Accounting for erroneous electrode data in Electrical Impedance Tomography. or, Salvaging EIT Data

Monitoring of the lung fluid movement and estimation of lung area using Electrical Impedance Tomography: A Simulation Study

Statistical and Computational Inverse Problems with Applications Part 2: Introduction to inverse problems and example applications

Inverse parameter identification problems

MRI of the airways and lungs Including hyperpolarized techniques

Velocity Images. Phase Contrast Technique. G. Reiter 1,2, U. Reiter 1, R. Rienmüller 1

Soft Nanocomposite Based Multi-point, Multidirectional Strain Mapping Sensor Using Anisotropic Electrical Impedance Tomography

CT-PET calibration : physical principles and operating procedures F.Bonutti. Faustino Bonutti Ph.D. Medical Physics, Udine University Hospital.

Two-dimensional Tissue Image Reconstruction Based on Magnetic Field Data

Numerical Reconstruction of Perfectly Conducting Inclusions from One Electrostatic Boundary Measurement

What are inverse problems?

arxiv: v1 [math.ap] 6 Mar 2017

Analysis of Sigmoidal Equations To Describe the Pulmonary Pressure- Volume Curve in Acute Respiratory Distress Syndrome

Amherst. University of Massachusetts Amherst. Tamara V. Sobers University of Massachusetts Amherst

EE 4BD4 Lecture 26. Plethysmography

Propagation of Measurement Noise Through Backprojection Reconstruction in Electrical Impedance Tomography

Quantification of the magnetic dust deposition in the magnetopneumographic diagnostic of the human respiratory tract

A uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance tomography

Electrical Impedance Tomography Based on Vibration Excitation in Magnetic Resonance System

Electromagnetic Tomography: Real-Time Imaging using Linear Approaches

A GENERIC FRAMEWORK FOR THE DEVELOPMENT OF THE SIGNAL SIMULATOR

A method of the forward problem for magneto-acousto-electrical tomography

Electromagnetic inverse problems in biomedical engineering Jens Haueisen

COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING QUESTION BANK SUBJECT CODE & NAME : EE1354-MODERN CONTROL SYSTEMS

Regularized Nonrigid Registration of Lung CT Images by Preserving Tissue Volume and Vesselness Measure

NONLINEAR DYNAMICS AND CHAOS. Facilitating medical diagnosis. Medical classifications

Electrical Impedance Tomography using EIDORS in a Closed Phantom


PRELIMINARY RESULTS ON BRAIN MONITORING OF MENINGITIS USING 16 CHANNELS MAGNETIC IN- DUCTION TOMOGRAPHY MEASUREMENT SYSTEM

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS IB PHYSICS

Optimum Design of the Structure of Electrode for Medical EIT System

Excerpt from the Proceedings of the COMSOL Conference 2009 Milan

A mathematical model and inversion procedure for Magneto-Acousto-Electric Tomography (MAET)

Regularizing the optimization-based solution of inverse coefficient problems with monotonicity constraints

Reconstructing conductivities with boundary corrected D-bar method

Inverse problems and medical imaging

Lectures on Medical Biophysics Dept. of Biophysics, Medical faculty, Masaryk University in Brno. Biosignals and their processing Thermometry

Introduction to Medical Imaging. Medical Imaging

Methods for Calculating the Electrode Position Jacobian for Impedance Imaging

Novel tomography techniques and parameter identification problems

Magnetoacoustic tomography with current injection

AN IMPLEMENTATION OF THE BACK-PROJECTION ALGORITHM ACCORDING TO SANTOSA AND VOGELIUS

Trends in hybrid data tomography Workshop at DTU Compute Wednesday January 24, 2018 Room 324/050

Electrical Capacitance Tomography

Field trip: Tuesday, Feb 5th

Electrical conductivity images of biological tissue phantoms in MREIT

Localization of Multiple Deep Epileptic Sources in a Realistic Head Model via Independent Component Analysis

Midterm II Sample Problems

M. DAROWSKI 1, W. KLONOWSKI 1,2, M. KOZARSKI 1, G. FERRARI 3, K. ZIELIŃSKI 1, R. STEPIEN 1,2

Final Exam: 11 May 2001, Friday

UNIVERSITY OF DUBLIN

Inverse problems and medical imaging

EE 5345 Biomedical Instrumentation Lecture 6: slides

NUMERICAL SIMULATION OF THE HIGH FREQUENCY OSCILLATORY VENTILATION IN GENERIC MODELS OF THE HUMAN AIRWAYS

Monotonicity-based Electrical Impedance Tomography for Lung Imaging

Development of Tetrapolar Conductivity Cell for Liquid Measurement Application

An impedance method for spatial sensing of 3D cell constructs towards applications in tissue engineering

A Parameter-Choice Method That Exploits Residual Information

Malvin Carl Teich. Boston University and Columbia University Workshop on New Themes & Techniques in Complex Systems 2005

Tomography is imaging by sections. 1

Introduction and Background

Put Paper Number Here

Regional Time Constants Determined by Electrical Impedance Tomography are Affected by Ventilatory Parameters

Development of High-resolution EEG Devices

Advanced Topics and Diffusion MRI

REFAAT E. GABR, PHD Fannin Street, MSE R102D, Houston, Texas 77030

Microwave Ablation and associated Dielectric Properties- Modelling, Measurements and Sensitivity Investigation

Winters s formula revisited

The effect of different number of diffusion gradients on SNR of diffusion tensor-derived measurement maps

MR Advance Techniques. Flow Phenomena. Class I

Chapter 1. The Human Organism 1-1

MRI based quantification of outflow boundary conditions for wall shear stress calculations in stenosed human carotid bifurcations

Use of the magnetic field generated by the internal distribution of injected currents for Electrical Impedance Tomography (MR-EIT)

COPD Phenotyping: Hyperpolarized Helium-3 Magnetic Resonance Functional Imaging

HST-582J/6.555J/16.456J Biomedical Signal and Image Processing Spring Laboratory Project 3 Blind Source Separation: Fetal & Maternal ECG

Reaction Rates and Chemical Equilibrium

and Chemical Equilibrium Reaction Rates

Reaction Rates and Chemical Equilibrium. Chapter 10

A Riemannian Framework for Denoising Diffusion Tensor Images

Forward and inverse wave problems:quantum billiards and brain imaging p.1

Slope Stability Monitoring. Impedance Imaging

EL-GY 6813/BE-GY 6203 Medical Imaging, Fall 2016 Final Exam

Identifying pressure dependent elastance in lung mechanics with reduced influence of unmodelled effects

Professor Stuart Bunt 217

Wearable Sensor Scanner using Electrical Impedance Tomography

CONDITIONAL JOINT TRANSFER ENTROPY OF CARDIOVASCULAR AND CEREBROVASCULAR CONTROL SYSTEMS IN SUBJECTS PRONE TO POSTURAL SYNCOPE

THE VENTILATION CYCLE IN OCTOPUS

Automated P-Wave Quality Assessment for Wearable Sensors

Application of Maxwell Equations to Human Body Modelling

Analysing data: regression and correlation S6 and S7

A Neurosurgeon s Perspectives of Diffusion Tensor Imaging(DTI) Diffusion Tensor MRI (DTI) Background and Relevant Physics.

Overture: Why is Physiologic Variability Important?

Longitudinal growth analysis of early childhood brain using deformation based morphometry

An Introduction to Anatomy and Physiology

Navigator Echoes. BioE 594 Advanced Topics in MRI Mauli. M. Modi. BioE /18/ What are Navigator Echoes?

Sound Touch Elastography

CHAPTER 5 EEG SIGNAL CLASSIFICATION BASED ON NN WITH ICA AND STFT

B4 Organising animals and plants. Student Book answers. B4.1 The blood. Question Answer Marks Guidance

Radioactive Decedents What is the risk?

Transcription:

Impedance Imaging of the Thorax: Why it s difficult, and what we are doing about it? Biomedical Engineering Research Centre (BIRC) Western Univerity, London, ON, 6 May 2015 Andy Adler Professor & Canada Research Chair in Biomedical Engineering Systems and Computer Engineering, Carleton University, Ottawa Impedance Imaging, Andy Adler, 2015/05/06 1 / 39

Electrical Impedance Tomography 10-day old healthy baby with EIT electrodes Source: eidors3d.sf.net/data contrib/ifneonate-spontaneous Impedance Imaging, Andy Adler, 2015/05/06 2 / 39

Electronics Block Diagram Impedance Imaging, Andy Adler, 2015/05/06 3 / 39

Current Propagation Healthy Adult Male CT slide at heart Source: eidors3d.sf.net/tutorial/netgen/extrusion Impedance Imaging, Andy Adler, 2015/05/06 4 / 39

Current Propagation Healthy Adult Male CT slide at heart Source: eidors3d.sf.net/tutorial/netgen/extrusion Impedance Imaging, Andy Adler, 2015/05/06 4 / 39

Current Propagation I Healthy Adult Male CT slide at heart Source: eidors3d.sf.net/tutorial/netgen/extrusion Impedance Imaging, Andy Adler, 2015/05/06 4 / 39

Current Propagation I Healthy Adult Male CT slide at heart Source: eidors3d.sf.net/tutorial/netgen/extrusion Impedance Imaging, Andy Adler, 2015/05/06 4 / 39

Current Propagation I Healthy Adult Male CT slide at heart Source: eidors3d.sf.net/tutorial/netgen/extrusion Impedance Imaging, Andy Adler, 2015/05/06 4 / 39

Finite Element Modelling Impedance Imaging, Andy Adler, 2015/05/06 5 / 39

Finite Element Modelling Simulated Voltages Voxel Currents Impedance Imaging, Andy Adler, 2015/05/06 6 / 39

Thorax Propagation I CT Slice with simulated current streamlines and voltage equipotentials Impedance Imaging, Andy Adler, 2015/05/06 / 39

Thorax Propagation I CT Slice with simulated current streamlines and voltage equipotentials Impedance Imaging, Andy Adler, 2015/05/06 / 39

Changing Conductivity I Heart receives blood (diastole) and is more conductive Impedance Imaging, Andy Adler, 2015/05/06 8 / 39

Changing Conductivity I Heart receives blood (diastole) and is more conductive Impedance Imaging, Andy Adler, 2015/05/06 9 / 39

Changing Conductivity Impedance Imaging, Andy Adler, 2015/05/06 10 / 39

Application: Breathing Chest images of tidal breathing in healthy adult Impedance Imaging, Andy Adler, 2015/05/06 11 / 39

Application: Heart EIT Signal in ROI around heart (and ECG) Impedance Imaging, Andy Adler, 2015/05/06 12 / 39

Mechanical Ventilation Mechanical Ventilator with EIT monitor Source: Swisstom.com Impedance Imaging, Andy Adler, 2015/05/06 13 / 39

Acute Respiratory Distress Syndrome (ARDS) Impedance Imaging, Andy Adler, 2015/05/06 14 / 39

EIT + Lung State A B C D E F G Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT + Lung State A B C D E F G Overdistension A Collapse Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT + Lung State A B C D E F G Overdistension A B Collapse Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT + Lung State A B C D E F G Overdistension A B C Collapse Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT + Lung State A B C D E F G Overdistension A B C D Collapse Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT + Lung State A B C D E F G Overdistension A B C D Collapse E Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT + Lung State A B C D E F G Overdistension A B C D Collapse E F Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT + Lung State A B C D E F G Overdistension A B C D Collapse E F G Impedance Imaging, Andy Adler, 2015/05/06 15 / 39

EIT for Non-Invasive Blood Pressure Pulse transit time from heart to descending aorta using EIT Source: Sola et al, Med. Biol. Eng. Comput., 2011 Impedance Imaging, Andy Adler, 2015/05/06 16 / 39

Neonatal Breathing Preterm newborns have complex, unstable physiology Ventilatory support is often essential Currently, no adequate monitors of breathing These data are from a lamb model of neonates Source: D. Tingay et al, In Press Impedance Imaging, Andy Adler, 2015/05/06 1 / 39

EIT for Brain Imaging Applications: Epileptic foci Stroke (Ischaemic vs. Haemoragic) Fast Neural Imaging Source: Holder, www.ucl.ac.uk/medphys/ research/eit/pubs/brain EIT overview.pdf Impedance Imaging, Andy Adler, 2015/05/06 18 / 39

EIT for Cancer Imaging: Breast/Prostate Cancerous tissue has different electrical properties Image tissue Image increased vascularization Source: Khan, Mahara, Halter et al, Conf. EIT, 2014 Impedance Imaging, Andy Adler, 2015/05/06 19 / 39

Non-medical applications Flow in pipes Mixing tanks Imaging metalic ores Hydro-geology Figure shows resistivity in a cross-section of La Soufrière de Guadaloupe volcano. Source: N. Lesparre et al, Conf. EIT, 2014 Impedance Imaging, Andy Adler, 2015/05/06 20 / 39

Why is EIT hard? Impedance Imaging, Andy Adler, 2015/05/06 21 / 39

Inverse Problems... Plato s cave Impedance Imaging, Andy Adler, 2015/05/06 22 / 39

Plato s cave... Shadows on the wall Source: iamcriselleeee.files.wordpress.com/2013/11/cave-2.jpg Impedance Imaging, Andy Adler, 2015/05/06 23 / 39

Inverse Problems Forward Problem: Forms Shadows Impedance Imaging, Andy Adler, 2015/05/06 24 / 39

Inverse Problems Forward Problem: Inverse Problem: Forms Shadows Shadows Forms Impedance Imaging, Andy Adler, 2015/05/06 24 / 39

Inverse Problems Forward Problem: Inverse Problem: Forms Shadows Shadows Forms Ill-conditioned Sensitivity to some movements is low Impedance Imaging, Andy Adler, 2015/05/06 24 / 39

Inverse Problems Forward Problem: Inverse Problem: Forms Shadows Shadows Forms Ill-conditioned Sensitivity to some movements is low Ill-posed Some movements don t change shadows Impedance Imaging, Andy Adler, 2015/05/06 24 / 39

Inverse Problems Forward Problem: Inverse Problem: Forms Shadows Shadows Forms Ill-conditioned Sensitivity to some movements is low Ill-posed Some movements don t change shadows Noisy Flickering light Impedance Imaging, Andy Adler, 2015/05/06 24 / 39

Inverse Problems Techniques: to calculate stable & meaningful parameters in the presence of inversion difficulties Examples Image deblurring / restoration Medical imaging Geophyical imaging Model parameter fitting Reconstruction with incomplete/noisy data Impedance Imaging, Andy Adler, 2015/05/06 25 / 39

Forward Model (linearized) Reconstruction in Pictures Forward Problem Measurements (difference) Image (difference) = + noise Jacobian Impedance Imaging, Andy Adler, 2015/05/06 26 / 39

Forward Model (linearized) Reconstruction in Pictures Forward Problem Measurements (difference) Image (difference) = + noise Regularized linear inverse model Jacobian Linear Solution: Minimize norm 2 + Penalty Function Norm weighted by measurement accuracy Impedance Imaging, Andy Adler, 2015/05/06 26 / 39

Electrode errors: Zero bad data Idea #1: Reconstruction with Data Errors Traditional Solution Error Here Replace With zero 2-2 σ 1-2 σ 2-2 σ 3-2 σ 4 Impedance Imaging, Andy Adler, 2015/05/06 2 / 39

Electrode errors: Zero bad data Idea #1: Reconstruction with Data Errors Traditional Solution Error Here Replace With zero 2 Error Model Solution -2 σ 1-2 σ 2-2 σ 3-2 σ 4 Electrode Error: Regularized imaging Error Here 2-2 σ 1-2 σ 2 Low SNR here Replace With zero -2 σ 3-2 σ 4 Impedance Imaging, Andy Adler, 2015/05/06 2 / 39

Electrode Error compensation Offline compensation using jack-knife approach (2005) A B EIT images in anaesthetised, ventilated dog A: uncompensated, B: compensated. Left: ventilation Centre: saline (right lung) Right: ventilation and saline Automatic detection (via reciprocity comparison) (2009) New work to speed online calculation & use data quality Impedance Imaging, Andy Adler, 2015/05/06 28 / 39

Imaging Electrode Movement Idea #2: Electrode movement Sensitivity to sensor movement = + noise Jacobian now includes measurement change due to movement image now includes x,y sensor movement Impedance Imaging, Andy Adler, 2015/05/06 29 / 39

Imaging Electrode Movement Idea #2: Electrode movement Sensitivity to sensor movement = + noise Adapted penalty function Jacobian now includes measurement change due to movement Penalty:Image and movement Expected image Expected movement 2 1 -½ -½ 1 -½ -½ 1 -½ image now includes x,y sensor movement -½ 1 -½ -½ 1 -½ -½ 1 Unlikelyhood of movement and image co-variance Unlikelyhood of movement Impedance Imaging, Andy Adler, 2015/05/06 29 / 39

Electrode movement compensation Source: Gómez-Laberge et al, Phys. Meas., 2006 Impedance Imaging, Andy Adler, 2015/05/06 30 / 39

Idea #3: Data Quality Impedance Imaging, Andy Adler, 2015/05/06 31 / 39

Idea #3: Data Quality Depth Sounder with analog and digital guages Impedance Imaging, Andy Adler, 2015/05/06 31 / 39

What s the problem? With strong priors and complex algorithms, algorithms give us pretty pictures, even when they are irrelevant. Question: how can we know when to trust a pretty picture? how can we know when the data are junk? Impedance Imaging, Andy Adler, 2015/05/06 32 / 39

Data Quality Measure: Concept Concept: High Quality Data are Consistent Idea: Use IP to predict each data point from all others Original d Remove i d (i) Solve ˆm (i) Predict ˆd (i) Calculate error ɛ i = d i ˆd (i) i Impedance Imaging, Andy Adler, 2015/05/06 33 / 39

Normalized data quality values Example: Data quality measures R1 R2 R3 R4 T1 T2 T3 T4 A Patient 1 Patient 0.95 0.9 B 0.85 0.8 0.5 0. 0.65 0.6 Patient1 Patient R1 R2 R3 R4 T1 T2 T3 T4 Protocol steps Clinical data and data quality metric for each stage of the protocol (R1 R4 recruitment: PEEP, T1 T4 titration: PEEP ). A: EIT images B: Calculated data quality. Impedance Imaging, Andy Adler, 2015/05/06 34 / 39

Perspectives Data analysis is hard powerful algorithms are useful we live in a world of big data complex systems fail in complex ways users like pretty pictures So... the situation will get worse Impedance Imaging, Andy Adler, 2015/05/06 35 / 39

Solutions? Impedance Imaging, Andy Adler, 2015/05/06 36 / 39

Solutions? Impedance Imaging, Andy Adler, 2015/05/06 36 / 39

Solutions? Impedance Imaging, Andy Adler, 2015/05/06 36 / 39

Solutions? Impedance Imaging, Andy Adler, 2015/05/06 36 / 39

Solutions? Thus, we need Impedance Imaging, Andy Adler, 2015/05/06 36 / 39

Solutions? Thus, we need Open Data Impedance Imaging, Andy Adler, 2015/05/06 36 / 39

Solutions? Thus, we need Open Data Open source analysis Impedance Imaging, Andy Adler, 2015/05/06 36 / 39

For EIT... Impedance Imaging, Andy Adler, 2015/05/06 3 / 39

For EIT... Impedance Imaging, Andy Adler, 2015/05/06 3 / 39

For EIT... Impedance Imaging, Andy Adler, 2015/05/06 3 / 39

Traffic jam on the way to Carleton Impedance Imaging, Andy Adler, 2015/05/06 38 / 39

Abstract for talk: Impedance Imaging of the Thorax: Why its difficult and what we are doing about it Abstract: Electrical impedance tomography (EIT) uses body surface electrical current stimulation and measurements to generate images of the internal tissue electrical impedance. Currently, the most successful application of EIT is for imaging the thorax, where the movement on conductivity contrasting air and blood can be imaged over time. The generation of EIT images requires solving an inverse problem, which is ill-conditioned because of the diffuse nature of current propagation. The technology is thus sensitive to electrode properties, data quality, and patient movement. To address these issues, several innovative strategies to analyze and interpret these data have been developed. This talk will explain our recent progress in imaging the chest with EIT, and the image generation and interpretation strategies that are required. Impedance Imaging, Andy Adler, 2015/05/06 39 / 39