Comparison of different calculation techniques for absorbed dose assessment in patient specific peptide receptor radionuclide therapy.

AIM: The present work concerns the comparison of the performances of three systems for dosimetry in RPT that use different techniques for absorbed dose calculation (organ-level dosimetry, voxel-level dose kernel convolution and Monte Carlo simulations). The aim was to assess the importance of the choice of the most adequate calculation modality, providing recommendations about the choice of the computation tool. METHODS: The performances were evaluated both on phantoms and patients in a multi-level approach. Different phantoms filled with a 177Lu-radioactive solution were used: a homogeneous cylindrical phantom, a phantom with organ-shaped inserts and two cylindrical phantoms with inserts different for shape and volume. A total of 70 patients with NETs treated by PRRT with 177Lu-DOTATOC were retrospectively analysed. RESULTS: The comparisons were performed mainly between the mean values of the absorbed dose in the regions of interest. A general better agreement was obtained between Dose kernel convolution and Monte Carlo simulations results rather than between either of these two and organ-level dosimetry, both for phantoms and patients. Phantoms measurements also showed the discrepancies mainly depend on the geometry of the inserts (e.g. shape and volume). For patients, differences were more pronounced than phantoms and higher inter/intra patient variability was observed. CONCLUSION: This study suggests that voxel-level techniques for dosimetry calculation are potentially more accurate and personalized than organ-level methods. In particular, a voxel-convolution method provides good results in a short time of calculation, while Monte Carlo based computation should be conducted with very fast calculation systems for a possible use in clinics, despite its intrinsic higher accuracy. Attention to the calculation modality is recommended in case of clinical regions of interest with irregular shape and far from spherical geometry, in which Monte Carlo seems to be more accurate than voxel-convolution methods.

Usefulness of an additional lead shielding device in reducing occupational radiation exposure during interventional endoscopic procedures: An observational study.

Adoption of interventional endoscopic procedures is increasing with increasing prevalence of diseases. However, medical radiation exposure is concerning; therefore, radiation protection for medical staff is important. However, there is limited information on the usefulness of an additional lead shielding device during interventional endoscopic procedures. Therefore, we aimed to determine whether an additional lead shielding device protects medical staff from radiation.An X-ray unit (CUREVISTA; Hitachi Medical Systems, Tokyo, Japan) with an over-couch X-ray system was used. Fluoroscopy-associated scattered radiation was measured using a water phantom placed at the locations of the endoscopist, assistant, nurse, and clinical engineer. For each location, measurements were performed at the gonad and thyroid gland/eye levels. Comparisons were performed between with and without the additional lead shielding device and with and without a gap in the shielding device. Additionally, a clinical study was performed with 27 endoscopic retrograde cholangiopancreatography procedures.The scattered radiation dose was lower with than without additional lead shielding at all medical staff locations and decreased by 84.7%, 82.8%, 78.2%, and 83.7%, respectively, at the gonad level and by 89.2%, 86.4%, 91.2%, and 87.0%, respectively, at the thyroid gland/eye level. Additionally, the scattered radiation dose was lower without than with a gap in the shielding device at all locations.An additional lead shielding device could protect medical staff from radiation during interventional endoscopic procedures. However, gaps in protective equipment reduce effectiveness and should be eliminated.

Towards guidelines to harmonize textural features in PET: Haralick textural features vary with image noise, but exposure-invariant domains enable comparable PET radiomics.

