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.

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.

Accurate determination of CT point-spread-function with high precision.

The measurement of modulation transfer functions (MTFs) in computed tomography (CT) is often performed by scanning a point source phantom such as a thin wire or a microbead. In these methods the region of interest (ROI) is generally placed on the scanned image to crop the point source response. The aim of the present study was to examine the effect of ROI size on MTF measurement, and to optimize the ROI size. Using a 4 multidetector-row CT, MTFs were measured by the wire and bead methods for three types of reconstruction kernels designated as 'smooth', 'standard', and 'edge-enhancement' kernels. The size of a square ROI was changed from 30 to 50 pixels (approximately 2.9 to 4.9 mm). The accuracies of the MTFs were evaluated using the verification method. The MTFs measured by the wire and bead methods were dependent on ROI size, particularly in MTF measurement for the 'edge-enhancement' kernel. MTF accuracy evaluated by the verification method changed with ROI size, and we were able to determine the optimum ROI size for each method (wire/bead) and for each kernel. Using these optimal ROI sizes, the MTF obtained by the wire method was in strong agreement with the MTF obtained by the bead method in each kernel. Our data demonstrate that the difficulties in obtaining accurate MTFs for some kernels such as edge-enhancement can be overcome by incorporating the verification method into the wire and bead methods, allowing optimization of the ROI size to accurately determine the MTF.

Validation of measurement-guided 3D VMAT dose reconstruction on a heterogeneous anthropomorphic phantom.

3DVH software (Sun Nuclear Corp., Melbourne, FL) is capable of generating a volumetric patient VMAT dose by applying a volumetric perturbation algorithm based on comparing measurement-guided dose reconstruction and TPS-calculated dose to a cylindrical phantom. The primary purpose of this paper is to validate this dose reconstruction on an anthropomorphic heterogeneous thoracic phantom by direct comparison to independent measurements. The dosimetric insert to the phantom is novel, and thus the secondary goal is to demonstrate how it can be used for the hidden target end-to-end testing of VMAT treatments in lung. A dosimetric insert contains a 4 cm diameter unit-density spherical target located inside the right lung (0.21 g/cm(3) density). It has 26 slots arranged in two orthogonal directions, milled to hold optically stimulated luminescent dosimeters (OSLDs). Dose profiles in three cardinal orthogonal directions were obtained for five VMAT plans with varying degrees of modulation. After appropriate OSLD corrections were applied, 3DVH measurement-guided VMAT dose reconstruction agreed 100% with the measurements in the unit density target sphere at 3%/3 mm level (composite analysis) for all profile points for the four less-modulated VMAT plans, and for 96% of the points in the highly modulated C-shape plan (from TG-119). For this latter plan, while 3DVH shows acceptable agreement with independent measurements in the unit density target, in the lung disagreement with experiment is relatively high for both the TPS calculation and 3DVH reconstruction. For the four plans excluding the C-shape, 3%/3 mm overall composite analysis passing rates for 3DVH against independent measurement ranged from 93% to 100%. The C-shape plan was deliberately chosen as a stress test of the algorithm. The dosimetric spatial alignment hidden target test demonstrated the average distance to agreement between the measured and TPS profiles in the steep dose gradient area at the edge of the 2 cm target to be 1.0 ± 0.7, 0.3 ± 0.3, and 0.3 ± 0.3 mm for the IEC X, Y, and Z directions, respectively.

Development of 5- and 10-year-old pediatric phantoms based on polygon mesh surfaces.

