H Ding

Quantification of breast density with spectral mammography based on a scanned multi-slit photon-counting detector: a feasibility study

A simple and accurate measurement of breast density is crucial for the understanding of its impact in breast cancer risk models. The feasibility to quantify volumetric breast density with a photon-counting spectral mammography system has been investigated using both computer simulations and physical phantom studies. A computer simulation model involved polyenergetic spectra from a tungsten anode x-ray tube and a Si-based photon-counting detector has been evaluated for breast density quantification. The figure-of-merit (FOM), which was defined as the signal-to-noise ratio of the dual energy image with respect to the square root of mean glandular dose, was chosen to optimize the imaging protocols, in terms of tube voltage and splitting energy. A scanning multi-slit photon-counting spectral mammography system has been employed in the experimental study to quantitatively measure breast density using dual energy decomposition with glandular and adipose equivalent phantoms of uniform thickness. Four different phantom studies were designed to evaluate the accuracy of the technique, each of which addressed one specific variable in the phantom configurations, including thickness, density, area and shape. In addition to the standard calibration fitting function used for dual energy decomposition, a modified fitting function has been proposed, which brought the tube voltages used in the imaging tasks as the third variable in dual energy decomposition. For an average sized 4.5 cm thick breast, the FOM was maximized with a tube voltage of 46 kVp and a splitting energy of 24 keV. To be consistent with the tube voltage used in current clinical screening exam (∼32 kVp), the optimal splitting energy was proposed to be 22 keV, which offered a FOM greater than 90% of the optimal value. In the experimental investigation, the root-mean-square (RMS) error in breast density quantification for all four phantom studies was estimated to be approximately 1.54% using standard calibration function. The results from the modified fitting function, which integrated the tube voltage as a variable in the calibration, indicated a RMS error of approximately 1.35% for all four studies. The results of the current study suggest that photon-counting spectral mammography systems may potentially be implemented for an accurate quantification of volumetric breast density, with an RMS error of less than 2%, using the proposed dual energy imaging technique.

Image-based spectral distortion correction for photon-counting x-ray detectors

PURPOSE To investigate the feasibility of using an image-based method to correct for distortions induced by various artifacts in the x-ray spectrum recorded with photon-counting detectors for their application in breast computed tomography (CT). METHODS The polyenergetic incident spectrum was simulated with the tungsten anode spectral model using the interpolating polynomials (TASMIP) code and carefully calibrated to match the x-ray tube in this study. Experiments were performed on a Cadmium-Zinc-Telluride (CZT) photon-counting detector with five energy thresholds. Energy bins were adjusted to evenly distribute the recorded counts above the noise floor. BR12 phantoms of various thicknesses were used for calibration. A nonlinear function was selected to fit the count correlation between the simulated and the measured spectra in the calibration process. To evaluate the proposed spectral distortion correction method, an empirical fitting derived from the calibration process was applied on the raw images recorded for polymethyl methacrylate (PMMA) phantoms of 8.7, 48.8, and 100.0 mm. Both the corrected counts and the effective attenuation coefficient were compared to the simulated values for each of the five energy bins. The feasibility of applying the proposed method to quantitative material decomposition was tested using a dual-energy imaging technique with a three-material phantom that consisted of water, lipid, and protein. The performance of the spectral distortion correction method was quantified using the relative root-mean-square (RMS) error with respect to the expected values from simulations or areal analysis of the decomposition phantom. RESULTS The implementation of the proposed method reduced the relative RMS error of the output counts in the five energy bins with respect to the simulated incident counts from 23.0%, 33.0%, and 54.0% to 1.2%, 1.8%, and 7.7% for 8.7, 48.8, and 100.0 mm PMMA phantoms, respectively. The accuracy of the effective attenuation coefficient of PMMA estimate was also improved with the proposed spectral distortion correction. Finally, the relative RMS error of water, lipid, and protein decompositions in dual-energy imaging was significantly reduced from 53.4% to 6.8% after correction was applied. CONCLUSIONS The study demonstrated that dramatic distortions in the recorded raw image yielded from a photon-counting detector could be expected, which presents great challenges for applying the quantitative material decomposition method in spectral CT. The proposed semi-empirical correction method can effectively reduce these errors caused by various artifacts, including pulse pileup and charge sharing effects. Furthermore, rather than detector-specific simulation packages, the method requires a relatively simple calibration process and knowledge about the incident spectrum. Therefore, it may be used as a generalized procedure for the spectral distortion correction of different photon-counting detectors in clinical breast CT systems.

