Electromechanical wave imaging for arrhythmias

An Ultrasound Matrix Transducer for High-Frame-Rate 3-D Intra-cardiac Echocardiography

OBJECTIVE

Described here is the development of an ultrasound matrix transducer prototype for high-frame-rate 3-D intra-cardiac echocardiography.

METHODS

The matrix array consists of 16 × 18 lead zirconate titanate elements with a pitch of 160 µm × 160 µm built on top of an application-specific integrated circuit that generates transmission signals and digitizes the received signals. To reduce the number of cables in the catheter to a feasible number, we implement subarray beamforming and digitization in receive and use a combination of time-division multiplexing and pulse amplitude modulation data transmission, achieving an 18-fold reduction. The proposed imaging scheme employs seven fan-shaped diverging transmit beams operating at a pulse repetition frequency of 7.7 kHz to obtain a high frame rate. The performance of the prototype is characterized, and its functionality is fully verified.

RESULTS

The transducer exhibits a transmit efficiency of 28 Pa/V at 5 cm per element and a bandwidth of 60% in transmission. In receive, a dynamic range of 80 dB is measured with a minimum detectable pressure of 10 Pa per element. The element yield of the prototype is 98%, indicating the efficacy of the manufacturing process. The transducer is capable of imaging at a frame rate of up to 1000 volumes/s and is intended to cover a volume of 70° × 70° × 10 cm.

CONCLUSION

These advanced imaging capabilities have the potential to support complex interventional procedures and enable full-volumetric flow, tissue, and electromechanical wave tracking in the heart.

Electromechanical Cycle Length Mapping for atrial arrhythmia detection and cardioversion success assessment

BACKGROUND

Direct current cardioversion (DCCV) is an established treatment to acutely convert atrial fibrillation (AF) to normal sinus rhythm. Yet, more than 70% of patients revert to AF shortly thereafter. Electromechanical Cycle Length Mapping (ECLM) is a high framerate, spectral analysis technique shown to non-invasively characterize electromechanical activation in paced canines and re-entrant flutter patients. This study assesses ECLM feasibility to map and quantify atrial arrhythmic electromechanical activation rates and inform on 1-day and 1-month DCCV response.

METHODS

Forty-five subjects (30 AF; 15 healthy sinus rhythm (SR) controls) underwent transthoracic ECLM in four standard apical 2D echocardiographic views. AF patients were imaged within 1 h pre- and post-DCCV. 3D-rendered atrial ECLM cycle length (CL) maps and spatial CL histograms were generated. CL dispersion and percentage of arrhythmic CLs≤333ms across the entire atrial myocardium were computed transmurally. ECLM results were subsequently used as indicators of DCCV success.

RESULTS

ECLM successfully confirmed the electrical atrial activation rates in 100% of healthy subjects (R2=0.96). In AF, ECLM maps localized the irregular activation rates pre-DCCV and confirmed successful post-DCCV with immediate reduction or elimination. ECLM metrics successfully distinguished DCCV 1-day and 1-month responders from non-responders, while pre-DCCV ECLM values independently predicted AF recurrence within 1-month post-DCCV.

CONCLUSIONS

ECLM can characterize electromechanical activation rates in AF, quantify their extent, and identify and predict short- and long-term AF recurrence. ELCM constitutes thus a noninvasive arrhythmia imaging modality that can aid clinicians in simultaneous AF severity quantification, prediction of AF DCCV response, and personalized treatment planning.

Personalized biomechanical insights in atrial fibrillation: opportunities & challenges

INTRODUCTION

Atrial fibrillation (AF) is an increasingly prevalent and significant worldwide health problem. Manifested as an irregular atrial electrophysiological activation, it is associated with many serious health complications. AF affects the biomechanical function of the heart as contraction follows the electrical activation, subsequently leading to reduced blood flow. The underlying mechanisms behind AF are not fully understood, but it is known that AF is highly correlated with the presence of atrial fibrosis, and with a manifold increase in risk of stroke.

AREAS COVERED

In this review, we focus on biomechanical aspects in atrial fibrillation, current and emerging use of clinical images, and personalized computational models. We also discuss how these can be used to provide patient-specific care.

EXPERT OPINION

Understanding the connection betweenatrial fibrillation and atrial remodeling might lead to valuable understanding of stroke and heart failure pathophysiology. Established and emerging imaging modalities can bring us closer to this understanding, especially with continued advancements in processing accuracy, reproducibility, and clinical relevance of the associated technologies. Computational models of cardiac electromechanics can be used to glean additional insights on the roles of AF and remodeling in heart function.

Imaging Scheme for 3-D High-Frame-Rate Intracardiac Echography: A Simulation Study

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is normally treated by RF ablation. Intracardiac echography (ICE) is widely employed during RF ablation procedures to guide the electrophysiologist in navigating the ablation catheter, although only 2-D probes are currently clinically used. A 3-D ICE catheter would not only improve visualization of the atrium and ablation catheter, but it might also provide the 3-D mapping of the electromechanical wave (EW) propagation pattern, which represents the mechanical response of cardiac tissue to electrical activity. The detection of this EW needs 3-D high-frame-rate imaging, which is generally only realizable in tradeoff with channel count and image quality. In this simulation-based study, we propose a high volume rate imaging scheme for a 3-D ICE probe design that employs 1-D micro-beamforming in the elevation direction. Such a probe can achieve a high frame rate while reducing the channel count sufficiently for realization in a 10-Fr catheter. To suppress the grating-lobe (GL) artifacts associated with micro-beamforming in the elevation direction, a limited number of fan-shaped beams with a wide azimuthal and narrow elevational opening angle are sequentially steered to insonify slices of the region of interest. An angular weighted averaging of reconstructed subvolumes further reduces the GL artifacts. We optimize the transmit beam divergence and central frequency based on the required image quality for EW imaging (EWI). Numerical simulation results show that a set of seven fan-shaped transmission beams can provide a frame rate of 1000 Hz and a sufficient spatial resolution to visualize the EW propagation on a large 3-D surface.

A new Electromechanical Wave Imaging dispersion metric for the characterization of ventricular activation in different Cardiac Resynchronization Therapy pacing schemes

Conventional biventricular (BiV) pacing cardiac resynchronization therapy (CRT) is an established treatment for heart failure patients. Recently, multiple novel CRT delivering technologies such as His-Bundle pacing have been investigated as alternative pacing strategies for optimal treatment benefit. Electromechanical Wave Imaging (EWI), a high frame-rate echocardiography-based modality, is capable of visualizing the change from dyssynchronous activation to resynchronized BiV-paced ventricles in 3D. This proof-of-concept study introduces a new EWI-based dispersion metric to further characterize ventricular activation. Patients with His-Bundle device implantation (n=4), left-bundle branch block (n=10), right-ventricular (RV) pacing (n=10), or BiV pacing (n=15) were imaged, as well as four volunteers in normal sinus rhythm (NSR). EWI successfully mapped the ventricular activation resulting from His-Bundle pacing. Additionally, very similar activation patterns were obtained in the NSR subjects, confirming recovery of physiological activation with His pacing. The dispersion metric was the most sensitive EWI-based metric that identified His pacing as the most efficient treatment (lowest activation time spread), followed by BiV and RV pacing. More specifically, the dispersion metric significantly (p 0.005) distinguished His pacing from the other two pacing schemes as well as LBBB. The initial findings presented herein indicate that EWI and its new dispersion metric may provide a useful resynchronization evaluation clinical tool in CRT patients under both novel His-Bundle pacing and more conventional BiV pacing strategies.