PURPOSE: The purpose of this study is the development of reference pediatric phantoms for 5- and 10-year-old children to be used for the calculation of organ and tissue equivalent doses in radiation protection. METHODS: The study proposes a method for developing anatomically highly sophisticated pediatric phantoms without using medical images. The 5- and 10-year-old male and female phantoms presented here were developed using 3D modeling software applied to anatomical information taken from atlases and textbooks. The method uses polygon mesh surfaces to model body contours, the shape of organs as well as their positions, and orientations in the human body. Organ and tissue masses comply with the corresponding data given by the International Commission on Radiological Protection (ICRP) for the 5- and 10-year-old reference children. Bones were segmented into cortical bone, spongiosa, medullary marrow, and cartilage to allow for the use of micro computer tomographic (microCT) images of trabecular bone for skeletal dosimetry. RESULTS: The four phantoms, a male and a female for each age, and their organs are presented in 3D images and their organ and tissue masses in tables which show the compliance of the ICRP reference values. Dosimetric data, calculated for the reference pediatric phantoms by Monte Carlo methods were compared with corresponding data from adult mesh phantoms and pediatric stylized phantoms. The comparisons show reasonable agreement if the anatomical differences between the phantoms are properly taken into account. CONCLUSIONS: Pediatric phantoms were developed without using medical images of patients or volunteers for the first time. The models are reference phantoms, suitable for regulatory dosimetry, however, the 3D modeling method can also be applied to medical images to develop patient-specific phantoms.

Assessing the clinical utility of quantitative computed tomography with a routinely used diagnostic computed tomography scanner in a cancer center.

The purpose of this study was to characterize quantitative computed tomography (QCT) in our multi-detector computed tomography (MDCT) scanner with regard to the influence of the QCT phantom on dose and the influence of varying mA and CIRS phantom size on bone mineral density (BMD) measurements. We accomplish this by scanning a commercially available QCT phantom and a corresponding quality assurance phantom. To assess the feasibility of having the QCT phantom in place while patients are being scanned, we measured radiation dose difference in a CT body phantom with and without the QCT phantom on the CT table and also, with and without the use of dose modulation programs. We also analyzed reconstructed QCT phantom images with the manufacturer's software to measure BMD. Although patient characteristics may be different, leading to different mA values, the influence of the QCT phantom on the dose to patients was minimal when compared with doses measured without the phantom in place. Average BMD measurements were not significantly affected by varying mA, for a fixed-size phantom. The average BMD exhibited a weak dependence on computerized imaging reference systems (CIRS) torso phantom size, with a propensity for decreasing BMD with increasing size. Measurement precision was unaffected by varying CIRS size. Having the ability to measure bone density as part of the routine management of cancer patients, with no added cost, time, or radiation dose, will allow for the prospective evaluation of bone mineral changes. We believe that this ability will facilitate the detection of abnormal bone loss and will lead to better management of this loss and, thus, reduce the complications and associated morbidity in these cancer survivors.

Graphics processing unit parallel accelerated solution of the discrete ordinates for photon transport in biological tissues.

As a widely used numerical solution for the radiation transport equation (RTE), the discrete ordinates can predict the propagation of photons through biological tissues more accurately relative to the diffusion equation. The discrete ordinates reduce the RTE to a serial of differential equations that can be solved by source iteration (SI). However, the tremendous time consumption of SI, which is partly caused by the expensive computation of each SI step, limits its applications. In this paper, we present a graphics processing unit (GPU) parallel accelerated SI method for discrete ordinates. Utilizing the calculation independence on the levels of the discrete ordinate equation and spatial element, the proposed method reduces the time cost of each SI step by parallel calculation. The photon reflection at the boundary was calculated based on the results of the last SI step to ensure the calculation independence on the level of the discrete ordinate equation. An element sweeping strategy was proposed to detect the calculation independence on the level of the spatial element. A GPU parallel frame called the compute unified device architecture was employed to carry out the parallel computation. The simulation experiments, which were carried out with a cylindrical phantom and numerical mouse, indicated that the time cost of each SI step can be reduced up to a factor of 228 by the proposed method with a GTX 260 graphics card.

Monte Carlo simulation of the effect of miniphantom on in-air output ratio.

