Quantitative evaluation of myocardial function by a volume-normalized map generated from relative blood flow

Our study aimed to quantitatively evaluate blood flow in the left ventricle (LV) of apical hypertrophic cardiomyopathy (APH) by combining wall thickness obtained from cardiac magnetic resonance imaging (MRI) and myocardial perfusion from single-photon emission computed tomography (SPECT). In this study, we considered paired MRI and myocardial perfusion SPECT from ten patients with APH and ten normals. Myocardial walls were detected using a level set method, and blood flow per unit myocardial volume was calculated using 3D surface-based registration between the MRI and SPECT images. We defined relative blood flow based on the maximum in the whole myocardial region. Accuracies of wall detection and registration were around 2.50 mm and 2.95 mm, respectively. We finally created a bull's-eye map to evaluate wall thickness, blood flow (cardiac perfusion) and blood flow per unit myocardial volume. In patients with APH, their wall thicknesses were over 10 mm. Decreased blood flow per unit myocardial volume was detected in the cardiac apex by calculation using wall thickness from MRI and blood flow from SPECT. The relative unit blood flow of the APH group was 1/7 times that of the normals in the apex. This normalization by myocardial volume distinguishes cases of APH whose SPECT images resemble the distributions of normal cases.     

Automated cardiac motion compensation in PET/CT for accurate reconstruction of PET myocardial perfusion images

Error-free reconstruction of PET data with a registered CT attenuation map is essential for accurate quantification and interpretation of cardiac perfusion. Misalignment of the CT and PET data can produce an erroneous attenuation map that projects lung attenuation parameters onto the heart wall, thereby underestimating the attenuation and creating artifactual areas of hypoperfusion that can be misinterpreted as myocardial ischemia or infarction. The major causes of misregistration between CT and PET images are the respiratory motion, cardiac motion and gross physical motion of the patient. The misalignment artifact problem is overcome with automated cardiac registration software that minimizes the alignment error between the two modalities. Results show that the automated registration process works equally well for any respiratory phase in which the CT scan is acquired. Further evaluation of this procedure on 50 patients demonstrates that the automated registration software consistently aligns the two modalities, eliminating artifactual hypoperfusion in reconstructed PET images due to PET/CT misregistration. With this registration software, only one CT scan is required for PET/CT imaging, which reduces the radiation dose required for CT-based attenuation correction and improves the clinical workflow for PET/CT.     

A continuous 4D motion model from multiple respiratory cycles for use in lung radiotherapy

Respiratory motion causes errors when planning and delivering radiotherapy treatment to lung cancer patients. To reduce these errors, methods of acquiring and using four-dimensional computed tomography (4DCT) datasets have been developed. We have developed a novel method of constructing computational motion models from 4DCT. The motion models attempt to describe an average respiratory cycle, which reduces the effects of variation between different cycles. They require substantially less memory than a 4DCT dataset, are continuous in space and time, and facilitate automatic target propagation and combining of doses over the respiratory cycle. The motion models are constructed from CT data acquired in cine mode while the patient is free breathing (free breathing CT - FBCT). A "slab" of data is acquired at each couch position, with 3-4 contiguous slabs being acquired per patient. For each slab a sequence of 20 or 30 volumes was acquired over 20 seconds. A respiratory signal is simultaneously recorded in order to calculate the position in the respiratory cycle for each FBCT. Additionally, a high quality reference CT volume is acquired at breath hold. The reference volume is nonrigidly registered to each of the FBCT volumes. A motion model is then constructed for each slab by temporally fitting the nonrigid registration results. The value of each of the registration parameters is related to the position in the respiratory cycle by fitting an approximating B spline to the registration results. As an approximating function is used, and the data is acquired over several respiratory cycles, the function should model an average respiratory cycle. This can then be used to calculate the value of each degree of freedom at any desired position in the respiratory cycle. The resulting nonrigid transformation will deform the reference volume to predict the contents of the slab at the desired position in the respiratory cycle. The slab model predictions are then concatenated to produce a combined prediction over the entire region of interest. We have performed a number of experiments to assess the accuracy of the nonrigid registration results and the motion model predictions. The individual slab models were evaluated by expert visual assessment and the tracking of easily identifiable anatomical points. The combined models were evaluated by calculating the discontinuities between the transformations at the slab boundaries. The experiments were performed on five patients with a total of 18 slabs between them. For the point tracking experiments, the mean distance between where a clinician manually identified a point and where the registration results located the point, the target registration error (TRE), was 1.3 mm. The mean distance between a manually identified point and the models prediction of the point's location, the target model error (TME), was 1.6 mm. The mean discontinuity between model predictions at the slab boundaries, the Continuity Error, was 2.2 mm. The results show that the motion models perform with a level of accuracy comparable to the slice thickness of 1.5 mm.     

