Comparison of breathing gated CT images generated using a 5DCT technique and a commercial clinical protocol in a porcine model

A Retrospective Analysis of the First Clinical 5DCT Workflow

BACKGROUND/OBJECTIVES

5DCT was first proposed in 2005 as a motion-compensated CT simulation approach for radiotherapy treatment planning to avoid sorting artifacts that arise in 4DCT when patients breathe irregularly. Since March 2019, 5DCT has been clinically implemented for routine use at our institution to leverage this technological advantage. The clinical workflow includes a quality assurance report that describes the output of primary workflow steps. This study reports on the challenges and quality of the clinical 5DCT workflow using these quality assurance reports.

METHODS

We evaluated all thoracic 5DCT simulation datasets consecutively acquired at our institution between March 2019 and December 2022 for thoracic radiotherapy treatment planning. The 5DCT datasets utilized motion models constructed from 25 fast-helical free-breathing computed tomography (FHFBCTs) with simultaneous respiratory bellows signal monitoring to reconstruct individual, user-specified breathing-phase images (termed 5DCT phase images) for internal target volume contouring. Each 5DCT dataset was accompanied by a structured quality assurance report composed of qualitative and quantitative measures of the breathing pattern, image quality, DIR quality, model fitting accuracy, and a validation process by which the original FHFBCT scans were regenerated with the 5DCT model. Measures of breathing irregularity, image quality, and DIR quality were retrospectively categorized on a grading scale from 1 (regular breathing and accurate registration/modeling) to 4 (irregular breathing and inaccurate registration/modeling). The validation process was graded according to the same scale, and this grade was termed the suitability-for-treatment-planning (STP) grade. We correlated the graded variables to the STP grade. In addition to the quality assurance reports, we reviewed the contour sessions to determine how often 5DCT phase images were used for treatment planning and delivery.

RESULTS

There were 169 5DCT simulation datasets available from 156 patients for analysis. The STP was moderately correlated with breathing irregularity, image quality, and DIR quality (Spearman coefficients: 0.26, 0.30, and 0.50, respectively). Multiple linear regression analysis demonstrated that STP was correlated with regular breathing patterns (p = 0.008), image quality (p < 0.001), and better DIR quality (p < 0.001). 5DCT datasets were used for treatment planning in 82% of cases, while in 12% of cases, a backup image process was used. In total, 6% of image datasets were not used for treatment planning due to factors unrelated to the 5DCT workflow quality.

CONCLUSIONS

The strongest association with STP was with DIR quality grades, as indicated by both Spearman and multiple linear regression analysis, implying that improvements to DIR accuracy and evaluation may be the best route for further improvement to 5DCT. The high rate of 5DCT phase image use for treatment planning showed that the workflow was reliable, and this has encouraged us to continue to develop and improve the workflow steps.

An analysis of the regional heterogeneity in tissue elasticity in lung cancer patients with COPD

Purpose

Recent advancements in obtaining image-based biomarkers from CT images have enabled lung function characterization, which could aid in lung interventional planning. However, the regional heterogeneity in these biomarkers has not been well documented, yet it is critical to several procedures for lung cancer and COPD. The purpose of this paper is to analyze the interlobar and intralobar heterogeneity of tissue elasticity and study their relationship with COPD severity.

Methods

We retrospectively analyzed a set of 23 lung cancer patients for this study, 14 of whom had COPD. For each patient, we employed a 5DCT scanning protocol to obtain end-exhalation and end-inhalation images and semi-automatically segmented the lobes. We calculated tissue elasticity using a biomechanical property estimation model. To obtain a measure of lobar elasticity, we calculated the mean of the voxel-wise elasticity values within each lobe. To analyze interlobar heterogeneity, we defined an index that represented the properties of the least elastic lobe as compared to the rest of the lobes, termed the Elasticity Heterogeneity Index (EHI). An index of 0 indicated total homogeneity, and higher indices indicated higher heterogeneity. Additionally, we measured intralobar heterogeneity by calculating the coefficient of variation of elasticity within each lobe.

Results

