Physical characteristics of scattered radiation in diagnostic radiology: Monte Carlo simulation studies

We applied Monte Carlo methods for the simulation of x-ray scattering in water phantoms. The phantom thickness was varied from 5 to 20 cm, and the monoenergetic incident x rays were varied from 15 to 100 keV. Eight screen pairs and a total absorption system were used as x-ray receptors. We determined the angular, spectral, and spatial distributions of the scattered radiation and the scatter fractions recorded in the image plane. The dependence of these properties on the incident x-ray energy, the phantom thickness, and the energy response of the recording system was examined. The results of this study provide useful information for the development of antiscatter techniques and for the evaluation of radiographic procedures.     

Studies of x-ray energy absorption and quantum noise properties of x-ray screens by use of Monte Carlo simulation

The imaging properties of the phosphor layer in fluorescent screens or image intensifiers are related to its x-ray absorption characteristics. In this study, we applied Monte Carlo methods for the simulation of x-ray photon diffusion in a phosphor layer. The K-reabsorption factor, absorbed x-ray energy, quantum absorption efficiency, statistical factor, and noise-equivalent absorption were determined as a function of the incident energy and angle of the x rays for eight commonly used phosphor layers. These basic physical quantities will be useful for the prediction of the information transfer properties of a phosphor layer.     

Optimum momentum transfer arguments for x-ray forward scatter imaging

In our research program we have shown through modeling, related numerical calculations, and experimental measurements that there exists a potential use of scattered radiation for medical x-ray imaging. Each incident photon of wavelength lambda which scatters at a small angle theta with respect to its initial direction of travel has a change in momentum characterized by the photon momentum transfer argument x = lambda(-1) sin(theta/2). In this work, we show that in order to maximize the signal-to-noise ratio (SNR) obtained with scattered x rays, one must detect photons with specific x values. Using a photon counting detector to distinguish 2-cm-thick polymethyl methacrylate and nylon targets situated within a 15-cm-diam spherical water phantom with an 80 kV beam yields experimentally SNR/square root(K(air)c) = 12.8 +/- 0.2 (mJ/kg)(-1/2) when using the photons between x = 0.5 and 0.7 nm(-1). Here K(air)c is the air collision kerma and the average momentum transfer argument, x, is calculated by weighting x by the incident photon fluence distribution. The model predicts a value of SNR/square root(K(air)c) = 12.9 (mJ/kg)(-1/2). If we choose to form the signal with the range in x extended to be from 0.5 to 1.0 nm(-1) then, despite the detection of more scattered photons, experimentally SNR/square root(K(air)c) decreases by 38% to 7.9 +/- 0.3 (mJ/kg)(-1/2). The model predicts a value of 9.46 (mJ/kg)(-1/2). Results for energy integrating detectors are in general similar to those for photon counters, but there exist cases where a significant decrease in SNR can occur. For example, for measurements in air with the two plastics at theta = 3 degrees the SNR for an energy integrator was found to be 52% that of a photon counter. Numerical calculations predict that the effects of spectral blur can be significant when a narrow angular range is used for detection. Preliminary numerical predictions for breast tissues suggest a potential use of x-ray scatter in the field of mammography.     

Photon scatter in portal images: physical characteristics of pencil beam kernels generated using the EGS Monte Carlo code

Pencil beam kernels describing scattered photon fluence behind homogeneous water slabs at various air gap distances were generated using the EGS Monte Carlo code. Photon scatter fluence was scored in separate bins based on the particle's history: singly scattered, multiply scattered, and bremsstrahlung and positron annihilation photons. Simultaneously, the mean energy and mean angle with respect to the incident photon pencil beam were tallied. Kernels were generated for incident photon pencil beams exhibiting monoenergetic spectra of 2.0 and 10.0 MeV, and polyenergetic spectra representative of 6 and 24 MV beams. Reciprocity was used to generate scatter fractions on the central axis for various field sizes, phantom thicknesses, and air gaps. The scatter kernels were further characterized by full width at half-maximum estimates. Modulation transfer functions were calculated, providing theoretical estimates of the limit of performance of portal imaging systems due to the intrinsic scattering of photon radiation through the patient.     

Compton-scattering measurement of diagnostic x-ray spectrum using high-resolution Schottky CdTe detector

