Displacement and velocity of the coronary arteries: cardiac and respiratory motion

Cardiac substructure delineation in radiation therapy - A state-of-the-art review

Delineation of cardiac substructures is crucial for a better understanding of radiation-related cardiotoxicities and to facilitate accurate and precise cardiac dose calculation for developing and applying risk models. This review examines recent advancements in cardiac substructure delineation in the radiation therapy (RT) context, aiming to provide a comprehensive overview of the current level of knowledge, challenges and future directions in this evolving field. Imaging used for RT planning presents challenges in reliably visualising cardiac anatomy. Although cardiac atlases and contouring guidelines aid in standardisation and reduction of variability, significant uncertainties remain in defining cardiac anatomy. Coupled with the inherent complexity of the heart, this necessitates auto-contouring for consistent large-scale data analysis and improved efficiency in prospective applications. Auto-contouring models, developed primarily for breast and lung cancer RT, have demonstrated performance comparable to manual contouring, marking a significant milestone in the evolution of cardiac delineation practices. Nevertheless, several key concerns require further investigation. There is an unmet need for expanding cardiac auto-contouring models to encompass a broader range of cancer sites. A shift in focus is needed from ensuring accuracy to enhancing the robustness and accessibility of auto-contouring models. Addressing these challenges is paramount for the integration of cardiac substructure delineation and associated risk models into routine clinical practice, thereby improving the safety of RT for future cancer patients.

Motion artifact correction in cardiac CT using cross-phase temporospatial information and synergistic attention gate and spatial transformer sub-networks

Objectives.To improve quality of coronary CT angiography (CCTA) images using a generalizable motion-correction algorithm.Approach. A neural network with attention gate and spatial transformer (ATOM) was developed to correct coronary motion. Phantom and patient CCTA images (39 males, 32 females, age range 19-92, scan date 02/2020 to 10/2021) retrospectively collected from dual-source CT were used to create training, development, and testing sets corresponding to 140- and 75 ms temporal resolution, with 75 ms images as labels. To test generalizability, ATOM was deployed for locally adaptive motion-correction in both 140- and 75 ms patient images. Objective metrics were used to assess motion-corrupted and corrected phantom and patient images, including structural-similarity-index (SSIM), dice-similarity-coefficient (DSC), peak-signal-noise-ratio (PSNR), and normalized root-mean-square-error (NRMSE). In objective quality assessment, ATOM was compared with several baseline networks, including U-net, U-net plus attention gate, U-net plus spatial transformer, VDSR, and ResNet. Two cardiac radiologists independently interpreted motion-corrupted and -corrected images at 75 and 140 ms in a blinded fashion and ranked diagnostic image quality (worst to best: 1-4, no ties).Main results. ATOM improved quality metrics (p< 0.05) before/after correction: in phantom, SSIM 0.87/0.95, DSC 0.85/0.93, PSNR 19.4/22.5, NRMSE 0.38/0.27; in patient images, SSIM 0.82/0.88, DSC 0.88/0.90, PSNR 30.0/32.0, NRMSE 0.16/0.12. ATOM provided more consistent improvement of objective image quality, compared to the presented baseline networks. The motion-corrected images received better ranks than un-corrected at the same temporal resolution (p< 0.05): 140 ms images 1.65/2.25, and 75 ms images 3.1/3.2. The motion-corrected 75 ms images received the best rank in 65% of testing cases. A fair-to-good inter-reader agreement was observed (Kappa score 0.58).Significance. ATOM reduces motion artifacts, improving visualization of coronary arteries. This algorithm can be used to virtually improve temporal resolution in both single- and dual-source CT.

Useful Respiratory Maneuvers Aiding Left Heart Cardiac Catheterization and Intervention. A Comprehensive Review

Left heart catheterizations, coronary angiography, and coronary interventions are important common cardiac procedures. Performing a successful cardiac catheterization and intervention and proper catheterization and device delivery is not always without difficulties, especially in the context of calcification or vessel tortuosity. Although there are some techniques to overcome these issues, performing respiratory maneuvers (inspiration or expiration) can be simply tried as the first step to increase successful procedures which is underreported and underutilized. The goal of this article is to review current literature regarding useful respiratory maneuvers that can aid left heart cardiac catheterization, coronary angiography, and intervention for a successful procedure.

Higher Order Singular Value Decomposition Filter for Contrast Echocardiography

Assessing the coronary circulation with contrast-enhanced echocardiography has high clinical relevance. However, it is not being routinely performed in clinical practice because the current clinical tools generally cannot provide adequate image quality. The contrast agent's visibility in the myocardium is generally poor, impaired by motion and nonlinear propagation artifacts. The established multipulse contrast schemes (MPCSs) and the more experimental singular value decomposition (SVD) filter also fall short to solve these issues. Here, we propose a scheme to process amplitude modulation/amplitude-modulated pulse inversion (AM/AMPI) echoes with higher order SVD (HOSVD) instead of conventionally summing the complementary pulses. The echoes from the complementary pulses form a separate dimension in the HOSVD algorithm. Then, removing the ranks in that dimension with dominant coherent signals coming from tissue scattering would provide the contrast detection. We performed both in vitro and in vivo experiments to assess the performance of our proposed method in comparison with the current standard methods. A flow phantom study shows that HOSVD on AM pulsing exceeds the contrast-to-background ratio (CBR) of conventional AM and an SVD filter by 10 and 14 dB, respectively. In vivo porcine heart results also demonstrate that, compared to AM, HOSVD improves CBR in open-chest acquisition (up to 19 dB) and contrast ratio (CR) in closed-chest acquisition (3 dB).

High-pitch, high temporal resolution, multi-energy cardiac imaging on a dual-source photon-counting-detector CT

OBJECTIVE

To measure the accuracy of material decomposition using a dual-source photon-counting-detector (DS-PCD) CT operated in the high-pitch helical scanning mode and compare the results against dual-source energy-integrating-detector (DS-EID) CT, which requires use of a low-pitch value in dual-energy mode.

METHODS

A DS-PCD CT and a DS-EID CT were used to scan a cardiac motion phantom consisting of a 3-mm diameter iodine cylinder. Iodine maps were reconstructed using DS-PCD in high-pitch mode and DS-EID in low-pitch mode. Image-based circularity, diameter, and iodine concentration of the iodine cylinder were calculated and compared between the two scanners. With institutional review board approval, in vivo exams were performed with the DS-PCD CT in high-pitch mode. Images were qualitatively compared against patients with similar heart rates that were scanned with DS-EID CT in low-pitch dual-energy mode.

RESULTS

On iodine maps, the mean circularity was 0.97 ± 0.02 with DS-PCD in high-pitch mode and 0.95 ± 0.06 with DS-EID in low-pitch mode. The mean diameter was 2.9 ± 0.2 mm with DS-PCD and 3.1 ± 0.2 mm with DS-EID, both of which are close to the 3 mm ground truth. For DS-PCD, the mean iodine concentration was 9.6 ± 0.8 mg/ml and this was consistent with the 9.4 ± 0.6 mg/ml value obtained with the cardiac motion disabled. For DS-EID, the concentration was 12.7 ± 1.2 mg/ml with motion enabled and 11.7 ± 0.5 mg/ml disabled. The background noise in the iodine maps was 15.1 HU with DS-PCD and 14.4 HU with DS-EID, whereas the volume CT dose index (CTDIvol ) was 3 mGy with DS-PCD and 11 mGy with DS-EID. On comparison of six patients (three on PCD, three on EID) with similar heart rates, DS-PCD provided iodine maps with well-defined coronaries even at a high heart rate of 86 beats per minute. Meanwhile, there were substantial motion artifacts in iodine maps obtained with DS-EID for patients with similar heart rates.

CONCLUSION

In a cardiac motion phantom, DS-PCD CT can perform accurate material decomposition in high-pitch mode, providing iodine maps with excellent geometric accuracy and robustness to motion at approximately 38% of the dose for similar noise as DS-EID CT.

Heart Rate Lowering Significantly Increases Feasibility in Doppler Recording Blood Flow Velocity in Coronaries during Transthoracic Doppler Echocardiography