Dynamic ventilation imaging from four-dimensional computed tomography

A novel method for dynamic ventilation imaging of the full respiratory cycle from four-dimensional computed tomography (4D CT) acquired without added contrast is presented. Three cases with 4D CT images obtained with respiratory gated acquisition for radiotherapy treatment planning were selected. Each of the 4D CT data sets was acquired during resting tidal breathing. A deformable image registration algorithm mapped each (voxel) corresponding tissue element across the 4D CT data set. From local average CT values, the change in fraction of air per voxel (i.e. local ventilation) was calculated. A 4D ventilation image set was calculated using pairs formed with the maximum expiration image volume, first the exhalation then the inhalation phases representing a complete breath cycle. A preliminary validation using manually determined lung volumes was performed. The calculated total ventilation was compared to the change in contoured lung volumes between the CT pairs (measured volume). A linear regression resulted in a slope of 1.01 and a correlation coefficient of 0.984 for the ventilation images. The spatial distribution of ventilation was found to be case specific and a 30% difference in mass-specific ventilation between the lower and upper lung halves was found. These images may be useful in radiotherapy planning.     

A new approach towards volumetric assessment of left ventricular function with MSCT

Cardiovascular CT is considered the diagnostic standard for establishing the presence of a functional and dynamic imaging system. It is difficult, however, to estimate the ventricular motion and volumes that are processed using hundreds and thousands of CT images, in a few moments.The main concept and design of our work are two fold - the development of effective semi-automatic tools for measuring the sequential left ventricular volumes from the hundreds or thousands of cardiac trans-axial images, and providing a simple interface with an interactive diagnostic tool for the volumetry of left ventricle and valuable cardiac 4D visualisation.We converted ten and more sequential volume data sets of the heart acquired from retrospective ECG-gating helical scan into 3D images by volume rendering. These sequential 3D images could be displayed as a movie (4D cardiac image) file. Furthermore, we developed a method for semi-automatic calculation of ejection fraction (EF) and cardiac cycle (%)-volume (ml) curve for estimation of the motion and the volume of the left ventricle. This method involved the use an interactive selection tool in the region of interest (ROI). All 3D processing methods, such as, cutting objects, segmentation, and image fusion were based on mask processing data. We now describe the software developed for cardiac 4D imaging and the estimation of ventricular volume.     

T1 measurement using a short acquisition period for quantitative cardiac applications

Myocardial signal intensity curves for myocardial perfusion studies may be made quantitative by the use of T1 measurements made after the first-pass of contrast agent. A short data acquisition method for T1 mapping is presented in which all data for each T1 map are acquired in a short breath hold, and the slice geometry and timing in the cardiac cycle exactly match that of the dynamic first-pass perfusion sequence. This allows accurate image registration of the T1 map with the first-pass series of images. The T1 method is based on varying the preparation-pulse delay time of a saturation recovery sequence, and in this implementation employs an ECG-triggered, single-shot, spoiled gradient echo technique with SENSE reconstruction. The method allows T1 estimates of three slices to be made in fifteen heartbeats. For a range of samples with T1 values equivalent to those found in the myocardium during the first-pass of contrast agent, T1 estimates were accurate to within 6%, and the variation between slices was 2% or less.     