PURPOSE: The aim of the study was to quantify the effect of miniphantoms on in-air output ratio measurements, i.e., to determine correction factors for in-air output ratio. METHODS: Monte Carlo (MC) simulations were performed to simulate in-air output ratio measurements by using miniphantoms made of various materials (PMMA, graphite, copper, brass, and lead) and with different longitudinal thicknesses or depths (2-30 g/cm2) in photon beams of 6 and 15 MV, respectively, and with collimator settings ranging from 3 x 3 to 40 x 40 cm2. EGSnrc and BEAMnrc (2007) software packages were used. Photon energy spectra corresponding to the collimator settings were obtained from BEAMnrc code simulations on a linear accelerator and were used to quantify the components of in-air output ratio correction factors, i.e., attenuation, mass energy absorption, and phantom scatter correction factors. In-air output ratio correction factors as functions of miniphantom material, miniphantom longitudinal thickness, and collimator setting were calculated and compared to a previous experimental study. RESULTS: The in-air output ratio correction factors increase with collimator opening and miniphantom longitudinal thickness for all the materials and for both energies. At small longitudinal thicknesses, the in-air output ratio correction factors for PMMA and graphite are close to 1. The maximum magnitudes of the in-air output ratio correction factors occur at the largest collimator setting (40 x 40 cm2) and the largest miniphantom longitudinal thickness (30 g/cm2): 1.008 +/- 0.001 for 6 MV and 1.012 +/- 0.001 for 15 MV, respectively. The MC simulations of the in-air output ratio correction factor confirm the previous experimental study. CONCLUSIONS: The study has verified that a correction factor for in-air output ratio can be obtained as a product of attenuation correction factor, mass energy absorption correction factor, and phantom scatter correction factor. The correction factors obtained in the present study can be used in studies involving in-air output ratio measurements using miniphantoms.

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.

Magnetic field effects on the energy deposition spectra of MV photon radiation.

Several groups worldwide have proposed various concepts for improving megavoltage (MV) radiotherapy that involve irradiating patients in the presence of a magnetic field-either for image guidance in the case of hybrid radiotherapy-MRI machines or for purposes of introducing tighter control over dose distributions. The presence of a magnetic field alters the trajectory of charged particles between interactions with the medium and thus has the potential to alter energy deposition patterns within a sub-cellular target volume. In this work, we use the MC radiation transport code PENELOPE with appropriate algorithms invoked to incorporate magnetic field deflections to investigate electron energy fluence in the presence of a uniform magnetic field and the energy deposition spectra within a 10 microm water sphere as a function of magnetic field strength. The simulations suggest only very minor changes to the electron fluence even for extremely strong magnetic fields. Further, calculations of the dose-averaged lineal energy indicate that a magnetic field strength of at least 70 T is required before beam quality will change by more than 2%.

Modelling optical scattering artefacts for varying pathlength in a gel dosimeter phantom.

A gelatin phantom containing an optically scattering funnel-shaped region of elevated optical density (OD) was used to examine light-scattering-induced artefacts in a cone-beam optical CT scanner used for gel dosimetry. To simulate polymer gel dosimeters, the opacity was introduced by adding a colloidal scatterer to the gelatin. Scatter results in an underestimate of OD (hence dose). In line profiles of OD taken from 3D reconstructions of the funnel, those profiles with a long pathlength through high OD regions exhibited a 'dishing' (or 'cupping') artefact, while those of short pathlength exhibited the opposite effect-'doming'. These phenomena are accounted for by a model that includes the effect of stray, scattered light.

Investigation of detection limits for diffuse optical tomography systems: II. Analysis of slab and cup geometry for breast imaging.

Using a statistical (chi-square) test on simulated data and a realistic noise model derived from the system's hardware we study the performance of diffuse optical tomography systems for fluorescence imaging. We compare the predicted smallest size of detectable lesions at various positions in slab and cup geometry and model how detection sensitivity depends on breast compression and lesion fluorescence contrast. Our investigation shows that lesion detection is limited by relative noise in slab geometry and by absolute noise in cup geometry.

Experimental evaluation of a simple lesion detection task with time-of-flight PET.

A new generation of high-performance, time-of-flight (TOF) PET scanners have recently been developed. In earlier works, the gain with TOF information was derived as a reduction of noise in the reconstructed image, or essentially a gain in scanner sensitivity. These derivations were applicable to analytical reconstruction techniques and 2D PET imaging. In this work, we evaluate the gain measured in the clinically relevant task of lesion detection with TOF information in fully 3D PET scanners using iterative reconstruction algorithms. We performed measurements in a fully 3D TOF PET scanner using spherical lesions in uniform, cylindrical phantom. Lesion detectability was estimated for 10 mm diameter lesions using a non-prewhitening matched filter signal-to-noise-ratio (NPW SNR) as the metric. Our results show that the use of TOF information leads to increased lesion detectability, which is achieved with less number of iterations of the reconstruction algorithm. These phantom results indicate that clinically, TOF PET will allow reduced scan times and improved lesion detectability, especially in large patients.