Myocardial Strain Imaging With Electrocardiogram-Gated and Coherent Compounding for Early Diagnosis of Coronary Artery Disease

Myocardial elastography (ME) is an ultrasound-based technique that uses radiofrequency signals for 2-D cardiac motion tracking and strain imaging at a high frame rate. Early diagnosis of coronary artery disease (CAD) is critical for timely treatment and improvement of patient outcome. The objective of this study was to assess the performance of ME radial and circumferential strains in the detection and characterization of CAD in patients. In this study, 86 patients suspected of CAD were imaged with ME prior to invasive coronary angiography (ICA). End-systolic radial and circumferential left ventricular strains were estimated in all patients in each of their perfusion territories: left anterior descending (LAD), left circumflex (LCX) and right coronary artery (RCA). ME radial strains were capable of differentiating the obstructive CAD group (55.3 ± 29.8%) from the non-obstructive CAD (72.5 ± 46.8%, p < 0.05) and no CAD groups (73.4 ± 30.4%, p < 0.05) in the RCA territory. ME circumferential strains were capable of differentiating the obstructive CAD group (-3.1 ± 7.5%) from the non-obstructive CAD (-7.2 ± 6.8%, p < 0.05) and normal (-6.9 ± 8.0%, p < 0.05) groups in the LAD territory and to differentiate the normal group (-17.1 ± 8.2%) from the obstructive (-12.8 ± 7.2%, p < 0.05) and non-obstructive CAD (-13.6 ± 8.5%, p < 0.05) groups in the RCA territory. ME circumferential strain performed better than ME radial strain in differentiating normal, non-obstructive and obstructive perfusion territories. In the LCX territory, both ME radial and circumferential strains decreased when the level of stenosis was higher. However, it was not statistically significant. The findings presented herein indicate that ME radial and circumferential estimation obtained from ECG-gated and compounded acquisitions is a promising tool for early, non-invasive and radiation-free detection of CAD in patients.

Cardiac Resynchronization Therapy Response Assessment with Electromechanical Activation Mapping within 24 Hours of Device Implantation: A Pilot Study

BACKGROUND

Cardiac resynchronization therapy (CRT) response assessment relies on the QRS complex narrowing criterion. Yet one third of patients do not improve despite narrowed QRS after implantation. Electromechanical wave imaging (EWI) is a quantitative echocardiography-based technique capable of noninvasively mapping cardiac electromechanical activation in three dimensions. The aim of this exploratory study was to investigate the EWI technique, sensitive to ventricular dyssynchrony, for informing CRT response on the day of implantation.

METHODS

Forty-four patients with heart failure with left bundle branch block or right ventricular (RV) paced rhythm and decreased left ventricular ejection fraction (LVEF; mean, 25.3 ± 9.6%) underwent EWI without and with CRT within 24 hours of device implantation. Of those, 16 were also scanned while in left ventricular (LV) pacing. Improvement in LVEF at 3-, 6-, or 9-month follow-up defined (1) super-responders (ΔLVEF ≥ 20%), (2) responders (10% ≤ ΔLVEF < 20%), and (3) nonresponders (ΔLVEF ≤ 5%). Three-dimensionally rendered electromechanical maps were obtained under RV, LV, and biventricular CRT pacing conditions. Mean RV free wall and LV lateral wall activation times were computed. The percentage of resynchronized myocardium was measured by quantifying the percentage of the left ventricle activated within 120 msec of QRS onset. Correlations between percentage of resynchronized myocardium and type of CRT response were assessed.

RESULTS

LV lateral wall activation time was significantly different (P ≤ .05) among all three pacing conditions in the 16 patients: LV lateral wall activation time with CRT in biventricular pacing (73.1 ± 17.6 msec) was lower compared with LV pacing (89.5 ± 21.5 msec) and RV pacing (120.3 ± 17.8 msec). Retrospective analysis showed that the percentage of resynchronized myocardium with CRT was a reliable response predictor within 24 hours of implantation for significantly (P ≤ .05) identifying super-responders (n = 7; 97.7 ± 1.9%) from nonresponders (n = 17; 89.9 ± 9.9%).

CONCLUSION

Electromechanical activation mapping constitutes a valuable three-dimensional visualization tool within 24 hours of implantation and could potentially aid in the timely assessment of CRT response rates, including during implantation for adjustment of lead placement and pacing outcomes.

Inverse mechano-electrical reconstruction of cardiac excitation wave patterns from mechanical deformation using deep learning

The inverse mechano-electrical problem in cardiac electrophysiology is the attempt to reconstruct electrical excitation or action potential wave patterns from the heart's mechanical deformation that occurs in response to electrical excitation. Because heart muscle cells contract upon electrical excitation due to the excitation-contraction coupling mechanism, the resulting deformation of the heart should reflect macroscopic action potential wave phenomena. However, whether the relationship between macroscopic electrical and mechanical phenomena is well-defined and unique enough to be utilized for an inverse imaging technique in which mechanical activation mapping is used as a surrogate for electrical mapping has yet to be determined. Here, we provide a numerical proof-of-principle that deep learning can be used to solve the inverse mechano-electrical problem in phenomenological two- and three-dimensional computer simulations of the contracting heart wall, or in elastic excitable media, with muscle fiber anisotropy. We trained a convolutional autoencoder neural network to learn the complex relationship between electrical excitation, active stress, and tissue deformation during both focal or reentrant chaotic wave activity and, consequently, used the network to successfully estimate or reconstruct electrical excitation wave patterns from mechanical deformation in sheets and bulk-shaped tissues, even in the presence of noise and at low spatial resolutions. We demonstrate that even complicated three-dimensional electrical excitation wave phenomena, such as scroll waves and their vortex filaments, can be computed with very high reconstruction accuracies of about 95% from mechanical deformation using autoencoder neural networks, and we provide a comparison with results that were obtained previously with a physics- or knowledge-based approach.

Noninvasive localization of cardiac arrhythmias using electromechanical wave imaging

Cardiac arrhythmias are a major cause of morbidity and mortality worldwide. The 12-lead electrocardiogram (ECG) is the current noninvasive clinical tool used to diagnose and localize cardiac arrhythmias. However, it has limited accuracy and is subject to operator bias. Here, we present electromechanical wave imaging (EWI), a high-frame rate ultrasound technique that can noninvasively map with high accuracy the electromechanical activation of atrial and ventricular arrhythmias in adult patients. This study evaluates the accuracy of EWI for localization of various arrhythmias in all four chambers of the heart before catheter ablation. Fifty-five patients with an accessory pathway (AP) with Wolff-Parkinson-White (WPW) syndrome, premature ventricular complexes (PVCs), atrial tachycardia (AT), or atrial flutter (AFL) underwent transthoracic EWI and 12-lead ECG. Three-dimensional (3D) rendered EWI isochrones and 12-lead ECG predictions by six electrophysiologists were applied to a standardized segmented cardiac model and subsequently compared to the region of successful ablation on 3D electroanatomical maps generated by invasive catheter mapping. There was significant interobserver variability among 12-lead ECG reads by expert electrophysiologists. EWI correctly predicted 96% of arrhythmia locations as compared with 71% for 12-lead ECG analyses [unadjusted for arrhythmia type: odds ratio (OR), 11.8; 95% confidence interval (CI), 2.2 to 63.2; P = 0.004; adjusted for arrhythmia type: OR, 12.1; 95% CI, 2.3 to 63.2; P = 0.003]. This double-blinded clinical study demonstrates that EWI can localize atrial and ventricular arrhythmias including WPW, PVC, AT, and AFL. EWI when used with ECG may allow for improved treatment for patients with arrhythmias.

