Search NASA⌕ Search

SEARCH · Search NASA

Results for “Image Processing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-J-40: Evaluation of Sensitivity of the Automatic Matching Between Cone-Beam CT Image and Simulation CT Image in TrueBeam 2.0 Imaging System 6DoF Considering Different Uncertainty Sources

Purpose: To estimate the sensitivity of TrueBeam 2.0 Imaging System 6DoF automatic matching tool through the acquisition of cone-beam CT images in different phantoms applying submillimeter translations and rotations of tenths of a degree and registered with image simulation CT. Methods: To evaluate overall system-wide image, we consider two uncertainties source; First, the uncertainty of the manual phantom displacement (ε-m). This uncertainty is calculated by a digital caliper (0.01 mm) for vertical (Vrt), lateral (Lat) and longitudinal (Lng). A digital inclinometer (0.01°) for the pitch and roll and the own phantom scale to evaluate the coordinate rotation (Rtn). The second uncertainty is the displacement detected by the algorithm system of matching (σ-d) that we obtain from the standard deviations of the different measurements. We use three different phantoms. The BrainLab Radiosurgery system for supporting masks with an anthropomorphic dummy adapted to allow displacements of 0.1 mm in Vrt, Lat and Lng dimensions and rotations of 0.1° in Pitch dimension. For the analysis of the Rtn and Roll dimensions we use two homemade phantoms (RinoRot and RinoRoll, La Fe Hospital, Valencia, Spain) that allow rotations of 0.3°. Results: In the case of manual displacement of 0.10 ± 0.03 mm in the translations, the system detect 0.10 ± 0.07 mm, 0.12 ± 0.07 mm and 0.13 ± 0.07 mm (mean ± SD) in Lat, Vrt and Lng respectively. In the case of rotational dimension, manual displacement of 0.3 ± 0.1° was detected with 0.19 ± 0.06°, 0.29 ± 0.03° and 0.27 ± 0.06° in Pitch, Roll and Rtn. Conclusion: We conclude that the sensitivity of the automatic matching system is within 0.10 mm in translations and 0.3° in rotations. These values are under the own sensitivity of the software.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-58: Evaluation of RayStation Hybrid Deformable Image Registration for Accurate Contour Propagation in Adaptive Planning

Purpose: Adaptive radiotherapy requires complete new sets of regions of interests (ROIs) delineation on the mid-treatment CT images. This work aims at evaluating the accuracy of the RayStation hybrid deformable image registration (DIR) algorithm for its overall integrity and accuracy in contour propagation for adaptive planning. Methods: The hybrid DIR is based on the combination of intensity-based algorithm and anatomical information provided by contours. Patients who received mid-treatment CT scans were identified for the study, including six lung patients (two mid-treatment CTs) and six head-and-neck (HN) patients (one mid-treatment CT). DIRpropagated ROIs were compared with physician-drawn ROIs for 8 ITVs and 7 critical organs (lungs, heart, esophagus, and etc.) for the lung patients, as well as 14 GTVs and 20 critical organs (mandible, eyes, parotids, and etc.) for the HN patients. Volume difference, center of mass (COM) difference, and Dice index were used for evaluation. Clinical-relevance of each propagated ROI was scored by two physicians, and correlated with the Dice index. Results: For critical organs, good agreement (Dice>0.9) were seen on all 7 for lung patients and 13 out of 20 for HN patients, with the rest requiring minimal edits. For targets, COM differences were within 5 mm on average for all patients. For Lung, 5 out of 8 ITVs required minimal edits (Dice 0.8–0.9), with the rest 2 needed re-drawn due to their small volumes (<10 cc). However, the propagated HN GTVs resulted in relatively low Dice values (0.5–0.8) due to their small volumes (3–40 cc) and high variability, among which 2 required re-drawn due to new nodal target identified on the mid-treatment CT scans. Conclusion: The hybrid DIR algorithm was found to be clinically useful and efficient for lung and HN patients, especially for propagated critical organ ROIs. It has potential to significantly improve the workflow in adaptive planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-76: Evaluation of the Performance of Different Deformable Image Registration Algorithms in Helical, Axial and Cone-Beam CT Images of a Mobile Phantom

