Engineering topics

Wang, J

Publications and source records attributed to Wang, J.

At least 73 records · Page 4

DOE OSTI2016

WE-AB-207A-06: Progressive Dose Control for Cone Beam CT with Deformation Assisted Temporal Nonlocal Means Method

Purpose: To reduce cone beam CT (CBCT) imaging dose, we previously proposed a progressive dose control (PDC) scheme to employ temporal correlation between CBCT images at different fractions for image quality enhancement. A temporal non-local means (TNLM) method was developed to enhance quality of a new low-dose CBCT using existing high-quality CBCT. To enhance a voxel value, the TNLM method searches for similar voxels in a window. Due to patient deformation among the two CBCTs, a large searching window was required, reducing image quality and computational efficiency. This abstract proposes a deformation-assisted TNLM (DA-TNLM) method to solve this problem. Methods: For a low-dose CBCT to be enhanced using a high-quality CBCT, we first performed deformable image registration between the low-dose CBCT and the high-quality CBCT to approximately establish voxel correspondence between the two. A searching window for a voxel was then set based on the deformation vector field. Specifically, the search window for each voxel was shifted by the deformation vector. A TNLM step was then applied using only voxels within this determined window to correct image intensity at the low-dose CBCT. Results: We have tested the proposed scheme on simulated CIRS phantom data and real patient data. The CITS phantom was scanned on Varian onboard imaging CBCT system with coach shifting and dose reducing for each time. The real patient data was acquired in four fractions with dose reduced from standard CBCT dose to 12.5% of standard dose. It was found that the DA-TNLM method can reduce total dose by over 75% on average in the first four fractions. Conclusion: We have developed a PDC scheme which can enhance the quality of image scanned at low dose using a DA-TNLM method. Tests in phantom and patient studies demonstrated promising results.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

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↗

DOE OSTI2016

WE-FG-202-12: Investigation of Longitudinal Salivary Gland DCE-MRI Changes

Purpose: To determine the correlation between dose and changes through treatment in dynamic contrast enhanced (DCE) MRI voxel parameters (Ktrans, kep, Ve, and Vp) within salivary glands of head and neck oropharyngeal squamous cell carcinoma (HNSCC) patients. Methods: 17 HNSCC patients treated with definitive radiation therapy completed DCE-MRI scans on a 3T scanner at pre-treatment, mid-treatment, and post-treatment time points. Mid-treatment and post-treatment DCE images were deformably registered to pre-treatment DCE images (Velocity software package). Pharmacokinetic analysis of the DCE images used a modified Tofts model to produce parameter maps with an arterial input function selected from each patient’s perivertebral space on the image (NordicICE software package). In-house software was developed for voxel-by-voxel longitudinal analysis of the salivary glands within the registered images. The planning CT was rigidly registered to the pre-treatment DCE image to obtain dose values in each voxel. Voxels within the lower and upper dose quartiles for each gland were averaged for each patient, then an average of the patients’ means for the two quartiles were compared. Dose-relationships were also assessed by Spearman correlations between dose and voxel parameter changes for each patient’s gland. Results: Changes in parameters’ means between time points were observed, but inter-patient variability was high. Ve of the parotid was the only parameter that had a consistently significant longitudinal difference between dose quartiles. The highest Spearman correlation was Vp of the sublingual gland for the change in the pre-treatment to mid-treatment values with only a ρ=0.29. Conclusion: In this preliminary study, there was large inter-patient variability in the changes of DCE voxel parameters with no clear relationship with dose. Additional patients may reduce the uncertainties and allow for the determination of the existence of parameter and dose relationships.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

WE-FG-207B-05: Iterative Reconstruction Via Prior Image Constrained Total Generalized Variation for Spectral CT

Purpose: To investigate iterative reconstruction via prior image constrained total generalized variation (PICTGV) for spectral computed tomography (CT) using fewer projections while achieving greater image quality. Methods: The proposed PICTGV method is formulated as an optimization problem, which balances the data fidelity and prior image constrained total generalized variation of reconstructed images in one framework. The PICTGV method is based on structure correlations among images in the energy domain and high-quality images to guide the reconstruction of energy-specific images. In PICTGV method, the high-quality image is reconstructed from all detector-collected X-ray signals and is referred as the broad-spectrum image. Distinct from the existing reconstruction methods applied on the images with first order derivative, the higher order derivative of the images is incorporated into the PICTGV method. An alternating optimization algorithm is used to minimize the PICTGV objective function. We evaluate the performance of PICTGV on noise and artifacts suppressing using phantom studies and compare the method with the conventional filtered back-projection method as well as TGV based method without prior image. Results: On the digital phantom, the proposed method outperforms the existing TGV method in terms of the noise reduction, artifacts suppression, and edge detail preservation. Compared to that obtained by the TGV based method without prior image, the relative root mean square error in the images reconstructed by the proposed method is reduced by over 20%. Conclusion: The authors propose an iterative reconstruction via prior image constrained total generalize variation for spectral CT. Also, we have developed an alternating optimization algorithm and numerically demonstrated the merits of our approach. Results show that the proposed PICTGV method outperforms the TGV method for spectral CT.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