A comparison between unfocused and focused transmit strategies in cardiac strain imaging

Unfocused ultrasound imaging, particularly coherent compounding with diverging waves, is a commonly employed high-frame rate transmit strategy in cardiac strain imaging. However, the accuracy and precision of diverging wave imaging compared to focused-beam transmit approaches in human subjects is unknown. Three transmit strategies-coherent compounding imaging, composite focused imaging with ECG gating and narrow-beams, and focused imaging with wide-beams-were compared in simulation and in transthoracic imaging of healthy human subjects (n  =  7). The focused narrow-beam sequence estimated radial end-systolic cumulative strains of a simulated left ventricular deformation with 26%  ±  1.5% and 34%  ±  1.5% greater accuracy compared with compounding and wide-beam imaging, respectively. Strain estimation precision in transthoracic imaging was then assessed with the Strain Filter on cumulative end-systolic radial strains. Within the strain values where statistically significant differences in precision (E(SNR_e|ε)) were found between transmit strategies, the narrow-beam sequence estimated radial strain 13%  ±  0.71% and 34%  ±  8.9% more precisely on average compared to compounding or wide-beam imaging, respectively.

Detection of Regional Mechanical Activation of the Left Ventricular Myocardium Using High Frame Rate Ultrasound Imaging

We have investigated the feasibility of noninvasive mapping of mechanical activation patterns in the left ventricular (LV) myocardium using high frame rate ultrasound imaging for the purpose of detecting conduction abnormalities. Five anesthetized, open-chest dogs with implanted combined sonomicrometry and electromyography (EMG) crystals were studied. The animals were paced from the specified locations of the heart, while crystal and ultrasound data were acquired. Isochrone maps of the mechanical activation patterns were generated from the ultrasound data using a novel signal processing method called clutter filter wave imaging (CFWI). The isochrone maps showed the same mechanical activation pattern as the sonomicrometry crystals in 90% of the cases. For electrical activation, the activation sequences from ultrasound were the same in 92% of the cases. The coefficient of determination between the activation delay measured with EMG and ultrasound was R ² = 0.79 , indicating a strong correlation. These results indicate that high frame rate ultrasound imaging processed with CFWI has the potential to be a valuable tool for mechanical activation detection.

Synchronization-based reconstruction of electromechanical wave dynamics in elastic excitable media

The heart is an elastic excitable medium, in which mechanical contraction is triggered by nonlinear waves of electrical excitation, which diffuse rapidly through the heart tissue and subsequently activate the cardiac muscle cells to contract. These highly dynamic excitation wave phenomena have yet to be fully observed within the depths of the heart muscle, as imaging technology is unable to penetrate the tissue and provide panoramic, three-dimensional visualizations necessary for adequate study. As a result, the electrophysiological mechanisms that are associated with the onset and progression of severe heart rhythm disorders such as atrial or ventricular fibrillation remain insufficiently understood. Here, we present a novel synchronization-based data assimilation approach with which it is possible to reconstruct excitation wave dynamics within the volume of elastic excitable media by observing spatiotemporal deformation patterns, which occur in response to excitation. The mechanical data are assimilated in a numerical replication of the measured elastic excitable system, and within this replication, the data drive the intrinsic excitable dynamics, which then coevolve and correspond to a reconstruction of the original dynamics. We provide a numerical proof-of-principle and demonstrate the performance of the approach by recovering even complicated three-dimensional scroll wave patterns, including vortex filaments of electrical excitation from within a deformable bulk tissue with fiber anisotropy. In the future, the reconstruction approach could be combined with high-speed imaging of the heart's mechanical contractions to estimate its electrophysiological activity for diagnostic purposes.

Fast qualitative two-dimensional mapping of ultrasound fields with acoustic cavitation-enhanced ultrasound imaging

Characterization of ultrasound fields is a routine procedure for both diagnostic and therapeutic ultrasound. Quantitative field mapping with a calibrated hydrophone and multi-axis positioning system can be difficult and time consuming. In this study, the use of acoustic cavitation field mapping as a qualitative surrogate to acoustic pressure field mapping, albeit without acoustic pressure values is demonstrated. This technique allows for fast qualitative mapping of ultrasound fields and thereby functionality of the corresponding transducers, in a matter of seconds. In addition, this technique could be used to rapidly image in vivo acoustic cavitation fields during therapeutic ultrasound applications.

Natural Shear Wave Imaging in the Human Heart: Normal Values, Feasibility, and Reproducibility

Left ventricular myocardial stiffness could offer superior quantification of cardiac systolic and diastolic function when compared to the current diagnostic tools. Shear wave elastography in combination with acoustic radiation force has been widely proposed to noninvasively assess tissue stiffness. Interestingly, shear waves can also result from intrinsic cardiac mechanical events (e.g., closure of valves) without the need for external excitation. However, it remains unknown whether these natural shear waves always occur, how reproducible they can be detected and what the normal range of shear wave propagation speed is. The present study, therefore, aimed at establishing the feasibility of detecting shear waves created after mitral valve closure (MVC) and aortic valve closure (AVC), the variability of the measurements, and at reporting the normal values of propagation velocity. Hereto, a group of 30 healthy volunteers was scanned with high-frame rate imaging (>1000 Hz) using an experimental ultrasound system transmitting a diverging wave sequence. Tissue Doppler velocity and acceleration were used to create septal color M-modes, on which the shear waves were tracked and their velocities measured. Overall, the methodology was capable of detecting the transient vibrations that spread throughout the intraventricular septum in response to the closure of the cardiac valves in 92% of the recordings. Reference velocities of 3.2±0.6 m/s at MVC and 3.5±0.6 m/s at AVC were obtained. Moreover, in order to show the diagnostic potential of this approach, two patients (one with cardiac amyloidosis and one undergoing a dobutamine stress echocardiography) were scanned with the same protocol and showed markedly higher propagation speeds: the former presented velocities of 6.6 and 5.6 m/s; the latter revealed normal propagation velocities at baseline, and largely increased during the dobutamine infusion (>15 m/s). Both cases showed values consistent with the expected changes in stiffness and cardiac loading conditions.

In Vivo Open- and Closed-chest Measurements of Left-Ventricular Myocardial Viscoelasticity using Lamb wave Dispersion Ultrasound Vibrometry (LDUV): A Feasibility Study