PURPOSE: Image texture is increasingly used to discriminate tissues and lesions in PET/CT. For quantification or in computer-aided diagnosis, textural feature analysis must produce robust and comparable values. Because statistical feature values depend on image count statistics, we investigated in depth the stability of Haralick features values as functions of acquisition duration, and for common image resolutions and reconstructions. METHODS: A homogeneous cylindrical phantom containing 9.6 kBq/ml Ge-68 was repeatedly imaged on a Siemens Biograph mCT, with acquisition durations ranging from three seconds to three hours. Images with 1.5, 2, and 4 mm isometrically spaced voxels were reconstructed with filtered back-projection (FBP), ordered subset expectation maximization (OSEM), and the Siemens TrueX algorithm. We analysed Haralick features derived from differently quantized (3 to 8-bit) grey level co-occurrence matrices (GLCMs) as functions of exposure E, which we defined as the product of activity concentration in a volume of interest (VOI) and acquisition duration. The VOI was a 50 mm wide cube at the centre of the phantom. Feature stability was defined for df/dE → 0. RESULTS: The most stable feature values occurred in low resolution FBPs, whereas some feature values from 1.5 mm TrueX reconstructions ranged over two orders of magnitude. Within the same reconstructions, most feature value-exposure curves reached stable plateaus at similar exposures, regardless of GLCM quantization. With 8-bit GLCM, median time to stability was 16 s and 22 s for FBPs, 18 s and 125 s for OSEM, and 23 s, 45 s, and 76 s for PSF reconstructions, with longer durations for higher resolutions. Stable exposures coincided in OSEM and TrueX reconstructions with image noise distributions converging to a Gaussian. In FBP, the occurrence of stable values coincided the disappearance of negatives image values in the VOI. CONCLUSIONS: Haralick feature values depend strongly on exposure, but invariance exists within defined domains of exposure. Here, we present an easily replicable procedure to identify said stable exposure domains, where image noise does not substantially add to textural feature values. Only by imaging at predetermined feature-invariant exposure levels and by adjusting exposure to expected activity concentrations, can textural features have a quantitative use in PET/CT. The necessary exposure levels are attainable by modern PET/CT systems in clinical routine.

Monte Carlo simulations and analysis of transmitted gamma ray spectra through various tissue phantoms.

Studies of radiation interactions with tissue equivalent material find importance in efforts that seek to avoid unjustifiable dose to patients, also in ensuring quality control of for instance nuclear medicine imaging equipment. Use of the Monte Carlo (MC) simulation tool in such characterization processes allows for the avoidance of costly experiments involving transmitted X- and γ-ray spectrometry. Present work investigates MC simulations of γ-ray transmission through tissue equivalent solid phantoms. Use has been made of a range of radionuclide gamma ray sources, 99mTc, 131I, 137Cs, 60Co (offering photons in the energy range from a few keV up to low MeV), popularly applied in medicine and in some cases for gauging in industry, obtaining the transmission spectra following their interaction with various phantom materials and thicknesses. In validation of the model, the simulated values of mass attenuation coefficients (µ/ρ) for different phantom materials and thicknesses were found to be in good agreement with reference values (NIST, 2004) to within 1.1% for all material compositions. For all of the primary photon energies and medium thicknesses of interest herein, results show that multiple scattering peaks are generally located at energies lower than 100 keV, although for the larger phantom thicknesses it is more difficult to distinguish single, double and multiple scattering in the gamma spectra. Transmitted photon spectra investigated for water, soft tissue, breast, brain and lung tissue slab phantoms are demonstrated to be practically independent of the phantom material, while a significant difference is observed for the spectra transmitted through bone that was proved to be due to the density effect and not material composition.

Wiener filter improves diagnostic accuracy of CAD SPECT images-comparison to angiography and CT angiography.

Many discrepancy in selection of proper filter and its parameters for individual cases exists. The authors investigate the impact of the most common filters on patient NM images with coronary artery disease (CAD), and compare the results with the computerized tomography (CT)-Angio and angiography for accuracy.The investigation initiated by performing various single photon emission computerized tomography (SPECT)/CT scan of the national electrical manufacturers association chest phantoms having hot and cold inserts. Data acquired on GE 670 PRO SPECT/CT; 360Ø, 64 frames, 60 seconds, low energy high resolution (LEHR) 128, low energy general purpose (LEGP) with CT attenuation (120 kV and 170 mA). The images reconstructed with filtered back projection and ITERATIVE ordered-subset expectation maximization utilizing filters; Hann, Butterworth, Metz, Hamming, and Wiener. The Image contrast was calculated to assess absolute nearness of the inserts. Based on the preliminary results, then scans of 92 patients with CAD; 64 males and 28 females, age 41 to 77 years old, who had been reported earlier reprocessed with the nominated filter and were reported by 2 NM expert. The results compared to the earlier reports and to the CT-Angio and angiography.The optimization suggested 3 filters; Wiener (Wi), Metz and Butterworth (But) provide the highest contrast (99- 66.4%) and (81- 32%) for the cold and hot inserts respectively, with the (Wi) filter to be the better option. The reprocessed patients scan with the (Wi) presented an elevated diagnostic accuracy, correlated well with the CT-Angio and angiography results (P < .001 and r = 0.79 for [Wi] and P = .004 and r = 0.39 for [But]). The percentage of the false negative for moderate to severe CAD cases reported using Wi filter reduced from 27% to 7% and similarly for mild CAD cases from 7% to 1%.It appears the Wiener filter could produce results with the highest contrast for phantom imaging of various cold and hot spheres and for the patient data which is more consistent with angiography results, with much-elevated accuracy in intermediate cases (r = 0.79 for Wiener and r = 0.39 for Butterworth vs angiography). However, the optimum parameters obtained for the filters have no relation with the resolution of the imaging system, but the details of the objects could be improved.

