A simple method for determining the modulation transfer function in digital radiography

PSF and MTF from a bar pattern in digital mammography

Background.The MTF has difficulties being determined (according to the provisions of the IEC standards) in the hospital setting due to the lack of resources.Purpose.The objective of this work is to propose a quantitative method for obtaining the point spread function (PSF) and the modulation transfer function (MTF) of a digital mammography system from an image of a bar pattern.Methods.The method is based on the measurement of the contrast transfer function (CTF) of the system over the image of the bar pattern. In addition, a theoretical model for thePSFis proposed, from which the theoreticalCTFof the system is obtained by means of convolution with a square wave (mathematical simulation of the bar pattern). Through an iterative process, the free parameters of thePSFmodel are varied until the experimentalCTFcoincides with the one calculated by convolution. Once thePSFof the system is obtained, we calculate theMTFby means of its Fourier transform. TheMTFcalculated from the modelPSFhave been compared with those calculated from an image of a 65μm diameter gold wire using an oversampling process.Results.TheCTFhas been calculated for three digital mammographic systems (DMS 1, DMS 2 and DMS 3), no differences of more than 5 % were found with the CTF obtained with the PSF model. The comparison of theMTFshows us the goodness of thePSFmodel.Conclusions.The proposed method for obtainingPSFandMTFis a simple and accessible method, which does not require a complex configuration or the use of phantoms that are difficult to access in the hospital world. In addition, it can be used to calculate other magnitudes of interest such as the normalized noise power spectrum (NNPS) and the detection quantum efficiency (DQE).

Development, validation, and simplification of a scanner-specific CT simulator

BACKGROUND

Simulated computed tomography (CT) images allow for knowledge of the underlying ground truth and for easy variation of imaging conditions, making them ideal for testing and optimization of new applications or algorithms. However, simulating all processes that affect CT images can result in simulations that are demanding in terms of processing time and computer memory. Therefore, it is of interest to determine how much the simulation can be simplified while still achieving realistic results.

PURPOSE

To develop a scanner-specific CT simulation using physics-based simulations for the position-dependent effects and shift-invariant image corruption methods for the detector effects. And to investigate the impact on image realism of introducing simplifications in the simulation process that lead to faster and less memory-demanding simulations.

METHODS

To make the simulator realistic and scanner-specific, the spatial resolution and noise characteristics, and the exposure-to-detector output relationship of a clinical CT system were determined. The simulator includes a finite focal spot size, raytracing of the digital phantom, gantry rotation during projection acquisition, and finite detector element size. Previously published spectral models were used to model the spectrum for the given tube voltage. The integrated energy at each element of the detector was calculated using the Beer-Lambert law. The resulting angular projections were subsequently corrupted by the detector modulation transfer function (MTF), and by addition of noise according to the noise power spectrum (NPS) and signal mean-variance relationship, which were measured for different scanner settings. The simulated sinograms were reconstructed on the clinical CT system and compared to real CT images in terms of CT numbers, noise magnitude using the standard deviation, noise frequency content using the NPS, and spatial resolution using the MTF throughout the field of view (FOV). The CT numbers were validated using a multi-energy CT phantom, the noise magnitude and frequency were validated with a water phantom, and the spatial resolution was validated with a tungsten wire. These metrics were compared at multiple scanner settings, and locations in the FOV. Once validated, the simulation was simplified by reducing the level of subsampling of the focal spot area, rotation and of detector pixel size, and the changes in MTFs were analyzed.

RESULTS

The average relative errors for spatial resolution within and across image slices, noise magnitude, and noise frequency content within and across slices were 3.4%, 3.3%, 4.9%, 3.9%, and 6.2%, respectively. The average absolute difference in CT numbers was 10.2 HU and the maximum was 22.5 HU. The simulation simplification showed that all subsampling can be avoided, except for angular, while the error in frequency at 10% MTF would be maximum 16.3%.

CONCLUSION

The simulation of a scanner-specific CT allows for the generation of realistic CT images by combining physics-based simulations for the position-dependent effects and image-corruption methods for the shift-invariant ones. Together with the available ground truth of the digital phantom, it results in a useful tool to perform quantitative analysis of reconstruction or post-processing algorithms. Some simulation simplifications allow for reduced time and computer power requirements with minimal loss of realism.

Integrated-mode proton radiography with 2D lateral projections

Integrated-mode proton radiography leading to water equivalent thickness (WET) maps is an avenue of interest for motion management, patient positioning, andin vivorange verification. Radiographs can be obtained using a pencil beam scanning setup with a large 3D monolithic scintillator coupled with optical cameras. Established reconstruction methods either (1) involve a camera at the distal end of the scintillator, or (2) use a lateral view camera as a range telescope. Both approaches lead to limited image quality. The purpose of this work is to propose a third, novel reconstruction framework that exploits the 2D information provided by two lateral view cameras, to improve image quality achievable using lateral views. The three methods are first compared in a simulated Geant4 Monte Carlo framework using an extended cardiac torso (XCAT) phantom and a slanted edge. The proposed method with 2D lateral views is also compared with the range telescope approach using experimental data acquired with a plastic volumetric scintillator. Scanned phantoms include a Las Vegas (contrast), 9 tissue-substitute inserts (WET accuracy), and a paediatric head phantom. Resolution increases from 0.24 (distal) to 0.33 lp mm-1(proposed method) on the simulated slanted edge phantom, and the mean absolute error on WET maps of the XCAT phantom is reduced from 3.4 to 2.7 mm with the same methods. Experimental data from the proposed 2D lateral views indicate a 36% increase in contrast relative to the range telescope method. High WET accuracy is obtained, with a mean absolute error of 0.4 mm over 9 inserts. Results are presented for various pencil beam spacing ranging from 2 to 6 mm. This work illustrates that high quality proton radiographs can be obtained with clinical beam settings and the proposed reconstruction framework with 2D lateral views, with potential applications in adaptive proton therapy.

Dynamic X-ray imaging with screen-printed perovskite CMOS array

High performance X-ray detector with ultra-high spatial and temporal resolution are crucial for biomedical imaging. This study reports a dynamic direct-conversion CMOS X-ray detector assembled with screen-printed CsPbBr3, whose mobility-lifetime product is 5.2 × 10-4 cm2 V-1 and X-ray sensitivity is 1.6 × 104 µC Gyair-1 cm-2. Samples larger than 5 cm[Formula: see text]10 cm can be rapidly imaged by scanning this detector at a speed of 300 frames per second along the vertical and horizontal directions. In comparison to traditional indirect-conversion CMOS X-ray detector, this perovskite CMOS detector offers high spatial resolution (5.0 lp mm-1) X-ray radiographic imaging capability at low radiation dose (260 nGy). Moreover, 3D tomographic images of a biological specimen are also successfully reconstructed. These results highlight the perovskite CMOS detector's potential in high-resolution, large-area, low-dose dynamic biomedical X-ray and CT imaging, as well as in non-destructive X-ray testing and security scanning.

Preparation and performance of a CsI scintillation screen with a double-period structure based on an oxidized silicon micropore array template

A structured double-period CsI scintillation screen was successfully developed to improve its detection efficiency based on an oxidized silicon micropore array template with a period value on the order of micro-scale. The structure comprises a main structure along with a sub-structure. The main structure with a period of 8 µm was arranged in a square array consisting of square columnar scintillator units. The micropore walls between the main structure units were purposely fabricated from a SiO2-Si-SiO2 layered structure. The pore walls in commonly used single-structure with a period of 4 µm use the same layered structure composition to obtain a fair comparison. The thickness of both Si and the SiO2 layers was around 0.4 µm. The unique feature of the double structure lies in the even separation of each unit within the main structure into four square columnar scintillator sub-units. These four sub-units within each sub-structure were isolated solely by SiO2 layers with a thickness of approximately 0.8 µm. As a result, the X-ray-induced optical luminescence intensity of the double-structure screen exhibited a 31% increase compared to the corresponding single-structure scintillation screen. In X-ray imaging, a spatial resolution of 109 lp/mm was achieved, which closely matched the results obtained with the single-structure CsI screen. Furthermore, the detective quantum efficiency also displayed a notable improvement.