Background: Coronary blood flow Doppler recording by Transthoracic Doppler in convergent mode (E-Doppler TTE) might be further improved by lowering heart rate (HRL) down to <60 bpm, since low HR < 60 b/m causes a disproportional lengthening of the diastole, so the coronaries are still for a longer time, very much improving the Doppler signal/noise ratio. Methods: A group of 26 patients underwent E-Doppler TTE before and after HR lowering in four branches of the coronary tree, namely, the left main (LMCA); left anterior descending (LAD), subdivided into three segments: proximal, mid and distal; proximal left circumflex (LCx); and obtuse marginal (OM). Color and PW coronary Doppler signal was judged by two expert observers as undetectable (SCORE 1), weak or with clutter artifacts (SCORE 2), or well delineated (SCORE 3). In addition, local accelerated stenotic flow (AsF) was measured in the LAD before and after HRL. Results: Beta-blockers significantly decreased the mean HR from 76 ± 5 to 57 ± 6 bpm (p < 0.001). Before HRL, the Doppler quality was very poor in the proximal and mid-LAD segments (median score value = 1 in both), while in the distal LAD, it was significantly better but still suboptimal (median score value = 1.5, p = 0.009 vs. proximal and mid-LAD score). After HRL, blood flow Doppler recording in the three LAD segments was strikingly improved (median score value = 3, 3 and 3, p = ns), so the effect of HRL was more efficacious in the two more proximal LAD segments. In 10 patients undergoing coronary angiography (CA), no AsF as expression of transtenotic velocity was detected at baseline. After HRL, thanks to the better quality and length of color flow, ASF was detected in five patients while in five others, it was not in perfect agreement with CA (Spearman correlation coefficient = 1, p < 0.01). The color flow in the proximal LCx and OM was extremely poor at baseline (color flow length 0 and 0, median (interquartile range) mm, respectively) and improved considerably after HRL (color flow length 23 [13.5] and 25 [12.0] mm, respectively, p < 0.001). Conclusions: HRL greatly improved the success rate of blood flow Doppler recording in coronaries, not only in the LAD, but also in the LCx. Therefore, AsF for stenosis detection and coronary flow reserve assessment can have wider clinical applications. However, further studies with larger samples are needed to confirm these results.

Evaluation of the feasibility of cardiac gating for SBRT of ventricular tachycardia based on real-time ECG signal acquisition

PURPOSE

To investigate the feasibility of cardiac synchronized gating in stereotactic body radiation therapy (SBRT) of ventricular tachycardia (VT) using a real-time electrocardiogram (ECG) signal acquisition.

METHODS AND MATERIALS

Stability of beam characteristics during simulated ECG gating was examined by developing a microcontroller interface to a Varian Clinac iX linear accelerator allowing gating at frequencies and duty cycles relevant to cardiac rhythm. Delivery accuracy was evaluated by measuring dose linearity with an ionization chamber, and flatness and symmetry with a two-dimensional detector array, for different gating windows within typical human cardiac cycle periods. To establish a practical method of gating based on actual ECG signals, an AD8232 Heart Monitor board was used to acquire the ECG signal and synchronize the beam delivery. Real-time cardiac gated delivery measurements were performed for a single 10 × 10 cm² field and for a VT-SBRT plan using intensity-modulated radiation therapy (IMRT).

RESULTS AND DISCUSSION

Dose per monitor unit (MU) values were found to be linear within most gating windows investigated with maximum differences relative to non-gated delivery of <2% for gating windows ≥200 ms and for >10 MUs. Beam profiles for both gated and non-gated modes were also found to agree with maximum differences of 0.5% relative to central axis dose for all sets of beam-on/beam-off combinations. Comparison of dose distributions for intensity-modulated SBRT plans between non-gating and cardiac gating modes provided a gamma passing rate of 97.2% for a 2% 2 mm tolerance.

CONCLUSIONS

Beam output is stable with respect to linearity, flatness, and symmetry for gating windows within cardiac cycle periods. Agreement between dose distributions for VT-SBRT using IMRT in non-gated and cardiac cycle gated delivery modes shows that the proposed methodology is feasible. Technically, gating for delivery of SBRT for VT is possible with regard to beam stability.

Linking the region-specific tissue microstructure to the biaxial mechanical properties of the porcine left anterior descending artery

Coronary atherosclerosis is the main cause of death worldwide. Advancing the understanding of coronary microstructure-based mechanics is fundamental for the development of therapeutic tools and surgical procedures. Although the passive biaxial properties of the coronary arteries have been extensively explored, their regional differences and the relationship between tissue microstructure and mechanics have not been fully characterized. In this study, we characterized the passive biaxial mechanical properties and microstructural properties of the proximal, medial, and distal regions of the porcine left anterior descending artery (LADA). We also attempted to relate the biaxial stress-stretch response of the LADA and its respective birefringent responses to the polarized light for obtaining information about the load-dependent microstructural variations. We found that the LADA extensibility is reduced in the proximal-to-distal direction and that the medial region exhibits more heterogeneous mechanical behavior than the other two regions. We have also observed highly dynamic microstructural behavior where fiber families realign themselves depending on loading. In addition, we found that the microstructure of the distal region exhibited highly aligned fibers along the longitudinal axis of the artery. To verify this microstructural feature, we imaged the LADA specimens with multi-photon microscopy and observed that the adventitia microstructure transitioned from a random fiber network in the proximal region to highly aligned fibers in the distal region. Our findings could offer new perspectives for understanding coronary mechanics and aid in the development of tissue-engineered vascular grafts, which are currently limited due to their mismatch with native tissue in terms of mechanical properties and microstructural features. STATEMENT OF SIGNIFICANCE: The tissue biomechanics of coronary arteries is fundamental for the development of revascularization techniques such as coronary artery bypass. These therapeutics require a deep understanding of arterial mechanics, microstructure, and mechanobiology to prevent graft failure and reoperation. The present study characterizes the unique regional mechanical and microstructural properties of the porcine left anterior descending artery using biaxial testing, polarized-light imaging, and confocal microscopy. This comprehensive characterization provides an improved understanding of the collagen/elastin architecture in response to mechanical loads using a region-specific approach. The unique tissue properties obtained from this study will provide guidance for the selection of anastomotic sites in coronary artery bypass grafting and for the design of tissue-engineered vascular grafts.

Independent Component Analysis Filter for Small Vessel Contrast Imaging During Fast Tissue Motion

Suppressing tissue clutter is an essential step in blood flow estimation and visualization, even when using ultrasound contrast agents. Blind source separation (BSS)-based clutter filter for high-framerate ultrasound imaging has been reported to perform better in tissue clutter suppression than the conventional frequency-based wall filter and nonlinear contrast pulsing schemes. The most notable BSS technique, singular value decomposition (SVD) has shown compelling results in cases of slow tissue motion. However, its performance degrades when the tissue motion is faster than the blood flow speed, conditions that are likely to occur when imaging the small vessels, such as in the myocardium. Independent component analysis (ICA) is another BSS technique that has been implemented as a clutter filter in the spatiotemporal domain. Instead, we propose to implement ICA in the spatial domain where motion should have less impact. In this work, we propose a clutter filter with the combination of SVD and ICA to improve the contrast-to-background ratio (CBR) in cases where tissue velocity is significantly faster than the flow speed. In an in vitro study, the range of fast tissue motion velocity was 5-25 mm/s and the range of flow speed was 1-12 mm/s. Our results show that the combination of ICA and SVD yields 7-10 dB higher CBR than SVD alone, especially in the tissue high-velocity range. The improvement is crucial for cardiac imaging where relatively fast myocardial motions are expected.

Respiration-averaged CT versus standard CT attenuation map for correction of 18F-sodium fluoride uptake in coronary atherosclerotic lesions on hybrid PET/CT

BACKGROUND

To evaluate the impact of respiratory-averaged computed tomography attenuation correction (RACTAC) compared to standard single-phase computed tomography attenuation correction (CTAC) map, on the quantitative measures of coronary atherosclerotic lesions of 18F-sodium fluoride (18F-NaF) uptake in hybrid positron emission tomography and computed tomography (PET/CT).

METHODS

This study comprised 23 patients who underwent 18F-NaF coronary PET in a hybrid PET/CT system. All patients had a standard single-phase CTAC obtained during free-breathing and a 4D cine-CT scan. From the cine-CT acquisition, RACTAC maps were obtained by averaging all images acquired over 5 seconds. PET reconstructions using either CTAC or RACTAC were compared. The quantitative impact of employing RACTAC was assessed using maximum target-to-background (TBRMAX) and coronary microcalcification activity (CMA). Statistical differences were analyzed using reproducibility coefficients and Bland-Altman plots.

RESULTS

In 23 patients, we evaluated 34 coronary lesions using CTAC and RACTAC reconstructions. There was good agreement between CTAC and RACTAC for TBRMAX (median [Interquartile range]): CTAC = 1.65 [1.23 to 2.38], RACTAC = 1.63 [1.23 to 2.33], p = 0.55), with coefficient of reproducibility of 0.18, and CMA: CTAC = 0.10 [0 to 1.0], RACTAC = 0.15 [0 to 1.03], p = 0.55 with coefficient of reproducibility of 0.17 CONCLUSION: Respiratory-averaged and standard single-phase attenuation correction maps provide similar and reproducible methods of quantifying coronary 18F-NaF uptake on PET/CT.

Quantitative assessment of motion effects in dual-source dual-energy CT and dual-source photon-counting detector CT

