Optimization of a dual-energy contrast-enhanced technique for a photon-counting digital breast tomosynthesis system: I. A theoretical model

Exploring the Potential of a Novel Iodine-Based Material as an Alternative Contrast Agent in X-ray Imaging Studies

BACKGROUND

Contrast-enhanced mammography is one of the new emerging imaging techniques used for detecting breast tissue lesions. Optimization of imaging protocols and reconstruction techniques for this modality, however, requires the involvement of physical phantoms. Their development is related to the use of radiocontrast agents. This study assesses the X-ray properties of a novel contrast material in clinical settings. This material is intended for experimental use with physical phantoms, offering an alternative to commonly available radiocontrast agents.

MATERIALS AND METHODS

The water-soluble sodium salt of the newly synthesized diiodine-substituted natural eudesmic acid, Sodium 2,6-DiIodo-3,4,5-TriMethoxyBenzoate [NaDITMB], has been investigated with respect to one of the most commonly applied radiocontrast medium in medical practice-Omnipaque®. For this purpose, simulation and experimental studies were carried out with a computational phantom and a physical counterpart, respectively. Synthetic and experimental X-ray images were subsequently produced under varying beam kilovoltage peaks (kVps), and the proposed contrast material was evaluated.

RESULTS AND DISCUSSION

Simulation results revealed equivalent absorptions between the two simulated radiocontrast agents. Experimental findings supported these simulations, showing a maximum deviation of 3.7% between the image gray values of contrast materials for NaDITMB and Omnipaque solutions for a 46 kVp X-ray beam. Higher kVp X-ray beams show even smaller deviations in the mean grey values of the imaged contrast agents, with the NaDITMB solution demonstrating less than a 2% deviation compared to Omnipaque.

CONCLUSION

The proposed contrast agent is a suitable candidate for use in experimental work related to contrast-enhanced imaging by utilizing phantoms. It boasts the advantages of easy synthesis and is recognized for its safety, ensuring a secure environment for both the experimenter and the environment.

Comparison of contrast-enhanced digital mammography and contrast-enhanced digital breast tomosynthesis for lesion assessment

Contrast-enhanced digital mammography (CEDM) reveals neovasculature of breast lesions in a two-dimensional contrast enhancement map. Contrast-enhanced digital breast tomosynthesis (CEDBT) provides contrast enhancement in three dimensions, which may improve lesion characterization and localization. We aim to compare CEDM and CEDBT for lesion assessment. Women with breast imaging-reporting and data system 4 or 5 suspicious breast lesion(s) were recruited in our study and were imaged with CEDM and CEDBT in succession under one breast compression. Two radiologists assessed CEDM and CEDBT with both images displayed side-by-side and compared (1) contrast enhancement of lesions and (2) lesion margin using a five-point scale ranging from - 2 (CEDM much better) to + 2 (CEDBT much better). Biopsy identified 19 malignant lesions with contrast enhancement. Our results show that CEDBT provides better lesion margins than CEDM with limited reduction in contrast enhancement. CEDBT delivers less radiation dose compared to CEDM + DBT. Synthetic CEDM can be generated from CEDBT data and provides lesion contrast enhancement comparable to CEDM. CEDBT has potential for clinical applications, such as treatment response monitoring and guidance for biopsy.

Non-stationary model of oblique x-ray incidence in amorphous selenium detectors: I. Point spread function

PURPOSE

In previous work, a theoretical model of the point spread function (PSF) for oblique x-ray incidence in amorphous selenium (a-Se) detectors was proposed. The purpose of this paper is to develop a complementary model that includes two additional features. First, the incidence angle and the directionality of ray incidence are calculated at each position, assuming a divergent x-ray beam geometry. This approach allows the non-stationarity of the PSF to be modeled. Second, this paper develops a framework that is applicable to a digital system, unlike previous work which did not model the presence of a thin-film transistor (TFT) array.

METHODS

At each point on the detector, the incidence angle and the ray incidence direction are determined using ray tracing. Based on these calculations, an existing model for the PSF of the x-ray converter (Med Phys. 1995;22:365-374) is generalized to a non-stationary model. The PSF is convolved with the product of two rectangle functions, which model the sampling of the TFT array. The rectangle functions match the detector element (del) size in two dimensions.

