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 1,585 records · Page 88

TU-FG-BRB-02: The Impact of Using Dual-Energy CT for Determining Proton Stopping Powers: Comparison Between Theory and Experiments

Purpose: To evaluate the clinical performance of dual-energy CT (DECT) in determining proton stopping power ratios (SPR) and demonstrate advantages over conventional single-energy CT (SECT). Methods: SECT and DECT scans of tissue-equivalent plastics as well as animal meat samples are performed with a Siemens SOMATOM Definition Flash. The methods of Schneider et al. (1996) and Bourque et al. (2014) are used to determine proton SPR on SECT and DECT images, respectively. Waterequivalent path length (WEPL) measurements of plastics and tissue samples are performed with a 195 MeV proton beam. WEPL values are determined experimentally using the depth-dose shift and dose extinction methods. Results: Comparison between CT-based and experimental WEPL is performed for 12 tissue-equivalent plastic as well as 6 meat boxes containing animal liver, kidney, heart, stomach, muscle and bones. For plastic materials, results show a systematic improvement in determining SPR with DECT, with a mean absolute error of 0.4% compared to 1.7% for SECT. For the meat samples, preliminary results show the ability for DECT to determine WEPL with a mean absolute value of 1.1% over all meat boxes. Conclusion: This work demonstrates the potential in using DECT for determining proton SPR with plastic materials in a clinical context. Further work is required to show the benefits of DECT for tissue samples. While experimental uncertainties could be a limiting factor to show the benefits of DECT over SECT for the meat samples, further work is required to adapt the DECT formalism in the context of clinical use, where noise and artifacts play an important role.

07 ISOTOPE AND RADIATION SOURCES↗

TU-H-BRC-08: Use and Validation of Flexible 3D Printed Tissue Compensators for Post-Mastectomy Radiation Therapy

Purpose: Patient-specific tissue equivalent compensators can be used for post-mastectomy radiation therapy (PMRT) to achieve homogenous dose distributions with single-field treatments. However, current fabrication methods are time consuming and expensive. 3D-printing technology could overcome these limitations. The purposes of this study were to [1] evaluate materials for 3D-printed compensators [2] design and print a compensator to achieve a uniform thickness to a clinical target volume (CTV), and [3] demonstrate that a single-field electron compensator plan is a clinically feasible treatment option for PMRT. Methods: Blocks were printed with three materials; print accuracy, density, Hounsfield units (HU), and percent depth doses (PDD) were evaluated. For a CT scan of an anthropomorphic phantom, we used a ray-tracing method to design a compensator that achieved uniform thickness from compensator surface to CTV. The compensator was printed with flexible tissue equivalent material whose physical and radiological properties were most similar to soft tissue. A single-field electron compensator plan was designed and compared with two standard-of-care techniques. The compensator plan was validated with thermoluminescent dosimeter (TLD) measurements. Results: We identified an appropriate material for 3D-printed compensators that had high print accuracy (99.6%) and was similar to soft tissue; density was 1.04, HU was - 45 ± 43, and PDD curves agreed with clinical curves within 3 mm. We designed and printed a compensator that conformed well to the phantom surface and created a uniform thickness to the CTV. In-house fabrication was simple and inexpensive (<$75). Compared with the two standard plans, the compensator plan resulted in overall more homogeneous dose distributions and performed similarly in terms of lung/heart doses and 90% isodose coverage of the CTV. TLD measurements agreed well with planned doses (within 5 %). Conclusions: We have demonstrated that 3D-printed compensators make single-field electron therapy a clinically feasible treatment option for PMRT.

60 APPLIED LIFE SCIENCES↗

TU-H-CAMPUS-JeP2-05: Can Automatic Delineation of Cardiac Substructures On Noncontrast CT Be Used for Cardiac Toxicity Analysis?