Conventional dual-source CT scanners can be used to either provide better temporal resolution or dual-energy imaging, but not both at the same time. This presents a dilemma in cardiac CT as both high temporal resolution and multi-energy imaging are desirable. The current study evaluated a dual-source photon-counting-detector (DS-PCD) CT which can acquire multi-energy images at high temporal resolution. A cardiac motion phantom with a 3-mm diameter iodinated rod, mimicking the right coronary artery, was scanned 25 times using a DS-PCD CT (66 ms resolution) and a dual-source dual-energy (DS-DE, 125 ms resolution) CT. Low/high energy images and iodine maps were reconstructed at 40% and 75% cardiac phases. To quantify the impact of motion on image quality, dice similarity coefficient was computed between the low/high energy images while the circularity and effective diameter of the iodinated rod were computed on the iodine maps. The dice coefficients were higher for DS-PCD with a mean of 0.89 and 0.91 at the 40% and 70% phases, while DS-DE had a lower mean of 0.20 and 0.78, respectively. The circularity was excellent for DS-PCD with a mean of 0.97 and 0.98 at the 40% and 75% phases, while DS-DE had a mean of 0.71 and 0.98, respectively. The effective diameter was accurate for DS-PCD with a mean of 2.9 mm (true size of 3 mm) at both phases, while DS-DE had a mean of 4.0 mm and 3.2 mm at the 40% and 75% phases, respectively. These results indicate that DS-PCD CT enables simultaneous high temporal resolution and multi-energy cardiac imaging with minimal motion artifacts.

Aortic valve imaging using 18F-sodium fluoride: impact of triple motion correction

BACKGROUND

Current 18F-NaF assessments of aortic valve microcalcification using 18F-NaF PET/CT are based on evaluations of end-diastolic or cardiac motion-corrected (ECG-MC) images, which are affected by both patient and respiratory motion. We aimed to test the impact of employing a triple motion correction technique (3 × MC), including cardiorespiratory and gross patient motion, on quantitative and qualitative measurements.

MATERIALS AND METHODS

Fourteen patients with aortic stenosis underwent two repeat 30-min PET aortic valve scans within (29 ± 24) days. We considered three different image reconstruction protocols; an end-diastolic reconstruction protocol (standard) utilizing 25% of the acquired data, an ECG-gated (four ECG gates) reconstruction (ECG-MC), and a triple motion-corrected (3 × MC) dataset which corrects for both cardiorespiratory and patient motion. All datasets were compared to aortic valve calcification scores (AVCS), using the Agatston method, obtained from CT scans using correlation plots. We report SUVmax values measured in the aortic valve and maximum target-to-background ratios (TBRmax) values after correcting for blood pool activity.

RESULTS

Compared to standard and ECG-MC reconstructions, increases in both SUVmax and TBRmax were observed following 3 × MC (SUVmax: Standard = 2.8 ± 0.7, ECG-MC = 2.6 ± 0.6, and 3 × MC = 3.3 ± 0.9; TBRmax: Standard = 2.7 ± 0.7, ECG-MC = 2.5 ± 0.6, and 3 × MC = 3.3 ± 1.2, all p values ≤ 0.05). 3 × MC had improved correlations (R2 value) to the AVCS when compared to the standard methods (SUVmax: Standard = 0.10, ECG-MC = 0.10, and 3 × MC = 0.20; TBRmax: Standard = 0.20, ECG-MC = 0.28, and 3 × MC = 0.46).

CONCLUSION

3 × MC improves the correlation between the AVCS and SUVmax and TBRmax and should be considered in PET studies of aortic valves using 18F-NaF.

A two-stage cardiac PET and late gadolinium enhancement MRI co-registration method for improved assessment of non-ischemic cardiomyopathies using integrated PET/MR

PURPOSE

Respiratory motion causes mismatches between PET images of the myocardium and the corresponding cardiac MR images in cardiac integrated PET/MR. The mismatch may affect the attenuation correction and the diagnosis of non-ischemic cardiomyopathies. In this study, we present a two-stage cardiac PET and MR late gadolinium enhancement (LGE) co-registration method, which seeks to improve diagnostic accuracy of non-ischemic cardiomyopathies via better image co-registration using an integrated whole-body PET/MR system.

METHODS

The proposed PET and LGE two-stage co-registration method was evaluated through comparison with one-stage direct co-registration and no-registration. One hundred and ninety-one slices of LGE and forty lesions were studied. Two trained nuclear medicine physicians independently assessed the displacement between LGE and PET to qualitatively evaluate the co-registration quality. The changes of the mean SUV in the normal myocardium and the LGE-enhanced lesions before and after image co-registration were measured to quantitatively evaluate the accuracy and value of image co-registration.

RESULTS

The two-stage method had an improved image registration score (4.93 ± 0.89) compared with the no-registration method (3.49 ± 0.84, p value < 0.001) and the single-stage method (4.23 ± 0.81, p value < 0.001). Furthermore, the two-stage method led to increased SUV value in the myocardium (3.87 ± 2.56) compared with the no-registration method (3.14 ± 1.92, p value < 0.001) and the single-stage method (3.32 ± 2.16, p value < 0.001). The mean SUV in the LGE lesion significantly increased from 2.51 ± 2.09 to 2.85 ± 2.35 (p value < 0.001) after the two-stage co-registration.

CONCLUSION

The proposed two-stage registration method significantly improved the co-registration between PET and LGE in integrated PET/MR imaging. The technique may improve diagnostic accuracy of non-ischemic cardiomyopathies via better image co-registration.

REGISTERED NO

DF-2020-085,2020.04.30.

Motion compensation for aortic valves using partial angle CT reconstructions motion compensation of cardiac valve CT

PURPOSE

A motion compensation method that is aimed at correcting motion artifacts of cardiac valves is proposed. The primary focus is the aortic valve.

METHODS

The method is based around partial angle reconstructions and a cost function including the image entropy. A motion model is applied to approximate the cardiac motion in the temporal and spatial domain. Based on characteristic values for velocities and strain during cardiac motion, penalties for the velocity and spatial derivatives are introduced to maintain anatomically realistic motion vector fields and avoid distortions. The model addresses global elastic deformation, but not the finer and more complicated motion of the valve leaflets.

RESULTS

The method is verified based on clinical data. Image quality was improved for most artifact impaired reconstructions. An image quality study with Likert scoring of the motion artifact severity on a scale from 1 (highest image quality) to 5 (lowest image quality/extreme artifact presence) was performed. The biggest improvements after applying motion compensation were achieved for strongly artifact impaired initial images scoring 4 and 5, resulting in an average change of the scores by -0.59 ± 0.06 and -1.33 ± 0.03, respectively. In case of artifact free images, a chance to introduce blurring was observed and their average score was raised by 0.42 ± 0.03.

CONCLUSION

Motion artifacts were consistently removed and image quality improved.

Effect of Hinge Motion on Stent Edge-Related Restenosis After Right Coronary Artery Treatment in the Current Drug-Eluting Stent Era

BACKGROUND

Stent edge-related restenosis (SER) remains a potential limitation of drug-eluting stent (DES). Hinge motion at the stent edge could lead to mechanical stress and contribute to incidents of SER. We investigated the effect of hinge motion on SER after implantation of current-generation DES in the right coronary artery (RCA), where excessive vessel movement is commonly observed.Methods and Results:Of 647 consecutive lesions in the RCA treated with second-generation or later DESs, 426 with follow-up angiography were included in this study. Intravascular imaging analysis was performed for 584 stent edges and reference segments. Binary restenosis occurred in 42 lesions (9.9%), and 55% were SERs. The hinge angle was significantly larger in the SER group than in the other restenosis or the no-restenosis group (17.9° vs. 11.6° and 10.6°, respectively; P<0.001). Lesions with an excessive hinge angle (>11.5°) had an increased rate of target lesion revascularization (19.1% vs. 7.2%; P<0.001) during the median follow-up period of 1,578 days. In per-edge analysis, hinge angle and residual plaque burden were independent predictors of SER. The coexistence of excessive hinge motion and residual plaque burden had a synergistic effect on stenotic progression in quantitative angiographic analysis (P_interaction<0.001) at follow-up angiography.

CONCLUSIONS

Substantial stress determined by angulation at a stent edge and its interaction with residual plaque can be considered as one plausible mechanism for SER.

Evaluation of Motion Compensation Methods for Noninvasive Cardiac Radioablation of Ventricular Tachycardia

PURPOSE

Noninvasive cardiac radioablation is increasingly used for treatment of refractory ventricular tachycardia. Attempts to limit normal tissue exposure are important, including managing motion of the target. An interplay between cardiac and respiratory motion exists for cardiac radioablation, which has not been studied in depth. The objectives of this study were to estimate target motion during abdominal compression free breathing (ACFB) and respiratory gated (RG) deliveries and to investigate the quality of either implanted cardioverter defibrillator lead tip or the diaphragm as a gating surrogate.

METHODS AND MATERIALS

Eleven patients underwent computed tomography (CT) simulation with an ACFB 4-dimensional CT (r4DCT) and an exhale breath-hold cardiac 4D-CT (c4DCT). The target, implanted cardioverter defibrillator lead tip and diaphragm trajectories were measured for each patient on the r4DCT and c4DCT using rigid registration of each 4D phase to the reference (0%) phase. Motion ranges for ACFB and exhale (40%-60%) RG delivery were estimated from the target trajectories. Surrogate quality was estimated as the correlation with the target motion magnitudes.

RESULTS