Technical Report: A Cost-Effective, Easily Available Tofu Model for Training Residents in Ultrasound-Guided Fine Needle Thyroid Nodule Targeting Punctures.

Objective: To establish a cost-effective and easily available phantom for training residents in ultrasound-guided fine needle thyroid nodule targeting punctures. Materials and Methods: Tofu, drinking straws filled with coupling gel, a urine tube, and 21-gauge needles were used to generate a phantom thyroid with nodules for training. Twelve radiology residents were involved in the study. The puncture success rates were recorded and compared before and after phantom training using the Wilcoxon signed-rank test. Results: On ultrasonography, tofu mimicked the texture of the thyroid. Drinking straws filled with coupling gel mimicked vessels. The urine tube filled with air mimicked the trachea, and 21-gauge needles mimicked small nodules in the transverse section. The entire phantom was similar to the structure of the thyroid and surrounding tissues. The puncture success rates of radiology residents were significantly increased from 34.4 ± 14.2% to 66.7 ± 19.5% after training (p = 0.003). The phantom was constructed in approximately 10 minutes and materials cost less than CNY 10 (approximately $ 1.5) at a local store. Conclusion: The tofu model was cost-effective, easily attainable, and effective for training residents in ultrasound-guided fine needle thyroid nodule targeting punctures in vitro.

Fusing Infrared and Visible Images of Different Resolutions via Total Variation Model.

In infrared and visible image fusion, existing methods typically have a prerequisite that the source images share the same resolution. However, due to limitations of hardware devices and application environments, infrared images constantly suffer from markedly lower resolution compared with the corresponding visible images. In this case, current fusion methods inevitably cause texture information loss in visible images or blur thermal radiation information in infrared images. Moreover, the principle of existing fusion rules typically focuses on preserving texture details in source images, which may be inappropriate for fusing infrared thermal radiation information because it is characterized by pixel intensities, possibly neglecting the prominence of targets in fused images. Faced with such difficulties and challenges, we propose a novel method to fuse infrared and visible images of different resolutions and generate high-resolution resulting images to obtain clear and accurate fused images. Specifically, the fusion problem is formulated as a total variation (TV) minimization problem. The data fidelity term constrains the pixel intensity similarity of the downsampled fused image with respect to the infrared image, and the regularization term compels the gradient similarity of the fused image with respect to the visible image. The fast iterative shrinkage-thresholding algorithm (FISTA) framework is applied to improve the convergence rate. Our resulting fused images are similar to super-resolved infrared images, which are sharpened by the texture information from visible images. Advantages and innovations of our method are demonstrated by the qualitative and quantitative comparisons with six state-of-the-art methods on publicly available datasets.

Measuring PET Spatial Resolution Using a Cylinder Phantom Positioned at an Oblique Angle.

A cylinder phantom positioned at a slightly oblique angle with respect to the z-axis of a PET scanner allows for fine sampling of the edge-spread function. We show how this technique can be used to measure the spatial resolution that can be expected with clinical PET protocols, potentially providing more relevant estimates than are typically obtained with established experimental procedures. Methods: A 20-cm-diameter water-filled cylinder phantom containing a uniform 18F solution was centrally positioned at a small angle with respect to the z-axis of a clinical PET/CT system. The oblique angle ensures that the phantom edge intersects the image matrix differently in different slices. Combining line profiles from multiple slices results in a composite profile with fine sampling. Spatial resolution was measured as the full width at half maximum (FWHM) by fitting a model to the finely sampled edge-spread functions in both radial and axial directions. The technique was validated by controlled modulation of image reconstruction parameters and by comparison with extended phantoms with fillable inserts. Separate experiments with uniform cylinders containing 18F, 11C, 13N, 68Ga, and 124I were used to further assess the proposed method. Results: Controlled adjustment of a gaussian postreconstruction filter was accurately reflected in the measured FWHM values. Recovery coefficients derived using the cylinder FWHM values agreed closely with recovery coefficients derived from physical phantoms over a range of insert-to-background ratios, phantom geometries, and reconstruction protocols. The effect of increasing positron energy was clearly reflected in the FWHM values measured with different isotopes. Conclusion: A method has been developed for measuring the spatial resolution that is achieved with clinical PET protocols, providing more relevant estimates than are typically obtained with established procedures. The proposed method requires no special equipment and is versatile, being capable of measuring resolution for different isotopes as well as for different reconstruction protocols. The new technique promises to aid standardization of PET data acquisition by allowing a more informed selection of reconstruction parameters.