Purpose: To investigate quantitatively the performance of different deformable-image-registration algorithms (DIR) with helical (HCT), axial (ACT) and cone-beam CT (CBCT) by evaluating the variations in the CT-numbers and lengths of targets moving with controlled motion-patterns. Methods: Four DIR-algorithms including demons, fast-demons, Horn-Schunk and Locas-Kanade from the DIRART-software are used to register CT-images of a mobile-phantom. A mobile-phantom is scanned with different imaging techniques that include helical, axial and cone-beam CT. The phantom includes three targets with different lengths that are made from water-equivalent material and inserted in low-density-foam which is moved with adjustable motion-amplitudes and frequencies. Results: Most of the DIR-algorithms are able to produce the lengths of the stationary-targets, however, they do not produce the CT-number values in CBCT. The image-artifacts induced by motion are more regular in CBCT imaging where the mobile-target elongation increases linearly with motion-amplitude. In ACT and HCT, the motion-artifacts are irregular where some mobile -targets are elongated or shrunk depending on the motion-phase during imaging. The DIR-algorithms are successful in deforming the images of the mobile-targets to the images of the stationary-targets producing the CT-number values and length of the target for motion-amplitudes < 20 mm. Similarly in ACT, all DIR-algorithms produced the actual CT-number and length of the stationary-targets for motion-amplitudes < 15 mm. As stronger motion-artifacts are induced in HCT and ACT, DIR-algorithms fail to produce CT-values and shape of the stationary-targets and fast-demons-algorithm has worst performance. Conclusion: Most of DIR-algorithms produce the CT-number values and lengths of the stationary-targets in HCT and ACT images that has motion-artifacts induced by small motion-amplitudes. As motion-amplitudes increase, the DIR-algorithms fail to deform mobile-target images to the stationary-images in HCT and ACT. In CBCT, DIR-algorithms are successful in producing length and shape of the stationary-targets, however, they fail to produce the accurate CT-number level.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-77: Variations in the Displacement Vector Fields Calculated by Different Deformable Image Registration Algorithms Used in Helical, Axial and Cone-Beam CT Images of a Mobile

Purpose: To investigate quantitatively the displacement-vector-fields (DVF) obtained from different deformable image registration algorithms (DIR) in helical (HCT), axial (ACT) and cone-beam CT (CBCT) to register CT images of a mobile phantom and its correlation with motion amplitudes and frequencies. Methods: HCT, ACT and CBCT are used to image a mobile phantom which includes three targets with different sizes that are manufactured from water-equivalent material and embedded in low density foam. The phantom is moved with controlled motion patterns where a range of motion amplitudes (0–40mm) and frequencies (0.125–0.5Hz) are used. The CT images obtained from scanning of the mobile phantom are registered with the stationary CT-images using four deformable image registration algorithms including demons, fast-demons, Horn-Schunk and Locas-Kanade from DIRART software. Results: The DVF calculated by the different algorithms correlate well with the motion amplitudes that are applied on the mobile phantom where maximal DVF increase linearly with the motion amplitudes of the mobile phantom in CBCT. Similarly in HCT, DVF increase linearly with motion amplitude, however, its correlation is weaker than CBCT. In ACT, the DVF’s do not correlate well with the motion amplitudes where motion induces strong image artifacts and DIR algorithms are not able to deform the ACT image of the mobile targets to the stationary targets. Three DIR-algorithms produce comparable values and patterns of the DVF for certain CT imaging modality. However, DVF from fast-demons deviated strongly from other algorithms at large motion amplitudes. Conclusion: In CBCT and HCT, the DVF correlate well with the motion amplitude of the mobile phantom. However, in ACT, DVF do not correlate with motion amplitudes. Correlations of DVF with motion amplitude as in CBCT and HCT imaging techniques can provide information about unknown motion parameters of the mobile organs in real patients as demonstrated in this phantom visibility study.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-84: Comparison of Quantitative Deformable Image Registration Evaluation Tools: Application to Prostate IGART