Diastolic dysfunction causes close to half of congestive heart failures and is associated with increased stiffness in left-ventricular myocardium. A clinical tool capable of measuring viscoelasticity of the myocardium could be beneficial in clinical settings. We used Lamb wave Dispersion Ultrasound Vibrometry (LDUV) for assessing the feasibility of making in vivo non-invasive measurements of myocardial elasticity and viscosity in pigs. In vivo open-chest measurements of myocardial elasticity and viscosity obtained using a Fourier space based analysis of Lamb wave dispersion are reported. The approach was used to perform ECG-gated transthoracic in vivo measurements of group velocity, elasticity and viscosity throughout a single heart cycle. Group velocity, elasticity and viscosity in the frequency range 50-500 Hz increased from diastole to systole, consistent with contraction and relaxation of the myocardium. Systolic group velocity, elasticity and viscosity were 5.0 m/s, 19.1 kPa, 6.8 Pa·s, respectively. In diastole, the measured group velocity, elasticity and viscosity were 1.5 m/s, 5.1 kPa and 3.2 Pa·s, respectively.

Non-invasive Characterization of Focal Arrhythmia with Electromechanical Wave Imaging in Vivo

There is currently no established method for the non-invasive characterization of arrhythmia and differentiation between endocardial and epicardial triggers at the point of care. Electromechanical wave imaging (EWI) is a novel ultrasound-based imaging technique based on time-domain transient strain estimation that can map and characterize electromechanical activation in the heart in vivo. The objectives of this initial feasibility study were to determine that EWI is capable of differentiating between endocardial and epicardial sources of focal rhythm and, as a proof-of-concept, that EWI could characterize focal arrhythmia in one patient with premature ventricular contractions (PVCs) before radiofrequency (RF) ablation treatment. First, validation of EWI for differentiation of surface of origin was performed in seven (n = 7) adult dogs using four epicardial and four endocardial pacing protocols. Second, one (n = 1) adult patient diagnosed with PVC was imaged with EWI before the scheduled RF ablation procedure, and EWI results were compared with mapping procedure results. In dogs, EWI was capable of detecting whether pacing was of endocardial or epicardial origin in six of seven cases (86% success rate). In the PVC patient, EWI correctly identified both regions and surface of origin, as confirmed by results from the electrical mapping obtained from the RF ablation procedure. These results reveal that EWI can map the electromechanical activation across the myocardium and indicate that EWI could serve as a valuable pre-treatment planning tool in the clinic.

Diverging wave compounding with spatio-temporal encoding using orthogonal Golay pairs for high frame rate imaging

Golay coded excitation for diverging wave compounding (DWC) has been demonstrated to increase the signal-to-noise ratio (SNR) and contrast for high frame rate cardiac imaging. However, the complementary codes need to be transmitted in two consecutive firings for decoding, which reduces the frame rate by 2 folds. This paper proposes an orthogonal Golay pairs coded (OGPs-coded) DWC sequence to overcome this problem, which implements spatio-temporal encoding for DWC. Two diverging waves (DWs) at different steering angles coded by an orthogonal Golay pair are transmitted simultaneously, thus compensating the frame rate reduction caused by transmissions of complementary codes. The two DWs can be separated based on the orthogonality of Golay pairs. To test the feasibility of the proposed sequence, we performed simulations of point targets and tissue phantoms in both static and moving states. Compared with non-coded DWC at the same frame rate, OGPs-coded DWC obtains comparable resolution, SNR gains of 7.5-10 dB and contrast gains of 3-5 dB. The OGPs-coded DWC sequence was also tested experimentally on a tissue-mimicking phantom. Compared with non-coded DWC, OGPs-coded DWC achieves improvements in the SNR (3-6 dB) and contrast (1-2 dB). Preliminary in vivo results show brighter myocardium and larger penetration depth with the proposed method. The proposed OGPs-coded DWC sequence has potential for high frame rate and high quality cardiac imaging.

Spatiotemporal matrix image formation for programmable ultrasound scanners

As programmable ultrasound scanners become more common in research laboratories, it is increasingly important to develop robust software-based image formation algorithms that can be obtained in a straightforward fashion for different types of probes and sequences with a small risk of error during implementation. In this work, we argue that as the computational power keeps increasing, it is becoming practical to directly implement an approximation to the matrix operator linking reflector point targets to the corresponding radiofrequency signals via thoroughly validated and widely available simulations software. Once such a spatiotemporal forward-problem matrix is constructed, standard and thus highly optimized inversion procedures can be leveraged to achieve very high quality images in real time. Specifically, we show that spatiotemporal matrix image formation produces images of similar or enhanced quality when compared against standard delay-and-sum approaches in phantoms and in vivo, and show that this approach can be used to form images even when using non-conventional probe designs for which adapted image formation algorithms are not readily available.

Electromechanical wave imaging and electromechanical wave velocity estimation in a large animal model of myocardial infarction

Echocardiography is often used in the clinic for detection and characterization of myocardial infarction. Electromechanical wave imaging (EWI) is a non-invasive ultrasound-based imaging technique based on time-domain incremental motion and strain estimation that can evaluate changes in contractility in the heart. In this study, electromechanical activation is assessed in infarcted heart to determine whether EWI is capable of detecting and monitoring infarct formation. Additionally, methods for estimating electromechanical wave (EW) velocity are presented, and changes in the EW propagation velocity after infarct formation are studied. Five (n  =  5) adult mongrels were used in this study. Successful infarct formation was achieved in three animals by ligation of the left anterior descending (LAD) coronary artery. Dogs were survived for a few days after LAD ligation and monitored daily with EWI. At the end of the survival period, dogs were sacrificed and TTC (tetrazolium chloride) staining confirmed the formation and location of the infarct. In all three dogs, as soon as day 1 EWI was capable of detecting late-activated and non-activated regions, which grew over the next few days. On final day images, the extent of these regions corresponded to the location of infarct as confirmed by staining. EW velocities in border zones of infarct were significantly lower post-infarct formation when compared to baseline, whereas velocities in healthy tissues were not. These results indicate that EWI and EW velocity might help with the detection of infarcts and their border zones, which may be useful for characterizing arrhythmogenic substrate.

Reproducibility and Angle Independence of Electromechanical Wave Imaging for the Measurement of Electromechanical Activation during Sinus Rhythm in Healthy Humans

Electromechanical wave imaging (EWI) is an ultrasound-based technique that can non-invasively map the transmural electromechanical activation in all four cardiac chambers in vivo. The objective of this study was to determine the reproducibility and angle independence of EWI for the assessment of electromechanical activation during normal sinus rhythm (NSR) in healthy humans. Acquisitions were performed transthoracically at 2000 frames/s on seven healthy human hearts in parasternal long-axis, apical four- and two-chamber views. EWI data was collected twice successively in each view in all subjects, while four successive acquisitions were obtained in one case. Activation maps were generated and compared (i) within the same acquisition across consecutive cardiac cycles; (ii) within same view across successive acquisitions; and (iii) within equivalent left-ventricular regions across different views. EWI was capable of characterizing electromechanical activation during NSR and of reliably obtaining similar patterns of activation. For consecutive heart cycles, the average 2-D correlation coefficient between the two isochrones across the seven subjects was 0.9893, with a mean average activation time fluctuation in LV wall segments across acquisitions of 6.19%. A mean activation time variability of 12% was obtained across different views with a measurement bias of only 3.2 ms. These findings indicate that EWI can map the electromechanical activation during NSR in human hearts in transthoracic echocardiography in vivo and results in reproducible and angle-independent activation maps.