Purpose: To evaluate the feasibility of using an automatic segmentation tool to delineate cardiac substructures from computed tomography (CT) images for cardiac toxicity analysis for non-small cell lung cancer (NSCLC) patients after radiotherapy. Methods: A multi-atlas segmentation tool developed in-house was used to delineate eleven cardiac substructures including the whole heart, four heart chambers, and six greater vessels automatically from the averaged 4DCT planning images for 49 NSCLC patients. The automatic segmented contours were edited appropriately by two experienced radiation oncologists. The modified contours were compared with the auto-segmented contours using Dice similarity coefficient (DSC) and mean surface distance (MSD) to evaluate how much modification was needed. In addition, the dose volume histogram (DVH) of the modified contours were compared with that of the auto-segmented contours to evaluate the dosimetric difference between modified and auto-segmented contours. Results: Of the eleven structures, the averaged DSC values ranged from 0.73 ± 0.08 to 0.95 ± 0.04 and the averaged MSD values ranged from 1.3 ± 0.6 mm to 2.9 ± 5.1mm for the 49 patients. Overall, the modification is small. The pulmonary vein (PV) and the inferior vena cava required the most modifications. The V30 (volume receiving 30 Gy or above) for the whole heart and the mean dose to the whole heart and four heart chambers did not show statistically significant difference between modified and auto-segmented contours. The maximum dose to the greater vessels did not show statistically significant difference except for the PV. Conclusion: The automatic segmentation of the cardiac substructures did not require substantial modification. The dosimetric evaluation showed no statistically significant difference between auto-segmented and modified contours except for the PV, which suggests that auto-segmented contours for the cardiac dose response study are feasible in the clinical practice with a minor modification to the PV vessel.

60 APPLIED LIFE SCIENCES↗

WE-AB-BRA-05: PET-Guided Delivery Quality Evaluation of Yttrium-90 Microsphere Radioembolizaton for Hepatocellular Carcinoma Patients: The Optimal Sequence of Radioembolizaton and Chemoembolization Treatments

Purpose: Selective-internal-radiation-therapy (SIRT) and transarterial-chemoembolization (TACE) are commonly used for treatment of liver tumors. The use of TACE, which is macroembolic, prior to SIRT may cause hemodynamic changes in tumor vasculature that impair yttrium-90 (90Y) microsphere delivery to the targeted lesions. This work aims to quantify dosimetric tumor coverage using 90Y positron emission tomography (PET) dosimetry after SIRT alone compared to TACE followed by SIRT. Methods: A total of 40 consecutive hepatocellular carcinoma (HCC) SIRT patients who had a post-SIRT 90Y PET/CT scan were evaluated. The patient-specific-3D-dose was reconstructed from the PET images. Patients were categorized into two groups: patients received TACE prior SIRT procedure (n=18) and patient received SIRT alone (n=22). The lesions and liver were delineated by a senior radiologist. We evaluated both the lesion-specific dose-volume-histogram (DVH) and the selectivity index (SI) defined as the ratio of the average dose inside the total lesion(s) and the average dose of the normal liver. The SI values of patients were compared based on whether TACE was previously used. Results: A wide spectrum was observed in the lesion-specific DVH-evaluation and SI appeared to be suitable of evaluating the quality of each SIRT infusion. The average SI of the entire patient group was 3.0, i.e. targeted lesion receiving three times higher dose than normal liver. The average SI was 1.8 for patients who had prior TACE and 3.9 for patients who did not have prior TACE (p=0.008). 85% of the patients with prior TACE demonstrated poor 90Y-microsphere delivery (SI <2) while none demonstrated excellent delivery (SI >4). On the other hand, the incidence SI >4 among patients with no prior TACE was 37%. Conclusion: 3D dose evaluation using post-SIRT PET suggests that 90Y microsphere delivery to liver tumors is impaired among patients who received prior TACE compared to those who receive SIRT alone.

60 APPLIED LIFE SCIENCES↗

WE-AB-202-09: Feasibility and Quantitative Analysis of 4DCT-Based High Precision Lung Elastography

Purpose: The purpose of this project is to derive high precision elastography measurements from 4DCT lung scans to facilitate the implementation of elastography in a radiotherapy context. Methods: 4DCT scans of the lungs were acquired, and breathing stages were subsequently registered to each other using an optical flow DIR algorithm. The displacement of each voxel gleaned from the registration was taken to be the ground-truth deformation. These vectors, along with the 4DCT source datasets, were used to generate a GPU-based biomechanical simulation that acted as a forward model to solve the inverse elasticity problem. The lung surface displacements were applied as boundary constraints for the model-guided lung tissue elastography, while the inner voxels were allowed to deform according to the linear elastic forces within the model. A biomechanically-based anisotropic convergence magnification technique was applied to the inner voxels in order to amplify the subtleties of the interior deformation. Solving the inverse elasticity problem was accomplished by modifying the tissue elasticity and iteratively deforming the biomechanical model. Convergence occurred when each voxel was within 0.5 mm of the ground-truth deformation and 1 kPa of the ground-truth elasticity distribution. To analyze the feasibility of the model-guided approach, we present the results for regions of low ventilation, specifically, the apex. Results: The maximum apical boundary expansion was observed to be between 2 and 6 mm. Simulating this expansion within an apical lung model, it was observed that 100% of voxels converged within 0.5 mm of ground-truth deformation, while 91.8% converged within 1 kPa of the ground-truth elasticity distribution. A mean elasticity error of 0.6 kPa illustrates the high precision of our technique. Conclusion: By utilizing 4DCT lung data coupled with a biomechanical model, high precision lung elastography can be accurately performed, even in low ventilation regions of the lungs. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE-1144087.