Mean (range) target motion across patients from r4DCT was as follows: left/right (LR), 3.9 (1.7-6.9); anteroposterior (AP), 4.1 (2.2-5.4); and superoinferior (SI), 4.7 (2.2-7.9) mm. Mean (range) target motion from c4DCT was as follows: LR, 3.4 (1.0-4.8); AP, 4.3 (2.6-6.5); and SI, 4.1 (1.4-8.0) mm. For an ACFB, treatment required mean (range) margins to be 4.5 (3.1-6.9) LR, 4.8 (3-6.5) AP, and 5.5 (2.3-8.0) mm SI. For RG, mean (range) internal target volume motion would be 3.6 (1.1-4.8) mm LR, 4.3 (2.6-6.5) mm AP, and 4.2 (2.2-8.0) mm SI. The motion correlations between the surrogates and target showed a high level of interpatient variability.

CONCLUSIONS

In ACFB patients, a simulated exhale-gated approach did not lead to large projected improvements in margin reduction. Furthermore, the variable correlation between readily available gating surrogates could mitigate any potential advantage to gating and should be evaluated on a patient-specific basis.

Dynamic cardiac PET motion correction using 3D normalized gradient fields in patients and phantom simulations

This work expands on the implementation of three-dimensional (3D) normalized gradient fields to correct for whole-body motion and cardiac creep in [N-13]-ammonia patient studies and evaluates its accuracy using a dynamic phantom simulation model.

METHODS

A full rigid-body algorithm was developed using 3D normalized gradient fields including a multi-resolution step and sampling off the voxel grid to reduce interpolation artifacts. Optimization was performed using a weighted similarity metric that accounts for opposing gradients between images of blood pool and perfused tissue without the need for segmentation. Forty-three retrospective dynamic [N-13]-ammonia PET/CT rest/adenosine-stress patient studies were motion corrected and the mean motion parameters plotted at each frame time point. Motion correction accuracy was assessed using a comprehensive dynamic XCAT simulation incorporating published physiologic parameters of the heart's trajectory following adenosine infusion as well as corrupted attenuation correction commonly observed in clinical studies. Accuracy of the algorithm was assessed objectively by comparing the errors between isosurfaces and centers of mass of the motion corrected XCAT simulations.

RESULTS

In the patient studies, the overall mean cranial-to-caudal translation was 7 mm at stress over the duration of the adenosine infusion. Noninvasive clinical measures of relative flow reserve and myocardial flow reserve were highly correlated with their invasive analogues. Motion correction accuracy assessed with the XCAT simulations showed an error of <1 mm in late perfusion frames that broadened gradually to <3 mm in earlier frames containing blood pool.

CONCLUSION

This work demonstrates that patients undergoing [N-13]-ammonia dynamic PET/CT exhibit a large cranial-to-caudal translation related to cardiac creep primarily at stress and to a lesser extent at rest, which can be accurately corrected by optimizing their 3D normalized gradient fields. Our approach provides a solution to the challenging condition where the image intensity and its gradients are opposed without the need for segmentation and remains robust in the presence of PET-CT mismatch.

Angiography-Based 4-Dimensional Superficial Wall Strain and Stress: A New Diagnostic Tool in the Catheterization Laboratory

A novel method for four-dimensional superficial wall strain and stress (4D-SWS) is derived from the arterial motion as pictured by invasive coronary angiography. Compared with the conventional finite element analysis of cardiovascular biomechanics using the estimated pulsatile pressure, the 4D-SWS approach can calculate the dynamic mechanical state of the superficial wall in vivo, which could be directly linked with plaque rupture or stent fracture. The validation of this approach using in silico models showed that the distribution and maximum values of superficial wall stress were similar to those calculated by conventional finite element analysis. The in vivo deformation was validated on 16 coronary arteries, from the comparison of centerlines predicted by the 4D-SWS approach against the actual centerlines reconstructed from angiograms at a randomly selected time-point, which demonstrated a good agreement of the centerline morphology between both approaches (scaling: 0.995 ± 0.018 and dissimilarity: 0.007 ± 0.014). The in silico vessel models with softer plaque and larger plaque burden presented more variation in mean lumen diameter and resulted in higher superficial wall stress. In more than half of the patients (n = 16), the maximum superficial wall stress was found at the proximal lesion shoulder. Additionally, in three patients who later suffered from acute coronary syndrome, the culprit plaque rupture sites co-localized with the site of highest superficial wall stress on their baseline angiography. These representative cases suggest that angiography-based superficial wall dynamics have the potential to identify coronary segments at high-risk of plaque rupture and fracture sites of implanted stents. Ongoing studies are focusing on identifying weak spots in coronary bypass grafts, and on exploring the biomechanical mechanisms of coronary arterial remodeling and aneurysm formation. Future developments involve integration of fast computational techniques to allow online availability of superficial wall strain and stress in the catheterization laboratory.

Task-dependent estimability index to assess the quality of cardiac computed tomography angiography for quantifying coronary stenosis

Purpose: Quantifying stenosis in cardiac computed tomography angiography (CTA) images remains a difficult task, as image noise and cardiac motion can degrade image quality and distort underlying anatomic information. The purpose of this study was to develop a computational framework to objectively assess the precision of quantifying coronary stenosis in cardiac CTA.

Approach: The framework used models of coronary vessels and plaques, asymmetric motion point spread functions, CT image blur (task-based modulation transfer functions) and noise (noise-power spectrums), and an automated maximum-likelihood estimator implemented as a matched template squared-difference operator. These factors were integrated into an estimability index (e′) as a task-based measure of image quality in cardiac CTA. The e′ index was applied to assess how well it can to predict the quality of 132 clinical cases selected from the Prospective Multicenter Imaging Study for Evaluation of Chest Pain trial. The cases were divided into two cohorts, high quality and low quality, based on clinical scores and the concordance of clinical evaluations of cases by experienced cardiac imagers. The framework was also used to ascertain protocol factors for CTA Biomarker initiative of the Quantitative Imaging Biomarker Alliance (QIBA).

Results: The e′ index categorized the patient datasets with an area under the curve of 0.985, an accuracy of 0.977, and an optimal e′ threshold of 25.58 corresponding to a stenosis estimation precision (standard deviation) of 3.91%. Data resampling and training–test validation methods demonstrated stable classifier thresholds and receiver operating curve performance. The framework was successfully applicable to the QIBA objective.

Conclusions: A computational framework to objectively quantify stenosis estimation task performance was successfully implemented and was reflective of clinical results in the context of a prominent clinical trial with diverse sites, readers, scanners, acquisition protocols, and patients. It also demonstrated the potential for prospective optimization of imaging protocols toward targeted precision and measurement consistency in cardiac CT images.

Non-Rigid Respiratory Motion Estimation of Whole-Heart Coronary MR Images Using Unsupervised Deep Learning

Non-rigid motion-corrected reconstruction has been proposed to account for the complex motion of the heart in free-breathing 3D coronary magnetic resonance angiography (CMRA). This reconstruction framework requires efficient and accurate estimation of non-rigid motion fields from undersampled images at different respiratory positions (or bins). However, state-of-the-art registration methods can be time-consuming. This article presents a novel unsupervised deep learning-based strategy for fast estimation of inter-bin 3D non-rigid respiratory motion fields for motion-corrected free-breathing CMRA. The proposed 3D respiratory motion estimation network (RespME-net) is trained as a deep encoder-decoder network, taking pairs of 3D image patches extracted from CMRA volumes as input and outputting the motion field between image patches. Using image warping by the estimated motion field, a loss function that imposes image similarity and motion smoothness is adopted to enable training without ground truth motion field. RespME-net is trained patch-wise to circumvent the challenges of training a 3D network volume-wise which requires large amounts of GPU memory and 3D datasets. We perform 5-fold cross-validation with 45 CMRA datasets and demonstrate that RespME-net can predict 3D non-rigid motion fields with subpixel accuracy (0.44 ± 0.38 mm) within ~10 seconds, being ~20 times faster than a GPU-implemented state-of-the-art non-rigid registration method. Moreover, we perform non-rigid motion-compensated CMRA reconstruction for 9 additional patients. The proposed RespME-net has achieved similar motion-corrected CMRA image quality to the conventional registration method regarding coronary artery length and sharpness.

Technical Note: Cardiac synchronized volumetric modulated arc therapy for stereotactic arrhythmia radioablation - Proof of principle

PURPOSE

Ventricular tachycardia (VT) is a rapid, abnormal heart rhythm that can lead to sudden cardiac death. Current treatment options include antiarrhythmic drug therapy and catheter ablation, both of which have only modest efficacy and have potential complications. Cardiac radiosurgery has the potential to be a noninvasive and efficient treatment option for VT. Cardiac motion, however, must be accounted for to ensure accurate dose delivery to the target region. Cardiac synchronized volumetric modulated arc therapy (CSVMAT) aims to minimize the dose delivered to normal tissues by synchronizing beam delivery with a cardiac signal, irradiating only during the quiescent intervals of the cardiac cycle (when heart motion is minimal) and adjusting the beam delivery speed in response to heart rate changes.

