Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linear interpolation”

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 127 records · Page 7

Theory for Glassy Behavior of Supercooled Liquid Mixtures

Here, we introduce a model for glassy dynamics in supercooled liquid mixtures. Given the relaxation behavior of individual supercooled liquids, the model predicts the relaxation times of their mixtures as temperature is decreased. The model is based on dynamical facilitation theory for glassy dynamics, which provides a physical basis for relaxation and vitrification of a supercooled liquid. This is in contrast to empirical linear interpolations such as the Gordon-Taylor equation typically used to predict glass transition temperatures of liquid mixtures. To understand the behavior of supercooled liquid mixtures we consider a multicomponent variant of the kinetically constrained East model in which components have a different energy scale and can also diffuse when locally mobile regions, i.e., excitations, are present. Using a variational approach we confirm an effective single component model with a single effective energy scale that best approximates a mixture. When scaled by this single effective energy, we show that experimental relaxation times of many liquid mixtures all collapse onto the "parabolic law" predicted by dynamical facilitation theory. The model can be used to predict transport properties and glass transition temperatures of mixtures of glassy materials, with implications in atmospheric chemistry, biology, and pharmaceuticals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Adaptive PID Gain Scheduling Control for Hydropower Turbine Using Neural CDE and Stochastic Distribution Shaping

This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.

13 HYDRO ENERGY↗

Effective Missing Value Imputation Methods for Building Monitoring Data

To understand behaviors of natural and man-made events, such as energy consumption of buildings, which accounts for 40% of energy uses in the US, we deploy automated monitoring devices to record periodic observations. However, such experimental and observation data often contains problems and irregularities that have to be cleaned up before analyses. Due to various conditions affecting sensor operations, the communication channels, recording steps, or the recording media, the recorded data might have missing values, errors, or anomalous values. An effective way to clean up these problems is to replace these missing values, errors and anomalous values with expected values, a process generally known as imputation. In this work, we survey commonly used missing value imputation techniques and compare their performance on a set of building monitoring data. To compare the different types of sensor measurements with widely varying characteristics, we use normalized root mean squared error (NRMSE) as the key metric for the effectiveness of the imputation methods. We additionally consider periodicity and run time when considering comparing methods. Through extensive testing, we find that for small gap sizes, up to 8 consecutive missing values, linear interpolation performs the best; for larger gaps stretching up to 48 consecutive missing values, K-nearest neighbors provides the most accurate imputations; for even larger gaps, more computational intensive methods, such as matrix factorization, achieve the smallest NRMSE. Additionally, we observe that these computationally intensive algorithms not only provide accurate imputations for large gaps, but are also more robust across all types of sensors.

Cho, B↗

Multi Time-scale Imputation aided State Estimation in Distribution System

With the transition to a smart grid, we are witnessing a significant growth in sensor deployments and smart metering infrastructure in the distribution system. However, information from these sensors and meters are typically unevenly sampled at different time-scales and are incomplete. It is critical to effectively aggregate these information sources for situational awareness. In order to reconcile the heterogeneous multi-scale time-series data, we present a multi-task Gaussian process framework. This framework exploits the spatio-temporal correlation across the time-series data to impute data at any desired timescale while providing confidence bounds on the imputations. The value of the imputed data for distribution system operation is illustrated via a matrix completion based state estimation strategy. Results on the IEEE 37 bus distribution system reveals the superior performance of the proposed approach relative to linear interpolation approaches.

Dahale, Shweta↗

GPU-Accelerated Analytic Simulation of Sparse Ionization Signal Formation in Pixelated Projection Detector

This paper presents a GPU-accelerated simulation package, TRED, for next-generation neutrino detectors with pixelated charge readout, leveraging community-driven software ecosystems to ensure adaptability and extensibility. We introduce two generic contributions: (i) an effective-charge representation based on Gaussian quadrature rules, in which the linear- interpolation factors for the field response inside each voxel are absorbed into the effective charge, and (ii) a sparse, block- binned tensor representation that enables efficient FFT-based computation of induced signals on readout electrodes for sparsely activated detector volumes. The former captures structure inside a voxel without dense sampling, while the latter achieves low memory usage and scalable runtime, as demonstrated in bench- mark studies. The underlying data representation is applicable to large-scale detectors and to other computational problems involving sparse activity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A cascaded model of spectral distortions due to spectral response effects and pulse pileup effects in a photon-counting x-ray detector for CT