RESULTS

It is shown that the PSF can be calculated in closed form. This solution is used to simulate a digital mammography (DM) system at two x-ray energies (20 and 40 keV). Based on the divergence of the x-ray beam, the direction of ray incidence varies with position. Along this direction, the PSF is broader than the reference rect function matching the del size. The broadening is more pronounced with increasing obliquity. At high energy, the PSF deviates more strongly from the reference rect function, indicating that there is more blurring. In addition, the PSF is calculated along the polar angle perpendicular to the ray incidence direction. For this polar angle, the shape of the PSF is dependent upon whether the ray incidence direction is parallel with the sides of the detector. If the ray incidence direction is parallel with either dimension, the PSF is a perfect rectangle function, matching the del size. However, if the ray incidence direction is at an oblique angle relative to the sides of the detector, the PSF is not rectangular. These results illustrate the non-stationarity of the PSF.

CONCLUSIONS

This paper demonstrates that an existing model of the PSF of a-Se detectors can be generalized to include the effects of non-stationarity and digitization. The PSF is determined in closed form. This solution offers the advantage of shorter computation time relative to approaches that use numerical methods. This model is a tool for simulating a-Se detectors in future work, such as in virtual clinical trials with computational phantoms.

An all-in-one nanoparticle (AION) contrast agent for breast cancer screening with DEM-CT-MRI-NIRF imaging

Conventional X-ray mammography has low diagnostic sensitivity for women with dense breasts. As a result, alternative contrast-enhanced screening tools such as dual energy mammography (DEM), computed tomography (CT), magnetic resonance imaging (MRI), and near-infrared fluorescence (NIRF) imaging are being used or investigated for these women. However, currently available contrast agents are non-ideal, have safety issues, and each imaging technique requires a different contrast agent. We therefore sought to develop a multimodal contrast agent that is functional for each breast imaging modality to simplify the diagnosis process and address the issues of existing contrast agents. Herein, we present a novel "all-in-one" nanoparticle (AION) multimodal imaging probe that has potent DEM, CT, MRI, and NIRF contrast properties and improved biocompatibility. AION were formed by co-encapsulation of a near-infrared fluorophore (DiR), silver sulfide nanoparticles (Ag2S-NP), and iron oxide nanoparticles (IO-NP) in PEGylated micelles. AION showed negligible cytotoxicity, which was in agreement with its minimal silver ion release profiles. AION generated strong contrast with all imaging modalities as demonstrated in phantom imaging. AION allowed in vivo tumor imaging as evidenced by the increase in contrast after injection. This study indicates the potential of AION as an effective multimodal contrast agent for breast cancer diagnosis with a range of imaging methods.

Susceptibility of iodine concentration map of dual-energy contrast-enhanced digital mammography for quantitative and tumor enhancement assessment

Background Dual-energy (DE) contrast-enhanced digital mammography (DE-CEDM) provides additional information on tumor angiogenesis. Purpose To investigate the susceptibility of reconstructing color-coded iodine concentration maps on the basis of quantitative calibrations of the iodine concentration and contrast-to-noise ratio (CNR) in DE-CEDM applications. Material and Methods A custom-made phantom filled with iodine concentrations in the range of 0.1-10 mg/cm² was used in calibrations. All DE images were acquired using the GE Senographe Essential system. From DE subtraction images, the image contrast and CNR were obtained, and the quantitative relationship between these two metrics and the iodine concentration at each phantom thickness was investigated. The quantitative CNR calibration curves were applied to reconstruct color-coded iodine maps on a pixel-by-pixel basis. Results Both the mean contrast and mean CNR increased linearly with the iodine concentration. The iodine concentration estimated from the iodine map reconstructed from quantitative CNR calibrations was highly consistent with the desired iodine concentration (R²= 0.989), and smaller relative errors (in the range of 3.0-19.5%) were observed with iodine concentrations not less than 1 mg/cm². Conclusion An iodine concentration map could be reconstructed based on the linear relationship between the CNR and iodine concentration. From the color-coded iodine concentration map, the contrast medium enhancement phenomenon could be further estimated quantitatively, and tumor enhancement patterns could be easily observed.

