Single pulse frequency compounding protocol for superharmonic imaging

Second harmonic imaging is currently accepted as the standard in commercial echographic systems. A new imaging technique, coined as superharmonic imaging (SHI), combines the third till the fifth harmonics, arising during nonlinear sound propagation. It could further enhance the resolution and quality of echographic images. To meet the bandwidth requirement for SHI a dedicated phased array has been developed: a low frequency subarray, intended for transmission, interleaved with a high frequency subarray, used in reception. As the bandwidth of the elements is limited, the spectral gaps in between the harmonics cause multiple (ghost) reflection artifacts. A dual-pulse frequency compounding method aims at suppressing those artifacts at a price of a reduced frame rate. In this study we explore a possibility of performing frequency compounding within a single transmission. The traditional frequency compounding method suppresses the ripples by consecutively emitting two short Gaussian bursts with a slightly different center frequency. In the newly proposed method, the transmit aperture is divided into two parts: the first half is used to send a pulse at the lower center frequency, while the other half simultaneously transmits at a slightly higher center frequency. The suitability of the protocol for medical imaging applications in terms of the steering capabilities was performed in a simulation study with INCS and the hydrophone measurements. Moreover, an experimental study was carried out to find the optimal parameters for the clinical imaging protocol. The latter was subsequently used to obtain the images of a tissue mimicking phantom containing strongly reflecting wires. Additionally, the images of a human heart in the parasternal projection were acquired. The scanning aperture with the developed protocol amounts to approximately 90°, which is sufficient to capture the cardiac structures in the standard anatomical projections. The theoretically estimated and experimentally measured grating lobe levels are equal to -28.3 dB and -35.9 dB, respectively. A considerable improvement in the axial resolution of the SHI component (0.73 mm) at -6 dB in comparison with the third harmonic (2.23 mm) was observed. A similar comparison in terms of the lateral resolution slightly favored the superharmonic component by 0.2 mm. Additionally, the images of the tissue mimicking phantom exhibited the absence of the multiple reflection artifacts. The in-vivo acquisition allows one to clearly observe the dynamic of the mitral valve leaflets. The new method is equally effective in eliminating the ripple artifacts associated with SHI as the dual-pulse technique, while the full frame rate is maintained.

Dual-pulse frequency compounded superharmonic imaging

Tissue second-harmonic imaging is currently the default mode in commercial diagnostic ultrasound systems. A new modality, superharmonic imaging (SHI), combines the third through fifth harmonics originating from nonlinear wave propagation through tissue. SHI could further improve the resolution and quality of echographic images. The superharmonics have gaps between the harmonics because the transducer has a limited bandwidth of about 70% to 80%. This causes ghost reflection artifacts in the superharmonic echo image. In this work, a new dual-pulse frequency compounding (DPFC) method to eliminate these artifacts is introduced. In the DPFC SHI method, each trace is constructed by summing two firings with slightly different center frequencies. The feasibility of the method was established using a single-element transducer. Its acoustic field was modeled in KZK simulations and compared with the corresponding measurements obtained with a hydrophone apparatus. Subsequently, the method was implemented on and optimized for a setup consisting of an interleaved phased-array transducer (44 elements at 1 MHz and 44 elements at 3.7 MHz, optimized for echocardiography) and a programmable ultrasound system. DPFC SHI effectively suppresses the ghost reflection artifacts associated with imaging using multiple harmonics. Moreover, compared with the single-pulse third harmonic, DPFC SHI improved the axial resolution by 3.1 and 1.6 times at the -6-dB and -20-dB levels, respectively. Hence, DPFC offers the possibility of generating harmonic images of a higher quality at a cost of a moderate frame rate reduction.

Comparison of fundamental, second harmonic, and superharmonic imaging: a simulation study