The analysis of x-ray spectra is important for quality assurance (QA) and quality control (QC) of radiographic systems. The aim of this study is to measure the diagnostic x-ray spectra under clinical conditions using a high-resolution Schottky CdTe detector. Under clinical conditions, the direct measurement of a diagnostic spectrum is difficult because of the high photon fluence rates that cause significant detector photon pile-up. An alternative way of measuring the output spectra from a tube is first to measure the 90 deg Compton scattered photons from a given sample. With this set-up detector, pile-up is not a problem. From the scattered spectrum one can then use an energy correction and the Klein-Nishina function to reconstruct the actual spectrum incident upon the scattering sample. The verification of whether our spectra measured by the Compton method are accurate was accomplished by comparing exposure rates calculated from the reconstructed spectra to those measured with an ionization chamber. We used aluminum (Al) filtration ranging in thickness from 0 to 6 mm. The half value layers (HVLs) obtained for a 70 kV beam were 2.78 mm via the ionization chamber measurements and 2.93 mm via the spectral measurements. For a 100 kV beam we obtained 3.98 and 4.32 mm. The small differences in HVLs obtained by both techniques suggest that Compton scatter spectroscopy with a Schottky CdTe detector is suitable for measuring the diagnostic x-ray spectra and useful for QA and QC of clinical x-ray equipment.     

A three-dimensional x-ray scattering system for multi-parameter imaging of the human head

This work examines the suitability of a non-rotating one-side 3D x-ray scatter system for imaging the human head. The system simultaneously produces images of the x-ray attenuation coefficients at two photon energies, as well as an image of the electron density. The system relies on measuring the scattered radiation at two directions orthogonal to an incident beam that scans the object from one side, in addition to the traditionally recorded transmitted radiation. Algorithms for this multi-parameter imaging process are presented, and their numerical viability is demonstrated, using both idealized detector responses and those independently estimated from Monte Carlo simulations. The absorbed radiation dose is also calculated, and was shown to be about one quarter of that of conventional CT systems, for 5 mm spatial-resolution images. The introduced system can therefore be useful in radiotherapy planning, and in post-treatment imaging.     

Design of synchrotron light source and its beamline dedicated to dual-energy x-ray computed tomography

A synchrotron light source dedicated to medical applications has been designed at National Institute of Radiological Sciences. The storage ring, with circumference of 80 m, is designed for acceleration of 2.3 GeV and a stored current of 420 mA. It is equipped with two multipole wigglers to produce sufficient photon flux in a hard x-ray region required for medical applications. The purposes of the synchrotron light source are clinical performance of medical diagnoses clinically and research and development relating with medical applications. One of the most interesting applications for us is dual-energy x-ray computed tomography (CT). It gives the information about electron density of human tissue. The information plays an important role in advancing heavy-ion radiotherapy of cancers. Electron density can be derived from attenuation coefficients measured by different energy x rays. In this paper, a practical method of the dual-energy x-ray CT with synchrotron radiation is proposed with the theoretical consideration. The primitive experiment using monochromatic x rays emitted from radioisotopes proved the procedure of analysis mentioned here effective to derive electron densities from linear attenuation coefficients for two x rays of a different energy. The beamline dedicated to dual-energy x-ray CT is also proposed. It has a multipole wiggler as a light source and it mainly consists of a dual crystal monochromator and a rotating filter for attenuating photon flux of x rays and two-dimensional detector.     

Fundamental x-ray interaction limits in diagnostic imaging detectors: spatial resolution

The practice of diagnostic x-ray imaging has been transformed with the emergence of digital detector technology. Although digital systems offer many practical advantages over conventional film-based systems, their spatial resolution performance can be a limitation. The authors present a Monte Carlo study to determine fundamental resolution limits caused by x-ray interactions in four converter materials: Amorphous silicon (a-Si), amorphous selenium, cesium iodide, and lead iodide. The "x-ray interaction" modulation transfer function (MTF) was determined for each material and compared in terms of the 50% MTF spatial frequency and Wagner's effective aperture for incident photon energies between 10 and 150 keV and various converter thicknesses. Several conclusions can be drawn from their Monte Carlo study. (i) In low-Z (a-Si) converters, reabsorption of Compton scatter x rays limits spatial resolution with a sharp MTF drop at very low spatial frequencies (< 0.3 cycles/mm), especially above 60 keV; while in high-Z materials, reabsorption of characteristic x rays plays a dominant role, resulting in a mid-frequency (1-5 cycles/mm) MTF drop. (ii) Coherent scatter plays a minor role in the x-ray interaction MTF. (iii) The spread of energy due to secondary electron (e.g., photoelectrons) transport is significant only at very high spatial frequencies. (iv) Unlike the spread of optical light in phosphors, the spread of absorbed energy from x-ray interactions does not significantly degrade spatial resolution as converter thickness is increased. (v) The effective aperture results reported here represent fundamental spatial resolution limits of the materials tested and serve as target benchmarks for the design and development of future digital x-ray detectors.     