Phantom-based study exploring the effects of different scatter correction approaches on the reconstructed images generated by contrast-enhanced stationary digital breast tomosynthesis

Stationary digital breast tomosynthesis (sDBT) is an emerging technology in which the single rotating x-ray tube is replaced by a fixed array of multiple carbon nanotube-enabled sources, providing a higher spatial and temporal resolution. As such, sDBT offers a promising platform for contrast-enhanced (CE) imaging. However, given the minimal enhancement above background with standard operational tube settings and iodine dosing, CE breast imaging requires additional acquisition steps to isolate the iodine signal, using either temporal or dual energy subtraction (TS or DES) protocols. Also, correcting for factors that limit contrast is critical, and scatter and noise pose unique challenges during tomosynthesis. This phantom-based study of CE sDBT compared different postacquisition scatter correction approaches on the quality of the reconstructed image slices. Beam-pass collimation was used to sample scatter indirectly, from which an interpolated scatter map was obtained for each projection image. Scatter-corrected projections provided the information for reconstruction. Comparison between the application of different scatter maps demonstrated the significant effect that processing has on the contrast-to-noise ratio and feature detectability ([Formula: see text]) in the final displayed images and emphasized the critical importance of scatter correction during DES.

Low-Dose Contrast-Enhanced Breast CT Using Spectral Shaping Filters: An Experimental Study

Iodinated contrast-enhanced X-ray imaging of the breast has been studied with various modalities, including full-field digital mammography (FFDM), digital breast tomosynthesis (DBT), and dedicated breast CT. Contrast imaging with breast CT has a number of advantages over FFDM and DBT, including the lack of breast compression, and generation of fully isotropic 3-D reconstructions. Nonetheless, for breast CT to be considered as a viable tool for routine clinical use, it would be desirable to reduce radiation dose. One approach for dose reduction in breast CT is spectral shaping using X-ray filters. In this paper, two high atomic number filter materials are studied, namely, gadolinium (Gd) and erbium (Er), and compared with Al and Cu filters currently used in breast CT systems. Task-based performance is assessed by imaging a cylindrical poly(methyl methacrylate) phantom with iodine inserts on a benchtop breast CT system that emulates clinical breast CT. To evaluate detectability, a channelized hoteling observer (CHO) is used with sums of Laguerre-Gauss channels. It was observed that spectral shaping using Er and Gd filters substantially increased the dose efficiency (defined as signal-to-noise ratio of the CHO divided by mean glandular dose) as compared with kilovolt peak and filter settings used in commercial and prototype breast CT systems. These experimental phantom study results are encouraging for reducing dose of breast CT, however, further evaluation involving patients is needed.

Model and reconstruction of a K-edge contrast agent distribution with an X-ray photon-counting detector

Contrast-enhanced computed tomography (CECT) helps enhance the visibility for tumor imaging. When a high-Z contrast agent interacts with X-rays across its K-edge, X-ray photoelectric absorption would experience a sudden increment, resulting in a significant difference of the X-ray transmission intensity between the left and right energy windows of the K-edge. Using photon-counting detectors, the X-ray intensity data in the left and right windows of the K-edge can be measured simultaneously. The differential information of the two kinds of intensity data reflects the contrast-agent concentration distribution. K-edge differences between various matters allow opportunities for the identification of contrast agents in biomedical applications. In this paper, a general radon transform is established to link the contrast-agent concentration to X-ray intensity measurement data. An iterative algorithm is proposed to reconstruct a contrast-agent distribution and tissue attenuation background simultaneously. Comprehensive numerical simulations are performed to demonstrate the merits of the proposed method over the existing K-edge imaging methods. Our results show that the proposed method accurately quantifies a distribution of a contrast agent, optimizing the contrast-to-noise ratio at a high dose efficiency.

Investigation of x-ray spectra for iodinated contrast-enhanced dedicated breast CT

