Protocol optimization for functional cardiac CT imaging using noise emulation in the raw data domain.

BACKGROUND: Four-dimensional (4D) wide coverage computed tomography (CT) is an effective imaging modality for measuring the mechanical function of the myocardium. However, repeated CT measurement across a number of heartbeats is still a concern. PURPOSE: A projection-domain noise emulation method is presented to generate accurate low-dose (mA modulated) 4D cardiac CT scans from high-dose scans, enabling protocol optimization to deliver sufficient image quality for functional cardiac analysis while using a dose level that is as low as reasonably achievable (ALARA). METHODS: Given a targeted low-dose mA modulation curve, the proposed noise emulation method injects both quantum and electronic noise of proper magnitude and correlation to the high-dose data in projection domain. A spatially varying (i.e., channel-dependent) detector gain term as well as its calibration method were proposed to further improve the noise emulation accuracy. To determine the ALARA dose threshold, a straightforward projection domain image quality (IQ) metric was proposed that is based on the number of projection rays that do not fall under the non-linear region of the detector response. Experiments were performed to validate the noise emulation method with both phantom and clinical data in terms of visual similarity, contrast-to-noise ratio (CNR), and noise-power spectrum (NPS). RESULTS: For both phantom and clinical data, the low-dose emulated images exhibited similar noise magnitude (CNR difference within 2%), artifacts, and texture to that of the real low-dose images. The proposed channel-dependent detector gain term resulted in additional increase in emulation accuracy. Using the proposed IQ metric, recommended kVp and mA settings were calculated for low dose 4D Cardiac CT acquisitions for patients of different sizes. CONCLUSIONS: A detailed method to estimate system-dependent parameters for a raw-data based low dose emulation framework was described. The method produced realistic noise levels, artifacts, and texture with phantom and clinical studies. The proposed low-dose emulation method can be used to prospectively select patient-specific minimal-dose protocols for functional cardiac CT.

Protocol Optimization for Functional Cardiac CT Imaging Using Noise Emulation in the Raw Data Domain.

Background: Four-dimensional (4D) wide coverage computed tomography (CT) is an effective imaging modality for measuring the mechanical function of the myocardium. However, repeated CT measurement across a number of heartbeats is still a concern. Purpose: A projection-domain noise emulation method is presented to generate accurate low-dose (mA modulated) 4D cardiac CT scans from high-dose scans, enabling protocol optimization to deliver sufficient image quality for functional cardiac analysis while using a dose level that is as low as reasonably achievable (ALARA). Methods: Given a targeted low-dose mA modulation curve, the proposed noise emulation method injects both quantum and electronic noise of proper magnitude and correlation to the high-dose data in projection domain. A spatially varying (i.e., channel-dependent) detector gain term as well as its calibration method were proposed to further improve the noise emulation accuracy. To determine the ALARA dose threshold, a straightforward projection domain image quality (IQ) metric was proposed that is based on the number of projection rays that do not fall under the non-linear region of the detector response. Experiments were performed to validate the noise emulation method with both phantom and clinical data in terms of visual similarity, contrast-to-noise ratio (CNR), and noise-power spectrum (NPS). Results: For both phantom and clinical data, the low-dose emulated images exhibited similar noise magnitude (CNR difference within 2%), artifacts, and texture to that of the real low-dose images. The proposed channel-dependent detector gain term resulted in additional increase in emulation accuracy. Using the proposed IQ metric, recommended kVp and mA settings were calculated for low dose 4D Cardiac CT acquisitions for patients of different sizes. Conclusions: A detailed method to estimate system-dependent parameters for a raw-data based low dose emulation framework was described. The method produced realistic noise levels, artifacts, and texture with phantom and clinical studies. The proposed low-dose emulation method can be used to prospectively select patient-specific minimal-dose protocols for functional cardiac CT.

Motion correction of wide-detector 4DCT images for cardiac resynchronization therapy planning.