Purpose: Errors in displacement vector fields (DVFs) generated by Deformable Image Registration (DIR) algorithms can give rise to significant uncertainties in contour propagation and dose accumulation in Image-Guided Adaptive Radiotherapy (IGART). The purpose of this work is to assess the accuracy of two DIR algorithms using a variety of quality metrics for prostate IGART. Methods: Pelvic CT images were selected from an anonymized database of nineteen prostate patients who underwent 8–12 serial scans during radiotherapy. Prostate, bladder, and rectum were contoured on 34 image-sets for three patients by the same physician. The planning CT was deformably-registered to daily CT using three variants of the Small deformation Inverse Consistent Linear Elastic (SICLE) algorithm: Grayscale-driven (G), Contour-driven (C, which utilizes segmented structures to drive DIR), combined (G+C); and also grayscale ITK demons (Gd). The accuracy of G, C, G+C SICLE and Gd registrations were evaluated using a new metric Edge Gradient Distance to Agreement (EGDTA) and other commonly-used metrics such as Pearson Correlation Coefficient (PCC), Dice Similarity Index (DSI) and Hausdorff Distance (HD). Results: C and G+C demonstrated much better performance at organ boundaries, revealing the lowest HD and highest DSI, in prostate, bladder and rectum. G+C demonstrated the lowest mean EGDTA (1.14 mm), which corresponds to highest registration quality, compared to G and C DVFs (1.16 and 2.34 mm). However, demons DIR showed the best overall performance, revealing lowest EGDTA (0.73 mm) and highest PCC (0.85). Conclusion: As expected, both C- and C+G SICLE more accurately reproduce manually-contoured target datasets than G-SICLE or Gd using HD and DSI metrics. In general, the Gd appears to have difficulty reproducing large daily position and shape changes in the rectum and bladder. However, Gd outperforms SICLE in terms of EGDTA and PCC metrics, possibly at the expense of topological quality of the estimated DVFs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-87: Impact Of The Dosimetric Consequences From Minimal Displacements Throughout The Treatment Time In APBI With SAVI Applicators

Purpose: To assess the variation in dose received by the organs at risk (OARs) due to inter-fractional motion by SAVI to determine the importance of providing proper immobilization Methods: An analysis of 15 patients treated with SAVI applicators were considered for this study. Treatment planning teams did not see significant changes in their CT scans through scout images and initial treatment plan was used for the entire treatment. These scans, taken before each treatment were imported in to the treatment planning system and were fused together with respective to the applicator, using landmark registration. Dosimetric evaluations were performed. Dose received by skin, ribs and PTV(Planning target volume) respect to the initial treatment plan were measured. Results: Contours of the OARs were not similar with the initial image. Deduction in volumes of PTV and cavity, small deviations in displacements from the applicator to the OARs, difference in doses received by the OARs between treatments were noticed. The maximum, minimum, average doses varied between 10% to 20% 5% to 8% and 15% to 20% in ribs and skin. The 0.1cc doses to OARs showed an average change of 10% of the prescribed dose. PTV was receiving a different dose than the estimated dose Conclusion: The variation in volumes and isodoses related to the OARs, PTV receiving a lesser dose than the prescribed dose indicate that the estimated doses are different from the received dose. This study reveals the urgent need of improving the immobilization methods. Taking a CT scan before each treatment and replanning is helpful to minimize the risk of delivering undesired high doses to the OARs. Patient positioning, motion, respiration, observer differences and time lap between the planning and treating can arise more complications. VacLock, Positioning cushions, Image guided brachytherapy and adjustable registration should be used for further improvements.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-94: Development of a Plug-in Based Image Analysis Tool for Integration Into Treatment Planning

Purpose: To extend the functionality of a commercial treatment planning system (TPS) to support (i) direct use of quantitative image-based metrics within treatment plan optimization and (ii) evaluation of dose-functional volume relationships to assist in functional image adaptive radiotherapy. Methods: A script was written that interfaces with a commercial TPS via an Application Programming Interface (API). The script executes a program that performs dose-functional volume analyses. Written in C#, the script reads the dose grid and correlates it with image data on a voxel-by-voxel basis through API extensions that can access registration transforms. A user interface was designed through WinForms to input parameters and display results. To test the performance of this program, image- and dose-based metrics computed from perfusion SPECT images aligned to the treatment planning CT were generated, validated, and compared. Results: The integration of image analysis information was successfully implemented as a plug-in to a commercial TPS. Perfusion SPECT images were used to validate the calculation and display of image-based metrics as well as dose-intensity metrics and histograms for defined structures on the treatment planning CT. Various biological dose correction models, custom image-based metrics, dose-intensity computations, and dose-intensity histograms were applied to analyze the image-dose profile. Conclusion: It is possible to add image analysis features to commercial TPSs through custom scripting applications. A tool was developed to enable the evaluation of image-intensity-based metrics in the context of functional targeting and avoidance. In addition to providing dose-intensity metrics and histograms that can be easily extracted from a plan database and correlated with outcomes, the system can also be extended to a plug-in optimization system, which can directly use the computed metrics for optimization of post-treatment tumor or normal tissue response models. Supported by NIH - P01 - CA059827.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-96: Comparison of Frame-Based and Mutual Information Registration Techniques for CT and MR Image Sets