Screening for breast cancer with mammography has been very successful, resulting in part to a reduction of breast cancer mortality by approximately 39% since 1990. However, mammography still has limitations in performance, especially for women with dense breast tissue. Iodinated contrast-enhanced, dedicated breast CT (BCT) has been proposed to improve lesion analysis and the accuracy of diagnostic workup for patients suspected of having breast cancer. A mathematical analysis to explore the use of various x-ray filters for iodinated contrast-enhanced BCT is presented. To assess task-based performance, the ideal linear observer signal-to-noise ratio (SNR) is used as a figure-of-merit under the assumptions of a linear, shift-invariant imaging system. To estimate signal and noise propagation through the BCT detector, a parallel-cascade model was used. The lesion model was embedded into a structured background and included a realistic level of iodine uptake. SNR was computed for 84,000 different exposure settings by varying the kV setting, x-ray filter materials and thickness, breast size, and composition and radiation dose. It is shown that some x-ray filter material/thickness combinations can provide up to 75% improvement in the linear ideal observer SNR over a conventionally used x-ray filter for BCT. This improvement in SNR can be traded off for substantial reductions in mean glandular dose.

Feasibility study of dual energy radiographic imaging for target localization in radiotherapy for lung tumors

Purpose

Dual-energy (DE) radiographic imaging improves tissue discrimination by separating soft from hard tissues in the acquired images. This study was to establish a mathematic model of DE imaging based on intrinsic properties of tissues and quantitatively evaluate the feasibility of applying the DE imaging technique to tumor localization in radiotherapy.

Methods

We investigated the dependence of DE image quality on the radiological equivalent path length (EPL) of tissues with two phantoms using a stereoscopic x-ray imaging unit. 10 lung cancer patients who underwent radiotherapy each with gold markers implanted in the tumor were enrolled in the study approved by the hospital's Ethics Committee. The displacements of the centroids of the delineated gross tumor volumes (GTVs) in the digitally reconstructed radiograph (DRR) and in the bone-canceled DE image were compared with the averaged displacements of the centroids of gold markers to evaluate the feasibility of using DE imaging for tumor localization.

Results

The results of the phantom study indicated that the contrast-to-noise ratio (CNR) was linearly dependent on the difference of EPL and a mathematical model was established. The objects and backgrounds corresponding to ΔEPL less than 0.08 are visually indistinguishable in the bone-canceled DE image. The analysis of patient data showed that the tumor contrast in the bone-canceled images was improved significantly as compared with that in the original radiographic images and the accuracy of tumor localization using the DE imaging technique was comparable with that of using fiducial makers.

Conclusion

It is feasible to apply the technique for tumor localization in radiotherapy.

Exploring silver as a contrast agent for contrast-enhanced dual-energy X-ray breast imaging

OBJECTIVE

Through prior monoenergetic modelling, we have identified silver as a potential alternative to iodine in dual-energy (DE) X-ray breast imaging. The purpose of this study was to compare the performance of silver and iodine contrast agents in a commercially available DE imaging system through a quantitative analysis of signal difference-to-noise ratio (SDNR).

METHODS

A polyenergetic simulation algorithm was developed to model the signal intensity and noise. The model identified the influence of various technique parameters on SDNR. The model was also used to identify the optimal imaging techniques for silver and iodine, so that the two contrast materials could be objectively compared.

RESULTS

The major influences on the SDNR were the low-energy dose fraction and breast thickness. An increase in the value of either of these parameters resulted in a decrease in SDNR. The SDNR for silver was on average 43% higher than that for iodine when imaged at their respective optimal conditions, and 40% higher when both were imaged at the optimal conditions for iodine.

CONCLUSION

A silver contrast agent should provide benefit over iodine, even when translated to the clinic without modification of imaging system or protocol. If the system were slightly modified to reflect the lower k-edge of silver, the difference in SDNR between the two materials would be increased.

ADVANCES IN KNOWLEDGE

These data are the first to demonstrate the suitability of silver as a contrast material in a clinical contrast-enhanced DE image acquisition system.

A unified statistical framework for material decomposition using multienergy photon counting x-ray detectors

PURPOSE

Material decomposition using multienergy photon counting x-ray detectors (PCXD) has been an active research area over the past few years. Even with some success, the problem of optimal energy selection and three material decomposition including malignant tissue is still on going research topic, and more systematic studies are required. This paper aims to address this in a unified statistical framework in a mammographic environment.

METHODS