BACKGROUND: Lead placement at the latest mechanically activated left ventricle (LV) segments is strongly correlated with response to cardiac resynchronization therapy (CRT). We demonstrate the feasibility of a cardiac 4DCT motion correction algorithm (ResyncCT) in estimating LV mechanical activation for guiding lead placement in CRT. METHODS: Subjects with full cardiac cycle 4DCT images acquired using a wide-detector CT scanner for CRT planning/upgrade were included. 4DCT images exhibited motion artifact-induced false-dyssynchrony, hindering LV mechanical activation time estimation. Motion-corrupted images were processed with ResyncCT to yield motion-corrected images. Time to onset of shortening (TOS) was estimated in each of 72 endocardial segments. A false-dyssynchrony index (FDI) was used to quantify the extent of motion artifacts in the uncorrected and the ResyncCT images. After motion correction, the change in classification of LV free-wall segments as optimal target sites for lead placement was investigated. RESULTS: Twenty subjects (70.7 â€‹± â€‹13.9 years, 6 female) were analyzed. Motion artifacts in the ResyncCT-processed images were significantly reduced (FDI: 28.9 â€‹± â€‹9.3 â€‹% vs 47.0 â€‹± â€‹6.0 â€‹%, p â€‹< â€‹0.001). In 10 (50 â€‹%) subjects, ResyncCT motion correction yielded statistically different TOS estimates (p â€‹< â€‹0.05). Additionally, 43 â€‹% of LV free-wall segments were reclassified as optimal target sites for lead placement after motion correction. CONCLUSIONS: ResyncCT significantly reduced motion artifacts in wide-detector cardiac 4DCT images, yielded statistically different time to onset of shortening estimates, and changed the location of optimal target sites for lead placement. These results highlight the potential utility of ResyncCT motion correction in CRT planning when using wide-detector 4DCT imaging.

Efficient multi-fidelity computation of blood coagulation under flow.

Clot formation is a crucial process that prevents bleeding, but can lead to severe disorders when imbalanced. This process is regulated by the coagulation cascade, a biochemical network that controls the enzyme thrombin, which converts soluble fibrinogen into the fibrin fibers that constitute clots. Coagulation cascade models are typically complex and involve dozens of partial differential equations (PDEs) representing various chemical species' transport, reaction kinetics, and diffusion. Solving these PDE systems computationally is challenging, due to their large size and multi-scale nature. We propose a multi-fidelity strategy to increase the efficiency of coagulation cascade simulations. Leveraging the slower dynamics of molecular diffusion, we transform the governing PDEs into ordinary differential equations (ODEs) representing the evolution of species concentrations versus blood residence time. We then Taylor-expand the ODE solution around the zero-diffusivity limit to obtain spatiotemporal maps of species concentrations in terms of the statistical moments of residence time, [Formula: see text], and provide the governing PDEs for [Formula: see text]. This strategy replaces a high-fidelity system of N PDEs representing the coagulation cascade of N chemical species by N ODEs and p PDEs governing the residence time statistical moments. The multi-fidelity order (p) allows balancing accuracy and computational cost providing a speedup of over N/p compared to high-fidelity models. Moreover, this cost becomes independent of the number of chemical species in the large computational meshes typical of the arterial and cardiac chamber simulations. Using a coagulation network with N = 9 and an idealized aneurysm geometry with a pulsatile flow as a benchmark, we demonstrate favorable accuracy for low-order models of p = 1 and p = 2. The thrombin concentration in these models departs from the high-fidelity solution by under 20% (p = 1) and 2% (p = 2) after 20 cardiac cycles. These multi-fidelity models could enable new coagulation analyses in complex flow scenarios and extensive reaction networks. Furthermore, it could be generalized to advance our understanding of other reacting systems affected by flow.

Efficient multi-fidelity computation of blood coagulation under flow.

Clot formation is a crucial process that prevents bleeding, but can lead to severe disorders when imbalanced. This process is regulated by the coagulation cascade, a biochemical network that controls the enzyme thrombin, which converts soluble fibrinogen into the fibrin fibers that constitute clots. Coagulation cascade models are typically complex and involve dozens of partial differential equations (PDEs) representing various chemical species' transport, reaction kinetics, and diffusion. Solving these PDE systems computationally is challenging, due to their large size and multi-scale nature. We propose a multi-fidelity strategy to increase the efficiency of coagulation cascade simulations. Leveraging the slower dynamics of molecular diffusion, we transform the governing PDEs into ordinary differential equations (ODEs) representing the evolution of species concentrations versus blood residence time. We then Taylor-expand the ODE solution around the zero-diffusivity limit to obtain spatiotemporal maps of species concentrations in terms of the statistical moments of residence time, , and provide the governing PDEs for . This strategy replaces a high-fidelity system of N PDEs representing the coagulation cascade of N chemical species by N ODEs and p PDEs governing the residence time statistical moments. The multi-fidelity order( p ) allows balancing accuracy and computational cost, providing a speedup of over N/p compared to high-fidelity models. Using a simplified coagulation network and an idealized aneurysm geometry with a pulsatile flow as a benchmark, we demonstrate favorable accuracy for low-order models of p = 1 and p = 2. These models depart from the high-fidelity solution by under 16% ( p = 1) and 5% ( p = 2) after 20 cardiac cycles. The favorable accuracy and low computational cost of multi-fidelity models could enable unprecedented coagulation analyses in complex flow scenarios and extensive reaction networks. Furthermore, it can be generalized to advance our understanding of other systems biology networks affected by blood flow.

