Purpose: Dual-energy computed tomography (DECT) using TwinBeam CT (TBCT) is a new option for radiation oncology simulators. TBCT scanning provides virtual monoenergetic images which are attractive in treatment planning since lower energies offer better contrast for soft tissues, and higher energies reduce noise. A protocol is needed to achieve optimal performance of this feature. In this study, we investigated the TBCT scan schema with the head-and-neck radiotherapy workflow at our clinic and selected the optimal energy with best contrast-noise-ratio (CNR) in organs-at-risks (OARs) delineation for head-and-neck treatment planning. Methods and materials: We synthesized monochromatic images from 40 keV to 190 keV at 5 keV increments from data acquired by TBCT. We collected the Hounsfield unit (HU) numbers of OARs (brainstem, mandible, spinal cord, and parotid glands), the HU numbers of marginal regions outside OARs, and the noise levels for each monochromatic image. We then calculated the CNR for the different OARs at each energy level to generate a serial of spectral curves for each OAR. Based on these spectral curves of CNR, the mono-energy corresponding to the max CNR was identified for each OAR of each patient. Results: Computed tomography scans of ten patients by TBCT were used to test the optimal monoenergetic image for the CNR of OAR. Based on the maximized CNR, the optimal energy values were 78.5 ± 5.3 keV for the brainstem, 78.0 ± 4.2 keV for the mandible, 78.5 ± 5.7 keV for the parotid glands, and 78.5 ± 5.3 keV for the spinal cord. Overall, the optimal energy for the maximum CNR of these OARs in head-and-neck cancer patients was 80 keV. Conclusion: We have proposed a clinically feasible protocol that selects the optimal energy level of the virtual monoenergetic image in TBCT for OAR delineation based on the CNR in head-and-neck OAR. This protocol can be applied in TBCT simulation.
To conduct a patient-specific computational modeling of the aortic valve, 3-D aortic valve anatomic geometries of an individual patient need to be reconstructed from clinical 3-D cardiac images. Currently, most of computational studies involve manual heart valve geometry reconstruction and manual finite element (FE) model generation, which is both time-consuming and prone to human errors. A seamless computational modeling framework, which can automate this process based on machine learning algorithms, is desirable, as it can not only eliminate human errors and ensure the consistency of the modeling results but also allow fast feedback to clinicians and permits a future population-based probabilistic analysis of large patient cohorts. In this study, we developed a novel computational modeling method to automatically reconstruct the 3-D geometries of the aortic valve from computed tomographic images. The reconstructed valve geometries have built-in mesh correspondence, which bridges harmonically for the consequent FE modeling. The proposed method was evaluated by comparing the reconstructed geometries from 10 patients with those manually created by human experts, and a mean discrepancy of 0.69 mm was obtained. Based on these reconstructed geometries, FE models of valve leaflets were developed, and aortic valve closure from end systole to middiastole was simulated for 7 patients and validated by comparing the deformed geometries with those manually created by human experts, and a mean discrepancy of 1.57 mm was obtained. The proposed method offers great potential to streamline the computational modeling process and enables the development of a preoperative planning system for aortic valve disease diagnosis and treatment.
We propose a simple method for reconstructing vascular trees from 3D images. Our algorithm extracts persistent maxima of the intensity on all axis-aligned 2D slices of the input image. The maxima concentrate along 1D intensity ridges, in particular along blood vessels. We build a forest connecting the persistent maxima with short edges. The forest tends to approximate the blood vessels present in the image, but also contains numerous spurious features and often fails to connect segments belonging to one vessel in low contrast areas. We improve the forest by applying simple geometric filters that trim short branches, fill gaps in blood vessels and remove spurious branches from the vascular tree to be extracted. Experiments show that our technique can be applied to extract coronary trees from heart CT scans.
Cardiovascular simulations have great potential as a clinical tool for planning and evaluating patient-specific treatment strategies for those suffering from congenital heart diseases, specifically Fontan patients. However, several bottlenecks have delayed wider deployment of the simulations for clinical use; the main obstacle is simulation cost. Currently, time-averaged clinical flow measurements are utilized as numerical boundary conditions (BCs) in order to reduce the computational power and time needed to offer surgical planning within a clinical time frame. Nevertheless, pulsatile blood flow is observed in vivo, and its significant impact on numerical simulations has been demonstrated. Therefore, it is imperative to carry out a comprehensive study analyzing the sensitivity of using time-averaged BCs. In this study, sensitivity is evaluated based on the discrepancies between hemodynamic metrics calculated using time-averaged and pulsatile BCs; smaller discrepancies indicate less sensitivity. The current study incorporates a comparison between 3D patient-specific CFD simulations using both the time-averaged and pulsatile BCs for 101 Fontan patients. The sensitivity analysis involves two clinically important hemodynamic metrics: hepatic flow distribution (HFD) and indexed power loss (iPL). Paired demographic group comparisons revealed that HFD sensitivity is significantly different between single and bilateral superior vena cava cohorts but no other demographic discrepancies were observed for HFD or iPL. Multivariate regression analyses show that the best predictors for sensitivity involve flow pulsatilities, time-averaged flow rates, and geometric characteristics of the Fontan connection. These predictors provide patient-specific guidelines to determine the effectiveness of analyzing patient-specific surgical options with time-averaged BCs within a clinical time frame.
Cardiac computed tomography (CT) angiography using prospective gating requires that data be acquired during intervals of minimal cardiac motion to obtain diagnostic images of the coronary vessels free of motion artifacts. This work is intended to assess B-mode echocardiography as a continuous-time indication of these quiescent periods to determine if echocardiography can be used as a cost-efficient, non-ionizing modality to develop new prospective gating techniques for cardiac CT. These new prospective gating approaches will not be based on echocardiography itself but on CT-compatible modalities derived from the mechanics of the heart (e.g. seismocardiography and impedance cardiography), unlike the current standard electrocardiogram. To this end, echocardiography and retrospectively-gated CT data were obtained from ten patients with varied cardiac conditions. CT reconstructions were made throughout the cardiac cycle. Motion of the interventricular septum (IVS) was calculated from both echocardiography and CT reconstructions using correlation-based, deviation techniques. The IVS was chosen because it (1) is visible in echocardiography images, whereas the coronary vessels generally are not, and (2) has been shown to be a suitable indicator of cardiac quiescence. Quiescent phases were calculated as the minima of IVS motion and CT volumes were reconstructed for these phases. The diagnostic quality of the CT reconstructions from phases calculated from echocardiography and CT data was graded on a four-point Likert scale by a board-certified radiologist fellowship-trained in cardiothoracic radiology. Using a Wilcoxon signed-rank test, no significant difference in the diagnostic quality of the coronary vessels was found between CT volumes reconstructed from echocardiography- and CT-selected phases. Additionally, there was a correlation of 0.956 between the echocardiography- and CT-selected phases. This initial work suggests that B-mode echocardiography can be used as a tool to develop CT-compatible gating techniques based on modalities derived from cardiac mechanics rather than relying on the ECG alone.
by
Kartik S. Sundareswaran;
Christopher M. Haggerty;
Diane de Zelicourt;
Lakshmi P. Dasi;
Kerem Pekkan;
David H. Frakes;
Andrew J. Powell;
Kirk R Kanter;
Mark A. Fogel;
Ajit Yoganathan
Objective: Our objective was to analyze 3-dimensional (3D) blood flow patterns within the total cavopulmonary connection (TCPC) using in vivo phase contrast magnetic resonance imaging (PC MRI).
Methods: Sixteen single-ventricle patients were prospectively recruited at 2 leading pediatric institutions for PC MRI evaluation of their Fontan pathway. Patients were divided into 2 groups. Group 1 comprised 8 patients with an extracardiac (EC) TCPC, and group 2 comprised 8 patients with a lateral tunnel (LT) TCPC. A coronal stack of 5 to 10 contiguous PC MRI slices with 3D velocity encoding (5-9 ms resolution) was acquired and a volumetric flow field was reconstructed.
Results: Analysis revealed large vortices in LT TCPCs and helical flow structures in EC TCPCs. On average, there was no difference between LT and EC TCPCs in the proportion of inferior vena cava flow going to the left pulmonary artery (43% ± 7% vs 46% ± 5%; P = .34). However, for EC TCPCs, the presence of a caval offset was a primary determinant of inferior vena caval flow distribution to the pulmonary arteries with a significant bias to the offset side.
Conclusions: 3D flow structures within LT and EC TCPCs were reconstructed and analyzed for the first time using PC MRI. TCPC flow patterns were shown to be different, not only on the basis of LT or EC considerations, but with significant influence from the superior vena cava connection as well. This work adds to the ongoing body of research demonstrating the impact of TCPC geometry on the overall hemodynamic profile.
Over the last few years, computed tomography (CT) has developed into a standard clinical test for a variety of cardiovascular conditions. The emergence of cardiovascular CT during a period of dramatic increase in radiation exposure to the population from medical procedures and heightened concern about the subsequent potential cancer risk has led to intense scrutiny of the radiation burden of this new technique. This has hastened the development and implementation of dose reduction tools and prompted closer monitoring of patient dose. In an effort to aid the cardiovascular CT community in incorporating patient-centered radiation dose optimization and monitoring strategies into standard practice, the Society of Cardiovascular Computed Tomography has produced a guideline document to review available data and provide recommendations regarding interpretation of radiation dose indices and predictors of risk, appropriate use of scanner acquisition modes and settings, development of algorithms for dose optimization, and establishment of procedures for dose monitoring.
Total cavopulmonary connection is the result of a series of palliative surgical repairs performed on patients with single ventricle heart defects. The resulting anatomy has complex and unsteady hemodynamics characterized by flow mixing and flow separation. Although varying degrees of flow pulsatility have been observed in vivo, non-pulsatile (time-averaged) boundary conditions have traditionally been assumed in hemodynamic modeling, and only recently have pulsatile conditions been incorporated without completely characterizing their effect or importance. In this study, 3D numerical simulations with both pulsatile and non-pulsatile boundary conditions were performed for 24 patients with different anatomies and flow boundary conditions from Georgia Tech database. Flow structures, energy dissipation rates and pressure drops were compared under rest and simulated exercise conditions. It was found that flow pulsatility is the primary factor in determining the appropriate choice of boundary conditions, whereas the anatomic configuration and cardiac output had secondary effects. Results show that the hemodynamics can be strongly influenced by the presence of pulsatile flow. However, there was a minimum pulsatility threshold, identified by defining a weighted pulsatility index (wPI), above which the influence was significant. It was shown that when wPI < 30%, the relative error in hemodynamic predictions using time-averaged boundary conditions was less than 10% compared to pulsatile simulations. In addition, when wPI < 50, the relative error was less than 20%. A correlation was introduced to relate wPI to the relative error in predicting the flow metrics with non-pulsatile flow conditions.