Fuzzy Object Skeletonization: Theory, Algorithms, and Applications.

Skeletonization offers a compact representation of an object while preserving important topological and geometrical features. Literature on skeletonization of binary objects is quite mature. However, challenges involved with skeletonization of fuzzy objects are mostly unanswered. This paper presents a new theory and algorithm of skeletonization for fuzzy objects, evaluates its performance, and demonstrates its applications. A formulation of fuzzy grassfire propagation is introduced; its relationships with fuzzy distance functions, level sets, and geodesics are discussed; and several new theoretical results are presented in the continuous space. A notion of collision-impact of fire-fronts at skeletal points is introduced, and its role in filtering noisy skeletal points is demonstrated. A fuzzy object skeletonization algorithm is developed using new notions of surface- and curve-skeletal voxels, digital collision-impact, filtering of noisy skeletal voxels, and continuity of skeletal surfaces. A skeletal noise pruning algorithm is presented using branch-level significance. Accuracy and robustness of the new algorithm are examined on computer-generated phantoms and micro- and conventional CT imaging of trabecular bone specimens. An application of fuzzy object skeletonization to compute structure-width at a low image resolution is demonstrated, and its ability to predict bone strength is examined. Finally, the performance of the new fuzzy object skeletonization algorithm is compared with two binary object skeletonization methods.

Treatment Planning System Calculation Errors Are Present in Most Imaging and Radiation Oncology Core-Houston Phantom Failures.

PURPOSE: The anthropomorphic phantom program at the Houston branch of the Imaging and Radiation Oncology Core (IROC-Houston) is an end-to-end test that can be used to determine whether an institution can accurately model, calculate, and deliver an intensity modulated radiation therapy dose distribution. Currently, institutions that do not meet IROC-Houston's criteria have no specific information with which to identify and correct problems. In the present study, an independent recalculation system was developed to identify treatment planning system (TPS) calculation errors. METHODS AND MATERIALS: A recalculation system was commissioned and customized using IROC-Houston measurement reference dosimetry data for common linear accelerator classes. Using this system, 259 head and neck phantom irradiations were recalculated. Both the recalculation and the institution's TPS calculation were compared with the delivered dose that was measured. In cases in which the recalculation was statistically more accurate by 2% on average or 3% at a single measurement location than was the institution's TPS, the irradiation was flagged as having a "considerable" institutional calculation error. The error rates were also examined according to the linear accelerator vendor and delivery technique. RESULTS: Surprisingly, on average, the reference recalculation system had better accuracy than the institution's TPS. Considerable TPS errors were found in 17% (n=45) of the head and neck irradiations. Also, 68% (n=13) of the irradiations that failed to meet the IROC-Houston criteria were found to have calculation errors. CONCLUSIONS: Nearly 1 in 5 institutions were found to have TPS errors in their intensity modulated radiation therapy calculations, highlighting the need for careful beam modeling and calculation in the TPS. An independent recalculation system can help identify the presence of TPS errors and pass on the knowledge to the institution.

Validation of a CT-guided intervention robot for biopsy and radiofrequency ablation: experimental study with an abdominal phantom.

