Ethan Bunting

Sparse matrix beamforming and image reconstruction for 2-D HIFU monitoring using harmonic motion imaging for focused ultrasound (HMIFU) with in vitro validation.

Harmonic motion imaging for focused ultrasound (HMIFU) utilizes an amplitude-modulated HIFU beam to induce a localized focal oscillatory motion simultaneously estimated. The objective of this study is to develop and show the feasibility of a novel fast beamforming algorithm for image reconstruction using GPU-based sparse-matrix operation with real-time feedback. In this study, the algorithm was implemented onto a fully integrated, clinically relevant HMIFU system. A single divergent transmit beam was used while fast beamforming was implemented using a GPU-based delay-and-sum method and a sparse-matrix operation. Axial HMI displacements were then estimated from the RF signals using a 1-D normalized cross-correlation method and streamed to a graphic user interface with frame rates up to 15 Hz, a 100-fold increase compared to conventional CPU-based processing. The real-time feedback rate does not require interrupting the HIFU treatment. Results in phantom experiments showed reproducible HMI images and monitoring of 22 in vitro HIFU treatments using the new 2-D system demonstrated reproducible displacement imaging, and monitoring of 22 in vitro HIFU treatments using the new 2-D system showed a consistent average focal displacement decrease of 46.7 ±14.6% during lesion formation. Complementary focal temperature monitoring also indicated an average rate of displacement increase and decrease with focal temperature at 0.84±1.15%/°C, and 2.03±0.93%/°C, respectively. These results reinforce the HMIFU capability of estimating and monitoring stiffness related changes in real time. Current ongoing studies include clinical translation of the presented system for monitoring of HIFU treatment for breast and pancreatic tumor applications.

Stochastic precision analysis of 2D cardiac strain estimation in vivo

Ultrasonic strain imaging has been applied to echocardiography and carries great potential to be used as a tool in the clinical setting. Two-dimensional (2D) strain estimation may be useful when studying the heart due to the complex, 3D deformation of the cardiac tissue. Increasing the framerate used for motion estimation, i.e. motion estimation rate (MER), has been shown to improve the precision of the strain estimation, although maintaining the spatial resolution necessary to view the entire heart structure in a single heartbeat remains challenging at high MERs. Two previously developed methods, the temporally unequispaced acquisition sequence (TUAS) and the diverging beam sequence (DBS), have been used in the past to successfully estimate in vivo axial strain at high MERs without compromising spatial resolution. In this study, a stochastic assessment of 2D strain estimation precision is performed in vivo for both sequences at varying MERs (65, 272, 544, 815 Hz for TUAS; 250, 500, 1000, 2000 Hz for DBS). 2D incremental strains were estimated during left ventricular contraction in five healthy volunteers using a normalized cross-correlation function and a least-squares strain estimator. Both sequences were shown capable of estimating 2D incremental strains in vivo. The conditional expected value of the elastographic signal-to-noise ratio (E(SNRe|ε)) was used to compare strain estimation precision of both sequences at multiple MERs over a wide range of clinical strain values. The results here indicate that axial strain estimation precision is much more dependent on MER than lateral strain estimation, while lateral estimation is more affected by strain magnitude. MER should be increased at least above 544 Hz to avoid suboptimal axial strain estimation. Radial and circumferential strain estimations were influenced by the axial and lateral strain in different ways. Furthermore, the TUAS and DBS were found to be of comparable precision at similar MERs.

3D quasi-static ultrasound elastography with plane wave in vivo.

In biological tissue, an increase in elasticity is often a marker of abnormalities. Techniques such as quasi-static ultrasound elastography have been developed to assess the strain distribution in soft tissues in two dimensions using a quasi-static compression. However, as abnormalities can exhibit very heterogeneous shapes, a three dimensional approach would be necessary to accurately measure their volume and remove operator dependency. Acquisition of volumes at high rates is also critical to performing real-time imaging with a simple freehand compression. In this study, we developed for the first time a 3D quasi-static ultrasound elastography method with plane waves that estimates axial strain distribution in vivo in entire volumes at high volume rate. Acquisitions were performed with a 2D matrix array probe of 2.5 MHz frequency and 256 elements. Plane waves were emitted at a volume rate of 100 volumes/s during a continuous motorized and freehand compression. 3D B-mode volumes and 3D cumulative axial strain volumes were successfully estimated in inclusion phantoms and in ex vivo canine liver before and after a high intensity focused ultrasound ablation. We also demonstrated the in vivo feasibility of the method using freehand compression on the calf muscle of a human volunteer and were able to retrieve 3D axial strain volume at a high volume rate depicting the differences in stiffness of the two muscles which compose the calf muscle. 3D ultrasound quasi-static elastography with plane waves could become an important technique for the imaging of the elasticity in human bodies in three dimensions using simple freehand scanning.

Electromechanical wave imaging (EWI) validation in all four cardiac chambers with 3D electroanatomic mapping in canines in vivo