Monte Carlo studies of x-ray scattering in transmission diagnostic radiology

Monte Carlo methods have been used to simulate the scattering of x rays in polystyrene and water phantoms. In particular, the ratio of the scattered to total x-ray fluence (scatter fraction) has been calculated for monoenergetic x-ray beams in the energy region relevant to diagnostic radiology and nuclear medicine (30-660 keV). Simulations have been made for representative values of the pertinent geometrical factors; phantom thickness from 5 to 21 cm, x-ray beam diameters of 10 and 25 cm, and scatterer-to-image-plane separations from 0 to 20 cm. As a function of x-ray energy, the scatter fraction was found to vary slowly between 30 and 100 keV, and to decrease between 100 and 660 keV. The present results were generated with a special transport code which included the effects of special geometries and the response of the x-ray detector. With the inclusion of these effects, the results resolved inconsistencies and showed good agreement with previous measured and calculated data.     

Analysis of spectral blur effects in x-ray scatter imaging.

Previous analysis in our research program investigating the potential use of scattered photons for medical x-ray imaging has been for monoenergetic beams. In practice, polyenergetic beams are almost always used due to their higher photon fluence rate. The effects of beam polychromaticity on x-ray scatter imaging are determined with the aid of our semianalytic model that images a target object against a background material of the same dimensions when both are situated within a water phantom. Our analysis involves four different photon beams with constant incident energy fluence: (1) a monoenergetic beam with photon energy E0, (2) a dual peak beam with two separate monoenergetic peaks of energies E1 and E2, (3) a clinical x-ray beam, and (4) a rectangular beam with uniform energy fluence between energies Emin and Emax. A comparison between the polyenergetic spectra is accomplished by matching the centroids and standard deviations of the dual peak and rectangular spectra to those of the clinical x-ray spectrum. For the task of imaging liver versus fat structures 1 cm thick in a 25-cm-diam spherical water phantom with the scattered photons between 2 degrees and 12 degrees, the predicted signal-to-noise ratio (SNR) obtained with a 100 kV beam is 87.5% of the SNR acquired with the optimum monoenergetic beam (SNRopt). The SNR for the corresponding dual peak beam is 84.4% of SNRopt and for the rectangular beam is 86.3%. Our analysis shows that monoenergetic x-ray beams are not necessary for x-ray scatter imaging.     

Monte Carlo simulation of diagnostic x-ray scatter

A Monte Carlo method was developed and implemented to simulate x-ray photon transport. Simulations consisted of a pencil beam of monoenergetic photons with energies from 50 to 110 keV incident on water and aluminum slabs. The dependence of scatter fraction and multiple scattering on x-ray energy, scatterer thickness, and material is reported in both number and energy fluence. The average energy of scattered photons reaching the detector plane is also reported. Comparisons are made to previous x-ray scatter computations.     

Changes in incident photon fluence of 6 and 18 MV x rays caused by blocks and block trays.

When a block and tray are placed in a x-ray beam the dose to a point in a phantom is changed by the following factors: (1) attenuation of photon and electron fluence from the head of the accelerator by the tray and the block, (2) decrease in the scatter in the phantom by a reduction in the phantom volume that receives radiation, and (3) generation of scatter off the tray and block. This third factor is generally ignored in dosimetry calculation but has been measured in this work. Measurements of incident photon fluence for 6 and 18 MV x rays were made with a columnar miniphantom of 10 cm depth. The tray factor for a 9 mm thick Lexan tray is found to be variable and to increase by 1.8% due to scatter off the tray when the field size is increased from a 3cm x 3 cm to 40cm x 40 cm field. Also, it was found that scatter off a block could increase the incident photon fluence by as much as 2%. The magnitude of this block scatter depends on the length of the inner edge of the opening in the block and on amount of block that is being irradiated, the overlap of the block by the radiation field. The total block-tray factor can be as much as 3% larger than the single-value tray factor measured with a 10cm x 10cm field that is traditionally used. An analytical equation is developed that accurately models the block-tray factor.     

Blurring artifacts in megavoltage radiography with a flat-panel imaging system: comparison of Monte Carlo simulations with measurements