Cardiac Strain Imaging With Coherent Compounding of Diverging Waves

Current methods of cardiac strain imaging at high frame rate suffer from motion matching artifacts or poor lateral resolution. Coherent compounding has been shown to improve echocardiographic image quality while maintaining a high frame rate, but has never been used to image cardiac strain. However, myocardial velocity can have an impact on coherent compounding due to displacements between frames. The objective of this paper was to investigate the feasibility and performance of coherent compounding for cardiac strain imaging at a low and a high myocardial velocity. Left-ventricular contraction in short-axis view was modeled as an annulus with radial thickening and circumferential rotation. Simulated radio-frequency channel data with a cardiac phased array were obtained using three different beamforming methods: single diverging wave, coherent compounding of diverging waves, and conventional focusing. Axial and lateral displacements and strains as well as radial strains were estimated and compared to their true value. In vivo feasibility of cardiac strain imaging with coherent compounding was performed and compared to single diverging wave imaging. At low myocardial velocities, the axial, lateral, and radial strain relative error for nine compounded waves (16.3%, 40.4%, and 18.9%) were significantly lower than those obtained with single diverging wave imaging (19.9%, 80.3%, and 30.6%) and closer to that obtained with conventional focusing (16.7%, 43.7%, and 16%). In vivo left-ventricular radial strains exhibited higher quality with nine compounded waves than with single diverging wave imaging. These results indicate that cardiac strain can be imaged using coherent compounding of diverging waves with a better performance than with single diverging wave imaging while maintaining a high frame rate, and therefore, has the potential to improve diagnosis of myocardial strain-based cardiac diseases.

Technical Note: A 3-D rendering algorithm for electromechanical wave imaging of a beating heart

PURPOSE

Arrhythmias can be treated by ablating the heart tissue in the regions of abnormal contraction. The current clinical standard provides electroanatomic 3-D maps to visualize the electrical activation and locate the arrhythmogenic sources. However, the procedure is time-consuming and invasive. Electromechanical wave imaging is an ultrasound-based noninvasive technique that can provide 2-D maps of the electromechanical activation of the heart. In order to fully visualize the complex 3-D pattern of activation, several 2-D views are acquired and processed separately. They are then manually registered with a 3-D rendering software to generate a pseudo-3-D map. However, this last step is operator-dependent and time-consuming.

METHODS

This paper presents a method to generate a full 3-D map of the electromechanical activation using multiple 2-D images. Two canine models were considered to illustrate the method: one in normal sinus rhythm and one paced from the lateral region of the heart. Four standard echographic views of each canine heart were acquired. Electromechanical wave imaging was applied to generate four 2-D activation maps of the left ventricle. The radial positions and activation timings of the walls were automatically extracted from those maps. In each slice, from apex to base, these values were interpolated around the circumference to generate a full 3-D map.

RESULTS

In both cases, a 3-D activation map and a cine-loop of the propagation of the electromechanical wave were automatically generated. The 3-D map showing the electromechanical activation timings overlaid on realistic anatomy assists with the visualization of the sources of earlier activation (which are potential arrhythmogenic sources). The earliest sources of activation corresponded to the expected ones: septum for the normal rhythm and lateral for the pacing case.

CONCLUSIONS

The proposed technique provides, automatically, a 3-D electromechanical activation map with a realistic anatomy. This represents a step towards a noninvasive tool to efficiently localize arrhythmias in 3-D.

Evaluation of Coronary Artery Disease Using Myocardial Elastography with Diverging Wave Imaging: Validation against Myocardial Perfusion Imaging and Coronary Angiography

Myocardial elastography (ME) is an ultrasound-based technique that can image 2-D myocardial strains. The objectives of this study were to illustrate that 2-D myocardial strains can be imaged with diverging wave imaging and differ, on average, between normal and coronary artery disease (CAD) patients. In this study, 66 patients with symptoms of CAD were imaged with myocardial elastography before a nuclear stress test or an invasive coronary angiography. Radial cumulative strains were estimated in all patients. The end-systolic radial strain in the total cross section of the myocardium was significantly higher in normal patients (17.9 ± 8.7%) than in patients with reversible perfusion defect (6.2 ± 9.3%, p < 0.001) and patients with significant (-0.9 ± 7.4%, p < 0.001) and non-significant (3.7 ± 5.7%, p < 0.01) lesions. End-systolic radial strain in the left anterior descending, left circumflex and right coronary artery territory was found to be significantly higher in normal patients than in CAD patients. These preliminary findings indicate that end-systolic radial strain measured with ME is higher on average in healthy persons than in CAD patients and that ME has the potential to be used for non-invasive, radiation-free early detection of CAD.

Coded excitation for diverging wave cardiac imaging: a feasibility study

Diverging wave (DW) based cardiac imaging has gained increasing interest in recent years given its capacity to achieve ultrahigh frame rate. However, the signal-to-noise ratio (SNR), contrast, and penetration depth of the resulting B-mode images are typically low as DWs spread energy over a large region. Coded excitation is known to be capable of increasing the SNR and penetration for ultrasound imaging. The aim of this study was therefore to test the feasibility of applying coded excitation in DW imaging to improve the corresponding SNR, contrast and penetration depth. To this end, two types of codes, i.e. a linear frequency modulated chirp code and a set of complementary Golay codes were tested in three different DW imaging schemes, i.e. 1 angle DW transmit without compounding, 3 and 5 angles DW transmits with coherent compounding. The performances (SNR, contrast ratio (CR), contrast-to-noise ratio (CNR), and penetration) of different imaging schemes were investigated by means of simulations and in vitro experiments. As for benchmark, corresponding DW imaging schemes with regular pulsed excitation as well as the conventional focused imaging scheme were also included. The results showed that the SNR was improved by about 10 dB using coded excitation while the penetration depth was increased by 2.5 cm and 1.8 cm using chirp code and Golay codes, respectively. The CNR and CR gains varied with the depth for different DW schemes using coded excitations. Specifically, for non-compounded DW imaging schemes, the gain in the CR was about 5 dB and 3 dB while the gain in the CNR was about 4.5 dB and 3.5 dB at larger depths using chirp code and Golay codes, respectively. For compounded imaging schemes, using coded excitation, the gain in the penetration and contrast were relatively smaller compared to non-compounded ones. Overall, these findings indicated the feasibility of coded excitation in improving the image quality of DW imaging. Preliminary in vivo cardiac images of a healthy volunteer were presented finally, and higher SNR and deeper penetration depth can be achieved by coded schemes.

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.

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.

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

BACKGROUND

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.

METHODS

EWI acquisitions were obtained in five healthy controls and 16 HF patients with and without CRT pacing. Responders (n = 8) and nonresponders (n = 8) to CRT were identified retrospectively on the basis of left ventricular (LV) reverse remodeling. Electromechanical activation maps were obtained in all patients and used to compute a quantitative parameter describing the mean LV lateral wall activation time (LWAT).

RESULTS