Purpose: Energy discriminating, photon-counting detectors (PCDs) are an emerging technology for computed tomography (CT) with various potential benefits for clinical CT. The photon energies measured by PCDs can be distorted due to the interactions of a photon with the detector and the interaction of multiple coincident photons. These effects result in distorted recorded x-ray spectra which may lead to artifacts in reconstructed CT images and inaccuracies in tissue identification. Model-based compensation techniques have the potential to account for the distortion effects. This approach requires only a small number of parameters and is applicable to a wide range of spectra and count rates, but it needs an accurate model of the spectral distortions occurring in PCDs. The purpose of this study was to develop a model of those spectral distortions and to evaluate the model using a PCD (model DXMCT-1; DxRay, Inc., Northridge, CA) and various x-ray spectra in a wide range of count rates. Methods: The authors hypothesize that the complex phenomena of spectral distortions can be modeled by: (1) separating them into count-rate independent factors that we call the spectral response effects (SRE), and count-rate dependent factors that we call the pulse pileup effects (PPE), (2) developing separate models for SRE and PPE, and (3) cascading the SRE and PPE models into a combined SRE+PPE model that describes PCD distortions at both low and high count rates. The SRE model describes the probability distribution of the recorded spectrum, with a photo peak and a continuum tail, given the incident photon energy. Model parameters were obtained from calibration measurements with three radioisotopes and then interpolated linearly for other energies. The PPE model used was developed in the authors’ previous work [K. Taguchi et al., “Modeling the performance of a photon counting x-ray detector for CT: Energy response and pulse pileup effects,” Med. Phys. 38(2), 1089–1102 (2011)]. The agreement between the x-ray spectra calculated by the cascaded SRE+PPE model and the measured spectra was evaluated for various levels of deadtime loss ratios (DLR) and incident spectral shapes, realized using different attenuators, in terms of the weighted coefficient of variation (COV{sub W}), i.e., the root mean square difference weighted by the statistical errors of the data and divided by the mean. Results: At low count rates, when DLR < 10%, the distorted spectra measured by the DXMCT-1 were in agreement with those calculated by SRE only, with COV{sub W}'s less than 4%. At higher count rates, the measured spectra were also in agreement with the ones calculated by the cascaded SRE+PPE model; with PMMA as attenuator, COV{sub W} was 5.6% at a DLR of 22% and as small as 6.7% for a DLR as high as 55%. Conclusions: The x-ray spectra calculated by the proposed model agreed with the measured spectra over a wide range of count rates and spectral shapes. The SRE model predicted the distorted, recorded spectra with low count rates over various types and thicknesses of attenuators. The study also validated the hypothesis that the complex spectral distortions in a PCD can be adequately modeled by cascading the count-rate independent SRE and the count-rate dependent PPE.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-D-19A-02: Electron and Photon Absorbed Fractions for Tumors of Varying Sizes and Compositions

Purpose: To calculate absorbed fractions for mono-energetic photons and electrons in tumors of varying compositions using Monte Carlo simulations in MCNPX. Although tumor dosimetry has been previously investigated, these studies are very limited as they only consider absorbed fractions for soft-tissue tumors. Methods: The tumors were modeled as spheres with radii ranging from 0.10 cm to 6.0 cm and with compositions varying from 100% soft tissue to 100% bone. The energies of both the photons and electrons were varied from 10 keV to 10 MeV and were homogenously distributed throughout the tumor volume. Furthermore, this investigation addresses the issue of spherical versus elliptical tumors. Both prolate and oblate spheroid tumors of different compositions were modeled, and absorbed fractions were calculated for various electron and photon energies. Results: The data clearly shows an absorbed fraction dependence on tumor composition. For example, a soft-tissue model for a 3 MeV electron emitted in a 1 cm diameter bone tumor would have an 83% error, and this same assumption for a 500 keV photon would yield a 74% error. Ultimately, empirical fits were created for each of the five material compositions in order to facilitate the absorbed fraction calculation, requiring only the tumor size and particle energy. Furthermore, the data shows that absorbed fractions for moderate spheroids can be well approximated by spherical tumors of equal mass to within 8%, but in the extreme cases where the spheroid resembles more of a disk, the errors can be as high as 30%. Conclusion: This comprehensive data set is most valuable for nuclear medicine dosimetry because it incorporates particle type, particle energy, tumor size, and tumor composition. Although mono-energetic particles were modeled, absorbed fractions and S-values may be calculated for any radionuclide via linear interpolation, as long as the particle energies or spectra are known.