High-Efficiency Narrow-Band Green-Emitting Manganese(II) Halide for Multifunctional Applications

Zero-dimensional (0D) Mn2+-based metal halides used as luminescent materials and scintillators have become a research hotspot in the field of photoelectric materials and devices due to their unique composition, structure, and fluorescence properties. It is of great value to explore new Mn2+-based metal halides to achieve multifunctional applications. Herein, the novel 0D Mn2+-based metal halide single crystal (BPTP)2MnBr4 is synthesized by a simple solvent-antisolvent recrystallization method. Under excitation at 468 nm, the (BPTP)2MnBr4 single crystal shows a pronounced narrow-band green luminescence centered at 515 nm derived from the d-d transition of the Mn2+ ion. This emission has a relatively narrow full width at half maximum of 43 nm and a high photoluminescence quantum yield (PLQY) of 82%. In addition, (BPTP)2MnBr4 exhibits good thermal stability at 393 K with a retention of 79% of the initial photoluminescence intensity at 298 K. Benefiting from its strong blue light excitation, high PLQY, and good thermal stability, we manufacture an ideal white light-emitting diode (LED) device using a 460 nm blue LED chip, green-emitting (BPTP)2MnBr4, and commercial K2SiF6:Mn4+ red phosphor. Under 20 mA drive current, the LED shows a high luminous efficiency of 112 lm/W and a wide color gamut of 110.8%, according to the National Television System Committee standard. In addition, (BPTP)2MnBr4 crystals show a strong X-ray absorption. Based on the commercial Lu3Al5O12:Ce3+ scintillator, the calculated light yield of (BPTP)2MnBr4 reaches up to about 136,000 photons/MeV and the detection limit reaches 0.282 μGyair s-1. Additionally, a melt quenching approach is used to construct a (BPTP)2MnBr4 clear glass scintillation screen, realizing a spatial resolution of 10.1 lp/mm. The proper performances of (BPTP)2MnBr4 as phosphor-converted LED materials and the X-ray scintillator with the addition of eco-friendly, low-cost solution processability make 0D Mn2+-based metal halides potential luminescent materials for multifunctional applications.

A simulation study to evaluate the effect of 2D antiscatter grid primary transmission on flat panel detector based CBCT image quality

Purpose. Two-dimensional antiscatter grids' (2D-ASGs) septal shadows and their impact on primary transmission play a critical role in cone-beam computed tomography (CBCT) image noise and artifact characteristics. Therefore, a numerical simulation platform was developed to evaluate the effect of 2D-ASG's primary transmission on image quality, as a function of grid geometry and CBCT system properties.Methods. To study the effect of 2D-ASG's septal shadows on primary transmission and CBCT image quality, two new methods were introduced; one to simulate projection signal gradients in septal shadows, and the other to simulate septal shadow variations due to gantry flex. Signal gradients in septal shadows were simulated by generating a system point spread function that was directly extracted from projection images of 2D-ASG prototypes in experiments. Variations in septal shadows due to gantry flex were simulated by generating oversampled shadow profiles extracted from experiments. Subsequently, the effect of 2D-ASG's septal shadows on primary transmission and image quality was evaluated.Results.For an apparent septal thickness of 0.15 mm, the primary transmission of 2D-ASG varied between 72%-90% for grid pitches 1-3 mm. In low-contrast phantoms, the effect of 2D-ASG's radiopaque footprint on information loss was subtle. At high spatial frequencies, information loss manifested itself as undersampling artifacts, however, its impact on image quality is subtle when compared to quantum noise. Effects of additive electronic noise and gantry flex induced ring artifacts on image quality varied as a function of grid pitch and septal thickness. Such artifacts were substantially less in lower resolution images.Conclusion. The proposed simulation platform allowed successful evaluation of CBCT image quality variations as a function of 2D-ASG primary transmission properties and CBCT system characteristics. This platform can be potentially used for optimizing 2D-ASG design properties based on the imaging task and properties of the CBCT system.

Sparse2Noise: Low-dose synchrotron X-ray tomography without high-quality reference data

BACKGROUND

Synchrotron radiation computed tomography (SR-CT) holds promise for high-resolution in vivo imaging. Notably, the reconstruction of SR-CT images necessitates a large set of data to be captured with sufficient photons from multiple angles, resulting in high radiation dose received by the object. Reducing the number of projections and/or photon flux is a straightforward means to lessen the radiation dose, however, compromises data completeness, thus introducing noises and artifacts. Deep learning (DL)-based supervised methods effectively denoise and remove artifacts, but they heavily depend on high-quality paired data acquired at high doses. Although algorithms exist for training without high-quality references, they struggle to effectively eliminate persistent artifacts present in real-world data.

METHODS

This work presents a novel low-dose imaging strategy namely Sparse2Noise, which combines the reconstruction data from paired sparse-view CT scan (normal-flux) and full-view CT scan (low-flux) using a convolutional neural network (CNN). Sparse2Noise does not require high-quality reconstructed data as references and allows for fresh training on data with very small size. Sparse2Noise was evaluated by both simulated and experimental data.

RESULTS

Sparse2Noise effectively reduces noise and ring artifacts while maintaining high image quality, outperforming state-of-the-art image denoising methods at same dose levels. Furthermore, Sparse2Noise produces impressive high image quality for ex vivo rat hindlimb imaging with the acceptable low radiation dose (i.e., 0.5 Gy with the isotropic voxel size of 26 μm).

CONCLUSIONS

This work represents a significant advance towards in vivo SR-CT imaging. It is noteworthy that Sparse2Noise can also be used for denoising in conventional CT and/or phase-contrast CT.

Seven general radiography x-ray detectors with pixel sizes ranging from 175 to 76μm: technical evaluation with the focus on orthopaedic imaging