Characterization and mapping of arrhythmias is currently performed through invasive insertion and manipulation of cardiac catheters. Electromechanical wave imaging (EWI) is a non-invasive ultrasound-based imaging technique, which tracks the electromechanical activation that immediately follows electrical activation. Electrical and electromechanical activations were previously found to be linearly correlated in the left ventricle, but the relationship has not yet been investigated in the three other chambers of the heart. The objective of this study was to investigate the relationship between electrical and electromechanical activations and validate EWI in all four chambers of the heart with conventional 3D electroanatomical mapping. Six (n  =  6) normal adult canines were used in this study. The electrical activation sequence was mapped in all four chambers of the heart, both endocardially and epicardially using the St Judes EnSite 3D mapping system (St. Jude Medical, Secaucus, NJ). EWI acquisitions were performed in all four chambers during normal sinus rhythm, and during pacing in the left ventricle. Isochrones of the electromechanical activation were generated from standard echocardiographic imaging views. Electrical and electromechanical activation maps were co-registered and compared, and electrical and electromechanical activation times were plotted against each other and linear regression was performed for each pair of activation maps. Electromechanical and electrical activations were found to be directly correlated with slopes of the correlation ranging from 0.77 to 1.83, electromechanical delays between 9 and 58 ms and R 2 values from 0.71 to 0.92. The linear correlation between electrical and electromechanical activations and the agreement between the activation maps indicate that the electromechanical activation follows the pattern of propagation of the electrical activation. This suggests that EWI may be used as a novel non-invasive method to accurately characterize and localize sources of arrhythmias.

Validation of electromechanical wave imaging in a canine model during pacing and sinus rhythm

BACKGROUND Accurate determination of regional areas of arrhythmic triggers is of key interest to diagnose arrhythmias and optimize their treatment. Electromechanical wave imaging (EWI) is an ultrasound technique that can image the transient deformation in the myocardium after electrical activation and therefore has the potential to detect and characterize location of triggers of arrhythmias. OBJECTIVES The objectives of this study were to investigate the relationship between the electromechanical and the electrical activation of the left ventricular (LV) endocardial surface during epicardial and endocardial pacing and during sinus rhythm as well as to map the distribution of electromechanical delays. METHODS In this study, 6 canines were investigated. Two external electrodes were sutured onto the epicardial surface of the LV. A 64-electrode basket catheter was inserted through the apex of the LV. Ultrasound channel data were acquired at 2000 frames/s during epicardial and endocardial pacing and during sinus rhythm. Electromechanical and electrical activation maps were synchronously obtained from the ultrasound data and the basket catheter, respectively. RESULTS The mean correlation coefficient between electromechanical and electrical activation was 0.81 for epicardial anterior pacing, 0.79 for epicardial lateral pacing, 0.69 for endocardial pacing, and 0.56 for sinus rhythm. CONCLUSION The electromechanical activation sequence determined by EWI follows the electrical activation sequence and more specifically in the case of pacing. This finding is of key interest in the role that EWI can play in the detection of the anatomical source of arrhythmias and the planning of pacing therapies such as cardiovascular resynchronization therapy.

Atrial electromechanical cycle length mapping in paced canine hearts in vivo

Atrial arrhythmias affect millions of people worldwide. Characterization and study of arrhythmias in the atria in the clinic is currently performed point by point using mapping catheters capable of generating maps of the electrical activation rate or cycle length. In this paper, we describe a new ultrasound-based mapping technique called electromechanical cycle length mapping (ECLM) capable of estimating the electromechanical activation rate, or cycle length, i.e., the rate of the mechanical activation of the myocardium which follows the electrical activation. ECLM relies on frequency analysis of the incremental strain within the atria and can be performed in a single acquisition. ECLM was validated in a canine model paced from the left atrial appendage, against pacing rates within the reported range of cycle lengths previously measured during atrial arrhythmias such as atrial fibrillation. Correlation between the global estimated electromechanical cycle lengths and pacing rates was shown to be excellent (slope = 0.983, intercept = 3.91, r2 = 0.9999). The effect of the number of cardiac cycles on the performance of ECLM was also investigated and the reproducibility of ECLM was demonstrated (error between consecutive acquisitions for all pacing rates: 6.3 ± 4.3%). These findings indicate the potential of ECLM for noninvasively characterizing atrial arrhythmias and provide feedback on the treatment planning of catheter ablation procedures in the clinic.

Pulse wave imaging using coherent compounding in a phantom and in vivo.