60 APPLIED LIFE SCIENCES↗

SU-E-J-23: An Accurate Algorithm to Match Imperfectly Matched Images for Lung Tumor Detection Without Markers

Purpose: In order to locate lung tumors on projection images without internal markers, digitally reconstructed radiograph (DRR) is created and compared with projection images. Since lung tumors always move and their locations change on projection images while they are static on DRRs, a special DRR (background DRR) is generated based on modified anatomy from which lung tumors are removed. In addition, global discrepancies exist between DRRs and projections due to their different image originations, scattering, and noises. This adversely affects comparison accuracy. A simple but efficient comparison algorithm is reported. Methods: This method divides global images into a matrix of small tiles and similarities will be evaluated by calculating normalized cross correlation (NCC) between corresponding tiles on projections and DRRs. The tile configuration (tile locations) will be automatically optimized to keep the tumor within a single tile which has bad matching with the corresponding DRR tile. A pixel based linear transformation will be determined by linear interpolations of tile transformation results obtained during tile matching. The DRR will be transformed to the projection image level and subtracted from it. The resulting subtracted image now contains only the tumor. A DRR of the tumor is registered to the subtracted image to locate the tumor. Results: This method has been successfully applied to kV fluoro images (about 1000 images) acquired on a Vero (Brainlab) for dynamic tumor tracking on phantom studies. Radiation opaque markers are implanted and used as ground truth for tumor positions. Although, other organs and bony structures introduce strong signals superimposed on tumors at some angles, this method accurately locates tumors on every projection over 12 gantry angles. The maximum error is less than 2.6 mm while the total average error is 1.0 mm. Conclusion: This algorithm is capable of detecting tumor without markers despite strong background signals.

60 APPLIED LIFE SCIENCES↗

SU-E-T-641: Development and Verification of Automatic Reading Dose of Interest From Eclipse's DVH

Purpose: According to clinical and research requirement, we develop a function of automatic reading dose of interest from dose volume histogram(DVH), to replace the traditional method with a mouse one by one point, and it's also verified. Methods: The DVH automatic reading function will be developed in an in-house developed radiotherapy information management system(RTIMS), which is based on Apache+PHP+MySQL. A DVH ASCII file is exported from Varian Eclipse V8.6, which includes the following contents: 1. basic information of patient; 2. dose information of plan; 3. dose information of structures, including basic information and dose volume data of target volume and organ at risk. And the default exported dose volume data also includes relative doses by 1% step and corresponding absolute doses and cumulative relative volumes, and the volumes are 4 decimal fraction. Clinically, we often need read the doses of some integer percent volumes, such as D50 and D30. So it couldn't be directly obtained from the above data, but we can use linear interpolation bye the near volumes and doses: Dx=D2−(V2−Vx)*(D2−D1)/(V2−V1), and program a function to search, read and calculate the corresponding data. And the doses of all preseted volume of interest of all structures can be automatically read one by one patient, and saved as a CSV file. To verify it, we select 24 IMRT plans for prostate cancer, and doses of interest are PTV D98/D95/D5/D2, bladder D30/D50, and rectum D25/D50. Two groups of data, using the automatic reading method(ARM) and pointed dose method(PDM), are analyzed with SPSS 16. The absolute difference=D-ARM-D-PDM, relative difference=absolute difference*100%/prescription dose(7600cGy). Results: The differences are as following: PTV D98/D95/D5/D2: −0.04%/− 0.04%/0.13%/0.19%, bladder D30/D50: −0.02%/0.01%, and rectum D25/D50: 0.03%/0.01%. Conclusion: Using this function, the error is very small, and can be neglected. It could greatly improve the efficiency of clinical work. Project supported by the National Natural Science Foundation of China (Grant No.81101694)

60 APPLIED LIFE SCIENCES↗

TU-F-17A-07: Real-Time Personalized Margins

