Search NASA⌕ Search

SEARCH · Search NASA

Results for “image simulation”

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,783 records · Page 99

Optimization of laser-driven cylindrical implosions on the OMEGA laser

Laser-driven cylindrical implosions were conducted on the OMEGA laser as part of the laser-driven mini-MagLIF (Magnetized Liner Inertial Fusion) Campaign. Gated x-ray images were analyzed to infer shell trajectories and study the energy coupling in these implosions. Two-dimensional and three-dimensional HYDRA simulations were performed and post-processed to produce synthetic x-ray self-emission images for comparison. An analysis technique, which could be applied to both experimental and simulated x-ray images, was developed to characterize the shape and uniformity of the implosion. The analysis leads to a measurement of the average implosion velocity and axial implosion length, which can then be used to optimize the beam pointing and energy balance for future experiments. Discrepancies between simulation results and experiments allude to important physical processes that are not accounted for in the simulations. In 2-D simulations, the laser beam’s azimuthal angle of incidence is not included because the $\phi$-direction is not simulated, thus energy absorption is over-predicted. The 3-D simulation results are more consistent with the experiments, but the simulations do not include a calculation of crossbeam energy transfer or non-local thermal transport, which affect the energy coupled to the implosion. Finally, by appropriately adjusting the simulated energy balance and flux limit, the simulations can accurately model the experiments, which have achieved uniform implosions over a 700-μm-long region at velocities of approximately 200 km=s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using [{sup 18}F]Fluorothymidine Imaged With Positron Emission Tomography to Quantify and Reduce Hematologic Toxicity Due to Chemoradiation Therapy for Pelvic Cancer Patients

Purpose: The purpose of the present prospective clinical trial was to determine the efficacy of [{sup 18}F]fluorothymidine (FLT)-identified active bone marrow sparing for pelvic cancer patients by correlating the FLT uptake change during and after chemoradiation therapy with hematologic toxicity. Methods and Materials: Simulation FLT positron emission tomography (PET) images were used to spare pelvic bone marrow using intensity modulated radiation therapy (IMRT BMS) for 32 patients with pelvic cancer. FLT PET scans taken during chemoradiation therapy after 1 and 2 weeks and 30 days and 1 year after completion of chemoradiation therapy were used to evaluate the acute and chronic dose response of pelvic bone marrow. Complete blood counts were recorded at each imaging point to correlate the FLT uptake change with systemic hematologic toxicity. Results: IMRT BMS plans significantly reduced the dose to the pelvic regions identified with FLT uptake compared with control IMRT plans (P<.001, paired t test). Radiation doses of 4 Gy caused an ∼50% decrease in FLT uptake in the pelvic bone marrow after either 1 or 2 weeks of chemoradiation therapy. Additionally, subjects with more FLT-identified bone marrow exposed to ≥4 Gy after 1 week developed grade 2 leukopenia sooner than subjects with less marrow exposed to ≥4 Gy (P<.05, Cox regression analysis). Apparent bone marrow recovery at 30 days after therapy was not maintained 1 year after chemotherapy. The FLT uptake in the pelvic bone marrow regions that received >35 Gy was 18.8% ± 1.8% greater at 30 days after therapy than at 1 year after therapy. The white blood cell, platelet, lymphocyte, and neutrophil counts at 1 year after therapy were all lower than the pretherapy levels (P<.05, paired t test). Conclusions: IMRT BMS plans reduced the dose to FLT-identified pelvic bone marrow for pelvic cancer patients. However, reducing hematologic toxicity is challenging owing to the acute radiation sensitivity (∼4 Gy) and chronic suppression of activity in bone marrow receiving radiation doses >35 Gy, as measured by the FLT uptake change correlated with the complete blood cell counts.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Mapping of deciduous forest cover using simulated Landsat-D TM data

An evaluation is presented of the supervised vegetation classification images of heavily forested terrain in the eastern US produced from simulated Landsat-D Thematic Mapper data as part of the Joint NASA/Geosat Test Case study of the Lost River, West Virginia, gas field. This investigation utilized data supplied by the NS-001 aircraft multispectral scanner (15 m IFOV) of NASA. The instruments recognized a total of 9 vegetation classes and one soil class. A hybrid parallelepiped-Bayesian maximum likelihood classifier was employed to produce the supervised vegetation classifications. The data were obtained at the height of the fall foliage display for this portion of West Virginia. Results show that both classifications are highly accurate, based on a comparison of assignment of vegetation classes on each of the images with actual ground conditions. It is noted that the anomalous distribution of certain woody species found in the resulting supervised classifications may be related to the microseepage of methane from the gas reservoir.