Pulse wave velocity (PWV) is a surrogate marker of arterial stiffness linked to cardiovascular morbidity. Pulse wave imaging (PWI) is a technique developed by our group for imaging the pulse wave propagation in vivo. PWI requires high temporal and spatial resolution, which conventional ultrasonic imaging is unable to simultaneously provide. Coherent compounding is known to address this tradeoff and provides full aperture images at high frame rates. This study aims to implement PWI using coherent compounding within a GPU-accelerated framework. The results of the implemented method were validated using a silicone phantom against static mechanical testing. Reproducibility of the measured PWVs was assessed in the right common carotid of six healthy subjects (n  =  6) approximately 10-15 mm before the bifurcation during two cardiac cycles over the course of 1-3 d. Good agreement of the measured PWVs (3.97  ±  1.21 m s-1, 4.08  ±  1.15 m s-1, p  =  0.74) was obtained. The effects of frame rate, transmission angle and number of compounded plane waves on PWI performance were investigated in the six healthy volunteers. Performance metrics such as the reproducibility of the PWVs, the coefficient of determination (r 2), the SNR of the PWI axial wall velocities ([Formula: see text]) and the percentage of lateral positions where the pulse wave appears to arrive at the same time-point, indicating inadequacy of the temporal resolution (i.e. temporal resolution misses) were used to evaluate the effect of each parameter. Compounding plane waves transmitted at 1° increments with a linear array yielded optimal performance, generating significantly higher r 2 and [Formula: see text] values (p  ⩽  0.05). Higher frame rates (⩾1667 Hz) produced improvements with significant gains in the r 2 coefficient (p  ⩽  0.05) and significant increase in both r 2 and [Formula: see text] from single plane wave imaging to 3-plane wave compounding (p  ⩽  0.05). Optimal performance was established at 2778 Hz with 3 plane waves and at 1667 Hz with 5 plane waves.

3D myocardial elastography in vivo

Strain evaluation is of major interest in clinical cardiology as it can quantify the cardiac function. Myocardial elastography, a radio-frequency (RF)-based cross-correlation method, has been developed to evaluate the local strain distribution in the heart in vivo. However, inhomogeneities such as RF ablation lesions or infarction require a three-dimensional approach to be measured accurately. In addition, acquisitions at high volume rate are essential to evaluate the cardiac strain in three dimensions. Conventional focused transmit schemes using 2D matrix arrays, trade off sufficient volume rate for beam density or sector size to image rapid moving structure such as the heart, which lowers accuracy and precision in the strain estimation. In this study, we developed 3D myocardial elastography at high volume rates using diverging wave transmits to evaluate the local axial strain distribution in three dimensions in three open-chest canines before and after radio-frequency ablation. Acquisitions were performed with a 2.5 MHz 2D matrix array fully programmable used to emit 2000 diverging waves at 2000 volumes/s. Incremental displacements and strains enabled the visualization of rapid events during the QRS complex along with the different phases of the cardiac cycle in entire volumes. Cumulative displacement and strain volumes depict high contrast between non-ablated and ablated myocardium at the lesion location, mapping the tissue coagulation. 3D myocardial strain elastography could thus become an important technique to measure the regional strain distribution in three dimensions in humans.

Imaging the Propagation of the Electromechanical Wave in Heart Failure Patients with Cardiac Resynchronization Therapy

Current electrocardiographic and echocardiographic measurements in heart failure (HF) do not take into account the complex interplay between electrical activation and local wall motion. The utilization of novel technologies to better characterize cardiac electromechanical behavior may lead to improved response rates with cardiac resynchronization therapy (CRT). Electromechanical wave imaging (EWI) is a noninvasive ultrasound‐based technique that uses the transient deformations of the myocardium to track the intrinsic EW that precedes myocardial contraction. In this paper, we investigate the performance and reproducibility of EWI in the assessment of HF patients and CRT.

Intracardiac myocardial elastography for lesion quantification in cardiac radiofrequency ablation

Radiofrequency ablation of the myocardium is used to treat various cardiac arrhythmias. The size, spacing, and transmurality of lesions have been shown to affect the success of the ablation procedure; however, there is currently no method to directly image the size and formation of ablation lesions in real time. Intracardiac myocardial elastrography has been used previously to image the reduction in end-systolic strain in the ablated region as a result of the lesion formation. However, the relationship between end-systolic strain change and lesion size has not been investigated. In this study, a large animal model is used to establish a relationship between the area affected by the strain reduction and lesion volume. Ablation lesions (n=10) were created in the left ventricular epicardium in five anesthetized canines. A clinical intracardiac echocardiography machine was programmed to emit a custom diverging beam sequence at 600 Hz and used to image the ablation site before and after the induction of a lesion. Cumulative strains were estimated over systole using a normalized cross-correlation displacement algorithm and a Savitzky-Golay strain kernel. The reduction in strain as a result of the ablation was computed by comparing cumulative end-systolic strains before and after ablation. Lesion volume was also measured ex vivo and compared to the area of significant strain change (>8% reduction) for each lesion. A good correlation was found between the area of significant strain change and lesion volume (r2 = 0.86). These results indicate that end-systolic strain measured using ME can be used to estimate the size of lesions induced during an RF ablation procedure, potentially assisting clinicians in lesion formation assessment during the procedure.