Mean LWAT was increased by 52.1 ms in HF patients in native rhythm compared to controls (P < 0.01). For all HF patients, CRT pacing initiated a different electromechanical activation sequence. Responders exhibited a 56.4-ms ± 28.9-ms reduction in LWAT with CRT pacing (P < 0.01), while nonresponders showed no significant change.

CONCLUSION

In this initial feasibility study, EWI was capable of characterizing local cardiac electromechanical behavior as it pertains to HF and CRT response. Activation sequences obtained with EWI allow for quantification of LV lateral wall electromechanical activation, thus providing a novel method for CRT assessment.

Extension of Fourier-Based Techniques for Ultrafast Imaging in Ultrasound With Diverging Waves

Ultrafast ultrasound imaging has become an intensive area of research thanks to its capability in reaching high frame rates. In this paper, we propose a scheme that allows the extension of the current Fourier-based techniques derived for planar acquisition to the reconstruction of sectorial scan with wide angle using diverging waves. The flexibility of the proposed formulation was assessed through two different Fourier-based techniques. The performance of the derived approaches was evaluated in terms of resolution and contrast from both simulations and in vitro experiments. The comparisons of the current state-of-the-art method with the conventional delay-and-sum technique illustrated the potential of the derived methods for producing competitive results with lower computational complexity.

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 Jude's 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 ² 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.

High-Frame-Rate Echocardiography Using Coherent Compounding With Doppler-Based Motion-Compensation

High-frame-rate ultrasonography based on coherent compounding of unfocused beams can potentially transform the assessment of cardiac function. As it requires successive waves to be combined coherently, this approach is sensitive to high-velocity tissue motion. We investigated coherent compounding of tilted diverging waves, emitted from a 2.5 MHz clinical phased array transducer. To cope with high myocardial velocities, a triangle transmit sequence of diverging waves is proposed, combined with tissue Doppler imaging to perform motion compensation (MoCo). The compound sequence with integrated MoCo was adjusted from simulations and was tested in vitro and in vivo. Realistic myocardial velocities were analyzed in an in vitro spinning disk with anechoic cysts. While a 8 dB decrease (no motion versus high motion) was observed without MoCo, the contrast-to-noise ratio of the cysts was preserved with the MoCo approach. With this method, we could provide high-quality in vivo B-mode cardiac images with tissue Doppler at 250 frames per second. Although the septum and the anterior mitral leaflet were poorly apparent without MoCo, they became well perceptible and well contrasted with MoCo. The septal and lateral mitral annulus velocities determined by tissue Doppler were concordant with those measured by pulsed-wave Doppler with a clinical scanner (r(2)=0.7,y=0.9 x+0.5,N=60) . To conclude, high-contrast echo cardiographic B-mode and tissue Doppler images can be obtained with diverging beams when motion compensation is integrated in the coherent compounding process.

Assessing the Stability of Aortic Aneurysms with Pulse Wave Imaging

Purpose To assess whether the stability of murine aortic aneurysms is associated with the homogeneity of pulse wave propagation within the saccular wall. Materials and Methods All animal procedures were approved by the institutional Animal Care and Use Committee. Apolipoprotein E and tissue inhibitor of metalloproteinases-1 knockout mice (n = 26) were infused with angiotensin II by using subcutaneously implanted osmotic pumps, with an additional control mouse used for histologic examination (n = 1). Pulse wave imaging (PWI) was performed just before infusion and 15 days after infusion by using 40-MHz ultrasonography at 8000 frames per second (with electrocardiographic gating). Aneurysm appearance on B-mode images was monitored every 2-3 days for 30 days. On the basis of B-mode images obtained after 30 days, aneurysms were deemed to have been unstable if they had ruptured; otherwise, they were deemed stable. Statistical significance was assessed by using two-tailed t tests. Results In normal aortas, the pulse waves propagated at relatively constant velocities (mean ± standard deviation, 2.8 m/sec ± 0.9). Fifteen days after infusion, all mice had developed aneurysms, with significant (P < .001/12) changes in maximum anterior-posterior diameter (increase of 54.9% ± 2.5) and pulse wave velocity (PWV) (decrease of 1.3 m/sec ± 0.8). While there was no significant difference in these parameters (P = .45 for diameter and P = .55 for PWV) between stable aneurysms (n = 12) and unstable aneurysms (n = 14), the standard deviation of the high-resolution PWV was significantly higher (P < .001/12) in unstable aneurysms (5.7 m/sec ± 1.6) than in stable ones (3.2 m/sec ± 0.9). Conclusion High-resolution PWI was used to measure the local homogeneity of pulse wave propagation within the saccular wall, which is lower in unstable aneurysms than in stable ones. Hence, if proven to add additional information beyond size and appearance in human studies, PWI could potentially be used to assess the stability of aneurysms by providing information that is complementary to the anatomic data obtained with conventional B-mode imaging. ^© RSNA, 2016 Online supplemental material is available for this article.

Assessing the atrial electromechanical coupling during atrial focal tachycardia, flutter, and fibrillation using electromechanical wave imaging in humans

Minimally-invasive treatments of cardiac arrhythmias such as radio-frequency ablation are gradually gaining importance in clinical practice but still lack a noninvasive imaging modality which provides insight into the source or focus of an arrhythmia. Cardiac deformations imaged at high temporal and spatial resolution can be used to elucidate the electrical activation sequence in normal and paced human subjects non-invasively and could potentially aid to better plan and monitor ablation-based arrhythmia treatments. In this study, a novel ultrasound-based method is presented that can be used to quantitatively characterize focal and reentrant arrhythmias. Spatio-temporal maps of the full-view of the atrial and ventricular mechanics were obtained in a single heartbeat, revealing with otherwise unobtainable detail the electromechanical patterns of atrial flutter, fibrillation, and tachycardia in humans. During focal arrhythmias such as premature ventricular complex and focal atrial tachycardia, the previously developed electromechanical wave imaging methodology is hereby shown capable of identifying the location of the focal zone and the subsequent propagation of cardiac activation. During reentrant arrhythmias such as atrial flutter and fibrillation, Fourier analysis of the strains revealed highly correlated mechanical and electrical cycle lengths and propagation patterns. High frame rate ultrasound imaging of the heart can be used non-invasively and in real time, to characterize the lesser-known mechanical aspects of atrial and ventricular arrhythmias, also potentially assisting treatment planning for intraoperative and longitudinal monitoring of arrhythmias.

3-D ultrafast Doppler imaging applied to the noninvasive mapping of blood vessels in vivo

Ultrafast Doppler imaging was introduced as a technique to quantify blood flow in an entire 2-D field of view, expanding the field of application of ultrasound imaging to the highly sensitive anatomical and functional mapping of blood vessels. We have recently developed 3-D ultrafast ultrasound imaging, a technique that can produce thousands of ultrasound volumes per second, based on a 3-D plane and diverging wave emissions, and demonstrated its clinical feasibility in human subjects in vivo. In this study, we show that noninvasive 3-D ultrafast power Doppler, pulsed Doppler, and color Doppler imaging can be used to perform imaging of blood vessels in humans when using coherent compounding of 3-D tilted plane waves. A customized, programmable, 1024-channel ultrasound system was designed to perform 3-D ultrafast imaging. Using a 32 × 32, 3-MHz matrix phased array (Vermon, Tours, France), volumes were beamformed by coherently compounding successive tilted plane wave emissions. Doppler processing was then applied in a voxel-wise fashion. The proof of principle of 3-D ultrafast power Doppler imaging was first performed by imaging Tygon tubes of various diameters, and in vivo feasibility was demonstrated by imaging small vessels in the human thyroid. Simultaneous 3-D color and pulsed Doppler imaging using compounded emissions were also applied in the carotid artery and the jugular vein in one healthy volunteer.

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, r(2) = 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.