In medical ultrasound, fundamental imaging (FI) uses the reflected echoes from the same spectral band as that of the emitted pulse. The transmission frequency determines the trade-off between penetration depth and spatial resolution. Tissue harmonic imaging (THI) employs the second harmonic of the emitted frequency band to construct images. Recently, superharmonic imaging (SHI) has been introduced, which uses the third to the fifth (super) harmonics. The harmonic level is determined by two competing phenomena: nonlinear propagation and frequency dependent attenuation. Thus, the transmission frequency yielding the optimal trade-off between the spatial resolution and the penetration depth differs for THI and SHI. This paper quantitatively compares the concepts of fundamental, second harmonic, and superharmonic echocardiography at their optimal transmission frequencies. Forward propagation is modeled using a 3D-KZK implementation and the iterative nonlinear contrast source (INCS) method. Backpropagation is assumed to be linear. Results show that the fundamental lateral beamwidth is the narrowest at focus, while the superharmonic one is narrower outside the focus. The lateral superharmonic roll-off exceeds the fundamental and second harmonic roll-off. Also, the axial resolution of SHI exceeds that of FI and THI. The far-field pulse-echo superharmonic pressure is lower than that of the fundamental and second harmonic. SHI appears suited for echocardiography and is expected to improve its image quality at the cost of a slight reduction in depth-of-field.

Left ventricular border tracking using cardiac motion models and optical flow

The use of automated methods is becoming increasingly important for assessing cardiac function quantitatively and objectively. In this study, we propose a method for tracking three-dimensional (3-D) left ventricular contours. The method consists of a local optical flow tracker and a global tracker, which uses a statistical model of cardiac motion in an optical-flow formulation. We propose a combination of local and global trackers using gradient-based weights. The algorithm was tested on 35 echocardiographic sequences, with good results (surface error: 1.35 ± 0.46 mm, absolute volume error: 5.4 ± 4.8 mL). This demonstrates the method's potential in automated tracking in clinical quality echocardiograms, facilitating the quantitative and objective assessment of cardiac function.

Super-harmonic imaging: development of an interleaved phased-array transducer

For several years, the standard in ultrasound imaging has been second-harmonic imaging. A new imaging technique dubbed "super-harmonic imaging" (SHI) was recently proposed. It takes advantage of the higher - third to fifth - harmonics arising from nonlinear propagation or ultrasound-contrast-agent (UCA) response. Next to its better suppression of near-field artifacts, tissue SHI is expected to improve axial and lateral resolutions resulting in clearer images than second-harmonic imaging. When SHI is used in combination with UCAs, a better contrast-to-tissue ratio can be obtained. The use of SHI implies a large dynamic range and requires a sufficiently sensitive array over a frequency range from the transmission frequency up to its fifth harmonic (bandwidth > 130%). In this paper, we present the characteristics and performance of a new interleaved dual frequency array built chiefly for SHI. We report the rationale behind the design choice, frequencies, aperture, and piezomaterials used. The array is efficient both in transmission and reception with well-behaved transfer functions and a combined -6-dB bandwidth of 144%. In addition, there is virtually no contamination of the harmonic components by spurious transducer transmission, due to low element-to-element crosstalk (< 30 dB) and a low transmission efficiency of the odd harmonics (< 46 dB). The interleaved array presented in this article possesses ideal characteristics for SHI and is suitable for other methods like second-harmonic, subharmonic, and second-order ultrasound field (SURF) imaging.

Reconstructive compounding for IVUS palpography

This study proposes a novel algorithm for luminal strain reconstruction from sparse irregularly sampled strain measurements. It is based on the normalized convolution (NC) algorithm. The novel extension comprises the multilevel scheme, which takes into account the variable sampling density of the available strain measurements during the cardiac cycle. The proposed algorithm was applied to restore luminal strain values in intravascular ultrasound (IVUS) palpography. The procedure of reconstructing and averaging the strain values acquired during one cardiac cycle forms a technique, coined as reconstructive compounding. The accuracy of strain reconstruction was initially tested on the luminal strain map, computed from 3 in vivo IVUS pullbacks. The high quality of strain restoration was observed after systematically removing up to 90% of the initial elastographic measurements. The restored distributions accurately reproduced the original strain patterns and the error did not exceed 5%. The experimental validation of the reconstructed compounding technique was performed on 8 in vivo IVUS pullbacks. It demonstrated that the relative decrease in number of invalid strain estimates amounts to 92.05 +/- 6.03% and 99.17 +/- 0.92% for the traditional and reconstructive strain compounding schemes, respectively. In conclusion, implementation of the reconstructive compounding scheme boosts the diagnostic value of IVUS palpography.