The mean EHI was 0.223 ± 0.183. The mean coefficient of variation of the elasticity distributions was 51.1% ± 16.6%. For mild COPD patients, the interlobar heterogeneity was low compared to the other categories. For moderate-to-severe COPD patients, the interlobar and intralobar heterogeneities were highest, showing significant differences from the other groups.

Conclusion

We observed a high level of lung tissue heterogeneity to occur between and within the lobes in all COPD severity cases, especially in moderate-to-severe cases. Heterogeneity results demonstrate the value of a regional, function-guided approach like elasticity for procedures such as surgical decision making and treatment planning.

Scalable quorum-based deep neural networks with adversarial learning for automated lung lobe segmentation in fast helical free-breathing CTs

PURPOSE

Fast helical free-breathing CT (FHFBCT) scans are widely used for 5DCT and 5D Cone Beam imaging protocols. For quantitative analysis of lung physiology and function, it is important to segment the lung lobes in these scans. Since the 5DCT protocols use up to 25 FHFBCT scans, it is important that this segmentation task be automated. In this paper, we present a deep neural network (DNN) framework for segmenting the lung lobes in near real time.

METHODS

A total of 22 patient datasets (550 3D CT scans) were used for the study. Each of the lung lobes was manually segmented and considered ground-truth. A supervised and constrained generative adversarial network (CGAN) was employed for learning each set of lobe segmentations for each patient with 12 patients designated for training data. The resulting generator DNNs represented the lobe segmentations for each training dataset. A quorum-based algorithm was then implemented to test validation data consisting of 10 separate patient datasets (250 3D CTs). Each of the DNNs predicted their corresponding lobes for the validation data, and equal weights were given to the 12 generator CGANs. The quorum process worked by selecting the weighted average result of all 12 CGAN results for each lobe.

RESULTS

When evaluated against ground-truth segmentations, the quorum-based lobe segmentation was observed to have average structural similarity index, normalized cross-correlation coefficient, and dice coefficient values of 0.929, 0.806, and 0.814, respectively, compared to values of 0.911, 0.698, and 0.696, respectively, using a conventional strategy.

CONCLUSION

The proposed quorum-based approach computed segmentations with clinically acceptable accuracy in near real time using a multi-GPU-based computing setup. This method is scalable as more patient-specific CGANs can be added to the quorum over time.

An adversarial machine learning framework and biomechanical model-guided approach for computing 3D lung tissue elasticity from end-expiration 3DCT

PURPOSE

Lung elastography aims at measuring the lung parenchymal tissue elasticity for applications ranging from diagnostic purposes to biomechanically guided deformations. Characterizing the lung tissue elasticity requires four-dimensional (4D) lung motion as an input, which is currently estimated by deformably registering 4D computed tomography (4DCT) datasets. Since 4DCT imaging is widely used only in a radiotherapy treatment setup, there is a need to predict the elasticity distribution in the absence of 4D imaging for applications within and outside of radiotherapy domain.

METHODS