Purpose: To maximize normal tissue sparing for treatments requiring motion encompassing margins. Motion mitigation techniques including DMLC or couch tracking can freeze tumor motion within the treatment aperture potentially allowing for smaller treatment margins and thus better sparing of normal tissue. To enable for a safe application of this concept in the clinic we propose adapting margins dynamically in real-time during radiotherapy delivery based on personalized tumor localization confidence. To demonstrate technical feasibility we present a phantom study. Methods: We utilize a realistic anthropomorphic dynamic thorax phantom with a lung tumor model embedded close to the spine. The tumor, a 3D-printout of a patient's GTV, is moved 15mm peak-to-peak by diaphragm compression and monitored by continuous EPID imaging in real-time. Two treatment apertures are created for each beam, one representing ITV -based and the other GTV-based margin expansion. A soft tissue localization (STiL) algorithm utilizing the continuous EPID images is employed to freeze tumor motion within the treatment aperture by means of DMLC tracking. Depending on a tracking confidence measure (TCM), the treatment aperture is adjusted between the ITV and the GTV leaf. Results: We successfully demonstrate real-time personalized margin adjustment in a phantom study. We measured a system latency of about 250 ms which we compensated by utilizing a respiratory motion prediction algorithm (ridge regression). With prediction in place we observe tracking accuracies better than 1mm. For TCM=0 (as during startup) an ITV-based treatment aperture is chosen, for TCM=1 a GTV-based aperture and for 0<TCM<1 a linear interpolation between the two apertures is chosen. Conclusions: We have proposed and demonstrated a concept to safely shrink treatment margins during radiotherapy and adapt to tracking confidence in real-time. Normal tissue sparing is maximized. The worst case scenario results in delivering a plan with standard margins used in the clinic today.

60 APPLIED LIFE SCIENCES↗

Correcting lateral response artifacts from flatbed scanners for radiochromic film dosimetry

Purpose: A known factor affecting the accuracy of radiochromic film dosimetry is the lateral response artifact (LRA) induced by nonuniform response of a flatbed scanner in the direction perpendicular to the scan direction. This work reports a practical solution to eliminate such artifacts for all forms of dose QA. Methods: EBT3 films from a single production lot (02181401) cut into rectangular 4 × 5 cm{sup 2} pieces, with the long dimension parallel to the long dimension of the original 20.3 × 25.4 cm{sup 2} sheets, were exposed at a depth of 5 cm on a Varian Trilogy at the center of a 20 × 20 cm{sup 2} open field at seven doses between 50 and 1600 cGy using 6 MV photons. These films together with an unexposed film from the same production lot were lined one next to the other on an Epson 10000XL or 11000XL scanner in portrait orientation with their long dimension parallel to the scan direction. Scanned images were then obtained with the line of films positioned at seven discrete lateral locations perpendicular to the scan direction. The process was repeated in landscape orientation and on three other Epson scanners. Data were also collected for three additional production lots of EBT3 film (11051302, 03031401, and 03171403). From measurements at the various lateral positions, the scanner response was determined as a function of the lateral position of the scanned film. For a given color channel X, the response at any lateral position L is related to the response at the center, C, of the scanner by Response(C, D, X) = A{sub L,X} + B{sub L,X}⋅Response(L, D, X), where D is dose and the coefficients A{sub L,X} and B{sub L,X} are determined from the film measurements at the center of the scanner and six other discrete lateral positions. The values at intermediate lateral positions were obtained by linear interpolation. The coefficients were determined for the red, green, and blue color channels, preserving the ability to apply triple-channel dosimetry once corrections were applied to compensate for the lateral position response artifact. To validate this method, corrections were applied to several films that were exposed to 15 × 15cm{sup 2} open fields and large IMRT and VMAT fields and scanned at the extreme edges of the scan window in addition to the central location. Calibration and response data were used to generate dose maps and perform gamma analysis using single- or triple-channel dosimetry with FilmQAPro 2014 software. Results: The authors’ study found that calibration curves at the different lateral positions could be correlated by a simple two-point rescaling using the response for unexposed film as well as the response of film exposed at high doses between 800 and 1600 cGy. The coefficients A{sub L,X} and B{sub L,X} for each color channel X were found to be independent of dose at each lateral location L. This made it possible to apply the relationship Response(C, D, X) = A{sub L,X} + B{sub L,X}⋅Response(L, D, X), to the raw film responses, permitting correction of the response values at any lateral position to an equivalent response, as if that part of the film was located at the center of the scanner. This correction method was validated for several films exposed to open as well as large IMRT and VMAT fields. Conclusions: The work reported elaborates on the process using the correction procedures to eliminate the lateral response artifact and demonstrates improvements in the accuracy of radiochromic film dosimetry for the radiation therapy quality assurance applications.