Improving IVUS palpography by incorporation of motion compensation based on block matching and optical flow

Intravascular ultrasound (IVUS) strain imaging of the luminal layer in coronary arteries, coined as IVUS palpography, utilizes conventional radio frequency (RF) signals acquired at 2 different levels of a compressional load. The signals are cross-correlated to obtain the microscopic tissue displacements, which can be directly translated into local strain of the vessel wall. However, (apparent) tissue motion and nonuniform deformation of the vessel wall, due to catheter wiggling, reduce signal correlation and result in invalid strain estimates. Implications of probe motion were studied on the tissue-mimicking phantom. The measured circumferential tissue displacement and level of the speckle decorrelation amounted to 12 degrees and 0.58, respectively, for the catheter displacement of 456 microm. To compensate for the motion artifacts in IVUS palpography, a novel method based on the feature-based scale-space optical flow (OF), and classical block matching (BM) algorithm, were employed. The computed OF vector and BM displacement fields quantify the amount of local tissue misalignment in consecutive frames. Subsequently, the extracted circumferential displacements are used to realign the signals before strain computation. Motion compensation reduces the RF signal decorrelation and increases the number of valid strain estimates. The advantage of applying the motion correction in IVUS palpography was demonstrated in a midscale validation study on 14 in vivo pullbacks. Both methods substantially increase the number of valid strain estimates in the partial and compounded palpograms. Mean relative improvement in the number of valid strain estimates with motion compensation in comparison to one without motion compensation amounts to 28% and 14%, respectively. Implementation of motion compensation methods boosts the diagnostic value of IVUS palpography.

An inverse method for imaging the local elasticity of atherosclerotic coronary plaques

The rupture of thin-cap fibroatheroma (TCFA) plaques is a major cause of acute coronary events. A TCFA has a trombogenic soft lipid core, shielded from the blood stream by a thin, possibly inflamed, stiff cap. The majority of atherosclerotic plaques resemble a TCFA in terms of overall structural composition, but have a more complex, heterogeneous morphology. An assessment of the material distribution is vital for quantifying the plaque's mechanical stability and for determining the effect of plaque-stabilizing pharmaceutical agents. We describe a new automated inverse elasticity method, intravascular ultrasound (IVUS) modulography, which is capable of reconstructing a heterogeneous Young's modulus distribution. The elastogram (i.e., spatial strain distribution) of the plaque is the input for the method, and is measured using the clinically available technique, IVUS elastography. Our method incorporates a novel divide-and-conquer strategy, allowing the reconstruction of TCFAs as well as heterogeneous plaques with localized regions of soft, weakened tissue. The method was applied to ex vivo elastograms, which were simulated from the cross sections of postmortem human coronary plaques. To demonstrate the clinical feasibility of the method, measured elastograms from human atherosclerotic coronary arteries were analyzed. One elastogram was measured in vitro; the other, in vivo. The method approximated the true Young's modulus distribution of all simulated plaques, while the in vitro reconstruction was in agreement with histology. In conclusion, the IVUS modulography in combination with the IVUS elastography has strong potential to become an all-encompassing modality for detecting plaques, for assessing the information related to their rupture-proneness, and for imaging their heterogeneous elastic material composition.

Accuracy in prediction of catheter rotation in IVUS with feature-based optical flow--a phantom study

The quantitative assessment of and compensation for catheter rotation in intravascular ultrasound images presents a fundamental problem for noninvasive characterization of the mechanical properties of the coronary arteries. A method based on the scale-space optical flow algorithm with a feature-based weighting scheme is proposed to account for the aforementioned artifact. The computed vector field, describing the misalignment between two consecutive frames, allows the quantitative assessment of the amount of vessel wall tissue motion, which is directly related to the catheter rotation. Algorithm accuracy and robustness were demonstrated on two tissue-mimicking phantoms, subjected to controlled amount of angular deviation. The proposed method shows a great reliability in the prediction of catheter rotational motion up to 4 degrees.