High-frame rate four dimensional optoacoustic tomography enables visualization of cardiovascular dynamics and mouse heart perfusion

Functional imaging of mouse models of cardiac health and disease provides a major contribution to our fundamental understanding of the mammalian heart. However, imaging murine hearts presents significant challenges due to their small size and rapid heart rate. Here we demonstrate the feasibility of high-frame-rate, noninvasive optoacoustic imaging of the murine heart. The temporal resolution of 50 three-dimensional frames per second provides functional information at important phases of the cardiac cycle without the use of gating or other motion-reduction methods. Differentiation of the blood oxygenation state in the heart chambers was enabled by exploiting the wavelength dependence of optoacoustic signals. Real-time volumetric tracking of blood perfusion in the cardiac chambers was also evaluated using indocyanine green. Taken together, the newly-discovered capacities offer a unique tool set for in-vivo structural and functional imaging of the whole heart with high spatio-temporal resolution in all three dimensions.

Intracardiac myocardial elastography in canines and humans in vivo

Intracardiac echocardiography (ICE) is a useful imaging modality which is used during RF ablation procedures to identify anatomical structures. Utilizing ICE in conjunction with myocardial elastography (ME) can provide additional information on the mechanical properties of cardiac tissue and provide information on mechanical changes caused by ablation. The objective of this study was to demonstrate that ICE can be used at high frame rate using a diverging beam transmit sequence to image myocardial strain and differentiate myocardial tissue properties before, during, and after ablation for a clinical ablation procedure. In this feasibility study, three normal canines and eight patients with atrial fibrillation (AF) were studied in vivo. A 5.8-MHz ICE transducer was used to image the heart with a diverging beam transmit method achieving 1200 frames per second (fps). Cumulative axial displacement estimation was performed using 1-D cross-correlation with a window size of 2.7 mm and 95% overlap. Axial cumulative strains were estimated in the left atrium (LA) and right atrium (RA) using a least-squares estimator with a kernel of 2 mm on the axial displacements. In the canine case, radial thickening was detected in the lateral wall and in the interatrial septum during LA emptying. For AF patients, the mean absolute strain in the ablated region was lower (6.7 ± 3.1%) than before the ablation (17.4 ± 9.3%) in LA at the end of the LA emptying phase. In the cavotricuspid isthmus (CTI) region, mean absolute strain magnitude at the end of the RA emptying phase was found to be higher during ablation (43.0 ± 18.1%) compared with after ablation (33.7 ± 15.8%). Myocardial strains in the LA of an AF patient were approximately 2.6 times lower in the ablated region than before ablation. This initial feasibility indicates that ME can be used as a new imaging modality in conjunction with ICE in RF ablation guidance and lesion monitoring.

Effect of subaperture beamforming on phase coherence imaging

High-frame-rate echocardiography using unfocused transmit beams and parallel receive beamforming is a promising method for evaluation of cardiac function, such as imaging of rapid propagation of vibration of the heart wall resulting from electrical stimulation of the myocardium. In this technique, high temporal resolution is realized at the expense of spatial resolution and contrast. The phase coherence factor has been developed to improve spatial resolution and contrast in ultrasonography. It evaluates the variance in phases of echo signals received by individual transducer elements after delay compensation, as in the conventional delay-andsum beamforming process. However, the phase coherence factor suppresses speckle echoes because phases of speckle echoes fluctuate as a result of interference of echoes. In the present study, the receiving aperture was divided into several subapertures, and conventional delay-and-sum beamforming was performed with respect to each subaperture to suppress echoes from scatterers except for that at a focal point. After subaperture beamforming, the phase coherence factor was obtained from beamformed RF signals from respective subapertures. By means of this procedure, undesirable echoes, which can interfere with the echo from a focal point, can be suppressed by subaperture beamforming, and the suppression of the phase coherence factor resulting from phase fluctuation caused by such interference can be avoided. In the present study, the effect of subaperture beamforming in high-frame-rate echocardiography with the phase coherence factor was evaluated using a phantom. By applying subaperture beamforming, the average intensity of speckle echoes from a diffuse scattering medium was significantly higher (-39.9 dB) than that obtained without subaperture beamforming (-48.7 dB). As for spatial resolution, the width at half-maximum of the lateral echo amplitude profile obtained without the phase coherence factor was 1.06 mm. By using the phase coherence factor, spatial resolution was improved significantly, and subaperture beamforming achieved a better spatial resolution of 0.75 mm than that of 0.78 mm obtained without subaperture beamforming.

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 ultrafast ultrasound imaging in vivo

Very high frame rate ultrasound imaging has recently allowed for the extension of the applications of echography to new fields of study such as the functional imaging of the brain, cardiac electrophysiology, and the quantitative imaging of the intrinsic mechanical properties of tumors, to name a few, non-invasively and in real time. In this study, we present the first implementation of Ultrafast Ultrasound Imaging in 3D based on the use of either diverging or plane waves emanating from a sparse virtual array located behind the probe. It achieves high contrast and resolution while maintaining imaging rates of thousands of volumes per second. A customized portable ultrasound system was developed to sample 1024 independent channels and to drive a 32  ×  32 matrix-array probe. Its ability to track in 3D transient phenomena occurring in the millisecond range within a single ultrafast acquisition was demonstrated for 3D Shear-Wave Imaging, 3D Ultrafast Doppler Imaging, and, finally, 3D Ultrafast combined Tissue and Flow Doppler Imaging. The propagation of shear waves was tracked in a phantom and used to characterize its stiffness. 3D Ultrafast Doppler was used to obtain 3D maps of Pulsed Doppler, Color Doppler, and Power Doppler quantities in a single acquisition and revealed, at thousands of volumes per second, the complex 3D flow patterns occurring in the ventricles of the human heart during an entire cardiac cycle, as well as the 3D in vivo interaction of blood flow and wall motion during the pulse wave in the carotid at the bifurcation. This study demonstrates the potential of 3D Ultrafast Ultrasound Imaging for the 3D mapping of stiffness, tissue motion, and flow in humans in vivo and promises new clinical applications of ultrasound with reduced intra--and inter-observer variability.

Transthoracic ultrafast Doppler imaging of human left ventricular hemodynamic function

Heart diseases can affect intraventricular blood flow patterns. Real-time imaging of blood flow patterns is challenging because it requires both a high frame rate and a large field of view. To date, standard Doppler techniques can only perform blood flow estimation with high temporal resolution within small regions of interest. In this work, we used ultrafast imaging to map in 2-D human left ventricular blood flow patterns during the whole cardiac cycle. Cylindrical waves were transmitted at 4800 Hz with a transthoracic phased-array probe to achieve ultrafast Doppler imaging of the left ventricle. The high spatio-temporal sampling of ultrafast imaging permits reliance on a much more effective wall filtering and increased sensitivity when mapping blood flow patterns during the pre-ejection, ejection, early diastole, diastasis, and late diastole phases of the heart cycle. The superior sensitivity and temporal resolution of ultrafast Doppler imaging makes it a promising tool for the noninvasive study of intraventricular hemodynamic function.