In this paper, we present a machine learning-based method that predicts the three-dimensional (3D) lung tissue elasticity distribution for a given end-expiration 3DCT. The method to predict the lung tissue elasticity from an end-expiration 3DCT employed a deep neural network that predicts the tissue elasticity for the given CT dataset. For training and validation purposes, we employed five-dimensional CT (5DCT) datasets and a finite element biomechanical lung model. The 5DCT model was first used to generate end-expiration lung geometry, which was taken as the source lung geometry for biomechanical modeling. The deformation vector field pointing from end expiration to end inhalation was computed from the 5DCT model and taken as input in order to solve for the lung tissue elasticity. An inverse elasticity estimation process was employed, where we iteratively solved for the lung elasticity distribution until the model reproduced the ground-truth deformation vector field. The machine learning process uses a specific type of learning process, namely a constrained generalized adversarial neural network (cGAN) that learned the lung tissue elasticity in a supervised manner. The biomechanically estimated tissue elasticity together with the end-exhalation CT was the input for the supervised learning. The trained cGAN generated the elasticity from a given breath-hold CT image. The elasticity estimated was validated in two approaches. In the first approach, a L2-norm-based direct comparison was employed between the estimated elasticity and the ground-truth elasticity. In the second approach, we generated a synthetic four-dimensional CT (4DCT0 using a lung biomechanical model and the estimated elasticity and compared the deformations with the ground-truth 4D deformations using three image similarity metrics: mutual Information (MI), structured similarity index (SSIM), and normalized cross correlation (NCC).

RESULTS

The results show that a cGAN-based machine learning approach was effective in computing the lung tissue elasticity given the end-expiration CT datasets. For the training data set, we obtained a learning accuracy of 0.44 ± 0.2 KPa. For the validation dataset, consisting of 13 4D datasets, we were able to obtain an accuracy of 0.87 ± 0.4 KPa. These results show that the cGAN-generated elasticity correlates well with that of the underlying ground-truth elasticity. We then integrated the estimated elasticity with the biomechanical model and applied the same boundary conditions in order to generate the end inhalation CT. The cGAN-generated images were very similar to that of the original end inhalation CT. The average value of the MI is 1.77 indicating the high local symmetricity between the ground truth and the cGAN elasticity-generated end inhalation CT data. The average value of the structural similarity for the 13 patients was observed to be 0.89 indicating the high structural integrity of the cGAN elasticity-generated end inhalation CT. Finally, the average NCC value of 0.97 indicates that potential variations in the contrast and brightness of the cGAN elasticity-generated end inhalation CT and the ground-truth end inhalation CT.

CONCLUSION

The cGAN-generated lung tissue elasticity given an end-expiration CT image can be computed in near real time. Using the lung tissue elasticity along with a biomechanical model, 4D lung deformations can be generated from a given end-expiration CT image within clinically acceptable numerical accuracy.

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.

Safety-oriented design of in-house software for new techniques: A case study using a model-based 4DCT protocol

PURPOSE

In-house software is commonly employed to implement new imaging and therapy techniques before commercial solutions are available. Risk analysis methods, as detailed in the TG-100 report of the American Association of Physicists in Medicine, provide a framework for quality management of processes but offer little guidance on software design. In this work, we examine a novel model-based four-dimensional computed tomography (4DCT) protocol using the TG-100 approach and describe two additional methods for promoting safety of the associated in-house software.

METHODS

To implement a previously published model-based 4DCT protocol, in-house software was necessary for tasks such as synchronizing a respiratory signal to computed tomography images, deformable image registration (DIR), model parameter fitting, and interfacing with a treatment planning system. A process map was generated detailing the workflow. Failure modes and effects analysis (FMEA) was performed to identify critical steps and guide quality interventions. Software system safety was addressed through writing "use cases," narratives that characterize the behavior of the software, for all major operations to elicit safety requirements. Safety requirements were codified using the easy approach to requirements syntax (EARS) to ensure testability and eliminate ambiguity.

RESULTS

Sixty-one failure modes were identified and assigned risk priority numbers using FMEA. Resultant quality management interventions include integration of a comprehensive reporting and logging system into the software, mandating daily and monthly equipment quality assurance procedures, and a checklist to be completed at image acquisition. Use cases and resulting safety requirements informed the design of needed in-house software as well as a suite of tests performed during the image generation process.

CONCLUSIONS

TG-100 methods were used to construct a process-level quality management program for a 4DCT imaging protocol. Two supplemental tools from the field of requirements engineering facilitated elicitation and codification of safety requirements that informed the design and testing of in-house software necessary to implement the protocol. These general tools can be applied to promote safety when in-house software is needed to bring new techniques to the clinic.

Multislice motion modeling for MRI-guided radiotherapy gating

PURPOSE

On-board magnetic resonance imaging (MRI) greatly enhances real-time target tracking capability during radiotherapy treatments. However, multislice and volumetric MRI techniques are frame rate limited and introduce unacceptable latency between the target moving out of position and the beam being turned off. We present a technique to estimate continuous volumetric tissue motion using motion models built from a repeated acquisition of a stack of MR slices. Applications including multislice target visualization and out-of-slice motion estimation during MRI-guided radiotherapy are demonstrated.

METHODS

Eight healthy volunteer studies were performed using a 0.35 T MRI-guided radiotherapy system. Images were acquired at three frames per second in an interleaved fashion across ten adjacent sagittal slice positions covering 4.5 cm using a balanced steady-state-free precession sequence. A previously published five-dimensional (5D) linear motion model used for MRI-guided radiotherapy gating was extended to include multiple slices. This model utilizes an external respiratory bellows signal recorded during imaging to simultaneously estimate motion across all imaged slices. For comparison to an image-based approach, the manifold learning technique local linear embedding (LLE) was used to derive a respiratory surrogate for motion modeling. Manifolds for every slice were aligned during LLE in a group-wise fashion, enabling motion estimation outside the current imaged slice using a motion model, a process which we denote as mSGA. Additionally, a method is developed to evaluate out-of-slice motion estimates. The multislice motion model was evaluated in a single slice with each newly acquired image using a leave-one-out approach. Model-generated gating decision accuracy and beam-on positive predictive value (PPV) are reported along with the median and 95th percentile distance between model and ground truth target centroids.

RESULTS

The average model gating decision accuracy and PPV across all volunteer studies was 93.7% and 92.8% using the 5D model, and 96.8% and 96.1% using the mSGA model, respectively. The median and 95th percentile distance between model and ground truth target centroids was 0.91 and 2.90 mm, respectively, using the 5D model and 0.58 and 1.49 mm using the mSGA model, averaged over all eight subjects. The mSGA motion model provided a statistically significant improvement across all evaluation metrics compared to the external surrogate-based 5D model.

CONCLUSION

The proposed techniques for out-of-slice target motion estimation demonstrated accuracy likely sufficient for clinical use. Results indicate the mSGA model may provide higher accuracy, however, the external surrogate-based model allows for unbiased in vivo accuracy evaluation.

Model-Interpolated Gating for Magnetic Resonance Image-Guided Radiation Therapy

PURPOSE

To develop and validate a technique for radiation therapy gating using slow (≤1 frame per second) magnetic resonance imaging (MRI) and a motion model. Proposed uses of the technique include radiation therapy gating using T2-weighted images and conducting additional imaging studies during gated treatments.

METHODS AND MATERIALS

The technique uses a physiologically guided breathing motion model to interpolate deformed target position between 2-dimensional (2D) MRI images acquired every 1 to 3 seconds. The model is parameterized by a 1-dimensional respiratory bellows surrogate and is continuously updated with the most recently acquired 2D images. A phantom and 8 volunteers were imaged with a 0.35T MRI-guided radiation therapy system. A balanced steady-state free precession sequence with a 2D frame rate of 3 frames per second was used to evaluate the technique. The accuracy and beam-on positive predictive value (PPV) of the model-based gating decisions were evaluated using the gating decisions derived from imaging as a ground truth. A T2-weighted gating offline proof-of-concept study using a half-Fourier, single-shot, turbo-spin echo sequence is reported.

RESULTS

Model-interpolated gating accuracy, beam-on PPV, and median absolute distances between model and image-tracked target centroids were, on average, 98.3%, 98.4%, and 0.33 mm, respectively, in the balanced steady-state free precession phantom studies and 93.7%, 92.1%, and 0.86 mm, respectively, in the volunteer studies. T2 model-interpolated gating in 6 volunteers yielded an average accuracy and PPV of 94.3% and 92.5%, respectively, and the mean absolute median distance between modeled and imaged target centroids was 0.86 mm.

CONCLUSIONS

This work demonstrates the concept of model-interpolated gating for MRI-guided radiation therapy. The technique was found to be potentially sufficiently accurate for clinical use. Further development is needed to accommodate out-of-plane motion and the use of an internal MR-based respiratory surrogate.

Dependence of subject-specific parameters for a fast helical CT respiratory motion model on breathing rate: an animal study

To determine if the parameters relating lung tissue displacement to a breathing surrogate signal in a previously published respiratory motion model vary with the rate of breathing during image acquisition. An anesthetized pig was imaged using multiple fast helical scans to sample the breathing cycle with simultaneous surrogate monitoring. Three datasets were collected while the animal was mechanically ventilated with different respiratory rates: 12 bpm (breaths per minute), 17 bpm, and 24 bpm. Three sets of motion model parameters describing the correspondences between surrogate signals and tissue displacements were determined. The model error was calculated individually for each dataset, as well asfor pairs of parameters and surrogate signals from different experiments. The values of one model parameter, a vector field denoted [Formula: see text] which related tissue displacement to surrogate amplitude, determined for each experiment were compared. The mean model error of the three datasets was 1.00  ±  0.36 mm with a 95th percentile value of 1.69 mm. The mean error computed from all combinations of parameters and surrogate signals from different datasets was 1.14  ±  0.42 mm with a 95th percentile of 1.95 mm. The mean difference in [Formula: see text] over all pairs of experiments was 4.7%  ±  5.4%, and the 95th percentile was 16.8%. The mean angle between pairs of [Formula: see text] was 5.0  ±  4.0 degrees, with a 95th percentile of 13.2 mm. The motion model parameters were largely unaffected by changes in the breathing rate during image acquisition. The mean error associated with mismatched sets of parameters and surrogate signals was 0.14 mm greater than the error achieved when using parameters and surrogate signals acquired with the same breathing rate, while maximum respiratory motion was 23.23 mm on average.

A prospective gating method to acquire a diverse set of free-breathing CT images for model-based 4DCT

Breathing motion modeling requires observation of tissues at sufficiently distinct respiratory states for proper 4D characterization. This work proposes a method to improve sampling of the breathing cycle with limited imaging dose. We designed and tested a prospective free-breathing acquisition protocol with a simulation using datasets from five patients imaged with a model-based 4DCT technique. Each dataset contained 25 free-breathing fast helical CT scans with simultaneous breathing surrogate measurements. Tissue displacements were measured using deformable image registration. A correspondence model related tissue displacement to the surrogate. Model residual was computed by comparing predicted displacements to image registration results. To determine a stopping criteria for the prospective protocol, i.e. when the breathing cycle had been sufficiently sampled, subsets of N scans where 5  ⩽  N  ⩽  9 were used to fit reduced models for each patient. A previously published metric was employed to describe the phase coverage, or 'spread', of the respiratory trajectories of each subset. Minimum phase coverage necessary to achieve mean model residual within 0.5 mm of the full 25-scan model was determined and used as the stopping criteria. Using the patient breathing traces, a prospective acquisition protocol was simulated. In all patients, phase coverage greater than the threshold necessary for model accuracy within 0.5 mm of the 25 scan model was achieved in six or fewer scans. The prospectively selected respiratory trajectories ranked in the (97.5  ±  4.2)th percentile among subsets of the originally sampled scans on average. Simulation results suggest that the proposed prospective method provides an effective means to sample the breathing cycle with limited free-breathing scans. One application of the method is to reduce the imaging dose of a previously published model-based 4DCT protocol to 25% of its original value while achieving mean model residual within 0.5 mm.

Comparison of lung tumor motion measured using a model-based 4DCT technique and a commercial protocol

PURPOSE

To compare lung tumor motion measured with a model-based technique to commercial 4-dimensional computed tomography (4DCT) scans and describe a workflow for using model-based 4DCT as a clinical simulation protocol.

METHODS AND MATERIALS

Twenty patients were imaged using a model-based technique and commercial 4DCT. Tumor motion was measured on each commercial 4DCT dataset and was calculated on model-based datasets for 3 breathing amplitude percentile intervals: 5th to 85th, 5th to 95th, and 0th to 100th. Internal target volumes (ITVs) were defined on the 4DCT and 5th to 85th interval datasets and compared using Dice similarity. Images were evaluated for noise and rated by 2 radiation oncologists for artifacts.

RESULTS

Mean differences in tumor motion magnitude between commercial and model-based images were 0.47 ± 3.0, 1.63 ± 3.17, and 5.16 ± 4.90 mm for the 5th to 85th, 5th to 95th, and 0th to 100th amplitude intervals, respectively. Dice coefficients between ITVs defined on commercial and 5th to 85th model-based images had a mean value of 0.77 ± 0.09. Single standard deviation image noise was 11.6 ± 9.6 HU in the liver and 6.8 ± 4.7 HU in the aorta for the model-based images compared with 57.7 ± 30 and 33.7 ± 15.4 for commercial 4DCT. Mean model error within the ITV regions was 1.71 ± 0.81 mm. Model-based images exhibited reduced presence of artifacts at the tumor compared with commercial images.

CONCLUSION

Tumor motion measured with the model-based technique using the 5th to 85th percentile breathing amplitude interval corresponded more closely to commercial 4DCT than the 5th to 95th or 0th to 100th intervals, which showed greater motion on average. The model-based technique tended to display increased tumor motion when breathing amplitude intervals wider than 5th to 85th were used because of the influence of unusually deep inhalations. These results suggest that care must be taken in selecting the appropriate interval during image generation when using model-based 4DCT methods.

A 4DCT imaging-based breathing lung model with relative hysteresis

To reproduce realistic airway motion and airflow, the authors developed a deforming lung computational fluid dynamics (CFD) model based on four-dimensional (4D, space and time) dynamic computed tomography (CT) images. A total of 13 time points within controlled tidal volume respiration were used to account for realistic and irregular lung motion in human volunteers. Because of the irregular motion of 4DCT-based airways, we identified an optimal interpolation method for airway surface deformation during respiration, and implemented a computational solid mechanics-based moving mesh algorithm to produce smooth deforming airway mesh. In addition, we developed physiologically realistic airflow boundary conditions for both models based on multiple images and a single image. Furthermore, we examined simplified models based on one or two dynamic or static images. By comparing these simplified models with the model based on 13 dynamic images, we investigated the effects of relative hysteresis of lung structure with respect to lung volume, lung deformation, and imaging methods, i.e., dynamic vs. static scans, on CFD-predicted pressure drop. The effect of imaging method on pressure drop was 24 percentage points due to the differences in airflow distribution and airway geometry.

A 4DCT imaging-based breathing lung model with relative hysteresis

To reproduce realistic airway motion and airflow, the authors developed a deforming lung computational fluid dynamics (CFD) model based on four-dimensional (4D, space and time) dynamic computed tomography (CT) images. A total of 13 time points within controlled tidal volume respiration were used to account for realistic and irregular lung motion in human volunteers. Because of the irregular motion of 4DCT-based airways, we identified an optimal interpolation method for airway surface deformation during respiration, and implemented a computational solid mechanics-based moving mesh algorithm to produce smooth deforming airway mesh. In addition, we developed physiologically realistic airflow boundary conditions for both models based on multiple images and a single image. Furthermore, we examined simplified models based on one or two dynamic or static images. By comparing these simplified models with the model based on 13 dynamic images, we investigated the effects of relative hysteresis of lung structure with respect to lung volume, lung deformation, and imaging methods, i.e., dynamic vs. static scans, on CFD-predicted pressure drop. The effect of imaging method on pressure drop was 24 percentage points due to the differences in airflow distribution and airway geometry.

Computed Tomography and Magnetic Resonance Imaging for Longitudinal Characterization of Lung Structure Changes in a Yucatan Miniature Pig Silicosis Model

Medical imaging is a rapidly advancing field enabling the repeated, noninvasive assessment of physiological structure and function. These beneficial characteristics can supplement studies in swine by mirroring the clinical functions of detection, diagnosis, and monitoring in humans. In addition, swine may serve as a human surrogate, facilitating the development and comparison of new imaging protocols for translation to humans. This study presents methods for pulmonary imaging developed for monitoring pulmonary disease initiation and progression in a pig exposure model with computed tomography and magnetic resonance imaging. In particular, a focus was placed on systematic processes, including positioning, image acquisition, and structured reporting to monitor longitudinal change. The image-based monitoring procedure was applied to 6 Yucatan miniature pigs. A subset of animals (n= 3) were injected with crystalline silica into the apical bronchial tree to induce silicosis. The methodology provided longitudinal monitoring and evidence of progressive lung disease while simultaneously allowing for a cross-modality comparative study highlighting the practical application of medical image data collection in swine. The integration of multimodality imaging with structured reporting allows for cross comparison of modalities, refinement of CT and MRI protocols, and consistently monitors potential areas of interest for guided biopsy and/or necropsy.