SPASM: A 3D-ASM for segmentation of sparse and arbitrarily oriented cardiac MRI data

A new technique (SPASM) based on a 3D-ASM is presented for automatic segmentation of cardiac MRI image data sets consisting of multiple planes with arbitrary orientations, and with large undersampled regions. Model landmark positions are updated in a two-stage iterative process. First, landmark positions close to intersections with images are updated. Second, the update information is propagated to the regions without image information, such that new locations for the whole set of the model landmarks are obtained. Feature point detection is performed by a fuzzy inference system, based on fuzzy C-means clustering. Model parameters were optimized on a computer cluster and the computational load distributed by grid computing. SPASM was applied to image data sets with an increasing sparsity (from 2 to 11 slices) comprising images with different orientations and stemming from different MRI acquisition protocols. Segmentation outcomes and calculated volumes were compared to manual segmentation on a dense short-axis data configuration in a 3D manner. For all data configurations, (sub-)pixel accuracy was achieved. Performance differences between data configurations were significantly different (p<0.05) for SA data sets with less than 6 slices, but not clinically relevant (volume differences<4 ml). Comparison to results from other 3D model-based methods showed that SPASM performs comparable to or better than these other methods, but SPASM uses considerably less image data. Sensitivity to initial model placement proved to be limited within a range of position perturbations of approximately 20 mm in all directions.

Influence of positional and angular variation of automatically planned short-axis stacks on quantification of left ventricular dimensions and function with cardiovascular magnetic resonance

PURPOSE

To theoretically and experimentally investigate the influence of the automated cardiovascular magnetic resonance (CMR) scan planning pitfalls, namely inaccurate positioning and tilting of short-axis (SA) imaging planes, on quantification of the left ventricular (LV) dimensions and function.

MATERIALS AND METHODS

Eleven healthy subjects and eight patients underwent CMR. Manually and automatically planned SA sets were acquired. To obtain the quantitative measurements of LV function, one observer performed image analysis twice. The agreement between planning methods, as well as the decomposition of the total variation into interstudy and intraobserver components was measured.

RESULTS

The decomposition of the total variation showed that the interstudy factor accounts for 70-85% of the total variation, while the rest is due to the intraobserver factor. Moreover, the relative contribution of the interstudy factor remains independent from errors in slice positioning and small angular deviation of SA stacks from the optimal orientation. Good agreement between the theoretical and measured variability factors was observed.

CONCLUSION

Global LV function derived from the automatically planned CMR acquisitions yield accurate quantification of the human cardiovascular system. Inaccurate positioning and tilting of SA images does not affect the quantitative measurements of LV function. The computer-aided system for automated CMR has proven clinical applicability.

Operator induced variability in cardiovascular MR: left ventricular measurements and their reproducibility

PURPOSE

To assess the intra- and inter-operator variability of the manual planning of cardiovascular magnetic resonance imaging and to evaluate the influence of these factors on the functional parameters of the left ventricle (LV).

METHOD

The study population consisted of 10 healthy volunteers. For each subject the manual planning of the short-axis cine acquisitions was carried out twice by one operator and once by a second operator. Left ventricular volume, mass, and function were manually evaluated twice by one experienced observer, resulting in an approximation of the intra-observer variability factor. The intra- and inter-operator variation factors were estimated as the difference between the total and intra-observer variation components.

RESULTS

LV end-diastolic volume varied by 3.3% and 4.16%, and LV end-systolic volume by 5.84% and 6.23% for intra- and inter-operator studies, respectively. The variability for LV mass at end-diastole was equal to 4.23% in both studies. For the ejection fraction the variability was 3.56% and 2.97% for intra- and inter-operator studies, respectively. Comparison of reproducibility between intra- and inter-operator studies resulted in insignificant statistical differences. Bland-Altman limits of agreements revealed no systematic bias in differences between measurements with respect to their means. Reliability of the planning expressed as the angular deviation of the short-axis imaging planes amounts to 2.67 -/+ 1.5 degrees and 4.99 +/- 2.17 degrees for the intra-operator and inter-operator studies, respectively. For EDV, ESV, and EF approximately 75-80% of the total variation can be explained by the within or between operator variation, while the same percentage is 60% for LVM.