PURPOSE: We aimed to evaluate the accuracy of a needle-placement robot for biopsy and radiofrequency ablation on an abdominal phantom. METHODS: A master-slave robotic system has been developed that includes a needle-path planning system and a needle-inserting robot arm with computed tomography (CT) and CT fluoroscopy guidance. For evaluation of its accuracy in needle placement, a commercially available abdominal phantom (Model 057A; CIRS Inc.) was used. The liver part of the phantom contains multiple spherical simulated tumors of three different size spheres. Various needle insertion trials were performed in the transverse plane and caudocranial plane two nodule sizes (10 mm and 20 mm in diameter) to test the reliability of this robot. To assess accuracy, a CT scan was performed after each trial with the needle in situ. RESULTS: The overall error was 2 mm (0-2.6 mm), which was calculated as the distance from the planned trajectory before insertion to the actual needle trajectory after insertion. The standard deviations of the insertions on two nodules (10 mm and 20 mm in diameter) were 0.5 mm and 0.2 mm, respectively. CONCLUSION: The CT-compatible needle placement robot for biopsy and radiofrequency ablation shows relatively acceptable accuracy and could be used for radiofrequency ablation of nodules ≥10 mm under CT fluoroscopy guidance.

The importance of BMI in dosimetry of 153Sm-EDTMP bone pain palliation therapy: A Monte Carlo study.

Using digital phantoms as an atlas compared to acquiring CT data for internal radionuclide dosimetry decreases patient overall radiation dose and reduces the required analysis effort and time for organ segmentation. The drawback is that the phantom may not match exactly with the patient. We assessed the effect of varying BMIs on dosimetry results for a bone pain palliation agent, 153Sm-EDTMP. The simulation was done using the GATE Monte Carlo code. Female XCAT phantoms with the following different BMIs were employed: 18.6, 20.8, 22.1, 26.8, 30.3 and 34.7kg/m2. S-factors (mGy/MBq.s) and SAFs (kg-1) were calculated for the dosimetry of the radiation from major source organs including spine, ribs, kidney and bladder into different target organs as well as whole body dosimetry from spine. The differences in dose estimates from different phantoms compared to those from the phantom with BMI of 26.8kg/m2 as the reference, were calculated for both gamma and beta radiations. The relative differences (RD) of the S-factors or SAFs from the values of reference phantom were calculated. RDs greater than 10% and 100% were frequent in radiations to organs for photon and beta particles, respectively. The relative differences in whole body SAFs from the reference phantom were 15.4%, 7%, 4.2%, -9.8% and -1.4% for BMIs of 18.6, 20.8, 22.1, 30.3 and 34.7kg/m2, respectively. The differences in whole body S-factors for the phantoms with BMIs of 18.6, 20.8, 22.1, 30.3 and 34.7kg/m2 were 39.5%, 19.4%, 8.8%, -7.9% and -4.3%, respectively. The dosimetry of the gamma photons and beta particles changes substantially with the use of phantoms with different BMIs. The change in S-factors is important for dose calculation and can change the prescribed therapeutic dose of 153Sm-EDTMP. Thus a phantom with BMI better matched to the patient is suggested for therapeutic purposes where dose estimates closer to those in the actual patient are required.

Fast and Quantitative T1ρ-weighted Dynamic Glucose Enhanced MRI.

Common medical imaging techniques usually employ contrast agents that are chemically labeled, e.g. with radioisotopes in the case of PET, iodine in the case of CT or paramagnetic metals in the case of MRI to visualize the heterogeneity of the tumor microenvironment. Recently, it was shown that natural unlabeled D-glucose can be used as a nontoxic biodegradable contrast agent in Chemical Exchange sensitive Spin-Lock (CESL) magnetic resonance imaging (MRI) to detect the glucose uptake and potentially the metabolism of tumors. As an important step to fulfill the clinical needs for practicability, reproducibility and imaging speed we present here a robust and quantitative T1ρ-weighted technique for dynamic glucose enhanced MRI (DGE-MRI) with a temporal resolution of less than 7 seconds. Applied to a brain tumor patient, the new technique provided a distinct DGE contrast between tumor and healthy brain tissue and showed the detailed dynamics of the glucose enhancement after intravenous injection. Development of this fast and quantitative DGE-MRI technique allows for a more detailed analysis of DGE correlations in the future and potentially enables non-invasive diagnosis, staging and monitoring of tumor response to therapy.

Activity concentration measurements using a conjugate gradient (Siemens xSPECT) reconstruction algorithm in SPECT/CT.