07 ISOTOPE AND RADIATION SOURCES↗

SU-E-J-257: A PCA Model to Predict Adaptive Changes for Head&neck Patients Based On Extraction of Geometric Features From Daily CBCT Datasets

Purpose: Using daily cone beam CTs (CBCTs) to develop principal component analysis (PCA) models of anatomical changes in head and neck (H&N) patients and to assess the possibility of using these prospectively in adaptive radiation therapy (ART). Methods: Planning CT (pCT) images of 4 H&N patients were deformed to model several different systematic changes in patient anatomy during the course of the radiation therapy (RT). A Pinnacle plugin was used to linearly interpolate the systematic change in patient for the 35 fraction RT course and to generate a set of 35 synthetic CBCTs. Each synthetic CBCT represents the systematic change in patient anatomy for each fraction. Deformation vector fields (DVFs) were acquired between the pCT and synthetic CBCTs with random fraction-to-fraction changes were superimposed on the DVFs. A patient-specific PCA model was built using these DVFs containing systematic plus random changes. It was hypothesized that resulting eigenDVFs (EDVFs) with largest eigenvalues represent the major anatomical deformations during the course of treatment. Results: For all 4 patients, the PCA model provided different results depending on the type and size of systematic change in patient’s body. PCA was more successful in capturing the systematic changes early in the treatment course when these were of a larger scale with respect to the random fraction-to-fraction changes in patient’s anatomy. For smaller scale systematic changes, random changes in patient could completely “hide” the systematic change. Conclusion: The leading EDVF from the patientspecific PCA models could tentatively be identified as a major systematic change during treatment if the systematic change is large enough with respect to random fraction-to-fraction changes. Otherwise, leading EDVF could not represent systematic changes reliably. This work is expected to facilitate development of population-based PCA models that can be used to prospectively identify significant anatomical changes early in treatment. This work is supported in part by a grant from Varian Medical Systems, Palo Alto, CA.

60 APPLIED LIFE SCIENCES↗

SU-E-T-150: Brachytherapy QA Employing a High Resolution Liquid Filled Ionisation Chamber Array: Initial Experience and Limitations

Purpose: Verifying a complex 3D brachytherapy dose distribution by measurement is non-trivial. Ideally a photon detector array should be independent of energy and angle, have high spatial resolution and be robust for routine clinical use. Methods: An iridium-192 source was used. A PMMA jig was constructed comprising an outer slab and a central insert with eight milled channels for 1.33mm (outer diameter) steel needles, see figure. All calculations were performed using an empty CT study reconstructing eight virtual needles (QA-CT), using the v2 source model (Elekta AG, Sweden). A high resolution liquid filled ionisation chamber array SRS1000, together with Verisoft software v6.0 (PTW Freiburg, Germany), was used to perform measurements of plans of increasing complexity to evaluate its suitability for device- and patient-specific QA. The dimension of backscatter material was investigated. The patient plan dwell time distribution was entered manually into the QA-CT and the dose distribution was calculated. Results: Our measurements indicate that the array is independent of energy and angle. The resulting measured dose values are linearly interpolated to 2025 values. Shifts of 1mm of the entire needle are readily detectable. Individual dwell position shifts (2.5mm) are also readily measurable. Moreover a dwell time increase of 1 second both in the edge and central region are detectable. Conclusion: The high resolution SRS1000 array is a powerful instrument for brachytherapy QA enabling 977 simultaneous measurements to be performed. Our measurements suggest 60mm of RW3 backscatter material upstream and downstream are sufficient. Local percentage difference analysis is useful for device based QA, normalized relative percentage difference is arguably better for patient specific QA. Automated transfer of patient plan dwell time distribution to the QA plan is required to enable a comprehensive patient QA study to be performed. Moreover the described measurement technique can be used for commissioning new brachytherapy applicators.

60 APPLIED LIFE SCIENCES↗

SU-E-T-485: In Vivo Dosimetry with EBT3 Radiochromic Films for TBI Treatments