METHODS

A CSVMAT plan was adapted from a conventional VMAT plan and delivered on a Varian TrueBeam linear accelerator. The original VMAT plan was divided into three interleaved CSVMAT phases, each consisting of alternating beam-on and beam-off segments synchronized to a sample heart rate. Trajectory log files were collected for the original VMAT and CSVMAT deliveries and the dose distributions were measured with Gafchromic EBT-XD film.

RESULTS

Analysis of the trajectory log files showed successful synchronization with the sample cardiac signal. Film analysis comparing the original VMAT and CSVMAT dose distributions returned a gamma passing rate of 99.14% (2%/2 mm tolerance).

CONCLUSIONS

The film results indicated excellent agreement between the dose distributions of the original and cardiac synchronized beam deliveries. This study demonstrates a proof of principle cardiac synchronization strategy for precise radiation treatment plan delivery and adjustment to a variable heart rate. The cardiac synchronized technique may be advantageous in radioablation for VT.

Combined T2 -preparation and multidimensional outer volume suppression for coronary artery imaging with 3D cones trajectories

PURPOSE

To develop a modular magnetization preparation sequence for combined T2 -preparation and multidimensional outer volume suppression (OVS) for coronary artery imaging.

METHODS

A combined T2 -prepared 1D OVS sequence with fat saturation was defined to contain a 90°-60 180°60 composite nonselective tip-down pulse, two 180°Y hard pulses for refocusing, and a -90° spectral-spatial sinc tip-up pulse. For 2D OVS, 2 modules were concatenated, selective in X and then Y. Bloch simulations predicted robustness of the sequence to B0 and B1 inhomogeneities. The proposed sequence was compared with a T2 -prepared 2D OVS sequence proposed by Luo et al, which uses a spatially selective 2D spiral tip-up. The 2 sequences were compared in phantom studies and in vivo coronary artery imaging studies with a 3D cones trajectory.

RESULTS

Phantom results demonstrated superior OVS for the proposed sequence compared with the Luo sequence. In studies on 15 healthy volunteers, the proposed sequence had superior image edge profile acutance values compared with the Luo sequence for the right (P < .05) and left (P < .05) coronary arteries, suggesting superior vessel sharpness. The proposed sequence also had superior signal-to-noise ratio (P < .05) and passband-to-stopband ratio (P < .05). Reader scores and reader preference indicated superior coronary image quality of the proposed sequence for both the right (P < .05) and left (P < .05) coronary arteries.

CONCLUSION

The proposed sequence with concatenated 1D spatially selective tip-ups and integrated fat saturation has superior image quality and suppression compared with the Luo sequence with 2D spatially selective tip-up.

Deep neuromuscular block improves angiographic image quality during endovascular coiling of unruptured cerebral aneurysm: a randomized clinical trial

BACKGROUND

Neuromuscular block (NMB) used during general anesthesia induces transient skeletal muscle paralysis, but patient movements during endovascular coiling still occur to some degree. Compared with moderate NMB, deep NMB may further improve the intervention condition during endovascular coiling for unruptured cerebral aneurysms; however, little research has focused on the angiographic image quality.

METHODS

This prospective, randomized, double-blind clinical trial included 58 patients treated for unruptured cerebral aneurysms with endovascular coiling under general anesthesia. Patients were randomly allocated to either the deep NMB group (post-tetanic count 1 or 2) or the moderate NMB group (train-of-four 1 or 2). The primary outcome was the proportion of patients with a satisfactory intervention condition assessed by surgeons after the procedure using a 5-point intervention condition rating scale (ICRS) from 1 (unable to obtain image) to 5 (optimal); ICRS 5 was defined as satisfactory.

RESULTS

There were significantly more cases of satisfactory intervention condition in the deep NMB group than in the moderate NMB group (82.1% vs 51.7%, p=0.015). The frequency of each ICRS score was significantly different between the groups (ICRS 5/4/3/2/1: 23/5/0/0/0 vs 15/9/2/3/0, p=0.035). The incidence of major patient movement requiring rescue muscle relaxant was 10.3% in the moderate NMB group and 0% in the deep NMB group (p=0.237). The drugs used to maintain hemodynamic stability were not significantly different between the two groups.

CONCLUSIONS

Deep NMB improves the intervention condition during endovascular coiling by improving the image quality.

Optimization of reconstruction and quantification of motion-corrected coronary PET-CT

BACKGROUND

Coronary PET shows promise in the detection of high-risk atherosclerosis, but there remains a need to optimize imaging and reconstruction techniques. We investigated the impact of reconstruction parameters and cardiac motion-correction in 18F Sodium Fluoride (18F-NaF) PET.

METHODS

Twenty-two patients underwent 18F-NaF PET within 22 days of an acute coronary syndrome. Optimal reconstruction parameters were determined in a subgroup of six patients. Motion-correction was performed on ECG-gated data of all patients with optimal reconstruction. Tracer uptake was quantified in culprit and reference lesions by computing signal-to-noise ratio (SNR) in diastolic, summed, and motion-corrected images.

RESULTS

Reconstruction using 24 subsets, 4 iterations, point-spread-function modelling, time of flight, and 5-mm post-filtering provided the highest median SNR (31.5) compared to 4 iterations 0-mm (22.5), 8 iterations 0-mm (21.1), and 8 iterations 5-mm (25.6; all P < .05). Motion-correction improved SNR of culprit lesions (n = 33) (24.5[19.9-31.5]) compared to diastolic (15.7[12.4-18.1]; P < .001) and summed data (22.1[18.9-29.2]; P < .001). Motion-correction increased the SNR difference between culprit and reference lesions (10.9[6.3-12.6]) compared to diastolic (6.2[3.6-10.3]; P = .001) and summed data (7.1 [4.8-11.6]; P = .001).

CONCLUSIONS

The number of iterations and extent of post-filtering has marked effects on coronary 18F-NaF PET quantification. Cardiac motion-correction improves discrimination between culprit and reference lesions.

Point-Cloud Method for Automated 3D Coronary Tree Reconstruction From Multiple Non-Simultaneous Angiographic Projections

X-ray angiography is the most commonly used imaging modality for the detection of coronary stenoses due to its high spatial and temporal resolution of lumen contour and its utility to guide coronary interventions in real time. However, the high inter- and intra-observer variability in interpreting the geometry of 3D vascular structure based on multiple 2D image projections is a limitation in the accurate determination of lesion severity. This could be addressed by the 3D reconstruction of the coronary arterial (CA) tree. The automated reconstruction of 3D CA tree from 2D projections is challenging due to the existence of several imaging artifacts, such as vessel overlap, foreshortening, and most importantly respiratory and cardiac motion. Along with these artifacts, the acquisition geometry introduces the possibility of generating false vessel segments in the reconstruction. Our approach aims to reduce the motion artifacts in angiographic projections by developing a new method for rigid and non-rigid motion correction. A novel point-cloud based approach is subsequently introduced for reconstruction of 3D vessel centerlines by iteratively minimizing the reconstruction error. The performance of the proposed 3D reconstruction is evaluated using angiographic projections from 45 patients, producing average reprojection errors of 0.092 ±0.055 mm and 0.910 ±0.352 mm for 3D centerlines reconstruction, when co-registered with the parent vessels on projection planes that were/were not used to derive the 3D reconstruction, respectively. A comparison of the reconstructed 3D lumen surface with optical coherence tomography (OCT) measurements has been performed, showing no statistically significant difference in the luminal cross-sections reconstructed with our method, compared to OCT.

Vulnerable plaque imaging using 18F-sodium fluoride positron emission tomography

Positron emission tomography (PET) with ¹⁸F-sodium fluoride (¹⁸F-NaF) has emerged as a promising non-invasive imaging modality to identify high-risk and ruptured atherosclerotic plaques. By visualizing microcalcification, ¹⁸F-NaF PET holds clinical promise in refining how we evaluate coronary artery disease, shifting our focus from assessing disease burden to atherosclerosis activity. In this review, we provide an overview of studies that have utilized ¹⁸F-NaF PET for imaging atherosclerosis. We discuss the associations between traditional coronary artery disease measures (risk factors) and ¹⁸F-NaF plaque activity. We also present the data on the histological validation as well as show how ¹⁸F-NaF uptake is associated with plaque morphology on intravascular and CT imaging. Finally, we discuss the technical challenges associated with ¹⁸F-NaF coronary PET highlighting recent advances in this area.

Precise measurement of coronary stenosis diameter with CCTA using CT number calibration

PURPOSE

Coronary x-ray computed tomography angiography (CCTA) continues to develop as a noninvasive method for the assessment of coronary vessel geometry and the identification of physiologically significant lesions. The uncertainty of quantitative lesion diameter measurement due to limited spatial resolution and vessel motion reduces the accuracy of CCTA diagnoses. In this paper, we introduce a new technique called computed tomography (CT)-number-Calibrated Diameter to improve the accuracy of the vessel and stenosis diameter measurements with CCTA.