Rock, B. N.↗

Effects of low earth orbit on the optical performance of multi-layer enhanced high reflectance mirrors

Two mirror designs developed for space applications were flown along with a standard mid-infrared design on the leading and trailing edges of the Long Duration Exposure Facility (LDEF). Preliminary observations of induced changes in optical performance of ZnS-coated mirrors and impact-related microstructural and microchemical effects are described in the proceedings of the First LDEF Post-Retrieval Symposium. In this paper, effects of the induced environment and meteoroid/debris impacts on mirror performance are described in more detail. Also, an analysis of reflectance spectra using the results of Auger and secondary ion mass spectroscopy (SIMS) profiling measurements are used to identify an optical-degradation mechanism for the ZnS-coated mirrors. Structural damage associated with a high-velocity impact on a (Si/Al2O3)-coated mirror was imaged optically and with scanning electron and atomic force microscopy (SEM and AFM). Scanning Auger and SIMS analysis provided chemical mapping of selected impact sites. The impact data suggest design and fabrication modifications for obtaining improved mechanical performance using a design variation identified in preflight laboratory simulations. Auger surface profile and SIMS imaging data verified the conclusion that secondary impacts are the source of contamination associated with the dendrites grown on the leading-edge ZnS-coated test samples. It was also found that dendrites can be grown in the laboratory by irradiating contaminated sites on a trailing-edge ZnS-coated sample with a rastered electron beam. These results suggest a mechanism for dendrite growth.

Donovan, Terence↗

Comparison of Satellite-Derived TOA Shortwave Clear-Sky Fluxes to Estimates from GCM Simulations Constrained by Satellite Observations of Land Surface Characteristics

Clear-sky, upwelling shortwave flux at the top of the atmosphere (S(sub TOA raised arrow)), simulated using the atmospheric and land model components of the Community Climate System Model 3 (CCSM3), is compared to corresponding observational estimates from the Clouds and Earth's Radiant Energy System (CERES) sensor. Improvements resulting from the use of land surface albedo derived from Moderate Resolution Imaging Spectroradiometer (MODIS) to constrain the simulations are also examined. Compared to CERES observations, CCSM3 overestimates global, annual averaged S(sub TOA raised arrow) over both land and oceans. However, regionally, CCSM3 overestimates S(sub TOA raised arrow) over some land and ocean areas while underestimating it over other sites. CCSM3 underestimates S(sub TOA raised arrow) over the Saharan and Arabian Deserts and substantial differences exist between CERES observations and CCSM3 over agricultural areas. Over selected sites, after using groundbased observations to remove systematic biases that exist in CCSM computation of S(sub TOA raised arrow), it is found that use of MODIS albedo improves the simulation of S(sub TOA raised arrow). Inability of coarse resolution CCSM3 simulation to resolve spatial heterogeneity of snowfall over high altitude sites such as the Tibetan Plateau causes overestimation of S(sub TOA raised arrow) in these areas. Discrepancies also exist in the simulation of S(sub TOA raised arrow) over ocean areas as CCSM3 does not account for the effect of wind speed on ocean surface albedo. This study shows that the radiative energy budget at the TOA is improved through the use of MODIS albedo in Global Climate Models.

Anantharaj, Valentine G.↗

Visualizing the impact of quenched disorder on 2D electron Wigner solids

Electron Wigner solids (WSs) provide an ideal system for understanding the competing effects of electron–electron and electron–disorder interactions, a central unsolved problem in condensed matter physics. Progress in this topic has been limited by a lack of single-defect-resolved experimental measurements as well as accurate theoretical tools to enable realistic experiment/theory comparison. Here, in this study, we overcome these limitations by combining atomically resolved scanning tunnelling microscopy (STM) with neural-quantum-state quantum Monte Carlo (NQS-QMC) simulation of disordered 2D electron WSs to discover new disorder-induced physical regimes of correlated electron behaviour. STM was used to image the electron density (n e )-dependent evolution of electron WSs in gate-tunable bilayer MoSe 2 (BL-MoSe 2 ) devices with varying long-range (n LR ) and short-range (n SR ) disorder densities. These images were compared with NQS-QMC simulations using realistic disorder maps extracted from experiment, thus allowing the roles of different disorder types to be disentangled. We identify two distinct physical regimes for disordered electron WSs that depend on n SR . For n SR ≲ n e , the WS behaviour is dominated by long-range disorder and features extensive mixed solid–liquid phases, a new type of local re-entrant melting/crystallization and prominent Friedel oscillations. By contrast, when n SR ≫ n e , these features are suppressed and a more robust amorphous WS phase emerges that persists to higher n e , highlighting the importance of short-range disorder in this regime. Our work establishes a powerful framework for studying disordered quantum solids through a combined experimental–theoretical approach.