Purpose: Total body irradiation (TBI) is a technique that requires special equipment to control “in vivo” the dose to the patient because it is a complex technique performed in extraordinary conditions. There are several devices to perform this task (diodes, TLDs, ionization chambers, MOSFET). In this paper we study the possibility of performing these measurements with radiochromic films EBT3 properly calibrated. This method has been compared to the PTW diodes system for TBI. Methods: Once made the TC to the patients, we measured different thicknesses of the relevant areas of the body (head, neck, chest with or without arms, umbilicus area, knees and ankles); for each of these thicknesses we measured dose rate (cGy / UM) in RW3 phantom, in TBI conditions, with ionization chamber in the center; in turn, the input diode and the output of each configuration is placed to assign dose to each set of diodes. Movie calibration is performed according to manufacturer’s recommendations but TBI conditions. The dose at the center of each thickness compared to a linear interpolation of the dose at the entrance and exit, resulting in an adequate approximation. Finally in each session for each patient put a piece of film (2×2 cm2) at the entrance and another at the exit in each area, obtaining these readings and interpolating the estimated center dose, as with the diodes. Results: These results show a greater homogeneity in the distribution for use with film and validate the use of the same for this task and, if necessary, to avoid purchasing diode group if they have not. Conclusion: By using radiochromic films for this technique gives us a proper calculation of the dose received by the patient in the absence of other methods, or gives us a second additional track that already used normally.

60 APPLIED LIFE SCIENCES↗

The development of a population of 4D pediatric XCAT phantoms for imaging research and optimization

Purpose: We previously developed a set of highly detailed 4D reference pediatric extended cardiac-torso (XCAT) phantoms at ages of newborn, 1, 5, 10, and 15 yr with organ and tissue masses matched to ICRP Publication 89 values. In this work, we extended this reference set to a series of 64 pediatric phantoms of varying age and height and body mass percentiles representative of the public at large. The models will provide a library of pediatric phantoms for optimizing pediatric imaging protocols. Methods: High resolution positron emission tomography-computed tomography data obtained from the Duke University database were reviewed by a practicing experienced radiologist for anatomic regularity. The CT portion of the data was then segmented with manual and semiautomatic methods to form a target model defined using nonuniform rational B-spline surfaces. A multichannel large deformation diffeomorphic metric mapping algorithm was used to calculate the transform from the best age matching pediatric XCAT reference phantom to the patient target. The transform was used to complete the target, filling in the nonsegmented structures and defining models for the cardiac and respiratory motions. The complete phantoms, consisting of thousands of structures, were then manually inspected for anatomical accuracy. The mass for each major tissue was calculated and compared to linearly interpolated ICRP values for different ages. Results: Sixty four new pediatric phantoms were created in this manner. Each model contains the same level of detail as the original XCAT reference phantoms and also includes parameterized models for the cardiac and respiratory motions. For the phantoms that were 10 yr old and younger, we included both sets of reproductive organs. This gave them the capability to simulate both male and female anatomy. With this, the population can be expanded to 92. Wide anatomical variation was clearly seen amongst the phantom models, both in organ shape and size, even for models of the same age and sex. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. Conclusions: This work provides a large cohort of highly detailed pediatric phantoms with 4D capabilities of varying age, height, and body mass. The population of phantoms will provide a vital tool with which to optimize 3D and 4D pediatric imaging devices and techniques in terms of image quality and radiation-absorbed dose.

60 APPLIED LIFE SCIENCES↗

Assessment of metal artifact reduction methods in pelvic CT

Purpose: Metal artifact reduction (MAR) produces images with improved quality potentially leading to confident and reliable clinical diagnosis and therapy planning. In this work, the authors evaluate the performance of five MAR techniques for the assessment of computed tomography images of patients with hip prostheses. Methods: Five MAR algorithms were evaluated using simulation and clinical studies. The algorithms included one-dimensional linear interpolation (LI) of the corrupted projection bins in the sinogram, two-dimensional interpolation (2D), a normalized metal artifact reduction (NMAR) technique, a metal deletion technique, and a maximum a posteriori completion (MAPC) approach. The algorithms were applied to ten simulated datasets as well as 30 clinical studies of patients with metallic hip implants. Qualitative evaluations were performed by two blinded experienced radiologists who ranked overall artifact severity and pelvic organ recognition for each algorithm by assigning scores from zero to five (zero indicating totally obscured organs with no structures identifiable and five indicating recognition with high confidence). Results: Simulation studies revealed that 2D, NMAR, and MAPC techniques performed almost equally well in all regions. LI falls behind the other approaches in terms of reducing dark streaking artifacts as well as preserving unaffected regions (p < 0.05). Visual assessment of clinical datasets revealed the superiority of NMAR and MAPC in the evaluated pelvic organs and in terms of overall image quality. Conclusions: Overall, all methods, except LI, performed equally well in artifact-free regions. Considering both clinical and simulation studies, 2D, NMAR, and MAPC seem to outperform the other techniques.