METHODS

A calibration phantom containing cylindrical holes (diameters spanning from 0.8 mm through 4.0 mm) capturing the range of diameters found in human coronary vessels was three-dimensional printed. We also printed a human stenosis phantom with 17 tubular channels having the geometry of lesions derived from patient data. We acquired CT scans of the two phantoms with seven different imaging protocols. Calibration curves relating vessel intraluminal maximum voxel value (maximum CT number of a voxel, described in Hounsfield Units, HU) to true diameter, and full-width-at-half maximum (FWHM) to true diameter were constructed for each CCTA protocol. In addition, we acquired scans with a small constant motion (15 mm/s) and used a motion correction reconstruction (Snapshot Freeze) algorithm to correct motion artifacts. We applied our technique to measure the lesion diameter in the 17 lesions in the stenosis phantom and compared the performance of CT-number-Calibrated Diameter to the ground truth diameter and a FWHM estimate.

RESULTS

In all cases, vessel intraluminal maximum voxel value vs diameter was found to have a simple functional form based on the two-dimensional point spread function yielding a constant maximum voxel value region above a cutoff diameter, and a decreasing maximum voxel value vs decreasing diameter below a cutoff diameter. After normalization, focal spot size and reconstruction kernel were the principal determinants of cutoff diameter and the rate of maximum voxel value reduction vs decreasing diameter. The small constant motion had a significant effect on the CT number calibration; however, the motion-correction algorithm returned the maximum voxel value vs diameter curve to that of stationary vessels. The CT number Calibration technique showed better performance than FWHM estimation of diameter, yielding a high accuracy in the tested range (0.8 mm through 2.5 mm). We found a strong linear correlation between the smallest diameter in each of 17 lesions measured by CT-number-Calibrated Diameter (DC ) and ground truth diameter (Dgt ), (DC  = 0.951 × Dgt  + 0.023 mm, r = 0.998 with a slope very close to 1.0 and intercept very close to 0 mm.

CONCLUSIONS

Computed tomography-number-Calibrated Diameter is an effective method to enhance the accuracy of the estimate of small vessel diameters and degree of coronary stenosis in CCTA.

Restricting motion effects in CT coronary angiography

OBJECTIVE

Evaluation of coronary CT image blur using multi segment reconstruction algorithm.

METHODS

Cardiac motion was simulated in a Catphan. CT coronary angiography was performed using 320 × 0.5 mm detector array and 275 ms gantry rotation. 1, 2 and 3 segment reconstruction algorithm, three heart rates (60, 80 and 100bpm), two peak displacements (4, 8 mm) and three cardiac phases (55, 35, 75%) were used. Wilcoxon test compared image blur from the different reconstruction algorithms.

RESULTS

Image blur for 1, 2 and 3 segments in: 60 bpm, 75% R-R interval and 8 mm peak displacement: 0.714, 0.588, 0.571 mm (1.18, 0.6, 0.4 mm displacement) 80 bpm, 35% R-R interval and 8 mm peak displacement: 0.869, 0.606, 0.606 mm (1.57, 0.79,0.52 mm displacement) 100 bpm, 35% R-R interval and 4 mm peak displacement: 0.645, 0.588, 0.571 mm (0.98, 0.49, 0.33 mm displacement). The median image blur overall for 1 and 2 segments was 0.714 mm and 0.588 mm respectively (p < 0.0001).

CONCLUSION

Two-segment reconstruction significantly reduces image blur.

ADVANCES IN KNOWLEDGE

Multisegment reconstruction algorithms during CT coronary angiography are a useful method to reduce image blur, improve visualization of the coronary artery wall and help the early detection of the plaque.

Influence of virtual monochromatic spectral image at different energy levels on motion artifact correction in dual-energy spectral coronary CT angiography

PURPOSE

To investigate the influence of virtual monochromatic spectral (VMS) CT images at different energy levels on the effectiveness of a motion correction technique (SSF) in dual-energy Spectral coronary CT angiography (CCTA).

MATERIALS AND METHODS

29 cases suspected of or diagnosed with coronary artery disease underwent Spectral CCTA using a prospective ECG triggering with 250 ms padding time. SSF was applied to the determined least-motion phase to generate 6 additional sets of VMS images with energy levels from 40 to 100 keV. CT value and standard deviation (SD) in the aortic root and epicardial adipose tissue were measured. Image quality of the RCA, LAD and LCX was evaluated on a per-vessel basis in each patient. Two reviewers evaluated the artery using the score of the segment.

RESULTS

The low energy VMS images increased CT value and image noise compared with higher-energy VMS images, except 90 keV and 100 keV. The CNR of 40-70 keV were higher than those of 80-100 keV (P < 0.05). The image quality scores for images at 50-80 keV were higher than those of 40, 90, and 100 keV (P < 0.05), and the VMS image quality at 50 keV and 60 keV with SSF was the highest.

CONCLUSION

SSF can effectively reduce the motion artifacts when coronary vessels have suitable contrast enhancement which can be achieved by adjusting energy levels of VMS images.

Deep feature descriptor based hierarchical dense matching for X-ray angiographic images

Backgroud and Objective: X-ray angiography, a powerful technique for blood vessel visualization, is widely used for interventional diagnosis of coronary artery disease because of its fast imaging speed and perspective inspection ability. Matching feature points in angiographic images is a considerably challenging task due to repetitive weak-textured regions.

METHODS

In this paper, we propose an angiographic image matching method based on the hierarchical dense matching framework, where a novel deep feature descriptor is designed to compute multilevel correlation maps. In particular, the deep feature descriptor is computed by a deep learning model specifically designed and trained for angiographic images, thereby making the correlation maps more distinctive for corresponding feature points in different angiographic images. Moreover, point correspondences are further hierarchically extracted from multilevel correlation maps with the highest similarity response(s), which is relatively robust and accurate. To overcome the problem regarding the lack of training samples, the convolutional neural network (designed for deep feature descriptor) is initially trained on samples from natural images and then fine-tuned on manually annotated angiographic images. Finally, a dense matching completion method, based on the distance between deep feature descriptors, is proposed to generate dense matches between images.

RESULTS

The proposed method has been evaluated on the number and accuracy of extracted matches and the performance of subtraction images. Experiments on a variety of angiographic images show promising matching accuracy, compared with state-of-the-art methods.

CONCLUSIONS

The proposed angiographic image matching method is shown to be accurate and effective for feature matching in angiographic images, and further achieves good performance in image subtraction.

Gating Approaches in Cardiac PET Imaging

Cardiac PET provides high sensitivity and high negative predictive value in the diagnosis of coronary artery disease and cardiomyopathies. Cardiac, respiratory as well as bulk patient motion have detrimental effects on thoracic PET imaging, in particular on cardiovascular PET imaging where the motion can affect the PET images quantitatively as well as qualitatively. Gating can ameliorate the unfavorable impact of motion additionally enabling evaluation of left ventricular systolic function. In this article, the authors review the recent advances in gating approaches and highlight the advances in data-driven approaches, which hold promise in motion detection without the need for complex hardware setup.

Intrafractional Displacement of Cardiac Substructures Among Patients With Mediastinal Lymphoma or Lung Cancer

Purpose

The radiation dose to specific substructures of the heart may be more critical than the dose to the whole heart. Yet, these substructures are sensitive to intrafractional motion from breathing and cardiac motion, which can affect their dose-volume histograms. We sought to investigate intrafractional motion of the heart and its substructures among free-breathing patients undergoing radiation for mediastinal lymphoma or lung cancer.

Methods and materials

After institutional review board approval, the medical records of 20 patients (12 with mediastinal lymphoma; 8 with lung cancer) were retrospectively reviewed. Patients underwent 4-dimensional computed tomography simulation and a contrasted scan for treatment planning. Using MIMVista software, the heart, coronary arteries, chambers, and valves were contoured on the 50% phase, and these contours were propagated to the other phases and edited. Each substructure was graded on the basis of its ease of contouring across all phases (1 = no difficulty; 2 = minor difficulty; 3 = moderate difficulty; and 4 = very difficult). The centroid position and volume of each substructure for all phases were exported to Excel to calculate basic statistics and the independent t test.

Results

The heart, 4 chambers, and atrioventricular valves were easily identified with a mean score of 1 to 1.2, and the pulmonic valve, left anterior descending artery, aortic valve, circumflex, and right coronary artery were minor-to-moderately difficult with a mean score of 2.1 to 3.2. The smallest centroid displacement was seen in the 4 chambers and mitral and pulmonic valves (0.7-1.1 cm). Greater displacement was seen in the coronary vessels and tricuspid and aortic valves (1.2-1.5 cm). The greatest displacement was in the Z direction (craniocaudal) for all substructures; however, the displacement was significantly greater among patients with lymphoma for the right ventricle, aortic valve, and left anterior descending artery (P < .05). However, patients with lung cancer had more displacement in the X and Y directions, which was statistically significant for the right atrium, tricuspid valve, right ventricle, and heart. When calculating overall displacement, no statistically significant difference was observed between patients with lymphoma and patients with lung cancer.

Conclusions

Intrafractional motion of the cardiac substructures ranged from 0.7 to 1.5 cm, mostly in the Z direction. Further investigation of the respiratory motion effect on the dose-volume histogram of the substructures is needed for patients treated with contemporary radiation techniques.

Nearest Neighbor Method to Estimate Internal Target for Real-Time Tumor Tracking

PURPOSE

This work proposed a nearest neighbor estimation method to track the respiration-induced tumor motion.

METHODS

Based on the simultaneously collected motion traces of external surrogate and internal target during the modeling phase prior to treatment, we first obtain the nearest neighbors of the current surrogate in external space. Subsequently, the concurrent targets in internal space are determined and used to estimate the current target position. The method was validated on 71 cases that were from 3 open access databases. In addition, to evaluate the method's estimation and prediction accuracy, we compared the method with other works.

RESULTS

Except for 2 cases, the nearest neighbor estimation achieved the root-mean-square error of <3 mm. The comparison indicated that the method had better estimation accuracy than polynomial model and good prediction performance.

DISCUSSION

The 2 exceptive cases were further analyzed for failure causes. We inferred that one was because of the lack of estimating new target in our method, and the other one was because of the mistake during data collection. Accordingly, the potential solutions were suggested. Besides, the method's estimation for surrogate outliers, effects of modeling length, calibration, and extension were discussed.

CONCLUSION

The results demonstrated nearest neighbor estimation's effectiveness. Except for this, the method imposes no restrictions on the modality of the pretreatment target images and does not assume a specific correspondence function between the surrogate and the target. With only 1 critical parameter, this nearest neighbor estimation method is easy to implement in clinical setting and thus has potential for broad applications.

Prior-Free Respiratory Motion Estimation in Rotational Angiography

Rotational coronary angiography using C-arm angiography systems enables intra-procedural 3-D imaging that is considered beneficial for diagnostic assessment and interventional guidance. Despite previous efforts, rotational angiography was not yet successfully established in clinical practice for coronary artery procedures due to challenges associated with substantial intra-scan respiratory and cardiac motion. While gating handles cardiac motion during reconstruction, respiratory motion requires compensation. State-of-the-art algorithms rely on 3-D / 2-D registration that requires an uncompensated reconstruction of sufficient quality. To overcome this limitation, we investigate two prior-free respiratory motion estimation methods based on the optimization of: 1) epipolar consistency conditions (ECCs) and 2) a task-based auto-focus measure (AFM). The methods assess redundancies in projection images or impose favorable properties of 3-D space, respectively, and are used to estimate the respiratory motion of the coronary arteries within rotational angiograms. We evaluate our algorithms on the publicly available CAVAREV benchmark and on clinical data. We quantify reductions in error due to respiratory motion compensation using a dedicated reconstruction domain metric. Moreover, we study the improvements in image quality when using an analytic and a novel temporal total variation regularized algebraic reconstruction algorithm. We observed substantial improvement in all figures of merit compared with the uncompensated case. Improvements in image quality presented as a reduction of double edges, blurring, and noise. Benefits of the proposed corrections were notable even in cases suffering little corruption from respiratory motion, translating to an improvement in the vessel sharpness of (6.08 ± 4.46)% and (14.7 ± 8.80)% when the ECC-based and the AFM-based compensation were applied. On the CAVAREV data, our motion compensation approach exhibits an improvement of (27.6 ± 7.5)% and (97.0 ± 17.7)% when the ECC and AFM were used, respectively. At the time of writing, our method based on AFM is leading the CAVAREV scoreboard. Both motion estimation strategies are purely image-based and accurately estimate the displacements of the coronary arteries due to respiration. While current evidence suggests the superior performance of AFM, future work will further investigate the use of ECC in the context of angiography as they solely rely on geometric calibration and projection-domain images.