60 APPLIED LIFE SCIENCES↗

WE-AB-207A-03: A CBCT Head Scanner for Point-Of-Care Imaging of Intracranial Hemorrhage

Purpose: This work reports the design, development, and first technical assessment of a cone-beam CT (CBCT) scanner developed specifically for imaging of acute intracranial hemorrhage (ICH) at the point of care, with target applications in diagnosis and monitoring of traumatic brain injury, stroke, and postsurgical hemorrhage. Methods: System design employed a task-based image quality model to quantify the influence of factors such as additive noise and high-gain (HG) detector readout on ICH detectability. Three bowtie filters with varying bare-beam attenuation strength and curvature were designed to enable HG readout without detector saturation, and a polyenergetic gain correction was developed to minimize artifacts from bowtie flood-field calibration. Image reconstruction used an iterative penalized weighted least squares (PWLS) method with artifact correction including Monte Carlo scatter estimation, Joseph-Spital beam hardening correction, and spatiotemporal deconvolution of detector glare and lag. Radiation dose was characterized for half-scan and full-scan protocols at various kV, and imaging performance was assessed in a head phantom presenting simulated ICH with diameter ranging 2–12 mm. Results: The image quality model guided system design and was validated by measurements on a CBCT imaging bench. Compared to low-gain readout without a bowtie filter, the combination of HG readout and a modest bowtie improved the contrast-to-noise ratio (CNR per unit square-root dose) by 20% in the center of the image but degraded noise performance near the periphery (20% reduction in CNR). Low-frequency bowtie artifacts (∼100 HU magnitude) were corrected by the polyenergetic gain correction. Image reconstructions on the prototype scanner demonstrate clear visibility of the smallest ICH insert (2 mm diameter) in both HG readout (with a bowtie) and dual-gain readout (without bowtie). Conclusion: Technical assessment of the prototype scanner suggests the capability for reliable visualization of small (2 mm), low-contrast (50 HU) ICH at <20 mGy dose and motivates translation to clinical studies, now underway. Xiaohui Wang and David Foos (co-authors) are employees of Carestream Health. Research funding support from Carestream Health.

60 APPLIED LIFE SCIENCES↗

WE-AB-207A-04: Random Undersampled Cone Beam CT: Theoretical Analysis and a Novel Reconstruction Method

Purpose: Reducing x-ray exposure and speeding up data acquisition motived studies on projection data undersampling. It is an important question that for a given undersampling ratio, what the optimal undersampling approach is. In this study, we propose a new undersampling scheme: random-ray undersampling. We will mathematically analyze its projection matrix properties and demonstrate its advantages. We will also propose a new reconstruction method that simultaneously performs CT image reconstruction and projection domain data restoration. Methods: By representing projection operator under the basis of singular vectors of full projection operator, matrix representations for an undersampling case can be generated and numerical singular value decomposition can be performed. We compared properties of matrices among three undersampling approaches: regular-view undersampling, regular-ray undersampling, and the proposed random-ray undersampling. To accomplish CT reconstruction for random undersampling, we developed a novel method that iteratively performs CT reconstruction and missing projection data restoration via regularization approaches. Results: For a given undersampling ratio, random-ray undersampling preserved mathematical properties of full projection operator better than the other two approaches. This translates to advantages of reconstructing CT images at lower errors. Different types of image artifacts were observed depending on undersampling strategies, which were ascribed to the unique singular vectors of the sampling operators in the image domain. We tested the proposed reconstruction algorithm on a Forbid phantom with only 30% of the projection data randomly acquired. Reconstructed image error was reduced from 9.4% in a TV method to 7.6% in the proposed method. Conclusion: The proposed random-ray undersampling is mathematically advantageous over other typical undersampling approaches. It may permit better image reconstruction at the same undersampling ratio. The novel algorithm suitable for this random-ray undersampling was able to reconstruct high-quality images.