Reconstruction of dynamic gated cardiac SPECT

In this paper we propose an image reconstruction procedure which aims to unify gated single photon emission computed tomography (SPECT) and dynamic SPECT into a single method. We divide the cardiac cycle into a number of gate intervals as in gated SPECT, but treat the tracer distribution for each gate as a time-varying signal. By using both dynamic and motion-compensated temporal regularization, our reconstruction procedure will produce an image sequence that shows both cardiac motion and time-varying tracer distribution simultaneously. To demonstrate the proposed reconstruction method, we simulated gated cardiac perfusion imaging using the gated mathematical cardiac-torso (gMCAT) phantom with Tc99m-Teboroxime as the imaging agent. Our results show that the proposed method can produce more accurate reconstruction of gated dynamic images than independent reconstruction of individual gate frames with spatial smoothness alone. In particular, our results show that the former could improve the contrast to noise ratio of a simulated perfusion defect by as much as 100% when compared to the latter.     

In vivo liver tracking with a high volume rate 4D ultrasound scanner and a 2D matrix array probe

The effectiveness of intensity-modulated radiation therapy (IMRT) is compromised by involuntary motion (e.g. respiration, cardiac activity). The feasibility of processing ultrasound echo data to automatically estimate 3D liver motion for real-time IMRT guidance was previously demonstrated, but performance was limited by an acquisition speed of 2 volumes per second due to hardware restrictions of a mechanical linear array probe. Utilizing a 2D matrix array probe with parallel receive beamforming offered increased acquisition speeds and an opportunity to investigate the benefits of higher volume rates. In?vivo livers of three volunteers were scanned with and without respiratory motion at volume rates of 24 and 48 Hz, respectively. Respiration was suspended via voluntary breath hold. Correlation-based, phase-sensitive 3D speckle tracking was applied to consecutively acquired volumes of echo data. Volumes were omitted at fixed intervals and 3D speckle tracking was re-applied to study the effect of lower scan rates. Results revealed periodic motion that corresponded with the heart rate or breathing cycle in the absence or presence of respiration, respectively. For cardiac-induced motion, volume rates for adequate tracking ranged from 8 to 12 Hz and was limited by frequency discrepancies between tracking estimates from higher and lower frequency scan rates. Thus, the scan rate of volume data acquired without respiration was limited by the need to sample the frequency induced by the beating heart. In respiratory-dominated motion, volume rate limits ranged from 4 to 12 Hz, interpretable from the root-mean-squared deviation (RMSD) from tracking estimates at 24 Hz. While higher volume rates yielded RMSD values less than 1 mm in most cases, lower volume rates yielded RMSD values of 2-6 mm.     

4D-CT motion estimation using deformable image registration and 5D respiratory motion modeling

Four-dimensional computed tomography (4D-CT) imaging technology has been developed for radiation therapy to provide tumor and organ images at the different breathing phases. In this work, a procedure is proposed for estimating and modeling the respiratory motion field from acquired 4D-CT imaging data and predicting tissue motion at the different breathing phases. The 4D-CT image data consist of series of multislice CT volume segments acquired in ciné mode. A modified optical flow deformable image registration algorithm is used to compute the image motion from the CT segments to a common full volume 3D-CT reference. This reference volume is reconstructed using the acquired 4D-CT data at the end-of-exhalation phase. The segments are optimally aligned to the reference volume according to a proposed a priori alignment procedure. The registration is applied using a multigrid approach and a feature-preserving image downsampling maxfilter to achieve better computational speed and higher registration accuracy. The registration accuracy is about 1.1 +/- 0.8 mm for the lung region according to our verification using manually selected landmarks and artificially deformed CT volumes. The estimated motion fields are fitted to two 5D (spatial 3D+tidal volume+airflow rate) motion models: forward model and inverse model. The forward model predicts tissue movements and the inverse model predicts CT density changes as a function of tidal volume and airflow rate. A leave-one-out procedure is used to validate these motion models. The estimated modeling prediction errors are about 0.3 mm for the forward model and 0.4 mm for the inverse model.     