A unified statistical framework for energy level optimization and decomposition of three materials is proposed. In particular, an energy level optimization algorithm is derived using the theory of the minimum variance unbiased estimator, and an iterative algorithm is proposed for material composition as well as system parameter estimation under the unified statistical estimation framework. To verify the performance of the proposed algorithm, the authors performed simulation studies as well as real experiments using physical breast phantom and ex vivo breast specimen. Quantitative comparisons using various performance measures were conducted, and qualitative performance evaluations for ex vivo breast specimen were also performed by comparing the ground-truth malignant tissue areas identified by radiologists.

RESULTS

Both simulation and real experiments confirmed that the optimized energy bins by the proposed method allow better material decomposition quality. Moreover, for the specimen thickness estimation errors up to 2 mm, the proposed method provides good reconstruction results in both simulation and real ex vivo breast phantom experiments compared to existing methods.

CONCLUSIONS

The proposed statistical framework of PCXD has been successfully applied for the energy optimization and decomposition of three material in a mammographic environment. Experimental results using the physical breast phantom and ex vivo specimen support the practicality of the proposed algorithm.

A review of breast tomosynthesis. Part II. Image reconstruction, processing and analysis, and advanced applications

Many important post-acquisition aspects of breast tomosynthesis imaging can impact its clinical performance. Chief among them is the reconstruction algorithm that generates the representation of the three-dimensional breast volume from the acquired projections. But even after reconstruction, additional processes, such as artifact reduction algorithms, computer aided detection and diagnosis, among others, can also impact the performance of breast tomosynthesis in the clinical realm. In this two part paper, a review of breast tomosynthesis research is performed, with an emphasis on its medical physics aspects. In the companion paper, the first part of this review, the research performed relevant to the image acquisition process is examined. This second part will review the research on the post-acquisition aspects, including reconstruction, image processing, and analysis, as well as the advanced applications being investigated for breast tomosynthesis.

X-ray spectral measurements for tungsten-anode from 20 to 49 kVp on a digital breast tomosynthesis system

PURPOSE

This paper presents new spectral measurements of a tungsten-target digital breast tomosynthesis (DBT) system, including spectra of 43-49 kVp.

METHODS

Raw x-ray spectra of 20-49 kVp were directly measured from the tube port of a Selenia Dimensions DBT system using a CdTe based spectrometer. Two configurations of collimation were employed: one with two tungsten pinholes of 25 μm and 200 μm diameters, and the other with a single pinhole of 25 μm diameter, for acquiring spectra from the focal spot and from the focal spot as well as its vicinity. Stripping correction was applied to the measured spectra to compensate distortions due to escape events. The measured spectra were compared with the existing mammographic spectra of the TASMIP model in terms of photon fluence per exposure, spectral components, and half-value layer (HVL). HVLs were calculated from the spectra with a numerical filtration of 0.7 mm aluminum and were compared against actual measurements on the DBT system using W/Al (target-filter) combination, without paddle in the beam.

RESULTS

The spectra from the double-pinhole configuration, in which the acceptance aperture pointed right at the focal spot, were harder than the single-pinhole spectra which include both primary and off-focus radiation. HVL calculated from the single-pinhole setup agreed with the measured HVL within 0.04 mm aluminum, while the HVL values from the double-pinhole setup were larger than the single-pinhole HVL by at most 0.1 mm aluminum. The spectra from single-pinhole setup agreed well with the TASMIP mammographic spectra, and are more relevant for clinical purpose.

CONCLUSIONS

The spectra data would be useful for future research on DBT system with tungsten targets.

Optimization of phosphor-based detector design for oblique x-ray incidence in digital breast tomosynthesis

PURPOSE

In digital breast tomosynthesis (DBT), a volumetric reconstruction of the breast is generated from a limited range of x-ray projections. One trade-off of DBT is resolution loss in the projections due to non-normal (i.e., oblique) x-ray incidence. Although degradation in image quality due to oblique incidence has been studied using empirical data and Monte Carlo simulations, a theoretical treatment has been lacking. The purpose of this work is to extend Swank's calculations of the transfer functions of turbid granular phosphors to oblique incidence. The model is ultimately used as a tool for optimizing the design of DBT detectors.

METHODS