60 APPLIED LIFE SCIENCES↗

WE-AB-207A-09: Optimization of the Design of a Moving Blocker for Cone-Beam CT Scatter Correction: Experimental Evaluation

Purpose: A moving blocker based strategy has shown promising results for scatter correction in cone-beam computed tomography (CBCT). Different geometry designs and moving speeds of the blocker affect its performance in image reconstruction accuracy. The goal of this work is to optimize the geometric design and moving speed of the moving blocker system through experimental evaluations. Methods: An Elekta Synergy XVI system and an anthropomorphic pelvis phantom CIRS 801-P were used for our experiment. A blocker consisting of lead strips was inserted between the x-ray source and the phantom moving back and forth along rotation axis to measure the scatter signal. Accoriding to our Monte Carlo simulation results, three blockers were used, which have the same lead strip width 3.2mm and different gap between neighboring lead strips, 3.2, 6.4 and 9.6mm. For each blocker, three moving speeds were evaluated, 10, 20 and 30 pixels per projection (on the detector plane). Scatter signal in the unblocked region was estimated by cubic B-spline based interpolation from the blocked region. CBCT image was reconstructed by a total variation (TV) based algebraic iterative reconstruction (ART) algorithm from the partially blocked projection data. Reconstruction accuracy in each condition is quantified as CT number error of region of interest (ROI) by comparing to a CBCT reconstructed image from analytically simulated unblocked and scatter free projection data. Results: Highest reconstruction accuracy is achieved when the blocker width is 3.2 mm, the gap between neighboring lead strips is 9.6 mm and the moving speed is 20 pixels per projection. RMSE of the CT number of ROIs can be reduced from 436 to 27. Conclusions: Image reconstruction accuracy is greatly affected by the geometry design of the blocker. The moving speed does not have a very strong effect on reconstruction result if it is over 20 pixels per projection.

60 APPLIED LIFE SCIENCES↗

WE-FG-202-04: Decomposition of FDG-PET Based Differential Uptake Volume Histograms in Rectal Cancer Patients

Purpose: The goal of this study is to test the possible use of the analytical decomposition of differential uptake volume histograms (dUVHs) obtained from FDG-PET/CT data to isolate sub-volumes within a tumor known as biological target volumes (BTVs). Methods: : A retrospective study was conducted on a cohort of 20 histo-pathologically confirmed rectal adenocarcinoma patients having PET/CT scans for staging. All patients (T3N0) underwent pre-operative endorectal brachytherapy. After surgery, patients were restaged: 10 patients were T0N0 and 10 were restaged as remaining T3N0. The extent of the disease was sampled in order to create dUVHs; subsequently decomposed into the fewest number of analytical Gaussian functions. Results: With the assumption that each function fit corresponded to a single sub-volume within the tumor, six sub-volumes were found to consistently emerge. The first two sub-volumes were influenced by contouring and were not considered in the analysis. For the T3N0 population, abundances for volumes V3-V6 were 63.6%±11.3%, 25.7%±8.4%, 6.1%±4.9%, and 4.7%±2.6%. For the T0N0 population, they were 50.2%±6.8%, 33.4%±4.3%, 11.8%±7.6%, and 4.7%±2.4%. The two populations were compared using two tailed T-tests: volumes 3 and 4 were statistically different with p values of 0.021 and 0.056 respectively. V6 was located at 8.63 ± 2.2 for T0N0 and 6.14 ± 0.78 for T3N0 group (p=0.016). Conclusion: We described a method for dUVH decomposition using FDG-PET images of rectal adenocarcinoma patients that subsequently went for pre-operative brachytherapy. In addition to extracting different sub-volumes corresponding to different FDG uptake levels, we observed different abundances of two sub-volumes as well as positions of the maximum uptake between the two patient groups. In addition to opening the door to further investigation into underlying physiological phenotypes of segmented subvolumes and their use for biological radiotherapy treatment planning, this method may also provide parameters that could correlate to clinical outcomes in radiotherapy patients.

60 APPLIED LIFE SCIENCES↗

WE-FG-207B-02: Material Reconstruction for Spectral Computed Tomography with Detector Response Function