The interest in quantitative single photon emission computer tomography (SPECT) shows potential in a number of clinical applications and now several vendors are providing software and hardware solutions to allow 'SUV-SPECT' to mirror metrics used in PET imaging. This brief technical report assesses the accuracy of activity concentration measurements using a new algorithm 'xSPECT' from Siemens Healthcare. SPECT/CT data were acquired from a uniform cylinder with 5, 10, 15 and 20 s/projection and NEMA image quality phantom with 25 s/projection. The NEMA phantom had hot spheres filled with an 8 : 1 activity concentration relative to the background compartment. Reconstructions were performed using parameters defined by manufacturer presets available with the algorithm. The accuracy of activity concentration measurements was assessed. A dose calibrator-camera cross-calibration factor (CCF) was derived from the uniform phantom data. In uniform phantom images, a positive bias was observed, ranging from ∼6% in the lower count images to ∼4% in the higher-count images. On the basis of the higher-count data, a CCF of 0.96 was derived. As expected, considerable negative bias was measured in the NEMA spheres using region mean values whereas positive bias was measured in the four largest NEMA spheres. Nonmonotonically increasing recovery curves for the hot spheres suggested the presence of Gibbs edge enhancement from resolution modelling. Sufficiently accurate activity concentration measurements can easily be measured on images reconstructed with the xSPECT algorithm without a CCF. However, the use of a CCF is likely to improve accuracy further. A manual conversion of voxel values into SUV should be possible, provided that the patient weight, injected activity and time between injection and imaging are all known accurately.

Characterization of a 0.35T MR system for phantom image quality stability and in vivo assessment of motion quantification.

ViewRay is a novel MR-guided radiotherapy system capable of imaging in near real-time at four frames per second during treatment using 0.35T field strength. It allows for improved gating techniques and adaptive radiotherapy. Three cobalt-60 sources (~ 15,000 Curies) permit multiple-beam, intensity-modulated radiation therapy. The primary aim of this study is to assess the imaging stability, accuracy, and automatic segmentation algorithm capability to track motion in simulated and in vivo targets. Magnetic resonance imaging (MRI) characteristics of the system were assessed using the American College of Radiology (ACR)-recommended phantom and accreditation protocol. Images of the ACR phantom were acquired using a head coil following the ACR scanning instructions. ACR recommended T1- and T2-weighted sequences were evaluated. Nine measurements were performed over a period of seven months, on just over a monthly basis, to establish consistency. A silicon dielectric gel target was attached to the motor via a rod. 40 mm total amplitude was used with cycles of 3 to 9 s in length in a sinusoidal trajectory. Trajectories of six moving clinical targets in four canine patients were quantified and tracked. ACR phantom images were analyzed, and the results were compared with the ACR acceptance levels. Measured slice thickness accuracies were within the acceptance limits. In the 0.35 T system, the image intensity uniformity was also within the ACR acceptance limit. Over the range of cycle lengths, representing a wide range of breathing rates in patients imaged at four frames/s, excellent agreement was observed between the expected and measured target trajectories. In vivo canine targets, including the gross target volume (GTV), as well as other abdominal soft tissue structures, were visualized with inherent MR contrast, allowing for preliminary results of target tracking.

Two Realistic Beagle Models for Dose Assessment.

Previously, the authors developed a series of eight realistic digital mouse and rat whole body phantoms based on NURBS technology to facilitate internal and external dose calculations in various species of rodents. In this paper, two body phantoms of adult beagles are described based on voxel images converted to NURBS models. Specific absorbed fractions for activity in 24 organs are presented in these models. CT images were acquired of an adult male and female beagle. The images were segmented, and the organs and structures were modeled using NURBS surfaces and polygon meshes. Each model was voxelized at a resolution of 0.75 × 0.75 × 2 mm. The voxel versions were implemented in GEANT4 radiation transport codes to calculate specific absorbed fractions (SAFs) using internal photon and electron sources. Photon and electron SAFs were then calculated for relevant organs in both models. The SAFs for photons and electrons were compatible with results observed by others. Absorbed fractions for electrons for organ self-irradiation were significantly less than 1.0 at energies above 0.5 MeV, as expected for many of these small-sized organs, and measurable cross irradiation was observed for many organ pairs for high-energy electrons (as would be emitted by nuclides like 32P, 90Y, or 188Re). The SAFs were used with standardized decay data to develop dose factors (DFs) for radiation dose calculations using the RADAR Method. These two new realistic models of male and female beagle dogs will be useful in radiation dosimetry calculations for external or internal simulated sources.