Ge, Zhehao [University of California, Berkeley, CA↗

SU-E-T-457: Impact of Interfractional Variations On Anterior Vs. Lateral-Field Proton Therapy of Prostate Cancer

Purpose: To investigate the effects of interfractional anatomy and setup variations on plans with anterior-oblique vs. lateral beams for prostate cancer pencil beam scanning (PBS) and passive scattered (PS) proton therapy. Methods: Six patients with low/intermediate risk prostate cancer treated with PS proton therapy at our institution were selected. All patients underwent weekly verification CT scans. Implanted fiducials were used for localization, and endorectal balloons for prostate immobilization. New PBS plans with lateral beams, as well as PBS and PS plans with anterior-oblique beams (±35 deg) were created. PBS plans used two different spot sizes: ∼10mm (large) and ∼5mm (medium) sigma at 25cm range and optimized as single-field-uniform-dose with ∼8% non-uniformity. No range uncertainty margins were applied in PBS plans to maximize rectal sparing. Field-specific apertures were used when planning with large spots to sharpen the penumbrae. The planned dose was recomputed on each weekly CT with fiducials aligned to the simulation CT, scaled and accumulated via deformable image registration. Results: The dose volume analysis showed that although difference between planned and accumulated dose remains negligible for plans with conventional lateral beams using both PS and PBS, this is not the case for plans with anterior beams. The target coverage in anterior plans was largely degraded due to the variations in the beam path length and the absence of range margins. The average prostate D95 was reduced by 7.5/15.9% (using PS/PBS) after accumulation for anterior plans, compared with 0/0.4% for lateral plans. The average mean dose in organs-at-risk decreased by 1% for lateral and 2% for anterior plans, similarly for PS and PBS. Spot size did not affect the dose changes. Conclusion: Prostate plans using anterior beams may undergo clinically relevant interfractional dose degradation. Corrective strategies guided by in-vivo range measurements should be studied before clinical application of this technique.

60 APPLIED LIFE SCIENCES↗

Downscaling-Based Segmentation for Unresolved Images of Highly Heterogeneous Granular Porous Samples

Numerical simulations of pore-scale flow and transport in natural sediments require the knowledge of pore-space topology. Limited resolution of X-ray tomography is often insufficient to fully characterize pore-space structure within fine-grained regions. Single and multilevel threshold-based segmentation approaches are customarily employed to identify solid, pore and porous-solid regions by means of grey intensity thresholds. While the choice of cutoff thresholds is often arbitrary, it dramatically affects the effective properties and the dynamical response of the reconstructed porous structure. As such, we propose an algorithm of downscaling, i.e., the process of increasing image resolution, followed by segmentation, i.e., the identification of different phases, to reconstruct the unresolved pore-space from XCT images of natural geological porous media. The method, applicable to moderately unresolved, chemically homogeneous granular media, is based on a map between local pixel porosity and pore size that does not rely on the definition of arbitrary thresholds and it allows to generate a high-resolution binary image of the porous medium from poorly resolved grey-scale images. First, we validate the method on synthetic unresolved images and compare their known pore-space distribution with the extracted one. Then, we consider a synthetic porous medium and compare the pore size distribution, conductivity, and tortuosity between the original and the reconstructed structures. Finally, we apply the method to extract the pore-space distribution from unresolved XCT images of two natural sediment columns and use it (i) to parametrize a capillary-bundle model and (ii) to estimate the hydraulic conductivity by matching breakthrough behavior of passive solute transport.

58 GEOSCIENCES↗

Online Adaptive Radiation Therapy

The current paradigm of radiation therapy has the treatment planned on a snapshot dataset of the patient's anatomy taken at the time of simulation. Throughout the course of treatment, this snapshot may vary from initial simulation. Although there is the ability to image patients within the treatment room with technologies such as cone beam computed tomography, the current state of the art is largely limited to rigid-body matching and not accounting for any geometric deformations in the patient's anatomy. A plan that was once attuned to the initial simulation can become suboptimal as the treatment progresses unless improved technologies are brought to bear. Adaptive radiation therapy (ART) is an evolving paradigm that seeks to address this deficiency by accounting for ongoing changes in the patient's anatomy and/or physiology during the course of treatment, affording an increasingly more accurate targeting of disease. ART relies on several components working in concert, namely in-room treatment image guidance, deformable image registration, automatic recontouring, plan evaluation and reoptimization, dose calculation, and quality assurance. Various studies have explored how a putative ART solution would improve the current state of the art of radiation therapy—some centers have even clinically implemented online adaptation. These explorations are reviewed here for a variety of sites.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

The Gemini NICI planet-finding campaign: The companion detection pipeline

We present high-contrast image processing techniques used by the Gemini NICI Planet-Finding Campaign to detect faint companions to bright stars. The Near-Infrared Coronographic Imager (NICI) is an adaptive optics instrument installed on the 8 m Gemini South telescope, capable of angular and spectral difference imaging and specifically designed to image exoplanets. The Campaign data pipeline achieves median contrasts of 12.6 mag at 0.''5 and 14.4 mag at 1'' separation, for a sample of 45 stars (V = 4.3-13.9 mag) from the early phase of the campaign. We also present a novel approach to calculating contrast curves for companion detection based on 95% completeness in the recovery of artificial companions injected into the raw data, while accounting for the false-positive rate. We use this technique to select the image processing algorithms that are more successful at recovering faint simulated point sources. We compare our pipeline to the performance of the Locally Optimized Combination of Images (LOCI) algorithm for NICI data and do not find significant improvement with LOCI.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of breast tissue composition with dual energy cone-beam computed tomography: A postmortem study

Purpose: To investigate the feasibility of a three-material compositional measurement of water, lipid, and protein content of breast tissue with dual kVp cone-beam computed tomography (CT) for diagnostic purposes. Methods: Simulations were performed on a flat panel-based computed tomography system with a dual kVp technique in order to guide the selection of experimental acquisition parameters. The expected errors induced by using the proposed calibration materials were also estimated by simulation. Twenty pairs of postmortem breast samples were imaged with a flat-panel based dual kVp cone-beam CT system, followed by image-based material decomposition using calibration data obtained from a three-material phantom consisting of water, vegetable oil, and polyoxymethylene plastic. The tissue samples were then chemically decomposed into their respective water, lipid, and protein contents after imaging to allow direct comparison with data from dual energy decomposition. Results: Guided by results from simulation, the beam energies for the dual kVp cone-beam CT system were selected to be 50 and 120 kVp with the mean glandular dose divided equally between each exposure. The simulation also suggested that the use of polyoxymethylene as the calibration material for the measurement of pure protein may introduce an error of -11.0%. However, the tissue decomposition experiments, which employed a calibration phantom made out of water, oil, and polyoxymethylene, exhibited strong correlation with data from the chemical analysis. The average root-mean-square percentage error for water, lipid, and protein contents was 3.58% as compared with chemical analysis. Conclusions: The results of this study suggest that the water, lipid, and protein contents can be accurately measured using dual kVp cone-beam CT. The tissue compositional information may improve the sensitivity and specificity for breast cancer diagnosis.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Evaluation of EPIC Oxygen Bands Stability With Radiative Transfer Simulations Over the South Pole

The Earth Polychromatic Imaging Camera (EPIC) onboard the Deep Space Climate Observatory (DSCOVR) satellite orbiting the Sun at the Lagrange-1 point was launched without onboard calibration systems. Vicarious calibration is conducted for 8 of the 10 UV/VIS/NIR channels using other low earth orbiting satellite instruments, while its two O 2 bands (688 nm and 764nm) rely on indirect moon-view calibrations because the same narrow-band O 2 bands are not readily available from other in-flight instruments. This study compares EPIC measurements from the four O 2 bands aiming at examining sensor stability over a uniquely suited location, i.e., the permanently snow-covered South Pole. The study utilizes radiative transfer model simulations with in-situ atmospheric soundings taken at South Pole during months of December and January from 2015 to 2022. The absolute discrepancy between the model simulations and observations is less than 1.0% for the two reference bands, but 5.75% and 15.63% for the 688nm, and 764 nm absorption bands, respectively. The simulated A-band and B-band ratios are 16.09% and 4.74% higher than that from the observations. Various sensitivities are conducted to estimate possible contributions to the discrepancies from input atmospheric profiles, spectral surface albedos and surface BRDF. While none of the input uncertainties is likely to account for the large discrepancies in the oxygen absorption bands, a small shift in the instrument response function could be the main reason for these biases. On the other hand, the model simulations are able to capture systematic variations with observed angular measurements and explain the multi-year trends found in observed O 2 band ratios due to satellite orbit shifting. When model simulated contributions from the angle variations are deducted from the observed O 2 band ratios, the residual O 2 band ratios are found to be stable since 2015.

EPIC↗

3D hybrid fluid-particle jet simulations and the importance of synchrotron radiative losses

Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization emission. Aims. Through the use of 3D hybrid fluid-particle jet simulations (with the PLUTO code), we overcome some of the commonly made assumptions in relativistic magnetohydrodynamic (RMHD) simulations by using non-thermal particle attributes to account for the resulting synchrotron radiation. Polarized radiative transfer and ray-tracing (via the RADMC-3D code) highlight the differences in total intensity maps when (i) the jet is simulated purely with the RMHD approach, (ii) a jet tracer is considered in the RMHD approach, and (iii) a hybrid fluid-particle approach is used. The resulting emission maps were compared to the example of the radio galaxy Centaurus A. Methods. We applied the Lagrangian particle module implemented in the latest version of the PLUTO code. This new module contains a state-of-the-art algorithm for modeling diffusive shock acceleration and for accounting for radiative losses in RMHD jet simulations. The module implements the physical postulates missing in RMHD jet simulations by accounting for a cooled ambient medium and strengthening the central jet emission. Results. We find a distinction between the innermost structure of the jet and the back-flowing material by mimicking the radio emission of the Seyfert II radio galaxy Centaurus A when considering an edge-brightened jet with an underlying purely toroidal magnetic field. We demonstrate the necessity of synchrotron cooling as well as the improvements gained when directly accounting for non-thermal synchrotron radiation via non-thermal particles.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

The LSST DESC data challenge 1: generation and analysis of synthetic images for next-generation surveys

Data Challenge 1 (DC1) is the first synthetic data set produced by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). DC1 is designed to develop and validate data reduction and analysis and to study the impact of systematic effects that will affect the LSST data set. DC1 is comprised of r -band observations of 40 deg 2 to 10 yr LSST depth. In this paper, we present each stage of the simulation and analysis process: (a) generation, by synthesizing sources from cosmological N -body simulations in individual sensor-visit images with different observing conditions; (b) reduction using a development version of the LSST Science Pipelines; and (c) matching to the input cosmological catalogue for validation and testing. We verify that testable LSST requirements pass within the fidelity of DC1. We establish a selection procedure that produces a sufficiently clean extragalactic sample for clustering analyses and we discuss residual sample contamination, including contributions from inefficiency in star–galaxy separation and imperfect deblending. We compute the galaxy power spectrum on the simulated field and conclude that: (i) survey properties have an impact of 50 per cent of the statistical uncertainty for the scales and models used in DC1; (ii) a selection to eliminate artefacts in the catalogues is necessary to avoid biases in the measured clustering; and (iii) the presence of bright objects has a significant impact (2σ–6σ) in the estimated power spectra at small scales (ℓ > 1200), highlighting the impact of blending in studies at small angular scales in LSST.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Solute atmospheres at dislocations