Purpose: To compare the accuracy of CT-MR registration using a mutual information method with registration using a frame-based localizer box. Methods: Ten patients having the Leksell head frame and scanned with a modality specific localizer box were imported into the treatment planning system. The fiducial rods of the localizer box were contoured on both the MR and CT scans. The skull was contoured on the CT images. The MR and CT images were registered by two methods. The frame-based method used the transformation that minimized the mean square distance of the centroids of the contours of the fiducial rods from a mathematical model of the localizer. The mutual information method used automated image registration tools in the TPS and was restricted to a volume-of-interest defined by the skull contours with a 5 mm margin. For each case, the two registrations were adjusted by two evaluation teams, each comprised of an experienced radiation oncologist and neurosurgeon, to optimize alignment in the region of the brainstem. The teams were blinded to the registration method. Results: The mean adjustment was 0.4 mm (range 0 to 2 mm) and 0.2 mm (range 0 to 1 mm) for the frame and mutual information methods, respectively. The median difference between the frame and mutual information registrations was 0.3 mm, but was not statistically significant using the Wilcoxon signed rank test (p=0.37). Conclusion: The difference between frame and mutual information registration techniques was neither statistically significant nor, for most applications, clinically important. These results suggest that mutual information is equivalent to frame-based image registration for radiosurgery. Work is ongoing to add additional evaluators and to assess the differences between evaluators.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-101: Stereotactic Body Radiation Therapy Planning for Primary Prostate Cancer with Selective Intraprostatic Boost Determined by 18F-Choline PET/CT

Purpose: To investigate the utility of {sup 18}F-choline positron emission tomography (PET) scans guidance for SBRT dose painting in patients with prostate cancer and its impact on tumor control probability (TCP) and normal tissue complication probability (NTCP). Methods: Twenty seven patients with localized prostate cancer who had {sup 18}F-choline PET/CT scan prior to treatment were included. A pair of nested intraprostatic dominant lesion (IDL) contours (IDL{sub suv60%} and IDL{sub suv70%}) were generated for each patient based on 60% and 70% of maximum prostate uptake on the {sup 18}F-choline PET images. GTV{sub reg} was delineated on prostate according to the gland boundary seen on CT images. The PTVs (PTV{sub suv60%} and PTV{sub suv70%}) were defined as respective IDLs with a 3-mm margin posteriorly and 5 mm in all other dimensions. Two 5-fraction SBRT plans using VMAT technique along with 10 MV FFF beams, plan{sub 36Gy} and plan{sub 50–55Gy}, were generated for each patient. All plans included a dose of 36.25 Gy prescribed to PTV{sub reg}. The Plan{sub 50–55Gy} also included a simultaneous boost dose of 50 Gy and 55 Gy prescribed to the PTV{sub suv60%} and PTV{sub suv70%}, respectively. The utility of {sup 18}F-Choline PET-guided SBRT dose escalation was evaluated by its ability to achieve the prescription dose objectives while adhering to organ-at-risk (OAR) dose constraints. The TCP and NTCP calculated by radiological models were also compared between two plans for each patient. Results: In all 54 SBRT plans generated, the planning objectives and dose constraints were met without exception. Plan{sub 50–55Gy} had a significantly higher dose in PTV{sub suv60%} and PTV{sub suv70%} than those in Plan{sub 36Gy} (p < 0.05), respectively, while still maintaining a safe OAR sparing profile. In addition, plan{sub 50–55Gy} had significantly higher TCP than plan{sub 36Gy}. Conclusion: Using VMAT with FFF beams to incorporate a simultaneous {sup 18}F-choline PET-guided radiation boost dose up to 55 Gy into a SBRT plan is technically feasible. This work was supported in part by Congressionally Directed Medical Research Programs Prostate Cancer Research Program grant PC04130, National Institutes of Health/National Cancer Institute grant R41CA110121, and the UNLV Lincy Endowed Assistant Professorship.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-113: Multi-Atlas Based Automatic Organ Segmentation for Lung Radiotherapy Planning