Towards dynamic cardiac scenes interpretation based on spatial-temporal knowledge

Cardiac motion analysis enables to identify pathologies related to myocardial anomalies or coronary arteries circulation deficiencies. Conventionally, bi-dimensional (2D) left ventricle contour images have been extensively used, to perform quantitative measurements and qualitative evaluations of the cardiac function. Nevertheless, there are other cardiac anatomical structures, the coronary arteries, imaged on routine procedures, upon which complementary motion interpretation can be conducted. This paper presents an experimental methodology to perform dynamic cardiac scenes interpretation, studying three-dimensional (3D) coronary arteries spatial-temporal behavior. Being an alternative way to approach computer assisted cardiac motion interpretation, it reveals a wide range of rarely explored spatial-temporal situations and proposes how to address them. Considering the challenges to achieve dynamic scene interpretation, it is explained how spatial and temporal knowledge, are connected to specialist knowledge and measured parameters, to obtain a dynamic scene interpretation. Global and local motion features are modeled according to cardiac motion and geometrical knowledge, before its transformation into symbols. Anatomical knowledge and spatial-temporal knowledge are applied, along with spatial-temporal reasoning schemes, to access symbols meaning. Experimental results obtained using real data are presented. Complexity of interpretation envisioning is discussed, taking the given results as an example.     

Dynamic model of the left ventricle for use in simulation of myocardial perfusion SPECT and gated SPECT

Simulation is a useful tool in cardiac SPECT to assess quantification algorithms. However, simple equation-based models are limited in their ability to simulate realistic heart motion and perfusion. We present a numerical dynamic model of the left ventricle, which allows us to simulate normal and anomalous cardiac cycles, as well as perfusion defects. Bicubic splines were fitted to a number of control points to represent endocardial and epicardial surfaces of the left ventricle. A transformation from each point on the surface to a template of activity was made to represent the myocardial perfusion. Geometry-based and patient-based simulations were performed to illustrate this model. Geometry-based simulations modeled (1) a normal patient, (2) a well-perfused patient with abnormal regional function, (3) an ischaemic patient with abnormal regional function, and (4) a patient study including tracer kinetics. Patient-based simulation consisted of a left ventricle including a realistic shape and motion obtained from a magnetic resonance study. We conclude that this model has the potential to study the influence of several physical parameters and the left ventricle contraction in myocardial perfusion SPECT and gated-SPECT studies.     

4D MR imaging of respiratory organ motion and its variability

This paper describes a method for 4D imaging, which is used to study respiratory organ motion, a key problem in various treatments. Whilst the commonly used imaging methods rely on simplified breathing patterns to acquire one breathing cycle, the proposed method was developed to study irregularities in organ motion during free breathing over tens of minutes. The method does not assume a constant breathing depth or even strict periodicity and does not depend on an external respiratory signal. Time-resolved 3D image sequences were reconstructed by retrospective stacking of dynamic 2D images using internal image-based sorting. The generic method is demonstrated for the liver and for the lung. Quantitative evaluations of the volume consistency show the advantages over one-dimensional measurements for image sorting. Dense deformation fields describing the respiratory motion were estimated from the reconstructed volumes using non-rigid 3D registration. All obtained motion fields showed variations in the range of minutes such as drifts and deformations, which changed both the exhalation position of the liver and the breathing pattern. The obtained motion data are used in proton therapy planning to evaluate dose delivery methodologies with respect to their motion sensitivity. Besides this application, the new possibilities of studying respiratory motion are valuable for other applications such as the evaluation of gating techniques with respect to residual motion.     

Non-invasive magnetic resonance imaging assessment of myocardial changes and the effects of angiotensin-converting enzyme inhibition in diabetic rats