In this study, a two-dimensional plane strain elastic solution is determined for the Cottrell solute atmosphere around an edge dislocation in an infinitely long cylinder of finite radius (the matrix), in which rows of solutes are represented by cylindrical rods with in-plane hydrostatic misfit (axial misfit is also considered). The periphery of the matrix is traction-free, thus introducing an image solute field which generates a solute-solute interaction energy that has not been considered previously. The relevant energy for the field of any distribution of solutes coexistent with a single edge dislocation along the (matrix) cylinder axis is determined, and coherency effects are discussed and studied. Monte Carlo simulations accounting for all pertinent interactions over a range of temperatures are found to yield solute distributions different from classical results, namely, (1) Fermi-Dirac condensations at low temperatures at the free surface, (2) the majority of the atmosphere lying within an unexpectedly large non-linear interaction region near the dislocation core, and (3) temperature-dependent asymmetrical solute arrangements that promote bending. The solute distributions at intermediate temperatures show a 1/r dependence in agreement with previous linearized approximations. With a standard state of solute corresponding to a mean concentration, c 0 , the relevant interaction energy expression presented in this work is valid when extended to large concentrations for which Henry's Law and Vegard's Law do not apply.

36 MATERIALS SCIENCE↗

A novel approach to assess the treatment response using Gaussian random field in PET