CONCLUSIONS

Our study confirms the excellent inter- and intra-operator reproducibility of the cardiovascular magnetic resonance measurements of the left ventricular volumes and mass in a group of healthy volunteers.

Optimizing the automatic segmentation of the left ventricle in magnetic resonance images

Automatic segmentation of the left ventricular (LV) myocardial borders in cardiovascular MR (CMR) images allows a significant speed-up of the procedure of quantifying LV function, and improves its reproducibility. The automated boundary delineation is usually based on a set of parameters that define the algorithms. Since the automatic segmentation algorithms are usually sensitive to the image quality and frequently depend heavily on the acquisition protocol, optimizing the parameters of the algorithm for such different protocols may be necessary to obtain optimal results. In other words, using a default set of parameters may be far from optimal for different scanners or protocols. For the MASS-software, for example, this means that a total of 14 parameters need to be optimized. This optimization is a difficult and labor-intensive process. To be able to more consistently and rapidly tune the parameters, an automated optimization system would be extremely desirable. In this paper we propose such an approach, which is based on genetic algorithms (GAs). The GA is an unsupervised iterative tool that generates new sets of parameters and converges toward an optimal set. We implemented and compared two different types of the genetic algorithms: a simple GA (SGA) and a steady state GA (2SGA). The difference between these two algorithms lies in the characteristics of the generated populations: "nonoverlapping populations" and "overlapping populations," respectively "nonoverlapping" population means that the two populations are disjoint, and "overlapping" means that the best parameters found in the previous generation are included in the present population. The performance of both algorithms was evaluated on twenty routinely obtained short-axis examinations (eleven examinations acquired with a steady-state free precession pulse sequence, and nine examinations with a gradient echo pulse sequence). The optimal parameters obtained with the GAs were used for the LV myocardial border delineation. Finally, the automatically outlined contours were compared to the gold standard--manually drawn contours by experts. The result of the comparison was expressed as a degree of similarity after a processing time of less than 72 h to a 59.5% of degree of similarity for SGA and a 66.7% of degree of similarity for 2SGA. In conclusion, genetic algorithms are very suitable to automatically tune the parameters of a border detection algorithm. Based on our data, the 2SGA was more suitable than the SGA method. This approach can be generalized to other optimization problems in medical image processing.

Accuracy of short-axis cardiac MRI automatically derived from scout acquisitions in free-breathing and breath-holding modes

To qualitatively assess the accuracy of automated cardiovascular magnetic resonance planning procedures devised from scout acquisitions in free-breathing and breath-holding modes, to quantitatively evaluate the accuracy of the derived left ventricular volumes, mass and function and compare these parameters with the ones obtained from the manually planned acquisitions. Ten healthy volunteers underwent cardiovascular MR (CMR) acquisitions for ventricular function assessment. Short-axis data sets of the left ventricle (LV) were manually planned and generated twice in an automatic fashion. Automated planning parameters were derived from gated scout acquisitions in free-breathing and breath-holding modes. End-diastolic volume (EDV), end-systolic volume (ESV), ejection fraction (EF), and left ventricular mass (LVM) were measured. The agreement between the manual and automatic planning methods, as well as the variability of the aforementioned measurements were assessed. The differences between two automated planning methods were also compared. The mean differences between the manual and automated CMR planning derived from gated scouts in free-breathing mode were 8.05 ml (EDV), 1.84 ml (ESV), 0.69% (EF), and 4.72 g (LVM). The comparison between manual and automated CMR planning derived from gated scouts in breath-holding mode yielded the following differences: 4.22 ml (EDV), 0.34 ml (ESV), 0.3% (EF), and -0.72 mg (LVM). The variability coefficients were 3.72 and 3.66 (EDV), 5.6 and 8.19 (ESV), 3.46 and 4.31 (EF), 6.49 and 5.20 (LVM) for the automated CMR planning methods derived from scouts in free-breathing and breath-holding modes, respectively. Automated CMR planning methods can provide accurate measurements of LV dimensions in normal subjects, and therefore may be utilized in the clinical environment to provide a cost-effective solution for functional assessment of the human cardiovascular system.