A non-invasive cine magnetic resonance imaging (MRI) technique was developed to allow, for the first time, detection and characterization of chronic changes in myocardial tissue volume and the effects upon these of treatment by the angiotensin-converting enzyme (ACE) inhibitor captopril in streptozotocin (STZ)-diabetic male Wistar rats. Animals that had been made diabetic at the ages of 7, 10 and 13 weeks and a captopril-treated group of animals made diabetic at the age of 7 weeks were scanned. The findings were compared with the results from age-matched controls. All animal groups (n = 4 animals in each) were consistently scanned at 16 weeks. Left and right ventricular myocardial volumes were reconstructed from complete data sets of left and right ventricular transverse sections which covered systole and most of diastole using twelve equally incremented time points through the cardiac cycle. The calculated volumes remained consistent through all twelve time points of the cardiac cycle in all five experimental groups and agreed with the corresponding post-mortem determinations. These gave consistent myocardial densities whose values could additionally be corroborated by previous reports, confirming the validity of the quantitative MRI results and analysis. The myocardial volumes were conserved in animals whose diabetes was induced at 13 weeks but were significantly increased relative to body weight in animals made diabetic at 7 and 10 weeks. Captopril treatment, which was started immediately after induction of diabetes, prevented the development of this relative hypertrophy in both the left and right ventricles. We have thus introduced and validated quantitative MRI methods in a demonstration, for the first time, of chronic myocardial changes in both the right and left ventricles of STZ-diabetic rats and their prevention by the ACE inhibitor captopril.     

Preservation of cardiac myocytes subjected to different preconditions: A comparative morphometric study of beating,fibrillating, and cardioplegically arrested canine hearts

This study compares the ultrastructure of beating canine hearts with that of hearts subjected to different clinically common forms of cardiac arrest. The contraction state per test field was ascertained according to a specially developed classification. The volume density of myofibrils and the surface to volume ratio of mitochondria were used as parameters for cellular and mitochondrial swelling. Contraction bands were not found in any of the differently pretreated hearts. Following immersion fixation, contractions as well as over- and hypercontractions in beating, fibrillating, and St. Thomas-arrested hearts are significantly more pronounced than in HTK-arrested hearts. Cellular and mitochondrial volumes were similar in beating and fibrillating hearts. St. Thomas-perfusion significantly decreased cellular and mitochondrial volume compared to beating hearts, but these values were in the same range as in fibrillating hearts. Only HTK-solution actually led to a strong reduction of these compartments. Compared to immersion, perfusion fixation after coronary perfusion with cardioplegic solutions led to comparable cellular volumes, but significantly elevated the percentage of relaxed sarcomeres and significantly reduced mitochondrial swelling. The best structural preservation of myocytes was found after HTK-perfusion and perfusion fixation. Such ultrastructural quantitative and morphometrical parameters are powerful tools since results confirm that the degree of myocardial preservation depends on the method of cardiac arrest. This forms the basis for the choice of preconditions for subsequent ischemia. Furthermore, significant alterations of myocardial ultrastructure depend on a combination of the functional state of the heart, the method of cardioplegia, and the technique of fixation. © 1993 Wiley-Liss, Inc.     

The accuracy and repeatability of an automatic 2D-3D fluoroscopic image-model registration technique for determining shoulder joint kinematics