A three-dimensional head-and-neck phantom for validation of multimodality deformable image registration for adaptive radiotherapy.

PURPOSE: To develop a three-dimensional (3D) deformable head-and-neck (H&N) phantom with realistic tissue contrast for both kilovoltage (kV) and megavoltage (MV) imaging modalities and use it to objectively evaluate deformable image registration (DIR) algorithms. METHODS: The phantom represents H&N patient anatomy. It is constructed from thermoplastic, which becomes pliable in boiling water, and hardened epoxy resin. Using a system of additives, the Hounsfield unit (HU) values of these materials were tuned to mimic anatomy for both kV and MV imaging. The phantom opens along a sagittal midsection to reveal radiotransparent markers, which were used to characterize the phantom deformation. The deformed and undeformed phantoms were scanned with kV and MV imaging modalities. Additionally, a calibration curve was created to change the HUs of the MV scans to be similar to kV HUs, (MC). The extracted ground-truth deformation was then compared to the results of two commercially available DIR algorithms, from Velocity Medical Solutions and mim software. RESULTS: The phantom produced a 3D deformation, representing neck flexion, with a magnitude of up to 8 mm and was able to represent tissue HUs for both kV and MV imaging modalities. The two tested deformation algorithms yielded vastly different results. For kV-kV registration, mim produced mean and maximum errors of 1.8 and 11.5 mm, respectively. These same numbers for Velocity were 2.4 and 7.1 mm, respectively. For MV-MV, kV-MV, and kV-MC Velocity produced similar mean and maximum error values. mim, however, produced gross errors for all three of these scenarios, with maximum errors ranging from 33.4 to 41.6 mm. CONCLUSIONS: The application of DIR across different imaging modalities is particularly difficult, due to differences in tissue HUs and the presence of imaging artifacts. For this reason, DIR algorithms must be validated specifically for this purpose. The developed H&N phantom is an effective tool for this purpose.

Image simulation and a model of noise power spectra across a range of mammographic beam qualities.

PURPOSE: The aim of this work is to create a model to predict the noise power spectra (NPS) for a range of mammographic radiographic factors. The noise model was necessary to degrade images acquired on one system to match the image quality of different systems for a range of beam qualities. METHODS: Five detectors and x-ray systems [Hologic Selenia (ASEh), Carestream computed radiography CR900 (CRc), GE Essential (CSI), Carestream NIP (NIPc), and Siemens Inspiration (ASEs)] were characterized for this study. The signal transfer property was measured as the pixel value against absorbed energy per unit area (E) at a reference beam quality of 28 kV, Mo/Mo or 29 kV, W/Rh with 45 mm polymethyl methacrylate (PMMA) at the tube head. The contributions of the three noise sources (electronic, quantum, and structure) to the NPS were calculated by fitting a quadratic at each spatial frequency of the NPS against E. A quantum noise correction factor which was dependent on beam quality was quantified using a set of images acquired over a range of radiographic factors with different thicknesses of PMMA. The noise model was tested for images acquired at 26 kV, Mo/Mo with 20 mm PMMA and 34 kV, Mo/Rh with 70 mm PMMA for three detectors (ASEh, CRc, and CSI) over a range of exposures. The NPS were modeled with and without the noise correction factor and compared with the measured NPS. A previous method for adapting an image to appear as if acquired on a different system was modified to allow the reference beam quality to be different from the beam quality of the image. The method was validated by adapting the ASEh flat field images with two thicknesses of PMMA (20 and 70 mm) to appear with the imaging characteristics of the CSI and CRc systems. RESULTS: The quantum noise correction factor rises with higher beam qualities, except for CR systems at high spatial frequencies, where a flat response was found against mean photon energy. This is due to the dominance of secondary quantum noise in CR. The use of the quantum noise correction factor reduced the difference from the model to the real NPS to generally within 4%. The use of the quantum noise correction improved the conversion of ASEh image to CRc image but had no difference for the conversion to CSI images. CONCLUSIONS: A practical method for estimating the NPS at any dose and over a range of beam qualities for mammography has been demonstrated. The noise model was incorporated into a methodology for converting an image to appear as if acquired on a different detector. The method can now be extended to work for a wide range of beam qualities and can be applied to the conversion of mammograms.