Five-minute whole-heart coronary MRA with sub-millimeter isotropic resolution, 100% respiratory scan efficiency, and 3D-PROST reconstruction

Purpose

To enable whole‐heart 3D coronary magnetic resonance angiography (CMRA) with isotropic sub‐millimeter resolution in a clinically feasible scan time by combining respiratory motion correction with highly accelerated variable density sampling in concert with a novel 3D patch‐based undersampled reconstruction (3D‐PROST).

Methods

An undersampled variable density spiral‐like Cartesian trajectory was combined with 2D image‐based navigators to achieve 100% respiratory efficiency and predictable scan time. 3D‐PROST reconstruction integrates structural information from 3D patch neighborhoods through sparse representation, thereby exploiting the redundancy of the 3D anatomy of the coronary arteries in an efficient low‐rank formulation. The proposed framework was evaluated in a static resolution phantom and in 10 healthy subjects with isotropic resolutions of 1.2 mm³ and 0.9 mm³ and undersampling factors of ×5 and ×9. 3D‐PROST was compared against fully sampled (1.2 mm³ only), conventional parallel imaging, and compressed sensing reconstructions.

Results

Phantom and in vivo (1.2 mm³) reconstructions were in excellent agreement with the reference fully sampled image. In vivo average acquisition times (min:s) were 7:57 ± 1:18 (×5) and 4:35 ± 0:44 (×9) for 0.9 mm³ resolution. Sub‐millimeter 3D‐PROST resulted in excellent depiction of the left and right coronary arteries including small branch vessels, leading to further improvements in vessel sharpness and visible vessel length in comparison with conventional reconstruction techniques. Image quality rated by 2 experts demonstrated that 3D‐PROST provides good image quality and is robust even at high acceleration factors.

Conclusion

The proposed approach enables free‐breathing whole‐heart 3D CMRA with isotropic sub‐millimeter resolution in <5 min and achieves improved coronary artery visualization in a short and predictable scan time.

Miniaturized Robotic End-Effector with Piezoelectric Actuation and Fiber Optic Sensing for Minimally Invasive Cardiac Procedures

Each year 35,000 cardiac ablation procedures are performed to treat atrial fibrillation through the use of catheter systems. The success rate of this treatment is highly dependent on the force which the catheter applies on the heart wall. If the magnitude of the applied force is much higher than a certain threshold the tissue perforates, whereas if the force is lower than this threshold the lesion size may be too large and is inconsistent. Furthermore, studies have shown large variability in the applied force from trained physicians during treatment, suggesting that physicians are unable to manually regulate the levels of the force at the site of treatment. Current catheter systems do not provide the physicians with active means for contact force control and are only at most aided by visual feedback of the forces measured in situ. This paper discusses a novel design of a robotic end-effector that integrates mechanisms of sensing and actively controlling of the applied forces into a miniaturized compact form. The required specifications for design and integration were derived from the current application under investigation. An off-the-shelf miniature piezoelectric motor was chosen for actuation, and a force sensing solution was developed to meet the specifications. Experimental characterization of the actuator and the force sensor within the integrated setup show compliance with the specifications and pave the way for future experimentation where closed-loop control of the system can be implemented according to the contact force control strategies for the application.

Definition of the margin of major coronary artery bifurcations during radiotherapy with electrocardiograph-gated 4D-CT

PURPOSE

The aim was to measure the cardiac motion-induced displacements of major coronary artery bifurcations utilizing electrocardiography (ECG)-gated four-dimensional computed tomography (4D-CT) and to determine the margin of coronary artery bifurcations.

METHODS

Thirty-seven female patients who underwent retrospective ECG-gated 4D-CT in inspiratory breath hold (IBH) were enrolled. The left main coronary artery bifurcation (LM), the obtuse marginal branch bifurcation (OM), the first diagonal branch bifurcation (D₁), the second diagonal branch bifurcation (D₂), the caudal portion of the left anterior descending branch (APX), the first right ventricular artery bifurcation (V) and the acute marginal branch bifurcation (AM) were contoured. The center of the contour of the coronary arterial bifurcations at end systole was defined as the standard, and the margin were then calculated.

RESULTS

The margin in the left-right (LR), cranio-caudal (CC), and anterior-posterior (AP) coordinates were as follows: LM 3, 3, and 3 mm; D₁ 6, 3, and 3 mm; D₂ 3, 3, and 3 mm; APX 4, 4, and 4 mm; OM 4, 6, and 5 mm; V 6, 8, and 7 mm; and AM 6, 8, and 7 mm, respectively.

CONCLUSION

Coronary artery bifurcations should be considered a separate organ at risk (OAR), and different margin should be provided due to the differences resulting from motion displacement. The maximum margin in the LR, CC, and AP coordinates of left coronary artery bifurcations were 6, 6, and 5 mm, and those of the right coronary artery bifurcations were 6, 8, and 7 mm, respectively.

Coronary artery rotation in native and stented porcine coronary arteries

INTRODUCTION

Coronary arteries are exposed to several complex biomechanical forces during the cardiac cycle. These biomechanical forces potentially contribute to both native coronary artery disease, development of atherosclerosis and eventual stent failure. The aim of the present study was to characterize and define coronary artery axial rotation and the effect of stent implantation on this biomechanical factor.

METHODS