Fluoroscopic imaging, using single plane or dual plane images, has grown in popularity to measure dynamic in vivo human shoulder joint kinematics. However, no study has quantified the difference in spatial positional accuracy between single and dual plane image-model registration applied to the shoulder joint. In this paper, an automatic 2D-3D image-model registration technique was validated for accuracy and repeatability with single and dual plane fluoroscopic images. Accuracy was assessed in a cadaver model, kinematics found using the automatic registration technique were compared to those found using radiostereometric analysis. The in vivo repeatability of the automatic registration technique was assessed during the dynamic abduction motion of four human subjects. The in vitro data indicated that the error in spatial positional accuracy of the humerus and the scapula was less than 0.30mm in translation and less than 0.58° in rotation using dual plane images. Single plane accuracy was satisfactory for in-plane motion variables, but out-of-plane motion variables on average were approximately 8 times less accurate. The in vivo test indicated that the repeatability of the automatic 2D-3D image-model registration was 0.50mm in translation and 1.04° in rotation using dual images. For a single plane technique, the repeatability was 3.31mm in translation and 2.46° in rotation for measuring shoulder joint kinematics. The data demonstrate that accurate and repeatable shoulder joint kinematics can be obtained using dual plane fluoroscopic images with an automatic 2D-3D image-model registration technique; and that out-of-plane motion variables are less accurate than in-plane motion variables using a single plane technique.     

Analysis and processing of laser Doppler perfusion monitoring signals recorded from the beating heart

Laser Doppler perfusion monitoring (LDPM) can be used for monitoring myocardial perfusion in the non-beating heart. However, the movement of the beating heart generates large artifacts. Therefore the aim of the study was to develop an LDPM system capable of correlating the laser Doppler signals to the cardiac cycle and to process the signals to reduce the movement artifacts. Measurements were performed on three calves, both on the normal beating heart and during occlusion of the left anterior descending coronary artery (LAD). The recorded LDPM signals were digitally processed and correlated to the sampled ECG. Large variations in the output (perfusion) and DC signals during the cardiac cycle were found, with average coefficients of variation of 0.36 and 0.14 (n-14), respectively. However, sections with a relatively low, stable output signal were found in late diastole, where the movement of the heart is at a minimum. Occlusion of the LAD showed the importance of recording the laser Doppler signals at an appropriate point in the cardiac cycle, in this case late systole, to minimise movement artifacts. It is possible to further reduce movement artifacts by increasing the lower cutoff frequency when calculating the output signal.     

Minimizing artifacts resulting from respiratory and cardiac motion by optimization of the transmission scan in cardiac PET/CT

The introduction of positron emission/computed tomography (PET/CT) systems coupled with multidetector CT arrays has greatly increased the amount of clinical information in myocardial perfusion studies. The CT acquisition serves the dual role of providing high spatial anatomical detail and attenuation correction for PET. However, the differences between the interaction of respiratory and cardiac cycles in the CT and PET acquisitions presents a challenge when using the CT to determine PET attenuation correction. Three CT attenuation correction protocols were tested for their ability to produce accurate emission images: gated, a step mode acquisition covering the diastolic heart phase; normal, a high-pitch helical CT; and slow, a low-pitch, low-temporal-resolution helical CT. The amount of cardiac tissue in the emission image that overlaid lung tissue in the transmission image was used as the measure of mismatch between acquisitions. Phantom studies simulating misalignment of the heart between the transmission and emission sequences were used to correlate the amount of mismatch with the artificial defect changes in the emission image. Consecutive patients were studied prospectively with either paired gated (diastolic phase, 120 kVp, 280 mA, 2.6 s) and slow CT (0.562:1 pitch, 120 kVp, Auto-mA, 16 s) or paired normal (0.938:1 pitch, 120 kVp, Auto-mA, 4.8 s) and slow CT protocols, prior to a Rb-82 perfusion study. To determine the amount of mismatch, the transmission and emission images were converted to binary representations of attenuating tissue and cardiac tissue and overlaid using their native registration. The number of cardiac tissue pixels from the emission image present in the CT lung field yielded the magnitude of misalignment represented in terms of volume, of where a small volume indicates better registration. Acquiring a slow CT improved registration between the transmission and emission acquisitions compared to the gated and normal CT protocols. The volume of PET cardiac tissue in the CT lung field was significantly lower (p < 0.03) for the slow CT protocol in both the rest and stress emission studies. Phantom studies showed that an overlaying volume greater than 2.6 mL would produce significant artificial defects as determined by a quantitative software package that employs a normal database. The percentage of patient studies with overlaying volume greater than 2.6 mL was reduced from 71% with the normal CT protocol to 28% with the slow CT protocol. The remaining 28% exhibited artifacts consistent with heart drift and patient motion that could not be corrected by adjusting the CT acquisition protocol. The low pitch of the slow CT protocol provided the best match to the emission study and is recommended for attenuation correction in cardiac PET/CT studies. Further reduction in artifacts arising from cardiac drift is required and warrants an image registration solution.     