A quantum-limited system and 20 keV x-rays are considered. Under these assumptions, the modulation transfer function (MTF) and noise power spectra (NPS) are derived using the diffusion approximation to the Boltzmann equation to model optical scatter within the phosphor. This approach is applicable to a nonstructured scintillator such as gadolinium oxysulfide doped with terbium (Gd(2)O(2)S:Tb), which is commonly used in breast imaging and which can reasonably approximate other detector materials. The detective quantum efficiency (DQE) is then determined from the Nishikawa formulation, where it is written as the product of the x-ray quantum detection efficiency, the Swank factor, and the Lubberts fraction. Transfer functions are calculated for both front- and back-screen configurations, which differ by positioning the photocathode at the exit or entrance point of the x-ray beam, respectively.

RESULTS

In the front-screen configuration, MTF and DQE are found to have considerable angular dependence, while NPS is shown to vary minimally with projection angle. As expected, the high frequency MTF and DQE are degraded substantially at large angles. By contrast, all transfer functions for the back-screen configuration have the advantage of significantly less angular dependence. Using these models, we investigated the possibility for optimizing the design of DBT detectors. As an example optimization strategy, the phosphor thickness which maximizes the DQE at a fixed frequency is analyzed. This work demonstrates that the optimal phosphor thickness for the front-screen is angularly dependent, shifting to lower thickness at higher angles. Conversely, the back-screen is not optimized by a single thickness but instead attains reasonably high DQE values over a large range of thicknesses. Although the back-screen configuration is not suited for current detectors using a glass substrate, it may prove to be preferred in future detectors using newly proposed plastic thin-film transistor (TFT) substrates.

CONCLUSIONS

Using the diffusion approximation to the Boltzmann equation to model the spread of light in a scintillator, this paper develops an analytical model of MTF, NPS, and DQE for a phosphor irradiated obliquely. The model is set apart from other studies on oblique incidence in being derived from first principles. This work has applications in the optimization of DBT detector design.

Optimization of a dual-energy contrast-enhanced technique for a photon-counting digital breast tomosynthesis system: II. An experimental validation

PURPOSE

Previously, the authors developed a dual-energy (DE) acquisition technique for a photon-counting digital breast tomosynthesis (DBT) imaging system. Low-energy (LE) and high-energy (HE) images are acquired in a single scan by covering alternate slits of a multislit prepatient collimator with Sn and Cu, respectively. A theoretical model was used to optimize the technique. In this article, an experimental validation of this technique is presented.

METHODS

Experiments were performed on a prototype DBT system. LE and HE projection images were acquired sequentially; either a Sn or a Cu filter was positioned in the filter holder at the exit window of the x-ray tube. Sn filters from 0.113 to 0.242 mm thick and Cu filters from 0.103 to 0.267 mm were used. The images were acquired with a W target at 49 kV. Tomographic images, hereafter referred to as DBT images, were reconstructed using a shift-and-add algorithm. DE-DBT images were obtained by weighted logarithmic subtraction of the LE and HE images. Weighting factors w(t) that optimally cancel breast tissues with two different glandularities were assessed for 20-80 mm thick phantoms with 0%, 50%, and 100% glandularity. The mean and standard deviation in the per-pixel signal intensity (SI) were calculated in the DBT images. These data were used to calculate signal-difference-to-noise ratios (SDNRs) between iodine enhanced and nonenhanced polymethyl methacrylate backgrounds. To illustrate the feasibility of the technique, DE-DBT images of a structured phantom containing iodine disks were assessed. The experimental results were compared against the values obtained from a theoretical model of the imaging system.

RESULTS

The average difference between theoretical and experimental w(t) was found to range from 8% to 21%. Experimental w(t) values increase with phantom thickness and Cu thickness, depend somewhat on Sn thickness, and vary more as a function of breast composition in thick breasts than in thin breasts. Theoretical and experimental mean and standard deviation in the per-pixel SI differ by -7% to 10% and by -3% to 4%. Theoretical and experimental SDNR values differ, on average, by 1.5%. Iodine concentrations can be predicted from SDNR; the relationship can be accurately fit to a quadratic. In the images of the structured phantom, iodine concentrations of 1 mg/cm2 and larger are discernable.

CONCLUSIONS

The strong agreement between experimental and theoretical results in this article indicates that the authors' computer model is accurate.