Search NASA⌕ Search

SEARCH · Search NASA

Results for “Digital image quality phantom”

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.

80 records · Page 5

Improved proton computed tomography by dual modality image reconstruction

Purpose: Proton computed tomography (CT) is a promising image modality for improving the stopping power estimates and dose calculations for particle therapy. However, the finite range of about 33 cm of water of most commercial proton therapy systems limits the sites that can be scanned from a full 360° rotation. In this paper the authors propose a method to overcome the problem using a dual modality reconstruction (DMR) combining the proton data with a cone-beam x-ray prior. Methods: A Catphan 600 phantom was scanned using a cone beam x-ray CT scanner. A digital replica of the phantom was created in the Monte Carlo code Geant4 and a 360° proton CT scan was simulated, storing the entrance and exit position and momentum vector of every proton. Proton CT images were reconstructed using a varying number of angles from the scan. The proton CT images were reconstructed using a constrained nonlinear conjugate gradient algorithm, minimizing total variation and the x-ray CT prior while remaining consistent with the proton projection data. The proton histories were reconstructed along curved cubic-spline paths. Results: The spatial resolution of the cone beam CT prior was retained for the fully sampled case and the 90° interval case, with the MTF = 0.5 (modulation transfer function) ranging from 5.22 to 5.65 linepairs/cm. In the 45° interval case, the MTF = 0.5 dropped to 3.91 linepairs/cm For the fully sampled DMR, the maximal root mean square (RMS) error was 0.006 in units of relative stopping power. For the limited angle cases the maximal RMS error was 0.18, an almost five-fold improvement over the cone beam CT estimate. Conclusions: Dual modality reconstruction yields the high spatial resolution of cone beam x-ray CT while maintaining the improved stopping power estimation of proton CT. In the case of limited angles, the use of prior image proton CT greatly improves the resolution and stopping power estimate, but does not fully achieve the quality of a 360° proton CT scan.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Tailoring automatic exposure control toward constant detectability in digital mammography

Purpose: The automatic exposure control (AEC) modes of most full field digital mammography (FFDM) systems are set up to hold pixel value (PV) constant as breast thickness changes. This paper proposes an alternative AEC mode, set up to maintain some minimum detectability level, with the ultimate goal of improving object detectability at larger breast thicknesses. Methods: The default “OPDOSE” AEC mode of a Siemens MAMMOMAT Inspiration FFDM system was assessed using poly(methyl methacrylate) (PMMA) of thickness 20, 30, 40, 50, 60, and 70 mm to find the tube voltage and anode/filter combination programmed for each thickness; these beam quality settings were used for the modified AEC mode. Detectability index (d′), in terms of a non-prewhitened model observer with eye filter, was then calculated as a function of tube current-time product (mAs) for each thickness. A modified AEC could then be designed in which detectability never fell below some minimum setting for any thickness in the operating range. In this study, the value was chosen such that the system met the achievable threshold gold thickness (T{sub t}) in the European guidelines for the 0.1 mm diameter disc (i.e., T{sub t} ≤ 1.10 μm gold). The default and modified AEC modes were compared in terms of contrast-detail performance (T{sub t}), calculated detectability (d′), signal-difference-to-noise ratio (SDNR), and mean glandular dose (MGD). The influence of a structured background on object detectability for both AEC modes was examined using a CIRS BR3D phantom. Computer-based CDMAM reading was used for the homogeneous case, while the images with the BR3D background were scored by human observers. Results: The default OPDOSE AEC mode maintained PV constant as PMMA thickness increased, leading to a reduction in SDNR for the homogeneous background 39% and d′ 37% in going from 20 to 70 mm; introduction of the structured BR3D plate changed these figures to 22% (SDNR) and 6% (d′), respectively. Threshold gold thickness (0.1 mm diameter disc) for the default AEC mode in the homogeneous background increased by 62% in going from 20 to 70 mm PMMA thickness; in the structured background, the increase was 39%. Implementation of the modified mode entailed an increase in mAs at PMMA thicknesses >40 mm; the modified AEC held threshold gold thickness constant above 40 mm PMMA with a maximum deviation of 5% in the homogeneous background and 3% in structured background. SDNR was also held constant with a maximum deviation of 4% and 2% for the homogeneous and the structured background, respectively. These results were obtained with an increase of MGD between 15% and 73% going from 40 to 70 mm PMMA thickness. Conclusions: This work has proposed and implemented a modified AEC mode, tailored toward constant detectability at larger breast thickness, i.e., above 40 mm PMMA equivalent. The desired improvement in object detectability could be obtained while maintaining MGD within the European guidelines achievable dose limit. (A study designed to verify the performance of the modified mode using more clinically realistic data is currently underway.)

60 APPLIED LIFE SCIENCES↗