Originally designed for use at medical-imaging x-ray energies, imaging systems comprising scintillating screens and amorphous Si detectors are also used at the megavoltage photon energies typical of portal imaging and industrial radiography. While image blur at medical-imaging x-ray energies is strongly influenced both by K-shell fluorescence and by the transport of optical photons within the scintillator layer, at higher photon energies the image blur is dominated by radiation scattered from the detector housing and internal support structures. We use Monte Carlo methods to study the blurring in a notional detector: a series of semi-infinite layers with material compositions, thicknesses, and densities similar to those of a commercially available flat-panel amorphous Si detector system comprising a protective housing, a gadolinium oxysulfide scintillator screen, and associated electronics. We find that the image blurring, as described by a point-spread function (PSF), has three length scales. The first component, with a submillimeter length scale, arises from electron scatter within the scintillator and detection electronics. The second component, with a millimeter-to-centimeter length scale, arises from electrons produced in the front cover of the detector. The third component, with a length scale of tens of centimeters, arises from photon scatter by the back cover of the detector. The relative contributions of each of these components to the overall PSF vary with incident photon energy. We present an algorithm that includes the energy-dependent sensitivity and energy-dependent PSF within a ray-tracing formalism. We find quantitative agreement (approximately 2%) between predicted radiographs with radiographs of copper step wedges, taken with a 9 MV bremsstrahlung source and a commercially available flat-panel system. The measured radiographs show the blurring artifacts expected from both the millimeter-scale electron transport and from the tens-of-centimeters length scale arising from the scattered photon transport. Calculations indicate that neglect of the energy-dependent blurring would lead to discrepancies in the apparent transmission of these wedges of the order of 9%.     

X-ray scatter background signals in transmission radiography.

With monoenergetic x-ray beams incident on polystyrene phantoms, the spectra of the tramsmitted x rays were measured with a scintillation spectrometer. The scattered and unscattered components of the transmitted x-ray fluence at a point on the beam axis were determined as a function of (i) the incident x-ray energy (18, 22, 32, 49, 58, 69, and 660 keV), (ii) the phantom thickness (5.3, 10, and 21 cm), (iii) the scatter solid angle determined by the exposed area of the phantom and the separation distance of the image plane (0.090, 0.31, 0.66, 1.8, 3.5 4.3, 4.8, and 5.1 sr), and (iv) the beam diameter at the image plane (25, 17, and 10 cm). The results indicate that, as the incident x-ray energy decreases from 660 to 30 keV, the contribution of the scattered component to the transmitted fluence increases from approximately 50% to 90% for the 21-cm phantom and from 21% to 50% for the 5.3-cm phantom. For typical cases, the data show the effect of the scatter component on the ratio of the image to the background signals. In addition, the examples show that optimum conditions for maximizing this signal ratio may be obtained by a careful selection of the incident x-ray energy for low-, medium, and high-contrast objects.     

Electronic portal imaging based on cerenkov radiation: a new approach and its feasibility

Most electronic portal imaging devices (EPIDs) developed so far use a Cu plate/phosphor screen to absorb x rays and convert their energies into light, and the light image is then read out. The main problem with this approach is that the Cu plate/phosphor screen must be thin (approximately 2 mm thick) in order to obtain a high spatial resolution, resulting in a low x-ray absorption or low quantum efficiency for megavoltage x rays (typically 2-4%). In addition, the phosphor screen contains high atomic number (high-Z) materials, resulting in an over-response of the detector to low-energy x rays in dosimetric verification. In this paper, we propose a new approach that uses Cerenkov radiation to convert x-ray energy absorbed by the detector into light for portal imaging applications. With our approach, a thick (approximately 10-30 cm) energy conversion layer made of a low-Z dielectric medium, such as a large-area, thick fiber-optic taper consisting of a matrix of optical fibers aligned with the incident x rays, is used to replace the thin Cu plate/phosphor screen. The feasibility of this approach has been investigated using a single optical fiber embedded in a solid material. The spatial resolution expressed by the modulation transfer function (MTF) and the sensitivity of the detector at low doses (approximately one Linac pulse) have been measured. It is predicted that, using this approach, a detective quantum efficiency of an order of magnitude higher at zero frequency can be obtained while maintaining a reasonable MTF, as compared to current EPIDs.     

Characteristics of low-angle x-ray scattering from some biological samples