60 APPLIED LIFE SCIENCES↗

An online replanning method using warm start optimization and aperture morphing for flattening-filter-free beams

Purpose: In a situation where a couch shift for patient positioning is not preferred or prohibited (e.g., MR-linac), segment aperture morphing (SAM) can address target dislocation and deformation. For IMRT/VMAT with flattening-filter-free (FFF) beams, however, SAM method would lead to an adverse translational dose effect due to the beam unflattening. Here the authors propose a new two-step process to address both the translational effect of FFF beams and the target deformation. Methods: The replanning method consists of an offline and an online step. The offline step is to create a series of preshifted-plans (PSPs) obtained by a so-called “warm start” optimization (starting optimization from the original plan, rather than from scratch) at a series of isocenter shifts. The PSPs all have the same number of segments with very similar shapes, since the warm start optimization only adjusts the MLC positions instead of regenerating them. In the online step, a new plan is obtained by picking the closest PSP or linearly interpolating the MLC positions and the monitor units of the closest PSPs for the shift determined from the image of the day. This two-step process is completely automated and almost instantaneous (no optimization or dose calculation needed). The previously developed SAM algorithm is then applied for daily deformation. The authors tested the method on sample prostate and pancreas cases. Results: The two-step interpolation method can account for the adverse dose effects from FFF beams, while SAM corrects for the target deformation. Plan interpolation method is effective in diminishing the unflat beam effect and may allow reducing the required number of PSPs. The whole process takes the same time as the previously reported SAM process (5–10 min). Conclusions: The new two-step method plus SAM can address both the translation effects of FFF beams and target deformation, and can be executed in full automation except the delineation of target contour required by the SAM process.

60 APPLIED LIFE SCIENCES↗

SU-F-I-75: Half-Value Layer Thicknesses and Homogeneity Coefficients for Fluoroscopic X-Ray Beam Spectra Incorporating Spectral Filtration

Purpose: The purpose of this investigation is to quantify various first half-value-layers (HVLs), second HVLs and homogeneity coefficients (HCs) for a state-of-the-art fluoroscope utilizing spectral (copper) filtration. Methods: A Radcal (Monrovia, Ca) AccuPro dosimeter with a 10×6-6 calibrated ionization chamber was used to measure air kerma for radiographic x-ray exposures made on a Siemens (Erlangen, Germany) Artis ZeeGo fluoroscope operated in the service mode. The ionization chamber was centered in the x-ray beam at 72 cm from the focal spot with a source-to-image-distance of 120 cm. The collimators were introduced to limit the x-ray field to approximately 5 cm × 5 cm at the ionization chamber plane. Type-1100 aluminum filters, in 0.5 mm increments, were used to determine the HVL. Two HVL calculation methods were used, log-linear interpolation and Lambert-W interpolation as described by Mathieu [Med Phys, 38(8), 4546 (2011)]. Multiple measurements were made at 60, 80, 100, 120 kVp at spectral filtration thicknesses of 0, 0.1, 0.3, 0.6 and 0.9 mm. Results: First HVL, second HVL, and HCs are presented for the fluoroscopic x-ray beam spectra indicated above, with nearly identical results from the two interpolation methods. Accuracy of the set kVp was also determined and deviated less than 2%. First HVLs for fluoroscopic x-ray beam spectra without spectral filtration determined in our study were 7%–16% greater than previously published data by Fetterly et al. [Med Phys, 28, 205 (2001)]. However, the FDA minimum HVL requirements changed since that publication, requiring larger HVLs as of 2006. Additionally, x-ray tube and generator architecture have substantially changed over the last 15 years providing different beam spectra. Conclusion: X-ray beam quality characteristics for state-of-the-art fluoroscopes with spectral filtration have not been published. This study provides reference data which will be useful for defining beam qualities encountered on fluoroscopes using spectral filtration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