Characterization of scatter in digital mammography from physical measurements

Purpose: That scattered radiation negatively impacts the quality of medical radiographic imaging is well known. In mammography, even slight amounts of scatter reduce the high contrast required for subtle soft-tissue imaging. In current clinical mammography, image contrast is partially improved by use of an antiscatter grid. This form of scatter rejection comes with a sizeable dose penalty related to the concomitant elimination of valuable primary radiation. Digital mammography allows the use of image processing as a method of scatter correction that might avoid effects that negatively impact primary radiation, while potentially providing more contrast improvement than is currently possible with a grid. For this approach to be feasible, a detailed characterization of the scatter is needed. Previous research has modeled scatter as a constant background that serves as a DC bias across the imaging surface. The goal of this study was to provide a more substantive data set for characterizing the spatially-variant features of scatter radiation at the image detector of modern mammography units. Methods: This data set was acquired from a model of the radiation beam as a matrix of very narrow rays or pencil beams. As each pencil beam penetrates tissue, the pencil widens in a predictable manner due to the production of scatter. The resultant spreading of the pencil beam at the detector surface can be characterized by two parameters: mean radial extent (MRE) and scatter fraction (SF). The SF and MRE were calculated from measurements obtained using the beam stop method. Two digital mammography units were utilized, and the SF and MRE were found as functions of target, filter, tube potential, phantom thickness, and presence or absence of a grid. These values were then used to generate general equations allowing the SF and MRE to be calculated for any combination of the above parameters. Results: With a grid, the SF ranged from a minimum of about 0.05 to a maximum of about 0.16, and the MRE ranged from about 3 to 13 mm. Without a grid, the SF ranged from a minimum of 0.25 to a maximum of 0.52, and the MRE ranged from about 20 to 45 mm. The SF with a grid demonstrated a mild dependence on target/filter combination and kV, whereas the SF without a grid was independent of these factors. The MRE demonstrated a complex relationship as a function of kV, with notable difference among target/filter combinations. The primary source of change in both the SF and MRE was phantom thickness. Conclusions: Because breast tissue varies spatially in physical density and elemental content, the effective thickness of breast tissue varies spatially across the imaging field, resulting in a spatially-variant scatter distribution in the imaging field. The data generated in this study can be used to characterize the scatter contribution on a point-by-point basis, for a variety of different techniques.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Gantry and isocenter displacements of a linear accelerator caused by an add-on micromultileaf collimator

Purpose: The delivery of high quality stereotactic radiosurgery (SRS) and stereotactic radiotherapy (SRT) treatments to the patient requires knowledge of the position of the isocenter to submillimeter accuracy. To meet the requirements the deviation between the radiation and mechanical isocenters must be less than 1 mm. The use of add-on micromultileaf collimators ({mu}MLCs) in SRS and SRT is an additional challenge to the anticipated high-level geometric and dosimetric accuracy of the treatment. The aim of this work was to quantify the gantry excursions during rotation with and without an add-on {mu}MLC attached to the gantry head. In addition, the shift in the position of the isocenter and its correlation to the kV beam center of the cone-beam CT system was included in the study. Methods: The quantification of the gantry rotational performance was done using a pointer supported by an in-house made rigid holder attached to the gantry head of the accelerator. The pointer positions were measured using a digital theodolite. To quantify the effect of an {mu}MLC of 50 kg, the measurements were repeated with the {mu}MLC attached to the gantry head. The displacement of the isocenter due to an add-on {mu}MLC of 50 kg was also investigated. In case of the pointer measurement the {mu}MLC was simulated by weights attached to the gantry head. A method of least squares was applied to determine the position and displacement of the mechanical isocenter. Additionally, the displacement of the radiation isocenter was measured using a ball-bearing phantom and the electronic portal image device system. These measurements were based on 8 MV photon beams irradiated onto the ball from the four cardinal angles and two opposed collimator angles. The measurements and analysis of the data were carried out automatically using software delivered by the manufacturer. Results: The displacement of the mechanical isocenter caused by a 50 kg heavy {mu}MLC was found to be (-0.01 {+-} 0.05, -0.10 {+-} 0.03, -0.26 {+-} 0.05) mm in lateral, longitudinal, and vertical direction, respectively. Similarly, the displacement of the radiation isocenter was found to be (0.00 {+-} 0.03, -0.08 {+-} 0.06, -0.32 {+-} 0.02) mm. Good agreement was found between the displacement of the two isocenters. A displacement of the kV cone-beam CT beam center due to the attached weight of 50 kg could not be detected. Conclusions: General characteristics of the gantry arm excursions and displacements caused by an add-on {mu}MLC have been reported. A 50 kg heavy add-on {mu}MLC results in a isocenter displacement downward of 0.26-0.32 mm. The authors recommend that the beam center of the kV cone-beam CT image system should be matched to the isocenter related to the weight of the {mu}MLC. Consequently, the imperfections in isocenter localizations are transferred to the conventional radiotherapy where the clinical consequences of uncertainties in the submillimeter regime are negligible.