Purpose: The assessment of early therapeutic response to anticancer therapy is vital for treatment planning and patient management in clinic. With the development of personal treatment plan, the early treatment response, especially before any anatomically apparent changes after treatment, becomes urgent need in clinic. Positron emission tomography (PET) imaging serves an important role in clinical oncology for tumor detection, staging, and therapy response assessment. Many studies on therapy response involve interpretation of differences between two PET images, usually in terms of standardized uptake values (SUVs). However, the quantitative accuracy of this measurement is limited. This work proposes a statistically robust approach for therapy response assessment based on Gaussian random field (GRF) to provide a statistically more meaningful scale to evaluate therapy effects. Methods: The authors propose a new criterion for therapeutic assessment by incorporating image noise into traditional SUV method. An analytical method based on the approximate expressions of the Fisher information matrix was applied to model the variance of individual pixels in reconstructed images. A zero mean unit variance GRF under the null hypothesis (no response to therapy) was obtained by normalizing each pixel of the post-therapy image with the mean and standard deviation of the pretherapy image. The performance of the proposed method was evaluated by Monte Carlo simulation, where XCAT phantoms (128{sup 2} pixels) with lesions of various diameters (2–6 mm), multiple tumor-to-background contrasts (3–10), and different changes in intensity (6.25%–30%) were used. The receiver operating characteristic curves and the corresponding areas under the curve were computed for both the proposed method and the traditional methods whose figure of merit is the percentage change of SUVs. The formula for the false positive rate (FPR) estimation was developed for the proposed therapy response assessment utilizing local average method based on random field. The accuracy of the estimation was validated in terms of Euler distance and correlation coefficient. Results: It is shown that the performance of therapy response assessment is significantly improved by the introduction of variance with a higher area under the curve (97.3%) than SUVmean (91.4%) and SUVmax (82.0%). In addition, the FPR estimation serves as a good prediction for the specificity of the proposed method, consistent with simulation outcome with ∼1 correlation coefficient. Conclusions: In this work, the authors developed a method to evaluate therapy response from PET images, which were modeled as Gaussian random field. The digital phantom simulations demonstrated that the proposed method achieved a large reduction in statistical variability through incorporating knowledge of the variance of the original Gaussian random field. The proposed method has the potential to enable prediction of early treatment response and shows promise for application to clinical practice. In future work, the authors will report on the robustness of the estimation theory for application to clinical practice of therapy response evaluation, which pertains to binary discrimination tasks at a fixed location in the image such as detection of small and weak lesion.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-F-T-645: To Test Spatial Anddosimetric Accuracy of Small Cranial Target Irradiation Based On 1.5 T MRIscans Using Static Arcs with MLCDefined Fields