Dual-dictionary learning-based iterative image reconstruction for spectral computed tomography application

In this study, we investigated the effectiveness of a novel iterative reconstruction (IR) method coupled with dual-dictionary learning (DDL) for image reconstruction in a dedicated breast computed tomography (CT) system based on a cadmium-zinc-telluride (CZT) photon-counting detector and compared it to the filtered-back-projection (FBP) method with the ultimate goal of reducing the number of projections necessary for reconstruction without sacrificing the image quality. Postmortem breast samples were scanned in a fan-beam CT system and were reconstructed from 100 to 600 projections with both IR and FBP methods. The contrast-to-noise ratio (CNR) between the glandular and adipose tissues of the postmortem breast samples was calculated to compare the quality of images reconstructed from IR and FBP. The spatial resolution of the two reconstruction techniques was evaluated using aluminum (Al) wires with diameters of 643, 813, 1020, 1290 and 1630 µm in a plastic epoxy resin phantom with a diameter of 13 cm. Both the spatial resolution and the CNR were improved with IR compared to FBP for the images reconstructed from the same number of projections. In comparison with FBP reconstruction, the CNR was improved from 3.4 to 7.5 by using the IR method with six-fold fewer projections while maintaining the same spatial resolution. The study demonstrated that the IR method coupled with DDL could significantly reduce the required number of projections for a CT reconstruction compared to the FBP method while achieving a much better CNR and maintaining the same spatial resolution. From this, the radiation dose and scanning time can potentially be reduced by a factor of approximately 6 by using this IR method for image reconstruction in a CZT-based breast CT system.

Tight-frame based iterative image reconstruction for spectral breast CT

PURPOSE To investigate tight-frame based iterative reconstruction (TFIR) technique for spectral breast computed tomography (CT) using fewer projections while achieving greater image quality. METHODS The experimental data were acquired with a fan-beam breast CT system based on a cadmium zinc telluride photon-counting detector. The images were reconstructed with a varying number of projections using the TFIR and filtered backprojection (FBP) techniques. The image quality between these two techniques was evaluated. The images spatial resolution was evaluated using a high-resolution phantom, and the contrast to noise ratio (CNR) was evaluated using a postmortem breast sample. The postmortem breast samples were decomposed into water, lipid, and protein contents based on images reconstructed from TFIR with 204 projections and FBP with 614 projections. The volumetric fractions of water, lipid, and protein from the image-based measurements in both TFIR and FBP were compared to the chemical analysis. RESULTS The spatial resolution and CNR were comparable for the images reconstructed by TFIR with 204 projections and FBP with 614 projections. Both reconstruction techniques provided accurate quantification of water, lipid, and protein composition of the breast tissue when compared with data from the reference standard chemical analysis. CONCLUSIONS Accurate breast tissue decomposition can be done with three fold fewer projection images by the TFIR technique without any reduction in image spatial resolution and CNR. This can result in a two-third reduction of the patient dose in a multislit and multislice spiral CT system in addition to the reduced scanning time in this system.

Breast composition measurement with a cadmium-zinc-telluride based spectral computed tomography system.