Purpose: Normal organ segmentation is one time-consuming and labor-intensive step for lung radiotherapy treatment planning. The aim of this study is to evaluate the performance of a multi-atlas based segmentation approach for automatic organs at risk (OAR) delineation. Methods: Fifteen Lung stereotactic body radiation therapy patients were randomly selected. Planning CT images and OAR contours of the heart - HT, aorta - AO, vena cava - VC, pulmonary trunk - PT, and esophagus – ES were exported and used as reference and atlas sets. For automatic organ delineation for a given target CT, 1) all atlas sets were deformably warped to the target CT, 2) the deformed sets were accumulated and normalized to produce organ probability density (OPD) maps, and 3) the OPD maps were converted to contours via image thresholding. Optimal threshold for each organ was empirically determined by comparing the auto-segmented contours against their respective reference contours. The delineated results were evaluated by measuring contour similarity metrics: DICE, mean distance (MD), and true detection rate (TD), where DICE=(intersection volume/sum of two volumes) and TD = {1.0 - (false positive + false negative)/2.0}. Diffeomorphic Demons algorithm was employed for CT-CT deformable image registrations. Results: Optimal thresholds were determined to be 0.53 for HT, 0.38 for AO, 0.28 for PT, 0.43 for VC, and 0.31 for ES. The mean similarity metrics (DICE[%], MD[mm], TD[%]) were (88, 3.2, 89) for HT, (79, 3.2, 82) for AO, (75, 2.7, 77) for PT, (68, 3.4, 73) for VC, and (51,2.7, 60) for ES. Conclusion: The investigated multi-atlas based approach produced reliable segmentations for the organs with large and relatively clear boundaries (HT and AO). However, the detection of small and narrow organs with diffused boundaries (ES) were challenging. Sophisticated atlas selection and multi-atlas fusion algorithms may further improve the quality of segmentations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-114: On-Treatment Imagereconstruction Using Transit Images of Treatment Beams Through Patient and Thosethrough Planning CT Images

Purpose: To reconstruct patient images at the time of radiation delivery using measured transit images of treatment beams through patient and calculated transit images through planning CT images. Methods: We hypothesize that the ratio of the measured transit images to the calculated images may provide changed amounts of the patient image between times of planning CT and treatment. To test, we have devised lung phantoms with a tumor object (3-cm diameter) placed at iso-center (simulating planning CT) and off-center by 1 cm (simulating treatment). CT images of the two phantoms were acquired; the image of the off-centered phantom, unavailable clinically, represents the reference on-treatment image in the image quality of planning CT. Cine-transit images through the two phantoms were also acquired in EPID from a non-modulated 6 MV beam when the gantry was rotated 360 degrees; the image through the centered phantom simulates calculated image. While the current study is a feasibility study, in reality our computational EPID model can be applicable in providing accurate transit image from MC simulation. Changed MV HU values were reconstructed from the ratio between two EPID projection data, converted to KV HU values, and added to the planning CT, thereby reconstructing the on-treatment image of the patient limited to the irradiated region of the phantom. Results: The reconstructed image was compared with the reference image. Except for local HU differences>200 as a maximum, excellent agreement was found. The average difference across the entire image was 16.2 HU. Conclusion: We have demonstrated the feasibility of a method of reconstructing on-treatment images of a patient using EPID image and planning CT images. Further studies will include resolving the local HU differences and investigation on the dosimetry impact of the reconstructed image.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-129: Verification of Geometric and Dosimetric Accuracy of Respiratory Management Systems Using Homemade Phantom

Purpose: Different placements of Infrared Cameras (IRC) in CT and treatment rooms can cause gating window level (GWL) variations leading to differences between GWL used for planning and treatments. Although, Varian Clinac DHX-OBI sytem and CT are equipped with the same kind of IRC, Truebeam STx (TB) has a different type of IRC known as banana type. In this study; geometric and dosimetric accuracy of respiratory management system (RPM) for different machines were investigated with a special homemade phantom. Methods: Special phantom was placed on the respiratory simulator machine and a CT data set was obtained at the end of the expirium phase (EOE). Conformal and IMRT plans were generated on the EOE CT image series for both DHX-OBI and TB LINACs while a VMAT plan was generated only for TB.The acquired respiratory graphs in the CT were directly sent to DHX-OBI system, and they were converted with software before sending to TB. EBT3 films were placed inside the phantom and were irradiated using RPM system with two machines for different plans. Planar dose distributions were compared with gamma analysis (GA) method (3mm, %3) to evaluate planned-measured dose differences. In addition, radio-opac marker was placed in the center of the phantom to evaluate the geometric accuracy of treatment field with gated flouroscopy (GF). Results: There were no shifts detected between planning and treeatment GWL for both DHX-OBI and TB. Difference on the GF image between digital graticule and radio-opac marker was <1mm for TB and 1mm for DHX-OBI. Although, GA agreement was 97% for conformal and IMRT techniques in TB, it was 96% for VMAT technique. While GA agreement was 98% for conformal technique in DHX-OBI, IMRT was 95%.ConclusionThis study showed that RPM can be used accurately in spite of different IRC placements or different types of ICR used.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-135: Tumor Displacement-Based Binning for Respiratory-Gated Time-Independent 5DCT Treatment Planning