WE-H-207A-02: Attenuation Correction in 4D-PET Using a Single-Phase Attenuation Map

Purpose: 4D-PET imaging has been proposed as a potential solution to the respiratory motion effect in thoracic region. CT-based attenuation correction (AC) is an essential step toward quantitative imaging for PET. However, due to the temporal difference of 4D-PET and a single breath-hold CT, motion artifacts are observed in the attenuation-corrected PET images that can lead to error in tumor shape and uptake. We introduce a practical method for aligning single-phase CT to all other 4D-PET phases using a penalized non-rigid demons registration. Methods: Individual 4D-PET frames were reconstructed without AC. Non-rigid Demons registration was used to derive deformation vector fields (DVFs) between the PET matched with CT phase and other 4D-PET images. While attenuated PET images provide enough useful data for organ borders such as lung and liver, tumors are not distinguishable from background due to loss of contrast. To preserve tumor shape in different phases, from CT image an ROI covering tumor was excluded from non-rigid transformation. Mean DVF of the central region of the tumor was assigned to all voxels in the ROI. This process mimics a rigid transformation of tumor along with a non-rigid transformation of other organs. 4D XCAT phantom with spherical tumors in lung with diameters ranging from 10 to 40 mm was used to evaluate the algorithm. Results: Motion related induced artifacts in attenuation-corrected 4D-PET images were significantly reduced. For tumors smaller than 20 mm, non-rigid transformation was capable to provide quantitative results. However, for larger tumors, where tumor self-attenuation is considerable, our combined method yields superior results. Conclusion: We introduced a practical method for deforming a single CT to match all 4D-PET images for accurate AC. Although 4D-PET data include insignificant anatomical information, we showed that they are still useful to estimate DVFs for aligning attenuation map and accurate AC.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-CD-BRA-04: Effect of a Strong Magnetic Field On TLDs, OSLDs, and Gafchromic Films Using An Electromagnet

Purpose: To study the effect of strong magnetic field on three types of dosimeters using an electromagnet inside a Linac vault. Materials and Methods: Three types of dosimeters, thermoluminescent Dosimeters (TLDs), optically stimulated luminescent Dosimeters (OSLDs), and EBT3 Film were used to measure radiation dose response inside an electromagnet that could produce a strong magnetic field (B>1.5 T). The dosimeters were placed inside a plastic phantom between the two poles of the magnet, at approximately 3 meters from the iso-center of an Elekta Versa HD Linac. The B field was calibrated with a Gauss meter (Model: GM-2, AlphaLab Inc). The dosimeters received ∼2 Gy with and without the presence of the 1.5 T magnetic field. The EBT3 films were scanned 24 hours before and 24 hours after irradiation. The TLD dosimeters were read 1 week after irradiation. The OSLDs were read two weeks after irradiation. The ratios of signals of dosimeters irradiated with the B field to the signals without the B field were calculated. Two experiments have been conducted so far. Results: The ratios (averaged over two experiments) of dosimeter signals with vs without B field were 0.994 for films, 0.994 for OSLDs, and 1.002 for TLDs. The statistical uncertainty was ∼3%. Conclusions: The three types of dosimeters (film, TLD, OSLD) seem not affected by the presence of a magnetic field (B=1.5 T) with the uncertainty of ∼3%. They may be suitable for QA purposes in a strong B field up to 1.5 T. More measurements will be conducted for reproducibility testing. We acknowledge research support from Elekta AB.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-CD-BRA-08: Novel Iron-Based Radiation Reporting Systems as 4D Dosimeters for MR-Guided Radiation Therapy