PURPOSE To investigate the feasibility of breast tissue composition in terms of water, lipid, and protein with a cadmium-zinc-telluride (CZT) based computed tomography (CT) system to help better characterize suspicious lesions. METHODS Simulations and experimental studies were performed using a spectral CT system equipped with a CZT-based photon-counting detector with energy resolution. Simulations of the figure-of-merit (FOM), the signal-to-noise ratio (SNR) of the dual energy image with respect to the square root of mean glandular dose (MGD), were performed to find the optimal configuration of the experimental acquisition parameters. A calibration phantom 3.175 cm in diameter was constructed from polyoxymethylene plastic with cylindrical holes that were filled with water and oil. Similarly, sized samples of pure adipose and pure lean bovine tissues were used for the three-material decomposition. Tissue composition results computed from the images were compared to the chemical analysis data of the tissue samples. RESULTS The beam energy was selected to be 100 kVp with a splitting energy of 40 keV. The tissue samples were successfully decomposed into water, lipid, and protein contents. The RMS error of the volumetric percentage for the three-material decomposition, as compared to data from the chemical analysis, was estimated to be approximately 5.7%. CONCLUSIONS The results of this study suggest that the CZT-based photon-counting detector may be employed in the CT system to quantify the water, lipid, and protein mass densities in tissue with a relatively good agreement.

Characterization of energy response for photon-counting detectors using x-ray fluorescence.

PURPOSE To investigate the feasibility of characterizing a Si strip photon-counting detector using x-ray fluorescence. METHODS X-ray fluorescence was generated by using a pencil beam from a tungsten anode x-ray tube with 2 mm Al filtration. Spectra were acquired at 90° from the primary beam direction with an energy-resolved photon-counting detector based on an edge illuminated Si strip detector. The distances from the source to target and the target to detector were approximately 19 and 11 cm, respectively. Four different materials, containing silver (Ag), iodine (I), barium (Ba), and gadolinium (Gd), were placed in small plastic containers with a diameter of approximately 0.7 cm for x-ray fluorescence measurements. Linear regression analysis was performed to derive the gain and offset values for the correlation between the measured fluorescence peak center and the known fluorescence energies. The energy resolutions and charge-sharing fractions were also obtained from analytical fittings of the recorded fluorescence spectra. An analytical model, which employed four parameters that can be determined from the fluorescence calibration, was used to estimate the detector response function. RESULTS Strong fluorescence signals of all four target materials were recorded with the investigated geometry for the Si strip detector. The average gain and offset of all pixels for detector energy calibration were determined to be 6.95 mV/keV and -66.33 mV, respectively. The detectors energy resolution remained at approximately 2.7 keV for low energies, and increased slightly at 45 keV. The average charge-sharing fraction was estimated to be 36% within the investigated energy range of 20-45 keV. The simulated detector output based on the proposed response function agreed well with the experimental measurement. CONCLUSIONS The performance of a spectral imaging system using energy-resolved photon-counting detectors is very dependent on the energy calibration of the detector. The proposed x-ray fluorescence technique offers an accurate and efficient way to calibrate the energy response of a photon-counting detector.

A high-resolution photon-counting breast CT system with tensor-framelet based iterative image reconstruction for radiation dose reduction