Design of medical imaging devices based on the detection of low-angle coherent scattering is a subject of increasing interest. The technique is based on the differences in the distribution of photons coherently scattered from different body tissues. Coherent scattering is also useful in monitoring changes that may occur in a healthy tissue (e.g. carcinoma). In this work, low angle scattering properties of some tissues and tissue-equivalent materials are studied. Special care is given to the possibility of distinguishing between tissues of similar water content (e.g. muscle and blood). For this purpose, a Monte Carlo simulation is updated, introducing molecular form factor data, which include molecular interference effects. This program is used to simulate the angular distribution of scattered photons from two tissue-equivalent materials (lucite and water) and three biological samples (muscle, fat and blood). Simulation results agree well with previously measured angular distributions of scattered photons at 59.54 keV. Scattering from water and lucite is also measured at 8.047 keV. The effects of scattering geometry, sample thickness, incident photon energy and tissue type on the angular distribution of scattered photons are investigated. Results reveal the potential of measuring the full width at half maximum (FWHM) of the scattered photon distribution for tissue characterization. Energies up to 13 keV and sample thickness of 0.3 cm reported maximum differences between investigated samples. These conditions are expected to maximize the potential of using coherent scattering set-ups to monitor changes in biological samples even if their water contents are similar. Present results may act as a guide for the optimization of coherent scattering imaging systems.     

An energy- and depth-dependent model for x-ray imaging

In this paper, we model an x-ray imaging system, paying special attention to the energy- and depth-dependent characteristics of the inputs and interactions: x rays are polychromatic, interaction depth and conversion to optical photons is energy-dependent, optical scattering and the collection efficiency depend on the depth of interaction. The model we construct is a random function of the point process that begins with the distribution of x rays incident on the phosphor and ends with optical photons being detected by the active area of detector pixels to form an image. We show how the point-process representation can be used to calculate the characteristic statistics of the model. We then simulate a Gd2O2S:Tb phosphor, estimate its characteristic statistics, and proceed with a signal-detection experiment to investigate the impact of the pixel fill factor on detecting spherical calcifications (the signal). The two extremes possible from this experiment are that SNR2 does not change with fill factor or changes in proportion to fill factor. In our results, the impact of fill factor is between these extremes, and depends on the diameter of the signal.     

Signal and noise transfer properties of photoelectric interactions in diagnostic x-ray imaging detectors

Image quality in diagnostic x-ray imaging is ultimately limited by the statistical properties governing how, and where, x-ray energy is deposited in a detector. This in turn depends on the physics of the underlying x-ray interactions. In the diagnostic energy range (10-100 keV), most of the energy deposited in a detector is through photoelectric interactions. We present a theoretical model of the photoelectric effect that specifically addresses the statistical nature of energy absorption by photoelectrons, K and L characteristic x rays, and Auger electrons. A cascaded-systems approach is used that employs a complex structure of parallel cascades to describe signal and noise transfer through the photoelectric effect in terms of the modulation transfer function, Wiener noise power spectrum, and detective quantum efficiency (DQE). The model was evaluated by comparing results with Monte Carlo calculations for x-ray converters based on amorphous selenium (a-Se) and lead (Pb), representing both low and high-Z materials. When electron transport considerations can be neglected, excellent agreement (within 3%) is obtained for each metric over the entire diagnostic energy range in both a-Se and Pb detectors up to 30 cycles/mm, the highest frequency tested. The cascaded model overstates the DQE when the electron range cannot be ignored. This occurs at approximately two cycles/mm in a-Se at an incident photon energy of 80 keV, whereas in Pb, excellent agreement is obtained for the DQE over the entire diagnostic energy range. However, within the context of mammography (20 keV) and micro-computed tomography (40 keV), the effects of electron transport on the DQE are negligible compared to fluorescence reabsorption, which can lead to decreases of up to 30% and 20% in a-Se and Pb, respectively, at 20 keV; and 10% and 5%, respectively, at 40 keV. It is shown that when Swank noise is identified in a Fourier model, the Swank factor must be frequency dependent. This factor decreases quickly with frequency, and in the case of a-Se and Pb, decreases by up to a factor of 3 at five cycles/mm immediately above the K edge. The frequency-dependent Swank factor is also equivalent to what we call the "photoelectric DQE," which describes signal and noise transfer through photoelectric interactions.     

Absorption and noise in cesium iodide x-ray image intensifiers

The measured and theoretically predicted values of detective quantum efficiency (DQE) for a CsI x-ray image intensifier are compared for nine monoenergetic beams of x rays. The agreement between measurement and theory of better than +/- 5% indicates that we have a sound understanding of the physical parameters controlling the DQE. It is shown that the fraction of K-fluorescent x rays escaping from the input phosphor is independent of incident energy. The number of electrons released within the x-ray image intensifier (XRII) by an incident x ray has been measured. The mechanism for energy broadening within the XRII is shown to be predominantly the limited number of electrons and not light absorption.