Purpose: To compare novel radiation reporting systems utilizing ferric ion (Fe{sup 3+}) reduction versus ferrous ion (Fe{sup 2+}) oxidation in gelatin matrixes for 3D and 4D (3D+time) MR-guided radiation therapy dosimetry. Methods: Dosimeters were irradiated using an integrated 1.5T MRI and 7MV linear accelerator (MR-Linac). Dosimeters were read-out with both a spectrophotometer and the MRI component of the MR-Linac immediately after irradiation. Changes in optical density (OD) were measured using a spectrophotometer; changes in MR signal intensity due to the paramagnetic differences in the iron ions were measured using the MR-Linac in real-time during irradiation (balanced-FFE sequences) and immediately after irradiation (T{sub 1}-weighted and inversion recovery sequences). Results: Irradiation of Fe{sup 3+} reduction dosimeters resulted in a stable red color with an absorbance peak at 512 nm. The change in OD relative to dose exhibited a linear response up to 100 Gy (R{sup 2}=1.00). T{sub 1}-weighted-MR signal intensity (SI) changed minimally after irradiation with increases of 8.0% for 17 Gy and 9.7% after escalation to 35 Gy compared to the un-irradiated region. Irradiation of Fe{sup 2+} oxidation dosimeters resulted in a stable purple color with absorbance peaks at 440 and 585 nm. The changes in OD, T{sub 1}-weighted-MR SI, and R{sub 1} relative to dose exhibited a linear response up to at least 8 Gy (R{sup 2}=1.00, 0.98, and 0.99) with OD saturation above 40 Gy. The T{sub 1}-weighted-MR SI increased 50.3% for 17 Gy compared to the un-irradiated region. The change in SI was observed in both 2D+time and 4D (3D+time) acquisitions post-irradiation and in real-time during irradiation with a linear increase with respect to dose (R{sup 2}>0.93). Conclusion: The Fe{sup 2+} oxidation-based system was superior as 4D dosimeters for MR-guided radiation therapy due to its higher sensitivity in both optical and MR signal readout and feasibility for real-time 4D dose readout. The Fe{sup 3+} reduction system is recommended for high dose applications. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant No. LH-102SPS.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-CD-207A-05: Lung Surface Deformation Vector Fields Prediction by Monitoring Respiratory Surrogate Signals

Purpose: In this study, we developed and evaluated a method for predicting lung surface deformation vector fields (SDVFs) based on surrogate signals such as chest and abdomen motion at selected locations and spirometry measurements. Methods: A Patient-specific 3D triangular surface mesh of the lung region at end-expiration (EE) phase was obtained by threshold-based segmentation method. For each patient, a spirometer recorded the flow volume changes of the lungs; and 192 selected points at a regular spacing of 2cm X 2cm matrix points over a total area of 34cm X 24cm on the surface of chest and abdomen was used to detect chest wall motions. Preprocessing techniques such as QR factorization with column pivoting (QRCP) were employed to remove redundant observations of the chest and abdominal area. To create a statistical model between the lung surface and the corresponding surrogate signals, we developed a predictive model based on canonical ridge regression (CRR). Two unique weighting vectors were selected for each vertex on the surface of the lung, and they were optimized during the training process using the all other phases of 4D-CT except the end-inspiration (EI) phase. These parameters were employed to predict the vertices locations of a testing data set, which was the EI phase of 4D-CT. Results: For ten lung cancer patients, the deformation vector field of each vertex of lung surface mesh was estimated from the external motion at selected positions on the chest wall surface plus spirometry measurements. The average estimation of 98th percentile of error was less than 1 mm (AP= 0.85, RL= 0.61, and SI= 0.82). Conclusion: The developed predictive model provides a non-invasive approach to derive lung boundary condition. Together with personalized biomechanical respiration modelling, the proposed model can be used to derive the lung tumor motion during radiation therapy accurately from non-invasive measurements.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TH-EF-BRA-04: Individually Optimized Contrast-Enhanced 4D-CT for Radiotherapy Simulation in Pancreatic Ductal Adenocarcinoma