Purpose: A published 5DCT breathing motion model enables image reconstruction at any user-selected breathing phase, defined by the model as a specific amplitude (v) and rate (f). Generation of reconstructed phase-specific CT scans will be required for time-independent radiation dose distribution simulations. This work answers the question: how many amplitude and rate bins are required to describe the tumor motion with a specific spatial resolution? Methods: 19 lung-cancer patients with 21 tumors were scanned using a free-breathing 5DCT protocol, employing an abdominally positioned pneumatic-bellows breathing surrogate and yielding voxel-specific motion model parameters α and β corresponding to motion as a function of amplitude and rate, respectively. Tumor GTVs were contoured on the first (reference) of 25 successive free-breathing fast helical CT image sets. The tumor displacements were binned into widths of 1mm to 5mm in 1mm steps and the total required number of bins recorded. The simulation evaluated the number of bins needed to encompass 100% of the breathing-amplitude and between the 5th and 95th percentile amplitudes to exclude breathing outliers. Results: The mean respiration-induced tumor motion was 9.90mm ± 7.86mm with a maximum of 25mm. The number of bins required was a strong function of the spatial resolution and varied widely between patients. For example, for 2mm bins, between 1–13 amplitude bins and 1–9 rate bins were required to encompass 100% of the breathing amplitude, while 1–6 amplitude bins and 1–3 rate bins were required to encompass 90% of the breathing amplitude. Conclusion: The strong relationship between number of bins and spatial resolution as well as the large variation between patients implies that time-independent radiation dose distribution simulations should be conducted using patient-specific data and that the breathing conditions will have to be carefully considered. This work will lead to the assessment of the dosimetric impact of binning resolution. This study is supported by Siemens Healthcare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-138: An Extension of PCA-Based Respiratory Deformation Modeling Via Multi-Linear Decomposition

Purpose: To address and lift the limited degree of freedom (DoF) of globally bilinear motion components such as those based on principal components analysis (PCA), for encoding and modeling volumetric deformation motion. Methods: We provide a systematic approach to obtaining a multi-linear decomposition (MLD) and associated motion model from deformation vector field (DVF) data. We had previously introduced MLD for capturing multi-way relationships between DVF variables, without being restricted by the bilinear component format of PCA-based models. PCA-based modeling is commonly used for encoding patient-specific deformation as per planning 4D-CT images, and aiding on-board motion estimation during radiotherapy. However, the bilinear space-time decomposition inherently limits the DoF of such models by the small number of respiratory phases. While this limit is not reached in model studies using analytical or digital phantoms with low-rank motion, it compromises modeling power in the presence of relative motion, asymmetries and hysteresis, etc, which are often observed in patient data. Specifically, a low-DoF model will spuriously couple incoherent motion components, compromising its adaptability to on-board deformation changes. By the multi-linear format of extracted motion components, MLD-based models can encode higher-DoF deformation structure. Results: We conduct mathematical and experimental comparisons between PCA- and MLD-based models. A set of temporally-sampled analytical trajectories provides a synthetic, high-rank DVF; trajectories correspond to respiratory and cardiac motion factors, including different relative frequencies and spatial variations. Additionally, a digital XCAT phantom is used to simulate a lung lesion deforming incoherently with respect to the body, which adheres to a simple respiratory trend. In both cases, coupling of incoherent motion components due to a low model DoF is clearly demonstrated. Conclusion: Multi-linear decomposition can enable decoupling of distinct motion factors in high-rank DVF measurements. This may improve motion model expressiveness and adaptability to on-board deformation, aiding model-based image reconstruction for target verification. NIH Grant No. R01-184173.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-153: ACR MRI Phantom Test On a 1.5T MR-Simulator with Flexible Coils Setting and Performance Comparison to Volumetric Head Coil Setting