Purpose: Different from the conventional computed tomography (CT), spectral CT based on energy-resolved photon-counting detectors is able to provide the unprecedented material composition. However, an important missing piece for accurate spectral CT is to incorporate the detector response function (DRF), which is distorted by factors such as pulse pileup and charge-sharing. In this work, we propose material reconstruction methods for spectral CT with DRF. Methods: The polyenergetic X-ray forward model takes the DRF into account for accurate material reconstruction. Two image reconstruction methods are proposed: a direct method based on the nonlinear data fidelity from DRF-based forward model; a linear-data-fidelity based method that relies on the spectral rebinning so that the corresponding DRF matrix is invertible. Then the image reconstruction problem is regularized with the isotropic TV term and solved by alternating direction method of multipliers. Results: The simulation results suggest that the proposed methods provided more accurate material compositions than the standard method without DRF. Moreover, the proposed method with linear data fidelity had improved reconstruction quality from the proposed method with nonlinear data fidelity. Conclusion: We have proposed material reconstruction methods for spectral CT with DRF, whichprovided more accurate material compositions than the standard methods without DRF. Moreover, the proposed method with linear data fidelity had improved reconstruction quality from the proposed method with nonlinear data fidelity. Jiulong Liu and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000), and the Shanghai Pujiang Talent Program (#14PJ1404500).

60 APPLIED LIFE SCIENCES↗

TH-AB-207A-05: A Fully-Automated Pipeline for Generating CT Images Across a Range of Doses and Reconstruction Methods

Purpose: To evaluate the robustness of CAD or Quantitative Imaging methods, they should be tested on a variety of cases and under a variety of image acquisition and reconstruction conditions that represent the heterogeneity encountered in clinical practice. The purpose of this work was to develop a fully-automated pipeline for generating CT images that represent a wide range of dose and reconstruction conditions. Methods: The pipeline consists of three main modules: reduced-dose simulation, image reconstruction, and quantitative analysis. The first two modules of the pipeline can be operated in a completely automated fashion, using configuration files and running the modules in a batch queue. The input to the pipeline is raw projection CT data; this data is used to simulate different levels of dose reduction using a previously-published algorithm. Filtered-backprojection reconstructions are then performed using FreeCT-wFBP, a freely-available reconstruction software for helical CT. We also added support for an in-house, model-based iterative reconstruction algorithm using iterative coordinate-descent optimization, which may be run in tandem with the more conventional recon methods. The reduced-dose simulations and image reconstructions are controlled automatically by a single script, and they can be run in parallel on our research cluster. The pipeline was tested on phantom and lung screening datasets from a clinical scanner (Definition AS, Siemens Healthcare). Results: The images generated from our test datasets appeared to represent a realistic range of acquisition and reconstruction conditions that we would expect to find clinically. The time to generate images was approximately 30 minutes per dose/reconstruction combination on a hybrid CPU/GPU architecture. Conclusion: The automated research pipeline promises to be a useful tool for either training or evaluating performance of quantitative imaging software such as classifiers and CAD algorithms across the range of acquisition and reconstruction parameters present in the clinical environment. Funding support: NIH U01 CA181156; Disclosures (McNitt-Gray): Institutional research agreement, Siemens Healthcare; Past recipient, research grant support, Siemens Healthcare; Consultant, Toshiba America Medical Systems; Consultant, Samsung Electronics.

60 APPLIED LIFE SCIENCES↗

TH-AB-209-03: Overcoming Resolution Limitations of Diffuse Optical Signals in X-Ray Induced Luminescence (XIL) Imaging Via Selective Plane Illumination and 2D Deconvolution

Purpose: To evaluate the resolution and sensitivity of XIL imaging using a surface radiance simulation based on optical diffusion and maximum likelihood expectation maximization (MLEM) image reconstruction. XIL imaging seeks to determine the distribution of luminescent nanophosphors, which could be used as nanodosimeters or radiosensitizers. Methods: The XIL simulation generated a homogeneous slab with optical properties similar to tissue. X-ray activated nanophosphors were placed at 1.0 cm depth in the tissue in concentrations of 10{sup −4} g/mL in two volumes of 10 mm{sup 3} with varying separations between each other. An analytical optical diffusion model determined the surface radiance from the photon distributions generated at depth in the tissue by the nanophosphors. The simulation then determined the detected luminescent signal collected with a f/1.0 aperture lens and back-illuminated EMCCD camera. The surface radiance was deconvolved using a MLEM algorithm to estimate the nanophosphors distribution and the resolution. To account for both Poisson and Gaussian noise, a shifted Poisson imaging model was used in the deconvolution. The deconvolved distributions were fitted to a Gaussian after radial averaging to measure the full width at half maximum (FWHM) and the peak to peak distance between distributions was measured to determine the resolving power. Results: Simulated surface radiances for doses from 1mGy to 100 cGy were computed. Each image was deconvolved using 1000 iterations. At 1mGy, deconvolution reduced the FWHM of the nanophosphors distribution by 65% and had a resolving power is 3.84 mm. Decreasing the dose from 100 cGy to 1 mGy increased the FWHM by 22% but allowed for a dose reduction of a factor of 1000. Conclusion: Deconvolving the detected surface radiance allows for dose reduction while maintaining the resolution of the nanophosphors. It proves to be a useful technique in overcoming the resolution limitations of diffuse optical imaging in tissue. C. S. acknowledges support from the NIH National Institute of General Medical Sciences (Award number R25GM109439, Project Title: University of Chicago Initiative for Maximizing Student Development, IMSD). B. Q. and P. L. acknowledge support from NIH grant R01EB017293.

60 APPLIED LIFE SCIENCES↗

TH-AB-209-09: Quantitative Imaging of Electrical Conductivity by VHF-Induced Thermoacoustics

Purpose: To demonstrate that very high frequency (VHF) induced thermoacoustics has the potential to provide quantitative images of electrical conductivity in Siemens/meter, much as shear wave elastography provides tissue stiffness in kPa. Quantitatively imaging a large organ requires exciting thermoacoustic pulses throughout the volume and broadband detection of those pulses because tomographic image reconstruction preserves frequency content. Applying the half-wavelength limit to a 200-micron inclusion inside a 7.5 cm diameter organ requires measurement sensitivity to frequencies ranging from 4 MHz down to 10 kHz, respectively. VHF irradiation provides superior depth penetration over near infrared used in photoacoustics. Additionally, VHF signal production is proportional to electrical conductivity, and prostate cancer is known to suppress electrical conductivity of prostatic fluid. Methods: A dual-transducer system utilizing a P4-1 array connected to a Verasonics V1 system augmented by a lower frequency focused single element transducer was developed. Simultaneous acquisition of VHF-induced thermoacoustic pulses by both transducers enabled comparison of transducer performance. Data from the clinical array generated a stack of 96-images with separation of 0.3 mm, whereas the single element transducer imaged only in a single plane. In-plane resolution and quantitative accuracy were measured at isocenter. Results: The array provided volumetric imaging capability with superior resolution whereas the single element transducer provided superior quantitative accuracy. Combining axial images from both transducers preserved resolution of the P4-1 array and improved image contrast. Neither transducer was sensitive to frequencies below 50 kHz, resulting in a DC offset and low-frequency shading over fields of view exceeding 15 mm. Fresh human prostates were imaged ex vivo and volumetric reconstructions reveal structures rarely seen in diagnostic images. Conclusion: Quantitative whole-organ thermoacoustic tomography will be feasible by sparsely interspersing transducer elements sensitive to the low end of the ultrasonic range.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

TH-CD-202-03: Enhancing Soft-Tissue CT Contrast for Radiation Therapy Using Mono-Energetic Decompositions of Dual Energy CT

Purpose: It is desirable to increase CT soft-tissue contrast to improve delineation of tumor target and/or surrounding organs at risk (OAR) in RT planning and delivery guidance. The purpose of this work is to investigate the use of monoenergetic decompositions obtained from dual energy (DE) CT to improve soft-tissue contrast. Methods: CT data were acquired for 5 prostate and 5 pancreas patients and a phantom with a CT Scanner (Definition AS Open, Siemens) using both sequential DE protocols and standard protocols. For the DE protocols, the scanner rapidly performs two acquisitions at 80 kVp and 140 kVp. The CT numbers of soft tissue inserts in the phantom (CTED/Gammex) were measured across the spectrum of available monoenergetic decompositions (40 to 140 keV) and compared to the standard protocol (120 kVp, 0.6 pitch, 18 mGy CTDIvol). Contrast, defined as the difference in the average CT number between target and OAR, was measured for all subjects and compared between the DE and standard protocols. Results: Mono-energetic decompositions of the phantom demonstrate an enhancement of soft-tissue contrast as the energy is decreased. For instance, relative to the 120 kVp scans the Liver ED insert increased in CT number by 25 HU while the adipose ED insert decreased by 50 HU. The lowest energy decompositions featured the highest contrast between target and OAR. For every patient, the contrast increased by decomposing at 40 keV. The average increase in contrast relative to a 120 kVp scan for prostate patients at 40 keV was 25.05±17.28 HU while for pancreas patients it was 19.21±17.39 HU. Conclusion: Low energy monoenergetic decompositions from dual-energy CT substantially increase soft-tissue contrast. At the lowest achievable monoenergetic decompositions the maximum soft-tissue contrast is achieved and the delineation of target and OAR is improved. Thus it is beneficial to use DECT in radiation oncology. Supported by Siemens.

60 APPLIED LIFE SCIENCES↗

TH-CD-207A-03: A Surface Deformation Driven Respiratory Model for Organ Motion Tracking in Lung Cancer Radiotherapy

Purpose: To propose and validate a novel real-time surface-mesh-based internal organ-external surface motion and deformation tracking method for lung cancer radiotherapy. Methods: Deformation vector fields (DVFs) which characterizes the internal and external motion are obtained by registering the internal organ and tumor contours and external surface meshes to a reference phase in the 4D CT images using a recent developed local topology preserved non-rigid point matching algorithm (TOP). A composite matrix is constructed by combing the estimated internal and external DVFs. Principle component analysis (PCA) is then applied on the composite matrix to extract principal motion characteristics and finally yield the respiratory motion model parameters which correlates the internal and external motion and deformation. The accuracy of the respiratory motion model is evaluated using a 4D NURBS-based cardiac-torso (NCAT) synthetic phantom and three lung cancer cases. The center of mass (COM) difference is used to measure the tumor motion tracking accuracy, and the Dice’s coefficient (DC), percent error (PE) and Housdourf’s distance (HD) are used to measure the agreement between the predicted and ground truth tumor shape. Results: The mean COM is 0.84±0.49mm and 0.50±0.47mm for the phantom and patient data respectively. The mean DC, PE and HD are 0.93±0.01, 0.13±0.03 and 1.24±0.34 voxels for the phantom, and 0.91±0.04, 0.17±0.07 and 3.93±2.12 voxels for the three lung cancer patients, respectively. Conclusions: We have proposed and validate a real-time surface-mesh-based organ motion and deformation tracking method with an internal-external motion modeling. The preliminary results conducted on a synthetic 4D NCAT phantom and 4D CT images from three lung cancer cases show that the proposed method is reliable and accurate in tracking both the tumor motion trajectory and deformation, which can serve as a potential tool for real-time organ motion and deformation monitoring in lung cancer radiotherapy. This work is supported in part by grant from VARIAN MEDICAL SYSTEMS INC, the National Natural Science Foundation of China (no 81428019 and no 81301940), the Guangdong Natural Science Foundation (2015A030313302)and the 2015 Pearl River S&T Nova Program of Guangzhou (201506010096).

60 APPLIED LIFE SCIENCES↗

TH-CD-207B-03: How to Quantify Temporal Resolution in X-Ray MDCT Imaging?

Purpose: In modern CT scanners, a quantitative metric to assess temporal response, namely, to quantify the temporal resolution (TR), remains elusive. Rough surrogate metrics, such as half of the gantry rotation time for single source CT, a quarter of the gantry rotation time for dual source CT, or measurements of motion artifact’s size, shape, or intensity have previously been used. In this work, a rigorous framework which quantifies TR and a practical measurement method are developed. Methods: A motion phantom was simulated which consisted of a single rod that is in motion except during a static period at the temporal center of the scan, termed the TR window. If the image of the motion scan has negligible motion artifacts compared to an image from a totally static scan, then the system has a TR no worse than the TR window used. By repeating this comparison with varying TR windows, the TR of the system can be accurately determined. Motion artifacts were also visually assessed and the TR was measured across varying rod motion speeds, directions, and locations. Noiseless fan beam acquisitions were simulated and images were reconstructed with a short-scan image reconstruction algorithm. Results: The size, shape, and intensity of motion artifacts varied when the rod speed, direction, or location changed. TR measured using the proposed method, however, was consistent across rod speeds, directions, and locations. Conclusion: Since motion artifacts vary depending upon the motion speed, direction, and location, they are not suitable for measuring TR. In this work, a CT system with a specified TR is defined as having the ability to produce a static image with negligible motion artifacts, no matter what motion occurs outside of a static window of width TR. This framework allows for practical measurement of temporal resolution in clinical CT imaging systems. Funding support: GE Healthcare; Conflict of Interest: Employee, GE Healthcare.

60 APPLIED LIFE SCIENCES↗

TH-EF-BRA-03: Assessment of Data-Driven Respiratory Motion-Compensation Methods for 4D-CBCT Image Registration and Reconstruction Using Clinical Datasets

Purpose: To evaluate the performance of a 4D-CBCT registration and reconstruction method that corrects for respiratory motion and enhances image quality under clinically relevant conditions. Methods: Building on previous work, which tested feasibility of a motion-compensation workflow using image datasets superior to clinical acquisitions, this study assesses workflow performance under clinical conditions in terms of image quality improvement. Evaluated workflows utilized a combination of groupwise deformable image registration (DIR) and image reconstruction. Four-dimensional cone beam CT (4D-CBCT) FDK reconstructions were registered to either mean or respiratory phase reference frame images to model respiratory motion. The resulting 4D transformation was used to deform projection data during the FDK backprojection operation to create a motion-compensated reconstruction. To simulate clinically realistic conditions, superior quality projection datasets were sampled using a phase-binned striding method. Tissue interface sharpness (TIS) was defined as the slope of a sigmoid curve fit to the lung-diaphragm boundary or to the carina tissue-airway boundary when no diaphragm was discernable. Image quality improvement was assessed in 19 clinical cases by evaluating mitigation of view-aliasing artifacts, tissue interface sharpness recovery, and noise reduction. Results: For clinical datasets, evaluated average TIS recovery relative to base 4D-CBCT reconstructions was observed to be 87% using fixed-frame registration alone; 87% using fixed-frame with motion-compensated reconstruction; 92% using mean-frame registration alone; and 90% using mean-frame with motion-compensated reconstruction. Soft tissue noise was reduced on average by 43% and 44% for the fixed-frame registration and registration with motion-compensation methods, respectively, and by 40% and 42% for the corresponding mean-frame methods. Considerable reductions in view aliasing artifacts were observed for each method. Conclusion: Data-driven groupwise registration and motion-compensated reconstruction have the potential to improve the quality of 4D-CBCT images acquired under clinical conditions. For clinical image datasets, the addition of motion compensation after groupwise registration visibly reduced artifact impact. This work was supported by the National Cancer Institute of the National Institutes of Health under Award Number R01CA166119. Hugo and Weiss hold a research agreement with Philips Healthcare and license agreement with Varian Medical Systems. Weiss receives royalties from UpToDate. Christensen receives funds from Roger Koch to support research.

60 APPLIED LIFE SCIENCES↗

TH-EF-207A-01: Respiratory Motion Correction in 4D-Multi Pinhole Small Animal SPECT

Purpose: Respiration causes lesion smearing, image blurring and quality degradation, affecting lesion contrast and the ability to define correct lesion size. The spatial resolution of current multi pinhole SPECT (MPHS) scanners is sub-millimeter. Therefore, the effect of motion is more noticeable in comparison to conventional SPECT scanner. Gated imaging aims to reduce motion artifacts. A major issue in gating is the lack of statistics and individual reconstructed frames are noisy. The increased noise in each frame, deteriorates the quantitative accuracy of the MPHS Images. The objective of this work, is to enhance the image quality in 4D-MPHS imaging, by 4D image reconstruction. Methods: The new algorithm requires deformation vector fields (DVFs) that are calculated by non-rigid Demons registration. The algorithm is based on the motion-incorporated version of ordered subset expectation maximization (OSEM) algorithm. This iterative algorithm is capable to make full use of all projections to reconstruct each individual frame. To evaluate the performance of the proposed algorithm a simulation study was conducted. A fast ray tracing method was used to generate MPHS projections of a 4D digital mouse phantom with a small tumor in liver in eight different respiratory phases. To evaluate the 4D-OSEM algorithm potential, tumor to liver activity ratio was compared with other image reconstruction methods including 3D-MPHS and post reconstruction registered with Demons-derived DVFs. Results: Image quality of 4D-MPHS is greatly improved by the 4D-OSEM algorithm. When all projections are used to reconstruct a 3D-MPHS, motion blurring artifacts are present, leading to overestimation of the tumor size and 24% tumor contrast underestimation. This error reduced to 16% and 10% for post reconstruction registration methods and 4D-OSEM respectively. Conclusion: 4D-OSEM method can be used for motion correction in 4D-MPHS. The statistics and quantification are improved since all projection data are combined together to update the image.

60 APPLIED LIFE SCIENCES↗