Delivery of modulated electron beams with conventional photon multi-leaf collimators.

Electron beam radiotherapy is an accepted method to treat shallow tumors. However, modulation of electrons to customize dose distributions has not readily been achieved. Studies of bolus and tertiary collimation systems have been met with limitations. We pursue the use of photon multi-leaf collimators (MLC) for modulated electron radiotherapy (MERT) to achieve customized distributions for potential clinical use. As commercial planning systems do not support the use of MLC with electrons, planning was conducted using Monte Carlo calculations. Segmented and dynamic modulated delivery of multiple electron segments was configured, calculated and delivered for validation. Delivery of electrons with segmented or dynamic leaf motion was conducted. A phantom possessing an idealized stepped target was planned and optimized with subsequent validation by measurements. Finally, clinical treatment plans were conducted for post-mastectomy and cutaneous lymphoma of the scalp using forward optimization techniques. Comparison of calculations and measurements was successful with agreement of +/-2%/2 mm for the energies, segment sizes, depths tested for delivered segments for the dynamic and segmented delivery. Clinical treatment plans performed provided optimal dose coverage of the target while sparing distal organs at risk. Execution of plans using an anthropomorphic phantom to ensure safe and efficient delivery was conducted. Our study validates that MERT is not only possible using the photon MLC, but the efficient and safe delivery inherent with the dynamic delivery provides an ideal technique for shallow tumor treatment.

Nonlinear reconstruction of absorption and fluorescence contrast from measured diffuse transmittance and reflectance of a compressed-breast-simulating phantom.

We report on the nonlinear reconstruction of local absorption and fluorescence contrast in tissuelike scattering media from measured time-domain diffuse reflectance and transmittance of laser as well as laser-excited fluorescence radiation. Measurements were taken at selected source-detector offsets using slablike diffusely scattering and fluorescent phantoms containing fluorescent heterogeneities. Such measurements simulate in vivo data that would be obtained employing a scanning, time-domain fluorescence mammograph, where the breast is gently compressed between two parallel glass plates, and source and detector optical fibers scan synchronously at various source-detector offsets, allowing the recording of laser and fluorescence mammograms. The diffusion equations modeling the propagation of the laser and fluorescence radiation were solved in frequency domain by the finite element method simultaneously for several modulation frequencies using Fourier transformation and preprocessed experimental data. To reconstruct the concentration of the fluorescent contrast agent, the Born approximation including higher-order reconstructed photon densities at the excitation wavelength was used. Axial resolution was determined that can be achieved by various detection schemes. We show that remission measurements increase the depth resolution significantly.

A technique for simulating the effect of dose reduction on image quality in digital chest radiography.

PURPOSE: The purpose of this study is to provide a pragmatic tool for studying the relationship between dose and image quality in clinical chest images. To achieve this, we developed a technique for simulating the effect of dose reduction on image quality of digital chest images. MATERIALS AND METHODS: The technique was developed for a digital charge-coupled-device (CCD) chest unit with slot-scan acquisition. Raw pixel values were scaled to a lower dose level, and a random number representing noise to each specific pixel value was added. After adding noise, raw images were post processed in the standard way. Validation was performed by comparing pixel standard deviation, as a measure of noise, in simulated images with images acquired at actual lower doses. To achieve this, a uniform test object and an anthropomorphic phantom were used. Additionally, noise power spectra of simulated and actual images were compared. Also, detectability of simulated lesions was investigated using a model observer. RESULTS: The mean difference in noise values between simulated and real lower-dose phantom images was smaller than 5% for relevant clinical settings. Noise power spectra appeared to be comparable on average but simulated images showed slightly higher noise levels for higher spatial frequencies and slightly lower noise levels for lower spatial frequencies. Comparable detection performance was shown in simulated and actual images with slightly worse detectability for simulated lower dose images. CONCLUSION: We have developed and validated a method for simulating dose reduction. Our method seems an acceptable pragmatic tool for studying the relationship between dose and image quality.

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.

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.