Purpose: To accommodate head-and-neck (HN) RT simulation scan with immobilization, sub-optimized flexible radiofrequency coils have to be used on an MR-sim rather than a volumetric head coil for diagnostic radiological (DR) applications. In this study, we present and compare ACR MRI phantom test with flexible coils setting (RT-setting) and volumetric head coil setting (DR-setting). Methods: Two 8-channel array coils were wrapped close to an ACR MRI phantom placed on a flat couch-top in RT-setting for acquisition. ACR MRI phantom tests were conducted 10 times for each setting on a 1.5T MR-sim. Sagittal localizer (TE/TR=20/200ms), axial T1 (TE/TR=20/500ms) and T2 scans (2TEs/TR=20,80/2000ms) were acquired (NEX=1). Acqusition and image analysis were conducted following ACR guidelines. Extra 10 RT-setting acquisitions were obtained by enabling prescan normalization. Results: For RT-setting without prescan normalization, all tests passed the ACR recommended criteria but image intensity uniformity. Geometric accuracy was 147.8±0.3mm and 191.0±0.4mm. Spatial resolution was 0.9mm. Slice thickness was 5.2mm±0.1mm (T1) and 5.1±0.1mm (T2). Slice position accuracy was 1.4±0.6mm, −2.2±0.3mm (T1), and 1.5±0.2mm, −2.2±0.2mm (T2). Image uniformity was 80.6%±2.6% (T1) and 80.2%±3.1% (T2) (criteria: ≥87.5% for 1.5T). Ghosting ratio was 0.0006±0.0004. Low contrast detectability was 30.7±1.7 (T1) and 22.6±3.0 (T2). Image uniformity (91.4%±4.2% and 91.8%±4.3%) and low contrast detectability (36.5±1.4 and 31.6±2.2) of DR-setting were better than RT-setting. Prescan normalization substantially improved image uniformity to 93.4%±0.2% and 93.4%±0.3%, but slightly compromised ghosting (0.0033±0.0004) and spatial resolution (0.9mm or 1.0mm) for RT-setting. Conclusion: Flexible coils setting for HN-RT simulation scan could successfully pass ACR criteria (prescan normalization enabled) and generally achieve comparable performance to volumetric coil setting, although compromise on low contrast detectability was observed probably due to the lower signal-to-noise ratio of the flexible coils. Our results could also serve as a reference baseline of high-field MR-sim QA performance under ACR guideline.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-156: The Feasibility of MR-Only IMRT Planning for Prostate Anatomy

Purpose: For prostate anatomy, previous investigations have shown that simulated CT (sCT) generated from MR images can be used for accurate dose computation. In this study, we demonstrate the feasibility of MR-only IMRT planning for prostate case. Methods: Regular CT (rCT) and MR images of the same patient were acquired for prostate anatomy. Regions-of-interest (ROIs) i.e. target and risk structures are delineated on the rCT. A simulated CT (sCT) is generated from the MR image using the method described by Schadewaldt N et al. Their work establishes the clinical acceptability of dose calculation results on the sCT when compared to rCT. rCT and sCT are rigidly registered to ensure proper alignment between the two images. rCT and sCT are overlaid on each other and slice-wise visual inspection confirms excellent agreement between the two images. ROIs on the rCT are copied over to sCT. Philips AutoPlanning solution is used for generating treatment plans. The same treatment technique protocol (plan parameters and clinical goals) is used to generate AutoPlan-rCT and AutoPlan-sCT respectively for rCT and and sCT. DVH comparison on ROIs and slice-wise evaluation of dose is performed between AutoPlan-rCT and AutoPlan-sCT. Delivery parameters i.e. beam and corresponding segments from the AutoPlan-sCT are copied over to rCT and dose is computed to get AutoPlan-sCT-on-rCT. Results: Plan evaluation is done based on Dose Volume Histogram (DVH) of ROIs and manual slice-wise inspection of dose distribution. Both AutoPlan-rCT and AutoPlan-sCT provide a clinically acceptable plan. Also, AutoPlan-sCT-on-rCT shows excellent agreement with AutoPlan-sCT. Conclusion: The study demonstrates that it is feasible to do IMRT planning on the simulated CT image obtained from MR image for prostate anatomy. The research is supported by Philips India Ltd.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-161: Prostate Contouring in Patients with Bilateral Hip Prostheses: Impact of Using Artifact-Reduced CT Images and MRI