Both computer simulations and experimental phantom studies were carried out to investigate the radiation dose reduction with tensor framelet based iterative image reconstruction (TFIR) for a dedicated high-resolution spectral breast computed tomography (CT) based on a silicon strip photon-counting detector. The simulation was performed with a 10 cm-diameter water phantom including three contrast materials (polyethylene, 8 mg ml(-1) iodine and B-100 bone-equivalent plastic). In the experimental study, the data were acquired with a 1.3 cm-diameter polymethylmethacrylate (PMMA) phantom containing iodine in three concentrations (8, 16 and 32 mg ml(-1)) at various radiation doses (1.2, 2.4 and 3.6 mGy) and then CT images were reconstructed using the filtered-back-projection (FBP) technique and the TFIR technique, respectively. The image quality between these two techniques was evaluated by the quantitative analysis on contrast-to-noise ratio (CNR) and spatial resolution that was evaluated using the task-based modulation transfer function (MTF). Both the simulation and experimental results indicated that the task-based MTF obtained from TFIR reconstruction with one-third of the radiation dose was comparable to that from the FBP reconstruction for low contrast target. For high contrast target, the TFIR was substantially superior to the FBP reconstruction in terms of spatial resolution. In addition, TFIR was able to achieve a factor of 1.6-1.8 increase in CNR, depending on the target contrast level. This study demonstrates that the TFIR can reduce the required radiation dose by a factor of two-thirds for a CT image reconstruction compared to the FBP technique. It achieves much better CNR and spatial resolution for high contrast target in addition to retaining similar spatial resolution for low contrast target. This TFIR technique has been implemented with a graphic processing unit system and it takes approximately 10 s to reconstruct a single-slice CT image, which can potentially be used in a future multi-slit multi-slice spiral CT system.

Characteristic performance evaluation of a photon counting Si strip detector for low dose spectral breast CT imaging

PURPOSE The possible clinical applications which can be performed using a newly developed detector depend on the detectors characteristic performance in a number of metrics including the dynamic range, resolution, uniformity, and stability. The authors have evaluated a prototype energy resolved fast photon counting x-ray detector based on a silicon (Si) strip sensor used in an edge-on geometry with an application specific integrated circuit to record the number of x-rays and their energies at high flux and fast frame rates. The investigated detector was integrated with a dedicated breast spectral computed tomography (CT) system to make use of the detectors high spatial and energy resolution and low noise performance under conditions suitable for clinical breast imaging. The aim of this article is to investigate the intrinsic characteristics of the detector, in terms of maximum output count rate, spatial and energy resolution, and noise performance of the imaging system. METHODS The maximum output count rate was obtained with a 50 W x-ray tube with a maximum continuous output of 50 kVp at 1.0 mA. A109Cd source, with a characteristic x-ray peak at 22 keV from Ag, was used to measure the energy resolution of the detector. The axial plane modulation transfer function (MTF) was measured using a 67 μm diameter tungsten wire. The two-dimensional (2D) noise power spectrum (NPS) was measured using flat field images and noise equivalent quanta (NEQ) were calculated using the MTF and NPS results. The image quality parameters were studied as a function of various radiation doses and reconstruction filters. The one-dimensional (1D) NPS was used to investigate the effect of electronic noise elimination by varying the minimum energy threshold. RESULTS A maximum output count rate of 100 million counts per second per square millimeter (cps/mm2) has been obtained (1 million cps per 100×100 μm pixel). The electrical noise floor was less than 4 keV. The energy resolution measured with the 22 keV photons from a 109Cd source was less than 9%. A reduction of image noise was shown in all the spatial frequencies in 1D NPS as a result of the elimination of the electronic noise. The spatial resolution was measured just above 5 line pairs per mm (lp/mm) where 10% of MTF corresponded to 5.4 mm(-1). The 2D NPS and NEQ shows a low noise floor and a linear dependence on dose. The reconstruction filter choice affected both of the MTF and NPS results, but had a weak effect on the NEQ. CONCLUSIONS The prototype energy resolved photon counting Si strip detector can offer superior imaging performance for dedicated breast CT as compared to a conventional energy-integrating detector due to its high output count rate, high spatial and energy resolution, and low noise characteristics, which are essential characteristics for spectral breast CT imaging.

Measurement of breast tissue composition with dual energy cone-beam computed tomography: A postmortem study