New developments in paediatric cardiac functional ultrasound imaging

Ultrasound imaging can be used to estimate the morphology as well as the motion and deformation of tissues. If the interrogated tissue is actively deforming, this deformation is directly related to its function and quantification of this deformation is normally referred as 'strain imaging'. Tissue can also be deformed by applying an internal or external force and the resulting, induced deformation is a function of the mechanical tissue characteristics. In combination with the load applied, these strain maps can be used to estimate or reconstruct the mechanical properties of tissue. This technique was named 'elastography' by Ophir et al. in 1991. Elastography can be used for atherosclerotic plaque characterisation, while the contractility of the heart or skeletal muscles can be assessed with strain imaging. Rather than using the conventional video format (DICOM) image information, radio frequency (RF)-based ultrasound methods enable estimation of the deformation at higher resolution and with higher precision than commercial methods using Doppler (tissue Doppler imaging) or video image data (2D speckle tracking methods). However, the improvement in accuracy is mainly achieved when measuring strain along the ultrasound beam direction, so it has to be considered a 1D technique. Recently, this method has been extended to multiple directions and precision further improved by using spatial compounding of data acquired at multiple beam steered angles. Using similar techniques, the blood velocity and flow can be determined. RF-based techniques are also beneficial for automated segmentation of the ventricular cavities. In this paper, new developments in different techniques of quantifying cardiac function by strain imaging, automated segmentation, and methods of performing blood flow imaging are reviewed and their application in paediatric cardiology is discussed.

High-contrast ultrafast imaging of the heart

Noninvasive ultrafast imaging of intrinsic waves such as electromechanical waves or remotely induced shear waves in elastography imaging techniques for human cardiac applications remains challenging. In this paper, we propose ultrafast imaging of the heart with adapted sector size by coherently compounding diverging waves emitted from a standard transthoracic cardiac phased-array probe. As in ultrafast imaging with plane wave coherent compounding, diverging waves can be summed coherently to obtain high-quality images of the entire heart at high frame rate in a full field of view. To image the propagation of shear waves with a large SNR, the field of view can be adapted by changing the angular aperture of the transmitted wave. Backscattered echoes from successive circular wave acquisitions are coherently summed at every location in the image to improve the image quality while maintaining very high frame rates. The transmitted diverging waves, angular apertures, and subaperture sizes were tested in simulation, and ultrafast coherent compounding was implemented in a commercial scanner. The improvement of the imaging quality was quantified in phantoms and in one human heart, in vivo. Imaging shear wave propagation at 2500 frames/s using 5 diverging waves provided a large increase of the SNR of the tissue velocity estimates while maintaining a high frame rate. Finally, ultrafast imaging with 1 to 5 diverging waves was used to image the human heart at a frame rate of 4500 to 900 frames/s over an entire cardiac cycle. Spatial coherent compounding provided a strong improvement of the imaging quality, even with a small number of transmitted diverging waves and a high frame rate, which allows imaging of the propagation of electromechanical and shear waves with good image quality.

Electromechanical wave imaging of biologically and electrically paced canine hearts in vivo

Electromechanical wave imaging (EWI) has been show capable of directly and entirely non-invasively mapping the trans mural electromechanical activation in all four cardiac chambers in vivo. In this study, we assessed EWI repeatability and reproducibility, as well as its capability of localizing electronic and, for the first time, biological pacing locations in closed-chest, conscious canines. Electromechanical activation was obtained in six conscious animals during normal sinus rhythm (NSR) and idioventricular rhythms occurring in dogs with complete heart block instrumented with electronic and biologic pacemakers (EPM and BPM respectively). After atrioventricular node ablation, dogs were implanted with an EPM in the right ventricular (RV) endocardial apex (n = 4) and two additionally received a BPM at the left ventricular (LV) epicardial base (n = 2). EWI was performed trans thoracically during NSR, BPM and EPM pacing, in conscious dogs, using an unfocused transmit sequence at 2000 frames/s. During NSR, the EW originated at the right atrium (RA), propagated to the left atrium (LA) and emerged from multiple sources in both ventricles. During EPM, the EW originated at the RV apex and propagated throughout both ventricles. During BPM, the EW originated from the LV basal lateral wall and subsequently propagated throughout the ventricles. EWI differentiated BPM from EPM and NSR and identified the distinct pacing origins. Isochrone comparison indicated that EWI was repeatable and reliable. These findings thus indicate the potential for EWI to serve as a simple, non-invasive and direct imaging technology for mapping and characterizing arrhythmias as well as the treatments thereof.

Ultra-high frame rate tissue Doppler imaging

We describe a new tissue Doppler imaging (TDI) method, ultra-high frame rate tissue Doppler imaging (UFR-TDI). With two broad transmit beams covering only the ventricular walls, we achieve 1200 frames/s in a four-chamber apical view. We examined 10 healthy volunteers to study the feasibility of this method. Ultra-high-frame-rate TDI provided peak annular velocities and time to peak S' intervals in good agreement with those measured with conventional TDI. Moreover, UFR-TDI provided additional information in early and late systole: In all subjects, the method was able to separate the timing of electrical activation, start of mechanical contraction, mitral valve closure and start of ejection. The earliest mechanical activation was seen before mitral valve closure. The method was also able to measure the propagation speed of the mechanical wave created by aortic valve closure.

Stolt's f-k migration for plane wave ultrasound imaging

Ultrafast ultrasound is an emerging modality that offers new perspectives and opportunities in medical imaging. Plane wave imaging (PWI) allows one to attain very high frame rates by transmission of planar ultrasound wave-fronts. As a plane wave reaches a given scatterer, the latter becomes a secondary source emitting upward spherical waves and creating a diffraction hyperbola in the received RF signals. To produce an image of the scatterers, all the hyperbolas must be migrated back to their apexes. To perform beamforming of plane wave echo RFs and return high-quality images at high frame rates, we propose a new migration method carried out in the frequency-wavenumber (f-k) domain. The f-k migration for PWI has been adapted from the Stolt migration for seismic imaging. This migration technique is based on the exploding reflector model (ERM), which consists in assuming that all the scatterers explode in concert and become acoustic sources. The classical ERM model, however, is not appropriate for PWI. We showed that the ERM can be made suitable for PWI by a spatial transformation of the hyperbolic traces present in the RF data. In vitro experiments were performed to outline the advantages of PWI with Stolt's f-k migration over the conventional delay-and-sum (DAS) approach. The Stolt's f-k migration was also compared with the Fourier-based method developed by J.-Y. Lu. Our findings show that multi-angle compounded f-k migrated images are of quality similar to those obtained with a stateof- the-art dynamic focusing mode. This remained true even with a very small number of steering angles, thus ensuring a highly competitive frame rate. In addition, the new FFT-based f-k migration provides comparable or better contrast-to-noise ratio and lateral resolution than the Lu's and DAS migration schemes. Matlab codes for the Stolt's f-k migration for PWI are provided.