Purpose: To test spatial and dosimetric accuracy of small cranial target irradiation based on 1.5 T MRI scans using static arcs with MLC-defined fields Methods: A plastic (PMMA) phantom simulating a small brain lesion was mounted on a GammaKnife headframe equipped with MRI localizer. The lesion was a 3 mm long, 3.175 mm diameter cylindrical cavity filled with MRI contrast. Radiochromic film passing through the cavity was marked with pin pricks at the cavity center. The cavity was contoured on an MRI image and fused with CT to simulate treatment of a lesion not visible on CT. The transfer of the target to CT involved registering the MRI contrast cannels of the localizer that were visible on both modalities. Treatments were planned to deliver 800 cGy to the cavity center using multiple static arcs with 5.0×2.4 mm MLC-defined fields. The phantom was aligned on a STx accelerator by registering the conebeam CT with the planning CT. Films from coronal and sagittal planes were scanned and evaluated using ImageJ software Results: Geographic errors in treatment based on 1.5 T scans agreed within 0.33, −0.27 and 1.21 mm in the vertical, lateral and longitudinal dimensions, respectively. The doses delivered to the cavity center were 7.2% higher than planned. The dose distributions were similar to those of a GammaKnife. Conclusion: Radiation can be delivered with an accelerator at mm accuracy to small cranial targets based on 1.5 MRI scans fused to CTs using a standard GammaKnife headframe and MRI localizer. MLC-defined static arcs produce isodose lines very similar to the GammaKnife.

60 APPLIED LIFE SCIENCES↗