PURPOSE To investigate the feasibility of a three-material compositional measurement of water, lipid, and protein content of breast tissue with dual kVp cone-beam computed tomography (CT) for diagnostic purposes. METHODS Simulations were performed on a flat panel-based computed tomography system with a dual kVp technique in order to guide the selection of experimental acquisition parameters. The expected errors induced by using the proposed calibration materials were also estimated by simulation. Twenty pairs of postmortem breast samples were imaged with a flat-panel based dual kVp cone-beam CT system, followed by image-based material decomposition using calibration data obtained from a three-material phantom consisting of water, vegetable oil, and polyoxymethylene plastic. The tissue samples were then chemically decomposed into their respective water, lipid, and protein contents after imaging to allow direct comparison with data from dual energy decomposition. RESULTS Guided by results from simulation, the beam energies for the dual kVp cone-beam CT system were selected to be 50 and 120 kVp with the mean glandular dose divided equally between each exposure. The simulation also suggested that the use of polyoxymethylene as the calibration material for the measurement of pure protein may introduce an error of -11.0%. However, the tissue decomposition experiments, which employed a calibration phantom made out of water, oil, and polyoxymethylene, exhibited strong correlation with data from the chemical analysis. The average root-mean-square percentage error for water, lipid, and protein contents was 3.58% as compared with chemical analysis. CONCLUSIONS The results of this study suggest that the water, lipid, and protein contents can be accurately measured using dual kVp cone-beam CT. The tissue compositional information may improve the sensitivity and specificity for breast cancer diagnosis.PURPOSE To investigate the feasibility of a three-material compositional measurement of water, lipid, and protein content of breast tissue with dual kVp cone-beam computed tomography (CT) for diagnostic purposes. METHODS Simulations were performed on a flat panel-based computed tomography system with a dual kVp technique in order to guide the selection of experimental acquisition parameters. The expected errors induced by using the proposed calibration materials were also estimated by simulation. Twenty pairs of postmortem breast samples were imaged with a flat-panel based dual kVp cone-beam CT system, followed by image-based material decomposition using calibration data obtained from a three-material phantom consisting of water, vegetable oil, and polyoxymethylene plastic. The tissue samples were then chemically decomposed into their respective water, lipid, and protein contents after imaging to allow direct comparison with data from dual energy decomposition. RESULTS Guided by results from simulation, the beam energies for the dual kVp cone-beam CT system were selected to be 50 and 120 kVp with the mean glandular dose divided equally between each exposure. The simulation also suggested that the use of polyoxymethylene as the calibration material for the measurement of pure protein may introduce an error of -11.0%. However, the tissue decomposition experiments, which employed a calibration phantom made out of water, oil, and polyoxymethylene, exhibited strong correlation with data from the chemical analysis. The average root-mean-square percentage error for water, lipid, and protein contents was 3.58% as compared with chemical analysis. CONCLUSIONS The results of this study suggest that the water, lipid, and protein contents can be accurately measured using dual kVp cone-beam CT. The tissue compositional information may improve the sensitivity and specificity for breast cancer diagnosis.

Detecting Cardiovascular Disease from Mammograms With Deep Learning

Coronary artery disease is a major cause of death in women. Breast arterial calcifications (BACs), detected inmammograms, can be useful riskmarkers associated with the disease. We investigate the feasibility of automated and accurate detection ofBACsinmammograms for risk assessment of coronary artery disease. We develop a 12-layer convolutional neural network to discriminate BAC from non-BAC and apply a pixelwise, patch-based procedure for BAC detection. To assess the performance of the system, we conduct a reader study to provide ground-truth information using the consensus of human expert radiologists. We evaluate the performance using a set of 840 full-field digital mammograms from 210 cases, using both free-responsereceiveroperatingcharacteristic (FROC) analysis and calcium mass quantification analysis. The FROC analysis shows that the deep learning approach achieves a level of detection similar to the human experts. The calcium mass quantification analysis shows that the inferred calcium mass is close to the ground truth, with a linear regression between them yielding a coefficient of determination of 96.24%. Taken together, these results suggest that deep learning can be used effectively to develop an automated system for BAC detection inmammograms to help identify and assess patients with cardiovascular risks.