Purpose: To develop an individually optimized contrast-enhanced (CE) 4D-CT for radiotherapy simulation in pancreatic ductal adenocarcinomas (PDA). Methods: Ten PDA patients were enrolled. Each underwent 3 CT scans: a 4D-CT immediately following a CE 3D-CT and an individually optimized CE 4D-CT using test injection. Three physicians contoured the tumor and pancreatic tissues. We compared image quality scores, tumor volume, motion, tumor-to-pancreas contrast, and contrast-to-noise ratio (CNR) in the 3 CTs. We also evaluated interobserver variations in contouring the tumor using simultaneous truth and performance level estimation (STAPLE). Results: Average image quality scores for CE 3DCT and CE 4D-CT were comparable (4.0 and 3.8, respectively; P=0.47), and both were significantly better than that for 4D-CT (2.6, P<0.001). Tumor-to-pancreas contrast results were comparable in CE 3D-CT and CE 4D-CT (15.5 and 16.7 HU, respectively; P=0.71), and the latter was significantly higher than in 4D-CT (9.2 HU, P=0.03). Image noise in CE 3D-CT (12.5 HU) was significantly lower than in CE 4D-CT (22.1 HU, P<0.001) and 4D-CT (19.4 HU, P=0.005). CNRs were comparable in CE 3D-CT and CE 4DCT (1.4 and 0.8, respectively; P=0.23), and the former was significantly better than in 4D-CT (0.6, P = 0.04). Mean tumor volumes were smaller in CE 3D-CT (29.8 cm{sup 3}) and CE 4D-CT (22.8 cm{sup 3}) than in 4D-CT (42.0 cm{sup 3}), although these differences were not statistically significant. Mean tumor motion was comparable in 4D-CT and CE 4D-CT (7.2 and 6.2 mm, P=0.23). Interobserver variations were comparable in CE 3D-CT and CE 4D-CT (Jaccard index 66.0% and 61.9%, respectively) and were worse for 4D-CT (55.6%) than CE 3D-CT. Conclusion: CE 4D-CT demonstrated characteristics comparable to CE 3D-CT, with high potential for simultaneously delineating the tumor and quantifying tumor motion with a single scan. Supported in part by Philips Healthcare.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

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↗

DOE OSTI2016

TH-EF-207A-05: Feasibility of Applying SMEIR Method On Small Animal 4D Cone Beam CT Imaging

Purpose: Small animal cone beam CT imaging has been widely used in preclinical research. Due to the higher respiratory rate and heat beats of small animals, motion blurring is inevitable and needs to be corrected in the reconstruction. Simultaneous motion estimation and image reconstruction (SMEIR) method, which uses projection images of all phases, proved to be effective in motion model estimation and able to reconstruct motion-compensated images. We demonstrate the application of SMEIR for small animal 4D cone beam CT imaging by computer simulations on a digital rat model. Methods: The small animal CBCT imaging system was simulated with the source-to-detector distance of 300 mm and the source-to-object distance of 200 mm. A sequence of rat phantom were generated with 0.4 mm{sup 3} voxel size. The respiratory cycle was taken as 1.0 second and the motions were simulated with a diaphragm motion of 2.4mm and an anterior-posterior expansion of 1.6 mm. The projection images were calculated using a ray-tracing method, and 4D-CBCT were reconstructed using SMEIR and FDK methods. The SMEIR method iterates over two alternating steps: 1) motion-compensated iterative image reconstruction by using projections from all respiration phases and 2) motion model estimation from projections directly through a 2D-3D deformable registration of the image obtained in the first step to projection images of other phases. Results: The images reconstructed using SMEIR method reproduced the features in the original phantom. Projections from the same phase were also reconstructed using FDK method. Compared with the FDK results, the images from SMEIR method substantially improve the image quality with minimum artifacts. Conclusion: We demonstrate that it is viable to apply SMEIR method to reconstruct small animal 4D-CBCT images.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Observations and Modeling of the Green Ocean Amazon 2014/15. Nanoparticle Size Distribution (NPSD) Field Campaign Report

Aerosol nucleation and initial growth were investigated during the Green Ocean Amazon (GoAmazon) 2014/15 campaign. Aerosol sampling occurred during the wet and dry seasons of 2014, and took place at the T3 measurement site, downwind of the city of Manaus, Brazil. Characterization of the aerosol size distribution from 10 to 500 nm was accomplished through the deployment of a conventional Scanning Mobility Particle Spectrometer (SMPS) and a fine condensation particle counter (> 10 nm). In order to directly measure aerosol nucleation and initial growth, a Nano SMPS (1.5-20 nm) was also deployed, consisting of a condensation particle counter-based electrical mobility spectrometer that was modified for the detection of sub-3 nm aerosol. Measurements of the aerosol size distribution from 1.5 nm to 10 nm were obtained during the first observational period, and from 3 nm to 15 nm during the second observational period. Routine, stable measurement in this size range was complicated due to persistent water condensation in the Nano SMPS and diffusional transport losses

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016