Quantitative analysis of small coronary artery calcium detectability with an accurate simulation and validation on a clinical CT scanner.

BACKGROUND: The absence of coronary artery calcium (CAC) measured via CT is associated with very favorable prognosis, and current guidelines recommend low-density lipoprotein cholesterol (LDL-c) lowering therapy for individuals with any CAC. This motivates early detection of small granules of CAC; however, calcium scan sensitivity for detecting very low levels of calcium has not been quantified. PURPOSE: In this work, the size limit of detectability of small calcium hydroxyapatite (CaHA) granules with clinical CAC scanning was assessed using validated simulations. METHODS: CT projections of digital 3D mathematical phantoms containing small CaHA granules were simulated analytically; images were reconstructed using a filter designed to reproduce the point spread function of a specific commercial scanner, and a relationship of HU number versus diameter was derived. These simulation results were validated with experimental measurements of HU versus diameter from phantoms containing small granules of CaHA on a GE Revolution CT scanner in the clinic; ground truth measurements of the CaHA granule diameters were obtained using a Zeiss Xradia 510 Versa high-resolution 3D micro-CT imaging system. Using experimental measurements on the clinical CT scanner, detectability was quantified with a detectability index (d') using a non-prewhitened matched filter. The effect of changes to reconstruction slice thickness and reconstruction kernel on granule detectability was evaluated. RESULTS: Under typical clinical calcium scanning and reconstruction conditions, the minimum detectable diameter of a simulated spherical calcium granule with a clinically relevant CaHA density was 0.76 mm. The minimum detectable volume was 2.4 times smaller on images reconstructed at a slice thickness of 0.625 mm compared to 2.5 mm. The detectability index d' increased by a factor of 1.7 when images were reconstructed with 0.625 mm slices compared to 2.5 mm slices. d' did not change when images were reconstructed with the high-resolution BONE filter compared to the less sharp STANDARD resolution filter on the GE Revolution CT. CONCLUSIONS: We have quantified detectability versus size of small calcium granules at the resolution limit of a widely available clinical CT scanner. Detectability increased significantly with reduced slice thickness and did not change with a sharper reconstruction kernel. The simulation can be used to calculate the trade-off between dose and CAC detectability.

Pulmonary vein flow split effects in patient-specific simulations of left atrial flow.

Disruptions to left atrial (LA) blood flow, such as those caused by atrial fibrillation (AF), can lead to thrombosis in the left atrial appendage (LAA) and an increased risk of systemic embolism. LA hemodynamics are influenced by various factors, including LA anatomy and function, and pulmonary vein (PV) inflow conditions. In particular, the PV flow split can vary significantly among and within patients depending on multiple factors. In this study, we investigated how changes in PV flow split affect LA flow transport, focusing for the first time on blood stasis in the LAA, using a high-fidelity patient-specific computational fluid dynamics (CFD) model. We use an Immersed Boundary Method, simulating the flow in a fixed, uniform Cartesian mesh and imposing the movement of the LA walls with a moving Lagrangian mesh generated from 4D Computerized Tomography images. We analyzed LA anatomies from eight patients with varying atrial function, including three with AF and either a LAA thrombus or a history of Transient Ischemic Attacks (TIAs). Using four different flow splits (60/40% and 55/45% through right and left PVs, even flow rate, and same velocity through each PV), we found that flow patterns are sensitive to PV flow split variations, particularly in planes parallel to the mitral valve. Changes in PV flow split also had a significant impact on blood stasis and could contribute to increased risk for thrombosis inside the LAA, particularly in patients with AF and previous LAA thrombus or a history of TIAs. Our study highlights the importance of considering patient-specific PV flow split variations when assessing LA hemodynamics and identifying patients at increased risk for thrombosis and stroke. This knowledge is relevant to planning clinical procedures such as AF ablation or the implementation of LAA occluders.

Myocardial Regional Shortening from 4D Cardiac CT Angiography for the Detection of Left Ventricular Segmental Wall Motion Abnormality.

Purpose: To investigate whether endocardial regional shortening computed from four-dimensional (4D) CT angiography (RSCT) can be used as a decision classifier to detect the presence of left ventricular (LV) wall motion abnormalities (WMAs). Materials and Methods: One hundred electrocardiographically gated cardiac 4D CT studies (mean age, 59 years ± 14 [SD]; 61 male patients) conducted between April 2018 and December 2020 were retrospectively evaluated. Three experts labeled LV wall motion in each of the 16 American Heart Association (AHA) segments as normal or abnormal; they also measured peak RSCT across one heartbeat in each segment. The data set was split evenly into training and validation groups. During training, interchangeability of RSCT thresholding with experts to detect WMA was assessed using the individual equivalence index (γ), and an optimal threshold of the peak RSCT (RSCT*) that achieved maximum agreement was identified. RSCT* was then validated using the validation group, and the effect of AHA segment-specific thresholds was evaluated. Agreement was assessed using κ statistics. Results: The optimal threshold, RSCT* of -0.19, when applied to all AHA segments, led to high agreement (agreement rate = 92.17%, κ = 0.82) and interchangeability with experts (γ = -2.58%). The same RSCT* also achieved high agreement in the validation group (agreement rate = 90.29%, κ = 0.76, γ = -0.38%). The use of AHA segment-specific thresholds (range: 0.16 to -0.23 across AHA segments) slightly improved agreement (1.79% increase). Conclusion: RSCT thresholding was interchangeable with expert visual analysis in detecting segmental WMA from 4D CT and may be used as an objective decision classifier.Keywords: CT, Left Ventricle, Regional Endocardial Shortening, Wall Motion Abnormality Supplemental material is available for this article. © RSNA, 2023.

Detection of left ventricular wall motion abnormalities from volume rendering of 4DCT cardiac angiograms using deep learning.

Background: The presence of left ventricular (LV) wall motion abnormalities (WMA) is an independent indicator of adverse cardiovascular events in patients with cardiovascular diseases. We develop and evaluate the ability to detect cardiac wall motion abnormalities (WMA) from dynamic volume renderings (VR) of clinical 4D computed tomography (CT) angiograms using a deep learning (DL) framework. Methods: Three hundred forty-three ECG-gated cardiac 4DCT studies (age: 61 ± 15, 60.1% male) were retrospectively evaluated. Volume-rendering videos of the LV blood pool were generated from 6 different perspectives (i.e., six views corresponding to every 60-degree rotation around the LV long axis); resulting in 2058 unique videos. Ground-truth WMA classification for each video was performed by evaluating the extent of impaired regional shortening visible (measured in the original 4DCT data). DL classification of each video for the presence of WMA was performed by first extracting image features frame-by-frame using a pre-trained Inception network and then evaluating the set of features using a long short-term memory network. Data were split into 60% for 5-fold cross-validation and 40% for testing. Results: Volume rendering videos represent ~800-fold data compression of the 4DCT volumes. Per-video DL classification performance was high for both cross-validation (accuracy = 93.1%, sensitivity = 90.0% and specificity = 95.1%, κ: 0.86) and testing (90.9, 90.2, and 91.4% respectively, κ: 0.81). Per-study performance was also high (cross-validation: 93.7, 93.5, 93.8%, κ: 0.87; testing: 93.5, 91.9, 94.7%, κ: 0.87). By re-binning per-video results into the 6 regional views of the LV we showed DL was accurate (mean accuracy = 93.1 and 90.9% for cross-validation and testing cohort, respectively) for every region. DL classification strongly agreed (accuracy = 91.0%, κ: 0.81) with expert visual assessment. Conclusions: Dynamic volume rendering of the LV blood pool combined with DL classification can accurately detect regional WMA from cardiac CT.

Prediction of Cardiac Resynchronization Therapy Response Using a Lead Placement Score Derived From 4-Dimensional Computed Tomography.

BACKGROUND: Cardiac resynchronization therapy (CRT) is an effective treatment for patients with heart failure; however, 30% of patients do not respond to the treatment. We sought to derive patient-specific left ventricle maps of lead placement scores (LPS) that highlight target pacing lead sites for achieving a higher probability of CRT response. METHODS: Eighty-two subjects recruited for the ImagingCRT trial (Empiric Versus Imaging Guided Left Ventricular Lead Placement in Cardiac Resynchronization Therapy) were retrospectively analyzed. All 82 subjects had 2 contrast-enhanced full cardiac cycle 4-dimensional computed tomography scans: a baseline and a 6-month follow-up scan. CRT response was defined as a reduction in computed tomography-derived end-systolic volume ≥15%. Eight left ventricle features derived from the baseline scans were used to train a support vector machine via a bagging approach. An LPS map over the left ventricle was created for each subject as a linear combination of the support vector machine feature weights and the subject's own feature vector. Performance for distinguishing responders was performed on the original 82 subjects. RESULTS: Fifty-two (63%) subjects were responders. Subjects with an LPS≤Q1 (lower-quartile) had a posttest probability of responding of 14% (3/21), while subjects with an LPS≥ Q3 (upper-quartile) had a posttest probability of responding of 90% (19/21). Subjects with Q1

Regional left ventricular endocardial strains estimated from low-dose 4DCT: Comparison with cardiac magnetic resonance feature tracking.

BACKGROUND: Estimates of regional left ventricular (LV) strains provide additional information to global function parameters such as ejection fraction (EF) and global longitudinal strain (GLS) and are more sensitive in detecting abnormal regional cardiac function. The accurate and reproducible assessment of regional cardiac function has implications in the management of various cardiac diseases such as heart failure, myocardial ischemia, and dyssynchrony. PURPOSE: To develop a method that yields highly reproducible, high-resolution estimates of regional endocardial strains from 4DCT images. METHODS: A method for estimating regional LV endocardial circumferential ( ε c c ) $( {{\epsilon }_{cc}} )$ and longitudinal ( ε l l ${\epsilon }_{ll}$ ) strains from 4DCT was developed. Point clouds representing the LV endocardial surface were extracted for each time frame of the cardiac cycle from 4DCT images. 3D deformation fields across the cardiac cycle were obtained by registering the end diastolic point cloud to each subsequent point cloud in time across the cardiac cycle using a 3D point-set registration technique. From these deformation fields, ε c c and ε l l ${\epsilon }_{cc}\ {\rm{and\ }}{\epsilon }_{ll}$ were estimated over the entire LV endocardial surface by fitting an affine transformation with maximum likelihood estimation. The 4DCT-derived strains were compared with strains estimated in the same subjects by cardiac magnetic resonance (CMR); twenty-four subjects had CMR scans followed by 4DCT scans acquired within a few hours. Regional LV circumferential and longitudinal strains were estimated from the CMR images using a commercially available feature tracking software (cvi42). Global circumferential strain (GCS) and global longitudinal strain (GLS) were calculated as the mean of the regional strains across the entire LV for both modalities. Pearson correlation coefficients and Bland-Altman analyses were used for comparisons. Intraclass correlation coefficients (ICC) were used to assess the inter- and intraobserver reproducibility of the 4DCT-derived strains. RESULTS: The 4DCT-derived regional strains correlated well with the CMR-derived regional strains ( ε c c ${\epsilon }_{cc}$ : r = 0.76, p < 0.001; ε l l ${\epsilon }_{ll}$ : r = 0.64, p < 0.001). A very strong correlation was found between 4DCT-derived GCS and 4DCT-derived EF (r = -0.96; p < 0.001). The 4DCT-derived strains were also highly reproducible, with very low inter- and intraobserver variability (intraclass correlation coefficients in the range of [0.92, 0.99]). CONCLUSIONS: We have developed a novel method to estimate high-resolution regional LV endocardial circumferential and longitudinal strains from 4DCT images. Except for the definition of the mitral valve and LV outflow tract planes, the method is completely user independent, thus yielding highly reproducible estimates of endocardial strain. The 4DCT-derived strains correlated well with those estimated using a commercial CMR feature tracking software. The promising results reported in this study highlight the potential utility of 4DCT in the precise assessment of regional cardiac function for the management of cardiac disease.

Four-dimensional computed tomography of the left ventricle, Part I: Motion artifact reduction.

PURPOSE: Standard four-dimensional computed tomography (4DCT) cardiac reconstructions typically include spiraling artifacts that depend not only on the motion of the heart but also on the gantry angle range over which the data was acquired. We seek to reduce these motion artifacts and, thereby, improve the accuracy of left ventricular wall positions in 4DCT image series. METHODS: We use a motion artifact reduction approach (ResyncCT) that is based largely on conjugate pairs of partial angle reconstruction (PAR) images. After identifying the key locations where motion artifacts exist in the uncorrected images, paired subvolumes within the PAR images are analyzed with a modified cross-correlation function in order to estimate 3D velocity and acceleration vectors at these locations. A subsequent motion compensation process (also based on PAR images) includes the creation of a dense motion field, followed by a backproject-and-warp style compensation. The algorithm was tested on a 3D printed phantom, which represents the left ventricle (LV) and on challenging clinical cases corrupted by severe artifacts. RESULTS: The results from our preliminary phantom test as well as from clinical cardiac scans show crisp endocardial edges and resolved double-wall artifacts. When viewed as a temporal series, the corrected images exhibit a much smoother motion of the LV endocardial boundary as compared to the uncorrected images. In addition, quantitative results from our phantom studies show that ResyncCT processing reduces endocardial surface distance errors from 0.9 ± 0.8 to 0.2 ± 0.1 mm. CONCLUSIONS: The ResyncCT algorithm was shown to be effective in reducing motion artifacts and restoring accurate wall positions. Some perspectives on the use of conjugate-PAR images and on techniques for CT motion artifact reduction more generally are also given.

Non-Newtonian blood rheology impacts left atrial stasis in patient-specific simulations.

The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects are subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered when calculating patient-specific blood stasis indices by computational fluid dynamics.

Four-dimensional computed tomography of the left ventricle, Part II: Estimation of mechanical activation times.

PURPOSE: We demonstrate the viability of a four-dimensional X-ray computed tomography (4DCT) imaging system to accurately and precisely estimate mechanical activation times of left ventricular (LV) wall motion. Accurate and reproducible timing estimates of LV wall motion may be beneficial in the successful planning and management of cardiac resynchronization therapy (CRT). METHODS: We developed an anthropomorphically accurate in silico LV phantom based on human CT images with programmed septal-lateral wall dyssynchrony. Twenty-six temporal phases of the in silico phantom were used to sample the cardiac cycle of 1 s. For each of the 26 phases, 1 cm thick axial slabs emulating axial CT image volumes were extracted, 3D printed, and imaged using a commercially available CT scanner. A continuous dynamic sinogram was synthesized by blending sinograms from these static phases; the synthesized sinogram emulated the sinogram that would be acquired under true continuous phantom motion. Using the synthesized dynamic sinogram, images were reconstructed at 70 ms intervals spanning the full cardiac cycle; these images exhibited expected motion artifact characteristics seen in images reconstructed from real dynamic data. The motion corrupted images were then processed with a novel motion correction algorithm (ResyncCT) to yield motion corrected images. Five pairs of motion uncorrected and motion corrected images were generated, each corresponding to a different starting gantry angle (0 to 180 degrees in 45 degree increments). Two line profiles perpendicular to the endocardial surface were used to sample local myocardial motion trajectories at the septum and the lateral wall. The mechanical activation time of wall motion was defined as the time at which the endocardial boundary crossed a fixed position defined on either of the two line profiles while moving toward the center of the LV during systolic contraction. The mechanical activation times of these myocardial trajectories estimated from the motion uncorrected and the motion corrected images were then compared with those derived from the static images of the 3D printed phantoms (ground truth). The precision of the timing estimates was obtained from the five different starting gantry angle simulations. RESULTS: The range of estimated mechanical activation times observed across all starting gantry angles was significantly larger for the motion uncorrected images than for the motion corrected images (lateral wall: 58 ± 15 ms vs 12 ± 4 ms, p < 0.005; septal wall: 61 ± 13 ms vs 13 ± 9 ms, p < 0.005). CONCLUSIONS: 4DCT images processed with the ResyncCT motion correction algorithm yield estimates of mechanical activation times of LV wall motion with significantly improved accuracy and precision. The promising results reported in this study highlight the potential utility of 4DCT in estimating the timing of mechanical events of interest for CRT guidance.

Computational ECG mapping and respiratory gating to optimize stereotactic ablative radiotherapy workflow for refractory ventricular tachycardia.

BACKGROUND: Stereotactic ablative radiotherapy (SAbR) is an emerging therapy for refractory ventricular tachycardia (VT). However, the current workflow is complicated, and the precision and safety in patients with significant cardiorespiratory motion and VT targets near the stomach may be suboptimal. OBJECTIVE: We hypothesized that automated 12-lead electrocardiogram (ECG) mapping and respiratory-gated therapy may improve the ease and precision of SAbR planning and facilitate safe radiation delivery in patients with refractory VT. METHODS: Consecutive patients with refractory VT were studied at 2 hospitals. VT exit sites were localized using a 3-D computational ECG algorithm noninvasively and compared to available prior invasive mapping. Radiotherapy (25 Gy) was delivered at end-expiration when cardiac respiratory motion was ≥0.6 cm or targets were ≤2 cm from the stomach. RESULTS: In 6 patients (ejection fraction 29% ± 13%), 4.2 ± 2.3 VT morphologies per patient were mapped. Overall, 7 out of 7 computational ECG mappings (100%) colocalized to the identical cardiac segment when prior invasive electrophysiology study was available. Respiratory gating was associated with smaller planning target volumes compared to nongated volumes (71 ± 7 vs 153 ± 35 cc, P < .01). In 2 patients with inferior wall VT targets close to the stomach (6 mm proximity) or significant respiratory motion (22 mm excursion), no GI complications were observed at 9- and 12-month follow-up. Implantable cardioverter-defibrillator shocks decreased from 23 ± 12 shocks/patient to 0.67 ± 1.0 (P < .001) post-SAbR at 6.0 ± 4.9 months follow-up. CONCLUSIONS: A workflow including computational ECG mapping and protocol-guided respiratory gating is feasible, is safe, and may improve the ease of SAbR planning. Studies to validate this workflow in larger populations are required.

Novel 4DCT Method to Measure Regional Left Ventricular Endocardial Shortening Before and After Transcatheter Mitral Valve Implantation.

BACKGROUND: Regional left ventricular (LV) mechanics in mitral regurgitation (MR) patients, and local changes in function after transcatheter mitral valve implantation (TMVI) have yet to be evaluated. Herein, we introduce a method for creating high resolution maps of endocardial function from 4DCT images, leading to detailed characterization of changes in local LV function. These changes are particularly interesting when evaluating the effect of the Tendyne™ TMVI device in the region of the epicardial pad. METHODS: Regional endocardial shortening from CT (RSCT) was evaluated in Tendyne (Abbott Medical) TMVI patients with 4DCT exams pre- and post-implantation. Regional function was evaluated in 90 LV segments (5 longitudinal × 18 circumferential). LV volumes and ejection fraction (EF) were also computed. A reproducibility study was performed in a subset of patients to determine the precision of RSCT measurements in this population. RESULTS: Baseline and local changes in RSCT post TMVI were highly variable and extremely spatially heterogeneous. Both inter- and intra-observer variability were low and demonstrated the high precision of RSCT for evaluating regional LV function. CONCLUSION: RSCT is a reproducible metric which can be evaluated in patients with highly abnormal regional LV function and geometry. After TMVI, significant spatially heterogeneous changes in RSCT were observed in all subjects; therefore, it is unlikely that the functional state of TMVI patients can be fully described by changes in LV volume or EF. Measurement of RSCT provides precise characterization of the spatially heterogeneous effects of MR and TMVI on LV function and remodeling.

Demonstration of Patient-Specific Simulations to Assess Left Atrial Appendage Thrombogenesis Risk.

Atrial fibrillation (AF) alters left atrial (LA) hemodynamics, which can lead to thrombosis in the left atrial appendage (LAA), systemic embolism and stroke. A personalized risk-stratification of AF patients for stroke would permit improved balancing of preventive anticoagulation therapies against bleeding risk. We investigated how LA anatomy and function impact LA and LAA hemodynamics, and explored whether patient-specific analysis by computational fluid dynamics (CFD) can predict the risk of LAA thrombosis. We analyzed 4D-CT acquisitions of LA wall motion with an in-house immersed-boundary CFD solver. We considered six patients with diverse atrial function, three with either a LAA thrombus (removed digitally before running the simulations) or a history of transient ischemic attacks (LAAT/TIA-pos), and three without a LAA thrombus or TIA (LAAT/TIA-neg). We found that blood inside the left atrial appendage of LAAT/TIA-pos patients had marked alterations in residence time and kinetic energy when compared with LAAT/TIA-neg patients. In addition, we showed how the LA conduit, reservoir and booster functions distinctly affect LA and LAA hemodynamics. Finally, fixed-wall and moving-wall simulations produced different LA hemodynamics and residence time predictions for each patient. Consequently, fixed-wall simulations risk-stratified our small cohort for LAA thrombosis worse than moving-wall simulations, particularly patients with intermediate LAA residence time. Overall, these results suggest that both wall kinetics and LAA morphology contribute to LAA blood stasis and thrombosis.

Targeted Coronary Artery Calcium Screening in High-Risk Younger Individuals Using Consumer Genetic Screening Results.

OBJECTIVES: The goal of this study was to assess the utility of a genetic risk score (GRS) in targeted coronary artery calcium (CAC) screening among young individuals. BACKGROUND: Early CAC screening and preventive therapy may reduce long-term risk of a coronary heart disease (CHD) event. However, identifying younger individuals at increased risk remains a challenge. GRS for CHD are age independent and can stratify individuals on various risk trajectories. METHODS: Using 142 variants associated with CHD events, we calculated a GRS in 1,927 individuals in the CARDIA (Coronary Artery Risk Development in Young Adults) cohort (aged 32 to 47 years) and 6,600 individuals in the MESA (Multi-Ethnic Study of Atherosclerosis) cohort (aged 44 to 87 years). We assessed GRS utility to predict CAC presence in the CARDIA cohort and stratify individuals of varying risk for CAC presence over the lifetime in both cohorts. RESULTS: The GRS predicted CAC presence in CARDIA males. It was not predictive in CARDIA females, which had a CAC prevalence of 6.4%. In combined analysis of the CARDIA and MESA cohorts, the GRS was predictive of CAC in both males and females and was used to derive an equation for the age at which CAC probability crossed a predetermined threshold. When assessed in combination with traditional risk factors, the GRS further stratified individuals. For individuals with an equal number of traditional risk factors, probability of CAC reached 25% approximately 10 years earlier for those in the highest GRS quintile compared to the lowest. CONCLUSIONS: The GRS may be used to target high-risk younger individuals for early CAC screening.

Constructing a Human Atrial Fibre Atlas.

Atrial anisotropy affects electrical propagation patterns, anchor locations of atrial reentrant drivers, and atrial mechanics. However, patient-specific atrial fibre fields and anisotropy measurements are not currently available, and consequently assigning fibre fields to atrial models is challenging. We aimed to construct an atrial fibre atlas from a high-resolution DTMRI dataset that optimally reproduces electrophysiology simulation predictions corresponding to patient-specific fibre fields, and to develop a methodology for automatically assigning fibres to patient-specific anatomies. We extended an atrial coordinate system to map the pulmonary veins, vena cava and appendages to standardised positions in the coordinate system corresponding to the average location across the anatomies. We then expressed each fibre field in this atrial coordinate system and calculated an average fibre field. To assess the effects of fibre field on patient-specific modelling predictions, we calculated paced activation time maps and electrical driver locations during AF. In total, 756 activation time maps were calculated (7 anatomies with 9 fibre maps and 2 pacing locations, for the endocardial, epicardial and bilayer surface models of the LA and RA). Patient-specific fibre fields had a relatively small effect on average paced activation maps (range of mean local activation time difference for LA fields: 2.67-3.60 ms, and for RA fields: 2.29-3.44 ms), but had a larger effect on maximum LAT differences (range for LA 12.7-16.6%; range for RA 11.9-15.0%). A total of 126 phase singularity density maps were calculated (7 anatomies with 9 fibre maps for the LA and RA bilayer models). The fibre field corresponding to anatomy 1 had the highest median PS density map correlation coefficient for LA bilayer simulations (0.44 compared to the other correlations, ranging from 0.14 to 0.39), while the average fibre field had the highest correlation for the RA bilayer simulations (0.61 compared to the other correlations, ranging from 0.37 to 0.56). For sinus rhythm simulations, average activation time is robust to fibre field direction; however, maximum differences can still be significant. Patient specific fibres are more important for arrhythmia simulations, particularly in the left atrium. We propose using the fibre field corresponding to DTMRI dataset 1 for LA simulations, and the average fibre field for RA simulations as these optimally predicted arrhythmia properties.

M-SiSSR: Regional Endocardial Function Using Multilabel Simultaneous Subdivision Surface Registration.

Quantification of regional cardiac function is a central goal of cardiology. Multiple methods, such as Coherent Point Drift (CPD) and Simultaneous Subdivision Surface Registration (SiSSR), have been used to register meshes to the endocardial surface. However, these methods do not distinguish between cardiac chambers during registration, and consequently the mesh may "slip" across the interface between two structures during contraction, resulting in inaccurate regional functional measurements. Here, we present Multilabel-SiSSR (M-SiSSR), a novel method for registering a "labeled" cardiac mesh (with each triangle assigned to a cardiac structure). We compare our results to the original, label-agnostic version of SiSSR and find both a visual and quantitative improvement in tracking of the mitral valve plane.