60 APPLIED LIFE SCIENCES↗

Eddy current correction in volume-localized MR spectroscopy

The quality of volume-localized magnetic resonance spectroscopy is affected by eddy currents caused by gradient switching. Eddy currents can be reduced with improved gradient systems; however, it has been suggested that the distortion due to eddy currents can be compensated for during postprocessing with a single-frequency reference signal. The authors propose modifying current techniques for acquiring the single-frequency reference signal by using relaxation weighting to reduce interference from components that cannot be eliminated by digital filtering alone. Additional sequences with T1 or T2 weighting for reference signal acquisition are shown to have the same eddy current characteristics as the original signal without relaxation weighting. The authors also studied a new eddy current correction method that does not require a single-frequency reference signal. This method uses two free induction decays (FIDs) collected from the same volume with two sequences with opposite gradients. Phase errors caused by eddy currents are opposite in these two FIDs and can be canceled completely by combining the FIDs. These methods were tested in a phantom. Eddy current distortions were corrected, allowing quantitative measurement of structures such as the -CH = CH- component, which is otherwise undetectable.

NASA Discipline Musculoskeletal↗

Optimization of energy window for {sup 90}Y bremsstrahlung SPECT imaging for detection tasks using the ideal observer with model-mismatch

Purpose: In yttrium-90 ({sup 90}Y) microsphere brachytherapy (radioembolization) of unresectable liver cancer, posttherapy {sup 90}Y bremsstrahlung single photon emission computed tomography (SPECT) has been used to document the distribution of microspheres in the patient and to help predict potential side effects. The energy window used during projection acquisition can have a significant effect on image quality. Thus, using an optimal energy window is desirable. However, there has been great variability in the choice of energy window due to the continuous and broad energy distribution of {sup 90}Y bremsstrahlung photons. The area under the receiver operating characteristic curve (AUC) for the ideal observer (IO) is a widely used figure of merit (FOM) for optimizing the imaging system for detection tasks. The IO implicitly assumes a perfect model of the image formation process. However, for {sup 90}Y bremsstrahlung SPECT there can be substantial model-mismatch (i.e., difference between the actual image formation process and the model of it assumed in reconstruction), and the amount of the model-mismatch depends on the energy window. It is thus important to account for the degradation of the observer performance due to model-mismatch in the optimization of the energy window. The purpose of this paper is to optimize the energy window for {sup 90}Y bremsstrahlung SPECT for a detection task while taking into account the effects of the model-mismatch. Methods: An observer, termed the ideal observer with model-mismatch (IO-MM), has been proposed previously to account for the effects of the model-mismatch on IO performance. In this work, the AUC for the IO-MM was used as the FOM for the optimization. To provide a clinically realistic object model and imaging simulation, the authors used a background-known-statistically and signal-known-statistically task. The background was modeled as multiple compartments in the liver with activity parameters independently following a Gaussian distribution; the signal was modeled as a tumor with a Gaussian-distributed activity parameter located randomly with equal probability at one of three positions. The IO test statistics (i.e., likelihood ratios) were estimated using Markov-chain Monte Carlo methods. The authors realistically modeled human anatomy using a digital phantom code, and realistically simulated {sup 90}Y bremsstrahlung SPECT imaging with a clinical SPECT system and typical imaging parameters using a previously validated Monte Carlo bremsstrahlung simulation method. Model-mismatch was included by modeling image formation process in the calculation of IO test statistics using an analytic modeling method previously developed for quantitative {sup 90}Y bremsstrahlung SPECT. To demonstrate the effects of the model-mismatch on the detection task, the authors optimized the energy window both with and without model-mismatch included in the IO. Results: For all the energy windows, the AUC values for the IO-MM were smaller than that for the IO. The optimal windows for the IO-MM and the IO were 80-180 and 60-400 keV, respectively. Conclusions: The authors have demonstrated the degradation of the ideal performance due to model-mismatch and optimized the energy window for {sup 90}Y bremsstrahlung SPECT for detection tasks by accounting for the effects of the model-mismatch. The obtained optimal window was much narrower when taking into account the model-mismatch and similar to that obtained previously for estimation tasks.

61 RADIATION PROTECTION AND DOSIMETRY↗

Collimator design for a multipinhole brain SPECT insert for MRI