A method for the reconstruction of four-dimensional synchronized CT scans acquired during free breathing

Breathing motion is a significant source of error in radiotherapy treatment planning for the thorax and upper abdomen. Accounting for breathing motion has a profound effect on the size of conformal radiation portals employed in these sites. Breathing motion also causes artifacts and distortions in treatment planning computed tomography (CT) scans acquired during free breathing and also causes a breakdown of the assumption of the superposition of radiation portals in intensity-modulated radiation therapy, possibly leading to significant dose delivery errors. Proposed voluntary and involuntary breath-hold techniques have the potential for reducing or eliminating the effects of breathing motion, however, they are limited in practice, by the fact that many lung cancer patients cannot tolerate holding their breath. We present an alternative solution to accounting for breathing motion in radiotherapy treatment planning, where multislice CT scans are collected simultaneously with digital spirometry over many free breathing cycles to create a four-dimensional (4-D) image set, where tidal lung volume is the additional dimension. An analysis of this 4-D data leads to methods for digital-spirometry, based elimination or accounting of breathing motion artifacts in radiotherapy treatment planning for free breathing patients. The 4-D image set is generated by sorting free-breathing multislice CT scans according to user-defined tidal-volume bins. A multislice CT scanner is operated in the ciné mode, acquiring 15 scans per couch position, while the patient undergoes simultaneous digital-spirometry measurements. The spirometry is used to retrospectively sort the CT scans by their correlated tidal lung volume within the patient's normal breathing cycle. This method has been prototyped using data from three lung cancer patients. The actual tidal lung volumes agreed with the specified bin volumes within standard deviations ranging between 22 and 33 cm3. An analysis of sagittal and coronal images demonstrated relatively small (<1 cm) motion artifacts along the diaphragm, even for tidal volumes where the rate of breathing motion is greatest. While still under development, this technology has the potential for revolutionizing the radiotherapy treatment planning for the thorax and upper abdomen.     

Registration-based segmentation of murine 4D cardiac micro-CT data using symmetric normalization

Micro-CT can play an important role in preclinical studies of cardiovascular disease because of its high spatial and temporal resolution. Quantitative analysis of 4D cardiac images requires segmentation of the cardiac chambers at each time point, an extremely time consuming process if done manually. To improve throughput this study proposes a pipeline for registration-based segmentation and functional analysis of 4D cardiac micro-CT data in the mouse. Following optimization and validation using simulations, the pipeline was applied to in vivo cardiac micro-CT data corresponding to ten cardiac phases acquired in C57BL/6 mice (n = 5). After edge-preserving smoothing with a novel adaptation of 4D bilateral filtration, one phase within each cardiac sequence was manually segmented. Deformable registration was used to propagate these labels to all other cardiac phases for segmentation. The volumes of each cardiac chamber were calculated and used to derive stroke volume, ejection fraction, cardiac output, and cardiac index. Dice coefficients and volume accuracies were used to compare manual segmentations of two additional phases with their corresponding propagated labels. Both measures were, on average, >0.90 for the left ventricle and >0.80 for the myocardium, the right ventricle, and the right atrium, consistent with trends in inter- and intra-segmenter variability. Segmentation of the left atrium was less reliable. On average, the functional metrics of interest were underestimated by 6.76% or more due to systematic label propagation errors around atrioventricular valves; however, execution of the pipeline was 80% faster than performing analogous manual segmentation of each phase.