Purpose: To investigate the usefulness of metal artifact reduction in CT images of patients with bilateral hip prostheses (BHP) for contouring the prostate and determine if the inclusion of MR images provides additional benefits. Methods: Five patients with BHP were CT scanned using our clinical protocol (140kV, 300mAs, 3mm slices, 1.5mm increment, Philips Medical Systems, OH). Images were reconstructed with the orthopaedic metal artifact reduction (O-MAR) algorithm. MRI scanning was then performed (1.5T, GE Healthcare, WI) with a flat table-top (T{sub 2}-weighted, inherent body coil, FRFSE, 3mm slices with 0mm gap). All images were transferred to Pinnacle (Version 9.2, Philips Medical Systems). For each patient, two data sets were produced: one containing the O-MAR-corrected CT images and another containing fused MRI and O-MAR-corrected CT images. Four genito-urinary radiation oncologists contoured the prostate of each patient on the O-MAR-corrected CT data. Two weeks later, they contoured the prostate on the fused data set, blinded to all other contours. During each contouring session, the oncologists reported their confidence in the contours (1=very confident, 3=not confident) and the contouring difficulty that they experienced (1=really easy, 4=very challenging). Prostate volumes were computed from the contours and the conformity index was used to evaluate inter-observer variability. Results: Larger prostate volumes were found on the O-MAR-corrected CT set than on the fused set (p< 0.05, median=36.9cm{sup 3} vs. 26.63 cm{sup 3}). No significant differences were noted in the inter-observer variability between the two data sets (p=0.3). Contouring difficulty decreased with the addition of MRI (p<0.05) while the radiation oncologists reported more confidence in their contours when MRI was fused with the O-MAR-corrected CT data (p<0.05). Conclusion: This preliminary work demonstrated that, while O-MAR correction to CT images improves visualization of anatomy, the addition of MRI enhanced the oncologists’ confidence in contouring the prostate in patients with BHP.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-175: Evaluation of Metal Artifact Reduction Algorithms in Computed Tomography and Their Application to Radiation Therapy Treatment Planning

Purpose: High-Z (metal) implants in CT scans cause significant streak-like artifacts in the reconstructed dataset. This results in both inaccurate CT Hounsfield units for the tissue as well as obscuration of the target and organs at risk (OARs) for radiation therapy planning. Herein we analyze two metal artifact reduction algorithms: GE’s Smart MAR and a Metal Deletion Technique (MDT) for geometric and Hounsfield Unit (HU) accuracy. Methods: A CT-to-electron density phantom, with multiple inserts of various densities and a custom Cerrobend insert (Zeff=76.8), is utilized in this continuing study. The phantom is scanned without metal (baseline) and again with the metal insert. Using one set of projection data, reconstructed CT volumes are created with filtered-back-projection (FBP) and the MAR and the MDT algorithms. Regions-of-Interest (ROIs) are evaluated for each insert for HU accuracy; the metal insert’s Full-Width-Half-Maximum (FWHM) is used to evaluate the geometric accuracy. Streak severity is quantified with an HU error metric over the phantom volume. Results: The original FBP reconstruction has a Root-Mean-Square-Error (RMSE) of 57.55 HU (STD=29.19, range=−145.8 to +79.2) compared to baseline. The MAR reconstruction has a RMSE of 20.98 HU (STD=13.92, range=−18.3 to +61.7). The MDT reconstruction has a RMSE of 10.05 HU (STD=10.5, range=−14.8 to +18.6). FWHM for baseline=162.05; FBP=161.84 (−0.13%); MAR=162.36 (+0.19%); MDT=162.99 (+0.58%). Streak severity metric for FBP=19.73 (22.659% bad pixels); MAR=8.743 (9.538% bad); MDT=4.899 (5.303% bad). Conclusion: Image quality, in terms of HU accuracy, in the presence of high-Z metal objects in CT scans is improved by metal artifact reduction reconstruction algorithms. The MDT algorithm had the highest HU value accuracy (RMSE=10.05 HU) and best streak severity metric, but scored the worst in terms of geometric accuracy. Qualitatively, the MAR and MDT algorithms increased detectability of inserts, although there is a loss of in-plane resolution near the metallic insert.

60 APPLIED LIFE SCIENCES↗