Intravascular ultrasound (IVUS) images were obtained from porcine coronary arteries and analyzed in ultrasound analysis software used to evaluate myocardial strain and torsion in echocardiography. In this study the software was utilized for a novel application to evaluate coronary artery rotation and time-to-peak (TTP) rotation in porcine coronary arteries. Clockwise (CW) and counterclockwise (CCW) rotation of coronary arteries during the cardiac cycle and (TTP) rotation were measured.

RESULTS

A total of 11 (4 LAD, 4 LCX, 3 RCA) coronary artery segments were independently analyzed pre- and post-stent implantation for a total of 22 IVUS runs. CW and CCW rotation and TTP varied widely within coronary artery segments and between different coronary arteries. Stent implantation impacted degree, direction and TTP of coronary rotation. Measurement reliability was assessed and the intraclass correlation coefficient for maximum average CCW was 0.990 (95% confidence interval 0.980-0.996, P < 0.0001), indicating excellent agreement.

CONCLUSIONS

Coronary arteries display wide spectrum of CW and CCW rotation during the cardiac cycle. Coronary stents impact the degree and direction of coronary artery rotation. The implications of these findings on development of atherosclerosis and stent failure require further investigation.

Incorporation of the Living Heart Model into the 4D XCAT Phantom for Cardiac Imaging Research

The 4D extended cardiac-torso (XCAT) phantom has provided a valuable tool to study the effects of anatomy and motion on medical images, especially cardiac motion. One limitation of the XCAT was that it did not have a physiological basis which to realistically simulate variations in cardiac function. In this work, we incorporate into the XCAT anatomy the four-chamber FE Living Heart Model (LHM) developed by the Living Heart Project (LHP). The LHM represents the state of the art in cardiac FE simulation because of its ability to accurately replicate the biomechanical motion of the entire heart and its variations. We create a new series of 4D phantoms capable of simulating patients with varying body sizes and shapes; cardiac positions, orientations, and dynamics. While extendable to other imaging modalities and technologies, our goal is to use the FE-enhanced XCAT models to investigate the optimal use of computed tomography (CT) for the evaluation of coronary artery disease (CAD). With the ability to simulate realistic, predictive, patient quality 4D imaging data, the enhanced XCAT models will enable optimization studies to identify the most promising systems or system parameters for further clinical validation.

Motion Correction and Its Impact on Absolute Myocardial Blood Flow Measures with PET

PURPOSE OF REVIEW

Motion artifacts, due to cardiac and respiratory cycles, myocardial cardiac creep, or gross patient movements, have been extensively investigated in the context of relative myocardial perfusion imaging with SPECT and PET. These movements have been identified as a major source of errors in image quantification and diagnosis. Recently, as dynamic PET quantification for myocardial blood flow assessment has entered clinical practice, similar questions have arisen on the impact of motion on final blood flow values.

RECENT FINDINGS

While preliminary investigations have underlined the potential impact of these motions on MBF quantification, their correction on dynamic acquisition remains challenging and limited to research studies. Gross patient's body movements occur in a consistent number of cases, particularly during stress acquisition, typically involving a limited number of image frames. If undetected, these movements can lead to great differences in flow values and consequently misdiagnosis. Quality control routines can be applied to automatically inspect the shape of time activity curves and to help identify motion artifacts. Cyclic cardiac and respiratory motion may have a considerable impact on final flow values. Correction of gross body motion represents a priority in the context of optimizing absolute flow clinical routine utilization and protocol standardization.

Beating-heart registration for organ-mounted robots

BACKGROUND

Organ-mounted robots address the problem of beating-heart surgery by adhering to the heart, passively providing a platform that approaches zero relative motion. Because of the quasi-periodic deformation of the heart due to heartbeat and respiration, registration must address not only spatial registration but also temporal registration.

METHODS

Motion data were collected in the porcine model in vivo (N = 6). Fourier series models of heart motion were developed. By comparing registrations generated using an iterative closest-point approach at different phases of respiration, the phase corresponding to minimum registration distance is identified.

RESULTS

The spatiotemporal registration technique presented here reduces registration error by an average of 4.2 mm over the 6 trials, in comparison with a more simplistic static registration that merely averages out the physiological motion.

CONCLUSIONS

An empirical metric for spatiotemporal registration of organ-mounted robots is defined and demonstrated using data from animal models in vivo.

Positron emission tomography imaging of coronary atherosclerosis

Inflammation has a central role in the progression of coronary atherosclerosis. Recent developments in cardiovascular imaging with the advent of hybrid positron emission tomography have provided a window into the molecular pathophysiology underlying coronary plaque inflammation. Using novel radiotracers targeted at specific cellular pathways, the potential exists to observe inflammation, apoptosis, cellular hypoxia, microcalcification and angiogenesis in vivo. Several clinical studies are now underway assessing the ability of this hybrid imaging modality to inform about atherosclerotic disease activity and the prediction of future cardiovascular risk. A better understanding of the molecular mechanisms governing coronary atherosclerosis may be the first step toward offering patients a more stratified, personalized approach to treatment.

Quantification of uncertainty in the assessment of coronary plaque in CCTA through a dynamic cardiac phantom and 3D-printed plaque model

The purpose of this study was to develop a dynamic physical cardiac phantom with a realistic coronary plaque to investigate stenosis measurement accuracy under clinically relevant heart-rates. The coronary plaque model (5 mm diameter, 50% stenosis, and 32 mm long) was designed and 3D-printed with tissue equivalent materials (calcified plaque with iodine-enhanced lumen). Realistic cardiac motion was modeled by converting computational cardiac motion vectors into compression and rotation profiles executed by a commercial base cardiac phantom. The phantom was imaged on a dual-source CT system applying a retrospective gated coronary CT angiography (CCTA) protocol using synthesized motion-synchronized electrocardiogram (ECG) waveforms. Multiplanar reformatted images were reconstructed along vessel centerlines. Enhanced lumens were segmented by five independent operators. On average, stenosis measurement accuracy was 0.9% positively biased for the motion-free condition. Average measurement accuracy monotonically decreased from 0.9% positive bias for the motion-free condition to 18.5% negative bias at 90 beats per minute. Contrast-to-noise ratio, lumen circularity, and segmentation conformity also decreased monotonically with increasing heart-rate. These results demonstrate successful implementation of a base cardiac phantom with a 3D-printed coronary plaque model, relevant motion profile, and coordinated ECG waveform. They further show the utility of the model to ascertain metrics of CCTA accuracy and image quality under realistic plaque, motion, and acquisition conditions.

Registration of angiographic image on real-time fluoroscopic image for image-guided percutaneous coronary intervention

PURPOSE

In percutaneous coronary intervention (PCI), cardiologists must study two different X-ray image sources: a fluoroscopic image and an angiogram. Manipulating a guidewire while alternately monitoring the two separate images on separate screens requires a deep understanding of the anatomy of coronary vessels and substantial training. We propose 2D/2D spatiotemporal image registration of the two images in a single image in order to provide cardiologists with enhanced visual guidance in PCI.

METHODS

The proposed 2D/2D spatiotemporal registration method uses a cross-correlation of two ECG series in each image to temporally synchronize two separate images and register an angiographic image onto the fluoroscopic image. A guidewire centerline is then extracted from the fluoroscopic image in real time, and the alignment of the centerline with vessel outlines of the chosen angiographic image is optimized using the iterative closest point algorithm for spatial registration.

RESULTS

A proof-of-concept evaluation with a phantom coronary vessel model with engineering students showed an error reduction rate greater than 74% on wrong insertion to nontarget branches compared to the non-registration method and more than 47% reduction in the task completion time in performing guidewire manipulation for very difficult tasks. Evaluation with a small number of experienced doctors shows a potentially significant reduction in both task completion time and error rate for difficult tasks. The total registration time with real procedure X-ray (angiographic and fluoroscopic) images takes [Formula: see text] 60 ms, which is within the fluoroscopic image acquisition rate of 15 Hz.

CONCLUSIONS

By providing cardiologists with better visual guidance in PCI, the proposed spatiotemporal image registration method is shown to be useful in advancing the guidewire to the coronary vessel branches, especially those difficult to insert into.

Enhancing Cardiac PET by Motion Correction Techniques

PURPOSE OF REVIEW

Cardiac positron emission tomography (PET) images often contain errors due to cardiac, respiratory, and patient motion during relatively long image acquisition. Advanced motion compensation techniques may improve PET spatial resolution, eliminate potential artifacts, and ultimately improve the research and clinical capabilities of PET.

RECENT FINDINGS

Combined cardiac and respiratory gating has only recently been implemented in clinical PET systems. Considering that the gated image bins contain much lower counts than the original PET data, they need to be summed after correcting for motion, forming motion-corrected, high-count image volume. Furthermore, automated image registration techniques can be used to correct for motion between CT attenuation scan and PET acquisition. While motion correction methods are not yet widely used in clinical practice, approaches including dual-gated non-rigid motion correction and the incorporation of motion correction information into the reconstruction process have the potential to markedly improve cardiac PET imaging.