Purpose: Brain single photon emission computed tomography (SPECT) imaging is an important clinical tool, with unique tracers for studying neurological diseases. Nowadays, most commercial SPECT systems are combined with x-ray computed tomography (CT) in so-called SPECT/CT systems to obtain an anatomical background for the functional information. However, while CT images have a high spatial resolution, they have a low soft-tissue contrast, which is an important disadvantage for brain imaging. Magnetic resonance imaging (MRI), on the other hand, has a very high soft-tissue contrast and does not involve extra ionizing radiation. Therefore, the authors designed a brain SPECT insert that can operate inside a clinical MRI. Methods: The authors designed and simulated a compact stationary multipinhole SPECT insert based on digital silicon photomultiplier detector modules, which have shown to be MR-compatible and have an excellent intrinsic resolution (0.5 mm) when combined with a monolithic 2 mm thick LYSO crystal. First, the authors optimized the different parameters of the SPECT system to maximize sensitivity for a given target resolution of 7.2 mm in the center of the field-of-view, given the spatial constraints of the MR system. Second, the authors performed noiseless simulations of two multipinhole configurations to evaluate sampling and reconstructed resolution. Finally, the authors performed Monte Carlo simulations and compared the SPECT insert with a clinical system with ultrahigh-resolution (UHR) fan beam collimators, based on contrast-to-noise ratio and a visual comparison of a Hoffman phantom with a 9 mm cold lesion. Results: The optimization resulted in a stationary multipinhole system with a collimator radius of 150.2 mm and a detector radius of 172.67 mm, which corresponds to four rings of 34 diSPM detector modules. This allows the authors to include eight rings of 24 pinholes, which results in a system volume sensitivity of 395 cps/MBq. Noiseless simulations show sufficient axial sampling (in a Defrise phantom) and a reconstructed resolution of 5.0 mm (in a cold-rod phantom). The authors compared the 24-pinhole setup with a 34-pinhole system (with the same detector radius but a collimator radius of 156.63 mm) and found that 34 pinholes result in better uniformity but a worse reconstruction of the cold-rod phantom. The authors also compared the 24-pinhole system with a clinical triple-head UHR fan beam system based on contrast-to-noise ratio and found that the 24-pinhole setup performs better for the 6 mm hot and the 16 mm cold lesions and worse for the 8 and 10 mm hot lesions. Finally, the authors reconstructed noisy projection data of a Hoffman phantom with a 9 mm cold lesion and found that the lesion was slightly better visible on the multipinhole image compared to the fan beam image. Conclusions: The authors have optimized a stationary multipinhole SPECT insert for MRI and showed the feasibility of doing brain SPECT imaging inside a MRI with an image quality similar to the best clinical SPECT systems available.

60 APPLIED LIFE SCIENCES↗

TH-CD-BRA-05: First Water Calorimetric Dw Measurement and Direct Measurement of Magnetic Field Correction Factors, KQ,B, in a 1.5 T B-Field of An MRI Linac

Purpose: Reference dosimetry in MR-guided radiotherapy is performed in the presence of a B-field. As a consequence the response of ionization chambers changes considerably and depends on parameters not considered in traditional reference dosimetry. Therefore future Codes of Practices need ionization chamber correction factors to correct for both the change in beam quality and the presence of a B-field. The objective was to study the feasibility of water calorimetric absorbed-dose measurements in a 1.5 T B-field of an MRLinac and the direct measurement of kQ,B calibration of ionization chambers. Methods: Calorimetric absorbed dose to water Dw was measured with a new water calorimeter in the bore of an MRLinac (TPR20,10 of 0.702). Two waterproof ionization chambers (PTW 30013, IBA FC-65G) were calibrated inside the calorimeter phantom (ND,w,Q,B). Both measurements were normalized to a monitor ionization chamber. Ionization chamber measurements were corrected for conventional influence parameter. Based on the chambers’ Co-60 calibrations (ND,w,Q0), measured directly against the calorimeter. In this study the correction factors kQ,B was determined as the ratio of the calibration coefficients in the MRLinac and in Co-60. Additionally, kB was determined based on kQ values obtained with the IAEA TRS-398 Code of Practice. Results: The kQ,B factors of the ionization chambers mentioned above were respectively 0.9488(8) and 0.9445(8) with resulting kB factors of 0.961(13) and 0.952(13) with standard uncertainties on the least significant digit(s) between brackets. Conclusion: Calorimetric Dw measurements and calibration of waterproof ionization chambers were successfully carried out in the 1.5 T B-field of an MRLinac with a standard uncertainty of 0.7%. Preliminary kQ,B and kB factors were determined with standard uncertainties of respectively 0.8% and 1.3%. The kQ,B agrees with an alternative method within 0.4%. The feasibility of water calorimetry in the presence of B-fields was demonstrated by the direct determination of Dw and kQ,B. This work was supported by EMRP grant HLT06. The EMRP is jointly funded by the EMRP participating countries within EURAMET and the European Union.

60 APPLIED LIFE SCIENCES↗