3D model-based approach to lung registration and prediction of respiratory cardiac motion

This paper presents a new approach for lung registration and cardiac motion prediction, based on a 3D geometric model of the left lung. Feature points, describing a shape of this anatomical object, are automatically extracted from acquired tomographic images. The "goodness-of-fit" measure is assessed at each step in the iterative scheme until spatial alignment between the model and subject's specific data is achieved. We applied the proposed methods to register the 3D lung surfaces of 5 healthy volunteers of thoracic MRI acquired in different respiratory phases. We also utilized this approach to predict the spatial displacement of the human heart due to respiration. The obtained results demonstrate a promising registration performance.

Automated Short-Axis Cardiac Magnetic Resonance Image Acquisitions

RATIONALE AND OBJECTIVE

This study investigates the use of an automated observer-independent planning system for short-axis cardiovascular magnetic resonance (MR) acquisitions in the clinical environment. The capacity of the automated method to produce accurate measurements of left ventricular dimensions and function was quantitatively assessed in normal subjects and patients.

METHODS

Fourteen healthy volunteers and 8 patients underwent cardiovascular MR (CMR) acquisitions for ventricular function assessment. Short-axis datasets of the left ventricle (LV) were acquired in 2 ways: manually planned and generated in an automatic fashion. End-diastolic volume (EDV), end-systolic volume (ESV), ejection fraction (EF), and left ventricular mass (LVM) were derived from the 2 datasets. The agreement between the manual and automatic planning methods was assessed.

RESULTS

The mean differences between the manual and automated CMR planning methods for the normal subjects and patients were 5.89 mL and 1.93 mL (EDV), 1.14 mL and -0.41 mL (ESV), 0.81% and 0.89% (EF), and 4.35 g and 3.88 g (LVM), respectively. There was no significant difference in ESV and EF. LVM significantly differed in both groups, whereas EDV was significantly different in the normal subjects and insignificantly different in the patients. The variability coefficients were 2.8 and 3.59 (EDV), 3.3 and 5.03 (ESV), 1.79 and 2.65 (EF), and 4.36 and 2.27 (LVM) for the normal subjects and patients, respectively. The mean angular deviation of the LV axes turned out to be 8.58 +/- 5.76 degrees for the normal subjects and 8.35 +/- 5.15 degrees for the patients.

CONCLUSIONS

Automated CMR planning method can provide accurate measurements of LV dimensions in normal subjects and patients, and therefore, can be used in the clinical environment for functional assessment of the human cardiovascular system.

Mitral valve regurgitation: accurate blood flow quantification with MRI

BACKGROUND

The quantification of transvalvular blood flow through the mitral valve (MV) and regurgitant flow in particular is difficult with echocardiography, which is the method of choice to diagnose patients selected for valve repair or replacement. With magnetic resonance imaging, information on the intraventricular blood flow can be obtained. Several scanning techniques have attempted to assess the regurgitant flow. These techniques either do not directly assess the complete flow through the MV, or they do not measure the flow at the location of the valve.

AIM

To investigate the accuracy of a novel method using three-directional velocity-encoded MRI to acquire the transvalvular blood flow directly from the intraventricular blood flow field, also representing the regurgitant flow during systole.

METHODS

Ten volunteers without cardiac valvular disease were recruited. The transvalvular MV flow volume was measured with three-directional velocity-encoded MRI (3-dir MV flow).

RESULTS

The transvalvular flow measurements correlate very well with the flow measured in the aorta (r_p=0.92, p<0.01). The small differences (mean -5±7 ml) are insignificant (p=0.06) and demonstrate the high accuracy of the new method. Intra- and inter-observer studies showed non-significant mean differences of 0.9±5.1 ml and 1.3±5.6 ml, respectively, thereby proving the high reproducibility.

CONCLUSION

Three-directional velocity-encoded MRI is a patient-friendly and easy-to-use method suitable for quantifying the regurgitant MV flow in clinical practice.