Aim. Flat panel detectors with small pixel sizes general can potentially improve imaging performance in radiography applications requiring fine detail resolution. This study evaluated the imaging performance of seven detectors, covering a wide range of pixel sizes, in the frame of orthopaedic applications.Material and methods. Pixel sizes ranged from 175 (detector A175) to 76μm (detector G76). Modulation transfer function (MTF) and detective quantum efficiency (DQE) were measured using International Electrotechnical Commission (IEC) RQA3 beam quality. Threshold contrast (CT) and a detectability index (d') were measured at three air kerma/image levels. Rabbit shoulder images acquired at 60 kV, over five air kerma levels, were evaluated in a visual grading study for anatomical sharpness, image noise and overall diagnostic image quality by four radiologists. The detectors were compared to detector E124.Results. The 10% point of the MTF ranged from 3.21 to 4.80 mm-1, in going from detector A175to detector G76. DQE(0.5 mm-1) measured at 2.38μGy/image was 0.50 ± 0.05 for six detectors, but was higher for F100at 0.62. High frequency DQE was superior for the smaller pixel detectors, howeverCTfor 0.25 mm discs correlated best with DQE(0.5 mm-1). Correlation betweenCTand the detectability model was good (R2= 0.964).CTfor 0.25 mm diameter discs was significantly higher for D150and F100compared to E124. The visual grading data revealed higher image quality ratings for detectors D125and F100compared to E124. An increase in air kerma was associated with improved perceived sharpness and overall quality score, independent of detector. Detectors B150, D125, F100and G76, performed well in specific tests, however only F100consistently outperformed the reference detector.Conclusion. Pixel size alone was not a reliable predictor of small detail detectability or even perceived sharpness in a visual grading analysis study.

Optimization of SiO2 reflective layer thickness for improving the performance of structured CsI scintillation screen based on oxidized Si micropore array template in X-ray imaging

Structured scintillation screen based on oxidized Si micropore array template can effectively improve the spatial resolution of X-ray imaging. The purpose of this study is to investigate the effect of SiO₂ layer thickness on the light guide and X-ray imaging performance of CsI scintillation screen when the structural period is as small as microns. Cylindrical micropores with a period of 4.3 µm, an average diameter of 3.3 µm and a depth of about 40 µm were prepared in Si wafers. SiO₂ layer was formed on the pore walls after thermal oxidation. Increasing SiO₂ layer thickness would be beneficial to the propagation of scintillation light along the cylindrical channels. What was not previously anticipated was that the pore size gradually shrank as the SiO₂ layer thickened. The pore shrinkage would reduce the filling rate of CsI in the templates and thus would reduce the production of scintillation light. The structured CsI scintillation screens with different SiO₂ layer thicknesses were fabricated by filling CsI scintillator into the oxidized silicon micropore array template. The morphology, crystallinity, X-ray excited optical luminescence, and X-ray imaging performance of the screens were studied. The results show that the spatial resolutions of X-ray images measured using the structured CsI scintillation screens with different SiO₂ layer thicknesses are close to each other, and they are all about 110 lp/mm. However, the X-ray excited optical luminescence of the screen and detective quantum efficiency of X-ray imaging vary with the thickness of the SiO₂ layer. The optimal thickness is about 350 nm.

Deep learning-based computed tomographic image super-resolution via wavelet embedding

Effort to realize high-resolution medical images have been made steadily. In particular, super resolution technology based on deep learning is making excellent achievement in computer vision recently. In this study, we developed a model that can dramatically increase the spatial resolution of medical images using deep learning technology, and we try to demonstrate the superiority of proposed model by analyzing it quantitatively. We simulated the computed tomography images with various detector pixel size and tried to restore the low-resolution image to high resolution image. We set the pixel size to 0.5, 0.8 and 1 mm2 for low resolution image and the high-resolution image, which were used for ground truth, was simulated with 0.25 mm2 pixel size. The deep learning model that we used was a fully convolution neural network based on residual structure. The result image demonstrated that proposed super resolution convolution neural network improve image resolution significantly. We also confirmed that PSNR and MTF was improved up to 38 % and 65% respectively. The quality of the prediction image is not significantly different depending on the quality of the input image. In addition, the proposed technique not only increases image resolution but also has some effect on noise reduction. In conclusion, we developed deep learning architectures for improving image resolution of computed tomography images. We quantitatively confirmed that the proposed technique effectively improves image resolution without distorting the anatomical structures.

Low-density tissue scaffold imaging by synchrotron radiation propagation-based imaging computed tomography with helical acquisition mode

Visualization of low-density tissue scaffolds made from hydrogels is important yet challenging in tissue engineering and regenerative medicine (TERM). For this, synchrotron radiation propagation-based imaging computed tomography (SR-PBI-CT) has great potential, but is limited due to the ring artifacts commonly observed in SR-PBI-CT images. To address this issue, this study focuses on the integration of SR-PBI-CT and helical acquisition mode (i.e. SR-PBI-HCT) to visualize hydrogel scaffolds. The influence of key imaging parameters on the image quality of hydrogel scaffolds was investigated, including the helical pitch (p), photon energy (E) and the number of acquisition projections per rotation/revolution (Np), and, on this basis, those parameters were optimized to improve image quality and to reduce noise level and artifacts. The results illustrate that SR-PBI-HCT imaging shows impressive advantages in avoiding ring artifacts with p = 1.5, E = 30 keV and Np = 500 for the visualization of hydrogel scaffolds in vitro. Furthermore, the results also demonstrate that hydrogel scaffolds can be visualized using SR-PBI-HCT with good contrast while at a low radiation dose, i.e. 342 mGy (voxel size of 26 µm, suitable for in vivo imaging). This paper presents a systematic study on hydrogel scaffold imaging using SR-PBI-HCT and the results reveal that SR-PBI-HCT is a powerful tool for visualizing and characterizing low-density scaffolds with a high image quality in vitro. This work represents a significant advance toward the non-invasive in vivo visualization and characterization of hydrogel scaffolds at a suitable radiation dose.

Full modulation transfer functions of thick parallel- and focused-element scintillator arrays obtained by a Monte Carlo optical transport model

BACKGROUND

Arrays of thick segmented crystalline scintillators are useful x-ray converters for image-guided radiation therapy using electronic portal imaging (EPI) and megavoltage cone-beam computed tomography (MV-CBCT). Ionizing-radiation-only simulations previously showed relatively low modulation transfer function (MTF) in parallel-element arrays because of beam divergence. Hence, a focused-element geometry (matching the beam divergence) has been proposed. The "full" (ionizing and optical) MTF performance of such a focused geometry compared to its radiation-only MTF has, however, not been fully investigated.

PURPOSE

To study the full MTF performance of such arrays in a more realistic situation in which optical characteristics are also included using an in-house detector model that supports light transport, and quantify the errors in MTF estimation when the optical stage is ignored.

METHODS

First, radiation (x-ray and electron) transport was simulated. Then, transport of the generated optical photons was modeled using ScintSim2, an optical Monte Carlo (MC) code developed in MATLAB for simulation of two-dimensional (2D) parallel- and focused-element scintillator arrays. The full-MTF responses of focused- and parallel-element geometries, for a large array of 3 × 3 mm² CsI:Tl detector elements of 10, 40, and 60 mm thicknesses, were examined. For each configuration, a composite line spread function (LSF) was calculated to obtain the MTF.

RESULTS

At the Nyquist frequency, for 10 mm-thick central elements and 60 mm-thick peripheral parallel elements, full-MTF exhibited a drop of up to 15 and 79 times, respectively, compared with radiation-only MTF. This was found to be partly attributable to the angular distribution of the light emerging from the detector-element exit face and the dependence on its aspect ratio, since the light exiting thicker scintillators exhibited a more forward-directed distribution. Focused elements provided an increase of up to nine times in peripheral-area full MTF values.

CONCLUSIONS

Full MTF was up to 79 times lower than radiation-only MTF. Focused arrays preserved full MTF by up to nine times compared to parallel elements. The differences in the results obtained with and without inclusion of optical photons emphasize the need to include light transport when optimizing thick segmented scintillation detectors. Besides their application in detector optimization for radiotherapy megavoltage photon imaging, these findings can also be useful for other segmented-scintillator-based imaging systems, for example, in nuclear medicine, or in 2D detection systems for quality assurance of MR-linacs.

Simulation Analysis on the Performance of a Circular-Edge Technique in Measurements of the Modulation Transfer Function

The modulation transfer function (MTF) plays an important role in characterizing medical imaging systems. For such characterization, the circular-edge technique has become a prevalent task-based methodology. When determining the MTF with complicated task-based measurements, error factors must be well understood to properly interpret the results. In this context, the aim of this work was to study the changes in measurement performance in the analysis of the MTF using a circular edge. To eliminate the systematic error related to the measurement and suitably manage the error factors, images were generated by Monte Carlo simulation. Further, a performance comparison with the conventional method was conducted; in addition, the influence of the edge size and contrast and the setting error of the center coordinates were investigated. The difference from the true value and the standard deviation relative to the average value were applied to the index as the accuracy and precision, respectively. The results demonstrated that the smaller the circular object used and the lower the contrast, the grater the deterioration in the measurement performance. Furthermore, this study clarified the underestimating of the MTF in proportion to the square of the distance with respect to the setting error of the center position, which is important for the synthesis of the edge profile. Evaluations in the backgrounds wherein there are multiple factors affecting the results are complex, and the system users must properly judge the validity of the characterization results. These findings provide meaningful insight in the context of MTF measurement techniques.

Compact and versatile neutron imaging detector with sub-4μm spatial resolution based on a single-crystal thin-film scintillator

A large and increasing number of scientific domains pushes for high neutron imaging resolution achieved in reasonable times. Here we present the principle, design and performance of a detector based on infinity corrected optics combined with a crystalline Gd₃Ga₅O₁₂ : Eu scintillator, which provides an isotropic sub-4 µm true resolution. The exposure times are only of a few minutes per image. This is made possible also by the uniquely intense cold neutron flux available at the imaging beamline NeXT-Grenoble. These comparatively rapid acquisitions are compatible with multiple high quality tomographic acquisitions, opening new venues for in-operando testing, as briefly exemplified here.

Review of high energy x-ray computed tomography for non-destructive dimensional metrology of large metallic advanced manufactured components

Advanced manufacturing technologies, led by additive manufacturing, have undergone significant growth in recent years. These technologies enable engineers to design parts with reduced weight while maintaining structural and functional integrity. In particular, metal additive manufacturing parts are increasingly used in application areas such as aerospace, where a failure of a mission-critical part can have dire safety consequences. Therefore, the quality of these components is extremely important. A critical aspect of quality control is dimensional evaluation, where measurements provide quantitative results that are traceable to the standard unit of length, the metre. Dimensional measurements allow designers, manufacturers and users to check product conformity against engineering drawings and enable the same quality standard to be used across the supply chain nationally and internationally. However, there is a lack of development of measurement techniques that provide non-destructive dimensional measurements beyond common non-destructive evaluation focused on defect detection. X-ray computed tomography (XCT) technology has great potential to be used as a non-destructive dimensional evaluation technology. However, technology development is behind the demand and growth for advanced manufactured parts. Both the size and the value of advanced manufactured parts have grown significantly in recent years, leading to new requirements of dimensional measurement technologies. This paper is a cross-disciplinary review of state-of-the-art non-destructive dimensional measuring techniques relevant to advanced manufacturing of metallic parts at larger length scales, especially the use of high energy XCT with source energy of greater than 400 kV to address the need in measuring large advanced manufactured parts. Technologies considered as potential high energy x-ray generators include both conventional x-ray tubes, linear accelerators, and alternative technologies such as inverse Compton scattering sources, synchrotron sources and laser-driven plasma sources. Their technology advances and challenges are elaborated on. The paper also outlines the development of XCT for dimensional metrology and future needs.

Neutron flat-panel detector using In-Ga-Zn-O thin-film transistor

Neutron imaging is a powerful tool for observing the internal structure of an object without destroying the object. Neutron imaging (neutron radiography) is a prominent application of neutrons but still requires significant improvements, for example, in sensitivity, resolution, radiation hardness, and handling of neutron imaging detectors. This paper presents the development and the first neutron imaging results of a neutron flat-panel detector (nFPD) based on an In-Ga-Zn-O (IGZO) thin-film transistor (TFT)/photodiode array coupled with a LiF/ZnS scintillator sheet. Direct photo-coupling to the scintillator increases the light collection efficiency. Moreover, unlike lens-coupled neutron cameras, the proposed detector is compact and easy to handle. Owing to the high off-state resistance of IGZO TFTs, their leakage current is lower than that of conventional TFTs, enabling the IGZO TFTs to hold an accumulated charge for a longer period of time and allowing longer exposure times for imaging. This would be a powerful feature for imaging at compact neutron sources with limited flux. This paper reports on the first neutron imaging results with an IGZO nFPD, its performance evaluation, and a demonstration of three-dimensional computed tomography with neutrons.

Assessment of Uncertainty Depending on Various Conditions in Modulation Transfer Function Calculation Using the Edge Method

In medical X-ray imaging, to perform optimal operations, it is essential for the user to understand whether a required image quality level which depends on a diagnostic task can be achieved with the imaging system used. This study focuses on the effects of noise on the modulation transfer function (MTF) using the edge method, the most widely used to evaluate the task dependence property. The purpose is to verify the uncertainty of the MTF value at each spatial frequency and examine the conditions under which the accuracy is ensured. By using a Monte Carlo simulation, edge images with various contrast-to-noise ratio (CNR) are acquired. MTFs are then calculated with different edge spread function (ESF) lengths. The uncertainties for each spatial frequency are estimated based on the independent MTF calculations obtained from the five edge data. The uncertainty of the MTF is inversely proportional to the CNR. In the frequency range up to the Nyquist frequency, the uncertainty in five calculations is <0.01 when the CNR is more than 60. In addition, it is observed that the uncertainty increases as the ESF length increases. This relationship depends on the frequency range, but it is proportional to the 0.3–0.5 power of the ESF length. The results in which the uncertainty is most likely to be large in the MTF calculation are clearly shown. Therefore, it is expected to provide an important barometer and useful insights for a proper image quality measurement.

Determination of DQE as a quantitative assessment of detectors in digital mammography: Measurements and calculation in practice

Introduction: Advances in digital detector technology and methods of image presentation in digital mammography now offer the possibility of implementing mathematical assessment methods to quantitative image analysis. The aim of this work was to develop new software to simplify the application of the existing international standard for DQE in digital mammography and show in detail how it can be applied, using a Siemens Mammomat Inspiration as a model.

Material and methods: Consistent with the IEC standard a 2 mm Al filter at the tube exit and images in DICOM format as raw data, without applying any additional post-processing were used. Measurements were performed for W/Rh anode/filter combination and different tube voltage values (26 ÷ 34 kV) without any anti-scatter grid. To verify new software doses ranging from 20-600 µGy were used in measurements. Exposure (air kerma) was measured using a calibrated radiation meter (Piranha Black 457, RTI Electronics AB, Sweden). MTF was determined, using an edge test device constructed specifically for this work.

Results: It has been demonstrated that with the new software the DQE can be measured with the accuracy required by the international standard IEC 62220-1-2. DQE has been presented as a function of spatial frequency for W/Rh anode/filter combination and different tube voltage.

Conclusions: New software was used successfully to analyze image quality parameters for the Siemens Mammomat Inspiration detector. This was done on the basis of an internationally accepted methodology. In the next step, mammographs with different detector types can be compared.

Effect of beam quality and readout direction in the edge profile on the modulation transfer function of photostimulable phosphor systems via the edge method

Purpose: The objective of this study was to investigate the effects of (i) the difference in the beam qualities on the presampled modulation transfer function (MTF) using the edge method and (ii) the readout direction of the edge profile for the photostimulable phosphor (PSP) system. Approach: The International Electrotechnical Commission (IEC) defined a technique using the "radiation qualities based on a phantom made up of an aluminum added filter" (RQA). A general radiographic system with a tube voltage of 50 kV and 9.7 mm of additional aluminum filtration was used to conform the x-ray to the IEC-specified beam quality definition RQA3. Additionally, we employed two different beams with tube voltages of 60 and 70 kV using a dental x-ray unit. The MTF was measured in the readout direction from low-to-high exposure regions and vice versa with respect to the scanning and subscanning directions of the PSP system. Results: The difference in the 50%MTF value for all directions between 60 and 70 kV averaged less than 0.05 . The 50%MTF of RQA3 was on average 0.2 lower than the value for 60 kV for each direction. For all beam qualities, no difference was observed between the MTFs measured in the readout direction from the low-to-high exposure regions and vice versa. Conclusions: The MTFs, measured using the dental x-ray unit, were unaffected by the tube voltage, and they were slightly higher than those measured using the RQA3. Furthermore, the MTF was unaffected by the differences in the readout directions of the edge profile.

Abbreviated on-treatment CBCT using roughness penalized mono-energization of kV-MV data and a multi-layer MV imager

Simultaneous acquisition of cone beam CT (CBCT) projections using both the kV and MV imagers of an image guided radiotherapy system reduces set-up scan times-a benefit to lung cancer radiation oncology patients-but increases noise in the 3D reconstruction. In this article, we present a kV-MV scan time reduction technique that uses two noise-reducing measures to achieve superior performance. The first is a high-DQE multi-layer MV imager prototype. The second is a beam hardening correction algorithm which combines poly-energetic modeling with edge-preserving, regularized smoothing of the projections. Performance was tested in real acquisitions of the Catphan 604 and a thorax phantom. Percent noise was quantified from voxel values in a soft tissue volume of interest (VOI) while edge blur was quantified from a VOI straddling a boundary between air and soft material. Comparisons in noise/resolution performance trade-off were made between our proposed approach, a dose-equivalent kV-only scan, and a kV-MV reconstruction technique previously published by Yinet al(2005Med. Phys.329). The proposed technique demonstrated lower noise as a function of spatial resolution than the baseline kV-MV method, notably a 50% noise reduction at typical edge blur levels. Our proposed method also exhibited fainter non-uniformity artifacts and in some cases superior contrast. Overall, we find that the combination of a multi-layer MV imager, acquiring at a LINAC source energy of 2.5 MV, and a denoised beam hardening correction algorithm enables noise, resolution, and dose performance comparable to standard kV-imager only set-up CBCT, but with nearly half the gantry rotation time.

Frequency-dependent optimal weighting approach for megavoltage multilayer imagers

Multi-layer imaging (MLI) devices improve the detective quantum efficiency (DQE) while maintaining the spatial resolution of conventional mega-voltage (MV) x-ray detectors for applications in radiotherapy. To date, only MLIs with identical detector layers have been explored. However, it may be possible to instead use different scintillation materials in each layer to improve the final image quality. To this end, we developed and validated a method for optimally combining the individual images from each layer of MLI devices that are built with heterogeneous layers. Two configurations were modeled within the GATE Monte Carlo package by stacking different layers of a terbium doped gadolinium oxysulfide Gd2O2S:Tb (GOS) phosphor and a LKH-5 glass scintillator. Detector response was characterized in terms of the modulation transfer function (MTF), normalized noise power spectrum (NNPS) and DQE. Spatial frequency-dependent weighting factors were then analytically derived for each layer such that the total DQE of the summed combination image would be maximized across all spatial modes. The final image is obtained as the weighted sum of the sub-images from each layer. Optimal weighting factors that maximize the DQE were found to be the quotient of MTF and NNPS of each layer in the heterogeneous MLI detector. Results validated the improvement of the DQE across the entire frequency domain. For the LKH-5 slab configuration, DQE(0) increases between 2%-3% (absolute), while the corresponding improvement for the LKH-5 pixelated configuration was 7%. The performance of the weighting method was quantitatively evaluated with respect to spatial resolution, contrast-to-noise ratio (CNR) and signal-to-noise ratio (SNR) of simulated planar images of phantoms at 2.5 and 6 MV. The line pair phantom acquisition exhibited a twofold increase in CNR and SNR, however MTF was degraded at spatial frequencies greater than 0.2 lp mm^-1. For the Las Vegas phantom, the weighting improved the CNR by around 30% depending on the contrast region while the SNR values are higher by a factor of 2.5. These results indicate that the imaging performance of MLI systems can be enhanced using the proposed frequency-dependent weighting scheme. The CNR and SNR of the weighted combined image are improved across all spatial scales independent of the detector combination or photon beam energy.

Theory, method, and test tools for determination of 3D MTF characteristics in cone-beam CT

PURPOSE

The modulation transfer function (MTF) is widely used as an objective metric of spatial resolution of medical imaging systems. Despite advances in capability for three-dimensional (3D) isotropic spatial resolution in computed tomography (CT) and cone-beam CT (CBCT), MTF evaluation for such systems is typically reported only in the axial plane, and practical methodology for assessment of fully 3D spatial resolution characteristics is lacking. This work reviews fundamental theoretical relationships of two-dimensional (2D) and 3D spread functions and reports practical methods and test tools for analysis of 3D MTF in CBCT.

METHODS

Fundamental aspects of 2D and 3D MTF measurement are reviewed within a common notational framework, and three MTF test tools with analysis code are reported and made available online (https://istar.jhu.edu/downloads/): (a) a multi-wire tool for measurement of the axial plane MTF [denoted as M T F ( f r ; φ = 0 ∘ ) , where φ is the measurement angle out of the axial plane] as a function of position in the axial plane; (b) a wedge tool for measurement of the MTF in any direction in the 3D Fourier domain [e.g.,  φ  = 45°, denoted as M T F ( f r ; φ = 45 ∘ ) ]; and (c) a sphere tool for measurement of the MTF in any or all directions in the 3D Fourier domain. Experiments were performed on a mobile C-arm with CBCT capability, showing that M T F ( f r ; φ = 45 ∘ ) yields an informative one-dimensional (1D) representation of the overall 3D spatial resolution characteristics, capturing important characteristics of the 3D MTF that might be missed in conventional analysis. The effects of anisotropic filters and detector readout mode were investigated, and the extent to which a system can be said to provide "isotropic" resolution was evaluated by quantitative comparison of MTF at various φ .

RESULTS

All three test tools provided consistent measurement of M T F ( f r ; φ = 0 ∘ ) , and the wedge and sphere tools demonstrated how measurement of the MTF in directions outside the axial plane ( φ > 0 ∘ ) can reveal spatial resolution characteristics to which conventional axial MTF measurement is blind. The wedge tool was shown to reduce statistical measurement error compared to the sphere tool due to improved sampling, and the sphere tool was shown to provide a basis for measurement of the MTF in any or all directions (outside the null cone) from a single scan. The C-arm system exhibited non-isotropic spatial resolution with conventional non-isotropic 1D apodization filters (i.e., frequency cutoff filters) - which is common in CBCT - and implementation of isotropic 2D apodization yielded quantifiably isotropic MTF. Asymmetric pixel binning modes were similarly shown to impart non-isotropic effects on the 3D MTF, and the overall 3D MTF characteristics were evident in each case with a single, 1D measurement of the 1D M T F ( f r ; φ = 45 ∘ ).

CONCLUSION

Three test tools and corresponding MTF analysis methods were presented within a consistent framework for analysis of 3D spatial resolution characteristics in a manner amenable to routine, practical measurements. Experiments on a CBCT C-arm validated many intuitive aspects of 3D spatial resolution and quantified the extent to which a CBCT system may be considered to have isotropic resolution. Measurement of M T F ( f r ; φ = 45 ∘ ) provided a practical 1D measure of the underlying 3D MTF characteristics and is extensible to other CT or CBCT systems offering high, isotropic spatial resolution.

Anomalous edge response of cadmium telluride-based photon counting detectors jointly caused by high-flux radiation and inter-pixel communication

This work reports an edge enhancing effect experimentally observed in cadmium telluride (CdTe)-based photon counting detector (PCD) systems operated under the charge summing (CS) mode and irradiated by high-flux x-rays. Experimental measurements of the edge spread functions (ESFs) of a PCD system (100μm pixel size, 88 ns deadtime) were performed at different input flux levels from 4.5 × 105count per second (cps) mm-2to 1.5 × 109cps mm-2for the single pixel mode (SP) and the CS mode. A theoretical model that incorporates the impacts of inter-pixel communications and the arbitration process involved in the CS mode was developed to help explain the physical origin of the observed edge enhancing effect. Compared with the monotonically increasing ESF of the SP mode, the ESF of the CS mode measured at high-flux levels shows a peak at an intermediate location (50μm from the edge). The peak became more pronounced with increasing flux levels. The theoretically calculated ESFs agreed well with experimental results with relative errors less than 5% at all flux levels and tested. These results indicate that the anomalous edge enhancing effect is jointly caused by the pileup effect and the CS circuit that introduces negative correlations between adjacent pixels. When the input flux is high enough to deliver photons to multiple adjacent pixels within the same deadtime period, the CS mode may treat the coincident x-rays as shared charges and thus introduce count losses in addition to the well-known pileup count loss. When a high contrast object partially blocks certain pixels from x-rays, the adjacent unblocked pixels have an increased probability of registering counts as a result of the negative correlation. This leads to a peak on the ESF at a pixel-to-edge distance half of the pixel pitch.

Evaluation of cone-beam computed tomography over a small field of view in a water bath based on the modulation transfer function with repeating-edge oversampling

PURPOSE

The spatial resolution of cone-beam computed tomography (CBCT) in small fields of view (FOVs) is important for clinical applications. However, it is difficult to measure spatial resolution reliably due to error factors such as noise. The aim of this study was to obtain a modulation transfer function (MTF) more accurately.

METHODS

A CBCT apparatus was used with small FOV. An aluminum pipe slightly tilted at an inclination ratio of 77/3 (25.7) was used as the measurement phantom. The MTF was calculated from the edge image of the phantom. The actual oversampling ratio was determined by regression analysis. The experiment was repeated 16 times and the edge-spread function (ESF) was approximated by the least-square method. Furthermore, a low-pass filter (LPF) was applied to eliminate the component at frequencies above the Nyquist frequency. Finally, the MTF was calculated from the pre-processed ESF.

RESULTS

Results showed that pre-processing reduced the noise of the ESF. The MTFs at frequencies of 1.0 and 2.0 LP/mm were 0.59 and 0.18, respectively, in air and 0.52 and 0.16, respectively, in water.

CONCLUSION

The repeating-edge oversampling method combined with ESF pre-processing improved the accuracy of the MTF under clinically relevant conditions with a phantom.

Clinical translation of a new flat-panel detector for beam's-eye-view imaging

Electronic portal imaging devices (EPIDs) lend themselves to beams-eye view clinical applications, such as tumor tracking, but are limited by low contrast and detective quantum efficiency (DQE). We characterize a novel EPID prototype consisting of multiple layers and investigate its suitability for use under clinical conditions. A prototype multi-layer imager (MLI) was constructed utilizing four conventional EPID layers, each consisting of a copper plate, a Gd2O2S:Tb phosphor scintillator, and an amorphous silicon flat panel array detector. We measured the detector's response to a 6 MV photon beam with regards to modulation transfer function, noise power spectrum, DQE, contrast-to-noise ratio (CNR), signal-to-noise ratio (SNR), and the linearity of the detector's response to dose. Additionally, we compared MLI performance to the single top layer of the MLI and the standard Varian AS-1200 detector. Pre-clinical imaging was done on an anthropomorphic phantom, and the detector's CNR, SNR and spatial resolution were assessed in a clinical environment. Images obtained from spine and liver patient treatment deliveries were analyzed to verify CNR and SNR improvements. The MLI has a DQE(0) of 9.7%, about 5.7 times the reference AS-1200 detector. Improved noise performance largely drives the increase. CNR and SNR of clinical images improved three-fold compared to reference. A novel MLI was characterized and prepared for clinical translation. The MLI substantially improved DQE and CNR performance while maintaining the same resolution. Pre-clinical tests on an anthropomorphic phantom demonstrated improved performance as predicted theoretically. Preliminary patient data were analyzed, confirming improved CNR and SNR. Clinical applications are anticipated to include more accurate soft tissue tracking.

Visualization of blurring process due to analog components in a digital radiography system using a simple method

To obtain useful medical X-ray images for diagnosis, it is necessary to accurately evaluate the characteristics and fully understand the features of the system to determine the imaging target and conditions. This study describes the performance of a digital radiography (DR) system in detail and proposes a method for characterizing the blurring process due to the presence of analog components in an imaging system. Our method does not involve any specialized technique, such as a simulation that requires long computing time on a high-performance computer. The method is based on the concept of the modulation transfer function (MTF). The functions that corresponded to the MTF of the analog components of the system when the Fourier transform was performed were examined. Indirect conversion type flat-panel detectors (FPD) in the general radiography energy range and direct conversion type FPD in the mammography were used for verification. A Gaussian function and Lorentz function were synthesized for modelling the blurring process of the indirect type FPDs. The spreads and shapes of the modeled functions depended on the systems. For the direct conversion type FPD, a combination of a Lorentz function with a narrow width and an impulse function characterized the results. Furthermore, based on the structure of each detector and the signal transfer process, the visualized results were considered reliable. The results of these experiments will aid in the detailed understanding of the image quality characteristics of the DR system by the feasible approach.

Single-plate method for practical modulation transfer function measurement in magnetic resonance imaging

A new "single-plate method" is presented for measuring the modulation transfer function (MTF) in magnetic resonance imaging (MRI). This method sets a slice plane perpendicular to a single-plate phantom to eliminate contamination effects from the direction perpendicular to the measurement direction in the image plane, which occur with a conventional ramp method. As no practical method for measuring the MTF has been established for MRI, we examined whether the MTF can be measured practically using the single-plate method for various fast imaging sequences. Furthermore, the MTFs of T₁-weighted (T₁W) fast spin echo (FSE) and conventional spin echo (CSE) images obtained using the single-plate method and ramp method were compared. The measured MTFs of T₁W CSE images revealed rectangular shapes with a sharp decrease near the Nyquist cutoff frequency in both phase-encoding (PE) and frequency-encoding (readout, RO) directions. The measured MTFs of T₁W FSE images obtained with centric-order acquisition showed symmetric step-function shapes reflecting k-space segmentation determined by the echo train length (ETL). The measured MTFs of T₂-weighted (T₂W) FSE images showed asymmetric step-function shapes reflecting differences in T₂ decay of signals from samples. The MTFs obtained using the single-plate method significantly reduced the collapse caused by the contamination effect, which is observed in all the MTF measurements of the ramp method. The proposed "single-plate method" simplified the complicated MTF measurement procedure and eliminated the contamination effect. This method is expected to be useful for evaluating the resolution properties of MR fast imaging techniques with a complicated k-space trajectory.

Performance of a new star-bar phantom designed for MTF calculations in x-ray imaging systems

PURPOSE

A new phantom, designed and manufactured for modulation transfer function (MTF) calculations is presented in this work. The phantom has a star-bar pattern and is manufactured in stainless steel. Modulation transfer function determinations are carried out with the new phantom and with an edge phantom to compare their performance and to compare them with previous theoretical predictions.

METHODS

The phantoms are imaged in an x-ray imaging system using different beam qualities and different entrance air KERMA. Methods, previously developed for synthetic images and simulations, are adapted to real measurements, solving practical implementation issues.

RESULTS

In the case of the star-bar, in order to obtain optimal MTF determinations it is necessary to accurately determine the center of the pattern. Also, to avoid underestimates in MTF calculations, the length in pixels of each of the scanning circumferences must be an integer multiple of the number of cycles in the pattern. Both methods, star-bar and edge, give similar mean values of the MTF in all cases analyzed. Also, the dependence with frequency of the experimental MTF standard deviation (SD) agrees with the theoretical expressions presented in previous works. In this regard, the precision is better for the star-bar method than for the edge and differences in precision between both methods are higher for the lowest beam quality.

CONCLUSIONS

The star-bar phantom can be used for MTF determinations with the advantage of having an improved precision. However, precision is reduced when the radiation quality increases. This fact suggests that, for the highest beam qualities, materials with an attenuation coefficient greater than that of steel should be used to manufacture the phantom.

Flight model characterization of the wide-field off-axis telescope for the MATS satellite

We present optical characterization, calibration, and performance tests of the Mesospheric Airglow/Aerosol Tomography Spectroscopy (MATS) satellite, which for the first time, to the best of our knowledge, for a satellite, applies a linear-astigmatism-free confocal off-axis reflective optical design. Mechanical tolerances of the telescope were investigated using Monte Carlo methods and single-element perturbations. The sensitivity analysis results indicate that tilt errors of the tertiary mirror and a surface RMS error of the secondary mirror mainly degrade optical performance. From the Monte Carlo simulation, the tolerance limits were calculated to ±0.5mm, ±1mm, and ±0.15^∘ for decenter, despace, and tilt, respectively. We performed characterization measurements and optical tests with the flight model of the satellite. Multi-channel relative pointing, total optical system throughput, and distortion of each channel were characterized for end-users. Optical performance was evaluated by measuring the modulation transfer function (MTF) and point spread function (PSF). The final MTF performance was 0.25 MTF at 20 lp/mm for the ultraviolet channel (304.5 nm), and 0.25-0.54 MTF at 10 lp/mm for infrared channels. The salient fact of the PSF measurement of this system is that there is no noticeable linear astigmatism detected over a wide field of view (5.67^∘×0.91^∘). All things considered, the design method showed great advantages in wide field of view observations with satellite-level optical performance.

Evaluation of a hybrid direct-indirect active matrix flat-panel imager using Monte Carlo simulation

Purpose: Monte Carlo simulations were used to evaluate the imaging properties of a composite direct-indirect active matrix flat-panel imager (AMFPI) with potentially more favorable tradeoffs between x-ray quantum efficiency and spatial resolution than direct or indirect AMFPIs alone. This configuration, referred to as a hybrid AMFPI, comprises a scintillator that is optically coupled to an a-Se direct AMFPI through a transparent electrode and hole blocking layer, such that a-Se acts as both a direct x-ray converter and an optical sensor. Approach: GEANT4 was used to simulate x-ray energy deposition, optical transport, and charge signal generation processes in various hybrid AMPFI configurations under RQA5 and RQA9 x-ray beam conditions. The Fujita-Lubberts-Swank method was used to quantify the impact of irradiation geometry, x-ray converter thicknesses, conversion gain of each layer, and x-ray cross talk between layers on detective quantum efficiency (DQE). Results: Each hybrid configuration had a greater DQE than its direct AMFPI layer alone. The DQE improvement was largest at low spatial frequencies in both front- and back-irradiation (BI) geometries due to increased x-ray quantum efficiency provided by the scintillator. DQE improvements persisted at higher frequencies in BI geometry due to preferential x-ray absorption in a-Se. Matching the x-ray-to-charge conversion gains of a hybrid AMFPI's direct and indirect detection layers affects its Swank factor and, thus, DQE(0). X-ray cross talk has a negligible impact on the DQE ( f ) of hybrid AMFPIs with sufficiently high optical quantum efficiency. Conclusion: An optimized hybrid AMFPI can achieve greater DQE performance than current direct or indirect AMFPIs.

Development of a high resolution x-ray inspection system using a carbon nanotube based miniature x-ray tube

A new concept for a non-destructive testing device using a novel carbon nanotube (CNT) based miniature x-ray tube is proposed. The device can be used for small-scale internal inspection of objects. To investigate the effectiveness of the proposed concept, the device was fabricated and its performance was systematically analyzed. The non-destructive testing device consists of a CNT based miniature x-ray tube, a scintillator, an optical lens, and a detector. The size of the focal spot needed to identify objects as small as 5 µm was calculated through simulation. An electron optics simulation software, E-GUN, was used to optimize the geometries of both the focusing cup and the x-ray target to achieve the desired focal spot size of the x-ray tube. The CNT based miniature x-ray tube was fabricated using the brazing process, and an NdFeB focusing lens was used to further reduce the focal spot size. XR images were obtained using the fabricated device and the spatial resolutions of the images were evaluated using the modulation transfer function (MTF). The fields of view (FOVs) per probe are 7.1 mm² and 1.8 mm² when using a 5× optical lens and a 10× optical lens, respectively. The FOV can be increased by increasing the number of probes incorporated into the device. MTF10 values were determined to be 105 lp/mm and 230 lp/mm when using the 5× optical lens and 10× optical lens, respectively. By using an optical lens to enlarge the XR images, the effect of focal spot was minimized and clear XR images were obtained.

An On-Board Spectral-CT/CBCT/SPECT Imaging Configuration for Small-Animal Radiation Therapy Platform: A Monte Carlo Study

This study investigated the feasibility of a highly specific multiplexed image-guided small animal radiation therapy (SART) platform based on triple imaging from on-board single-photon emission computed tomography (SPECT), spectral-CT, and cone-beam CT (CBCT) guidance in radiotherapy treatment. As a proof-of-concept, the SART system was built with the capability of triple on-board image guidance by utilizing an x-ray tube and a single cadmium zinc telluride (CZT) semiconductor photon-counting imager via a Monte Carlo simulation study. The x-ray tube can be set at a low tube current for imaging mode and a high tube current for radiation therapy mode, respectively. In the imaging mode, both x-ray and gamma-ray projection data were collected by the imager to reconstruct CBCT, SPECT and spectral CT images of small animals being treated. The modulation transfer function (MTF) of the pixelated CZT imager measured was 8.6 lp/mm. The overall performances of the CBCT and SPECT imaging of the system were evaluated with sufficient spatial resolution and imaging quality to be fitted into the SART platform. The material differentiation and decomposition capacities of spectral CT within the system were verified using K-edge imaging, image-based optimal energy weighted imaging, and image-based linear material decomposition methods. The triple imaging capability of the system was demonstrated using a PMMA phantom containing gadolinium, iodine and radioisotope ^99mTc inserts. All the probes were clearly identified in the registered image. The results demonstrated that a novel SART platform with high-quality on-board CBCT, spectral-CT, SPECT image guidance is technically feasible by using a single semiconductor imager, thus affording comprehensive image guidance from anatomical, functional, and molecular levels for radiation treatment beam delivery.

Towards 4D dedicated breast CT perfusion imaging of cancer: development and validation of computer simulated images

Dedicated breast CT is a fully tomographic breast imaging modality with potential for various applications throughout breast cancer care. If implemented to perform dynamic contrast-enhanced (CE) imaging (4D breast CT), it could be useful to obtain functional information at high combined spatio-temporal resolution. Before developing a 4D dedicated breast CT system, a computer simulation method for breast CT perfusion imaging is proposed. The simulation uses previously developed patient-based 4D digital breast phantoms, and generates realistic images with the selected acquisition parameters, allowing to investigate the effect of different acquisition settings on image quality. The simulation pipeline includes all steps of the image generation process, from ray tracing and scatter map generation, to the addition of realistic resolution losses and noise models. The pipeline was validated against experimental measurements performed on physical phantoms with a dedicated breast CT system, in terms of average error compared to ground truth projections (6.0%  ±  1.65%), Hounsfield unit (HU) values in a homogeneous phantom (acquired: -149 HU  ±  2 HU; simulated:  -140 HU  ±  2 HU), signal-to-noise ratio (SNR) (average error 6.7%  ±  4.2%), noise power spectra (NPS) (average error 4.3%  ±  2.5%), modulation transfer function (MTF) (average error 8.4%  ±  7.2%), and attenuation of different adipose/glandular equivalent mixtures (average error 6.9%  ±  4.0%) and glandular plus iodinated contrast medium concentrations equivalent mixtures (average error of 9.1%  ±  9.0%). 4D patient images were then simulated for different 4D digital breast phantoms at different air kerma levels to determine the effect of noise on the extracted tumor perfusion curves. In conclusion, the proposed pipeline could simulate images with a good level of realism, resulting in a tool that can be used for the design, development, and optimization of a 4D dedicated breast CT system.

Task-driven optimization of the non-spectral mode of photon counting CT for intracranial hemorrhage assessment

Non-contrast CT (NCCT) is widely employed as the first-line imaging test to evaluate intracranial hemorrhage (ICH). Advances in mutidetector CT (MDCT) technology have greatly improved the image quality of NCCT for the detection of established, relatively large, and acute ICHs. Meanwhile, the reliability of MDCT in detecting microbleeds and chronic hemorrhage, and in predicting hemorrhagic transformation needs to be further improved. The purpose of this work was to investigate the potential use of non-spectral photon counting CT (PCCT) to address these challenges in ICH imaging. Towards this goal, the NCCT protocol of an experimental PCCT system that simulates the geometry of a general-purpose MDCT was optimized. The optimization was driven by three imaging tasks: detection of a 4.0 mm intraparenchymal hemorrhage, detection of a 1.5 mm subarachnoid hemorrhage, and discrimination of a sulcus in the insular cortex from the parenchymal background. These imaging tasks were custom-built into an anthropomorphic head phantom. Under the guidance of the frequency-dependent noise equivalent quanta and the ideal observer model detectability index [Formula: see text], the optimal PCD detection mode, energy threshold, and reconstruction kernel were found to be the anti-charge sharing mode, 15 keV, and an apodized ramp kernel, respectively. Compared with a clinical MDCT operated with an ICH protocol and at a matched dose level, the PCCT system provided at least 20% improvements in [Formula: see text] for all three ICH imaging tasks. These results demonstrated the potential benefits of non-spectral PCCT in ICH assessment.

High-numerical-aperture macroscope optics for time-resolved experiments

A novel high-quality custom-made macroscope optics, dedicated to high-resolution time-resolved X-ray tomographic microscopy at the TOMCAT beamline at the Swiss Light Source (Paul Scherrer Institut, Switzerland), is introduced. The macroscope offers 4× magnification, has a very high numerical aperture of 0.35 and it is modular and highly flexible. It can be mounted both in a horizontal and vertical configuration, enabling imaging of tall samples close to the scintillator, to avoid edge-enhancement artefacts. The macroscope performance was characterized and compared with two existing in-house imaging setups, one dedicated to high spatial and one to high temporal resolution. The novel macroscope shows superior performance for both imaging settings compared with the previous systems. For the time-resolved setup, the macroscope is 4 times more efficient than the previous system and, at the same time, the spatial resolution is also increased by a factor of 6. For the high-spatial-resolution setup, the macroscope is up to 8.5 times more efficient with a moderate spatial resolution improvement (factor of 1.5). This high efficiency, increased spatial resolution and very high image quality offered by the novel macroscope optics will make 10-20 Hz high-resolution tomographic studies routinely possible, unlocking unprecedented possibilities for the tomographic investigations of dynamic processes and radiation-sensitive samples.

[Assessment Tool of Image Quality for Medical Liquid-crystal Displays Using Commercially Available Digital Cameras]

We developed a simplified tool for measuring image quality of medical liquid-crystal displays (LCDs) using a commercially available color digital camera. This tool implemented as a plug-in software for ImageJ (open-source image processing program) was designed to compute modulation transfer functions (MTFs) and Wiener spectra (WS) of monochrome and color LCDs from LCD photographed images captured by a camera. The intensities of the red (R), green (G), and blue (B) signals of the unprocessed image data depend on the spectral sensitivity of the image sensor used in the camera. In order to evaluate image quality based on LCD luminance, the plug-in software calibrates the RGB signals from the camera using measured luminance of the LCD and converts them into grayscale signals that correspond to the luminance of the LCD. The MTFs and WS are determined based on the line response from a one-pixel line image and the one-dimensional noise profiles acquired by scanning the uniform image using numerically synthesized slit, respectively. With this plug-in software for ImageJ, we are able to readily compute MTFs and WS of both monochrome and color LCDs from unprocessed image data of cameras. Our simplified tool is helpful to evaluate and understand the physical performance of LCDs for a large number of display users in hospitals and medical centers.

Characterizing a novel scintillating glass for application to megavoltage cone-beam computed tomography

PURPOSE

The purpose of this study was to evaluate the performance of a prototype electric portal imaging device (EPID) with a high detective quantum efficiency (DQE) scintillator, LKH-5. Specifically, image quality in context of both planar and megavoltage (MV) cone-beam computed tomography (CBCT) is analyzed.

METHODS

Planar image quality in terms of modulation transfer function (MTF), noise power spectrum (NPS), and DQE are measured and compared to an existing EPID (AS-1200) using the 6 MV beamline for a Varian TrueBeam linac. Imager performance is contextualized for three-dimensional (3D), MV-CBCT performance by measuring imager lag and analyzing the expected degradation of the DQE as a function of dose. Finally, comparisons between reconstructed images of the Catphan phantom in terms of qualitative quality and signal-difference-to-noise ratio (SDNR) are made for 6 MV images using both conventional and LKH-5 EPIDs as well as for the kilovoltage (kV) on-board imager (OBI).

RESULTS

Analysis of the NPS reveals linearity at all measured doses using the prototype LKH-5 detector. While the first zero of the MTF is much lower for the LKH-5 detector than the conventional EPID (0.6 cycles/mm vs 1.6 cycles/mm), the normalized NPS (NNPS) multiplied by total quanta (qNNPS) of the LKH-5 detector is roughly a factor of seven to eight times lower, yielding a DQE(0) of approximately 8%. First, second, and third frame lag were measured at approximately 23%, 5%, and 1%, respectively, although no noticeable image artifacts were apparent in reconstructed volumes. Analysis of low-dose performance reveals that DQE(0) remains at 80% of its maximum value at a dose as low as 7.5 × 10^-6  MU. For a 400 projection technique, this represents a total scan dose of 0.0030 MU, suggesting that if imaging doses are increased to a value typical of kV-CBCT scans (~2.7 cGy), the LKH-5 detector will retain quantum noise limited performance. Finally, comparing Catphan scans, the prototype detector exhibits much lower image noise than the conventional EPID, resulting in improved small object representation. Furthermore, SDNR of H₂ O and polystyrene cylinders improved from -1.95 and 2.94 to -15 and 18.7, respectively.

CONCLUSIONS

Imaging performance of the prototype LKH-5 detector was measured and analyzed for both planar and 3D contexts. Improving noise transfer of the detector results in concurrent improvement of DQE(0). For 3D imaging, temporal characteristics were adequate for artifact-free performance and at relevant doses, the detector retained quantum noise limited performance. Although quantitative MTF measurements suggest poorer resolution, small object representation of the prototype imager is qualitatively improved over the conventional detector due to the measured reduction in noise.

Nonstationary model of oblique x-ray incidence in amorphous selenium detectors: II. Transfer functions

PURPOSE

One limitation of experimental techniques for quantifying resolution and noise in detectors is that the measurement is made in a region-of-interest (ROI). With theoretical modeling, these properties can be measured at a point, allowing for quantification of spatial anisotropy. This paper calculates nonstationary transfer functions for amorphous selenium (a-Se) detectors in breast imaging. We use this model to demonstrate the performance advantage of a "next-generation" tomosynthesis (NGT) system, which is capable of x-ray source motion with more degrees of freedom than a clinical tomosynthesis system.

METHODS

Using Swank's formulation, the optical transfer function (OTF) and presampled noise power spectra (NPS) are determined based on the point spread function derived in Part 1. The modulation transfer function (MTF) is found from the normalized modulus of the OTF. To take into account the presence of digitization, the presampled NPS is convolved with a two-dimensional comb function, for which the period along each direction is the reciprocal of the detector element size. The detective quantum efficiency (DQE) is then determined from combined knowledge of the OTF and NPS.

RESULTS

First, the model is used to demonstrate the loss of image quality due to oblique x-ray incidence. The MTF is calculated along various polar angles, corresponding to different orientations of the input frequency. The MTF is independent of the incidence angle if the polar angle is perpendicular to the ray incidence direction. However, along other polar angles, oblique incidence results in MTF degradation at high frequencies. The MTF degradation is most substantial along the ray incidence direction. Unlike the MTF, the normalized NPS (NNPS) is independent of the incidence angle. To measure the relative signal-to-noise, the DQE is also calculated. Oblique incidence yields high-frequency DQE degradation, which is more pronounced than the MTF degradation. This arises because the DQE is proportionate with the square of the MTF. Ultimately, this model is used to evaluate how the image quality varies over the detector area. For various projection images, we calculate the variation in the incidence angle over this area. With the NGT system, the source can be positioned in such a way that this variation is minimized, and hence the DQE exhibits less anisotropy. To achieve this improvement in the image quality, the source needs to have a component of motion in the posteroanterior (PA) direction, which is perpendicular to the conventional direction of source motion in tomosynthesis.

CONCLUSIONS

In a-Se detectors, the DQE at high frequencies is degraded due to oblique incidence. The DQE degradation is more pronounced than the MTF degradation. This model is used to quantify the spatial variation in DQE over the detector area. The use of PA source motion is a strategy for minimizing this variation and thus improving the image quality.