A post-reconstruction method to correct cupping artifacts in cone beam breast computed tomography

Effect of Mandible Phantom Inclination in the Axial Plane on Image Quality in the Presence of Implant Using Cone-Beam Computer Tomography

Purpose To assess the effect of 30° phantom inclination on image quality in the presence of an implant using cone-beam computed tomography (CBCT). Materials and methods Three series of eight scans were taken and categorized by a range of 87-90 kVp and 7.1 mA, and 8 mA. For the first CBCT series, the phantom was placed on a flat plane. For the second series, the phantom was inclined at 30° in the axial plane. For the third series, inclined scans were re-oriented and included for statistics. In total, 24 scans were used for statistics. i.e., eight scans at three different planes (flat plane, inclined plane, and re-oriented inclined plane). All the images were analyzed for artifact and contrast-to-noise ratio (CNR) on ImageJ software. Results The inclination of the dry human mandible phantom by 30° reduces the artifact (p <0.05). However, the CNR was not affected by the phantom inclination. Conclusion The appropriate inclination of the head can significantly reduce the metal artifact in the presence of implants and thus improve the CBCT image quality for post-operative follow-up.

Preliminary exploration of response the course of radiotherapy for stage III non-small cell lung cancer based on longitudinal CT radiomics features

Purpose

Explore the longitudinal CT-based radiomics to demonstrate the changing trend of radiotherapy response and to determine at which point after the onset of treatment radiomics exhibit the greatest change for stage III NSCLC patients.

Methods and materials

Ten stage III NSCLC patients in line with inclusion criteria were enrolled retrospectively, each of whom received radiotherapy or concurrent chemo-radiotherapy and performed eight series of follow-up CT imaging. Longitudinal radiomics were extracted on region of interest from the eight registered images, then two steps were conducted to select significant features as indicators of tumor change: 1) stable features were selected by Kendall rank correlation; 2) texture feature types with a steadily changing trend were retained and intensity features with stable change trends were selected to represent the large number of them. Next, the trend and rate of tumor change were analyzed using the Delta method and Curve-fitting method. Finally, the statistics in the distribution of stable features in patients were calculated.

Results

675 stable features were selected from a total number of 1371 radiomics features, then 12 texture features types were retained and three intensity features were chosen to represent their own category. Among the final selected feature types, it was found that the two time points were weeks 1 and 3 with the higher rate of change. One patient had very few stable tumor features out of a total of 101 features, and the rate of change of features of another patient was conspicuously higher than the average level with number of 301 features.

Conclusion

The longitudinal CT radiomics could demonstrate the change trend of tumor and at which point exhibit the greatest change during radiotherapy, and potentially be used for treatment decisions concerning adaptive radiotherapy.

A novel hardware duo of beam modulation and shielding to reduce scatter acquisition and dose in cone-beam breast CT

PURPOSE

In cone-beam breast CT, scattered photons form a large portion of the acquired signal, adversely impacting image quality throughout the frequency response of the imaging system. Prior simulation studies provided proof of concept for utilization of a hardware solution to prevent scatter acquisition. Here, we report the design, implementation, and characterization of an auxiliary apparatus of fluence modulation and scatter shielding that does indeed lead to projections with a reduced level of scatter.

METHODS

An apparatus was designed for permanent installation within an existing cone-beam CT system. The apparatus is composed of two primary assemblies: a "Fluence Modulator" (FM) and a "Scatter Shield" (SS). The design of the assemblies enables them to operate in synchrony during image acquisition, converting the sourced x-rays into a moving narrow beam. During a projection, this narrow beam sweeps the entire fan angle coverage of the imaging system. As the two assemblies are contingent on one another, their joint implementation is described in the singular as apparatus FM-SS. The FM and the SS assemblies are each comprised a metal housing, a sensory system, and a robotic system. A controller unit handles their relative movements. A series of comparative studies were conducted to evaluate the performance of a cone-beam CT system in two "modes" of operation: with and without FM-SS installed, and to compare the results of physical implementation with those previously simulated. The dynamic range requirements of the utilized detector in the cone-beam CT imaging system were first characterized, independent of the mode of operation. We then characterized and compared the spatial resolution of the imaging system with, and without, FM-SS. A physical breast phantom, representative of an average size breast, was developed and imaged. Actual differences in signal level obtained with, versus without, FM-SS were then compared to the expected level gains based on previously reported simulations. Following these initial assessments, the scatter acquisition in each projection in both modes of operation was investigated. Finally, as an initial study of the impact of FM-SS on radiation dose in an average size breast, a series of Monte Carlo simulations were coupled with physical measurements of air kerma, with and without FM-SS.

RESULTS

With implementation of FM-SS, the detector's required dynamic range was reduced by a factor of 5.5. Substantial reduction in the acquisition of the scattered rays, by a factor of 5.1 was achieved. With the implementation of FM-SS, deposited dose was reduced by 27% in the studied breast.

CONCLUSIONS

The disclosed implementation of FM-SS, within a cone-beam breast CT system, results in reduction of scatter-components in acquired projections, reduction of dose deposit to the breast, and relaxation of requirements for the detector's dynamic range. Controlling or correcting for patient motion occurring during image acquisition remains an open problem to be solved prior to practical clinical usage of FM-SS cone-beam breast CT.

Flat-panel conebeam CT in the clinic: history and current state

Research into conebeam CT concepts began as soon as the first clinical single-slice CT scanner was conceived. Early implementations of conebeam CT in the 1980s focused on high-contrast applications where concurrent high resolution ( < 200    μ m ), for visualization of small contrast-filled vessels, bones, or teeth, was an imaging requirement that could not be met by the contemporaneous CT scanners. However, the use of nonlinear imagers, e.g., x-ray image intensifiers, limited the clinical utility of the earliest diagnostic conebeam CT systems. The development of consumer-electronics large-area displays provided a technical foundation that was leveraged in the 1990s to first produce large-area digital x-ray detectors for use in radiography and then compact flat panels suitable for high-resolution and high-frame-rate conebeam CT. In this review, we show the concurrent evolution of digital flat panel (DFP) technology and clinical conebeam CT. We give a brief summary of conebeam CT reconstruction, followed by a brief review of the correction approaches for DFP-specific artifacts. The historical development and current status of flat-panel conebeam CT in four clinical areas-breast, fixed C-arm, image-guided radiation therapy, and extremity/head-is presented. Advances in DFP technology over the past two decades have led to improved visualization of high-contrast, high-resolution clinical tasks, and image quality now approaches the soft-tissue contrast resolution that is the standard in clinical CT. Future technical developments in DFPs will enable an even broader range of clinical applications; research in the arena of flat-panel CT shows no signs of slowing down.

Evaluation of scatter rejection and correction performance of 2D antiscatter grids in cone beam computed tomography

PURPOSE

We have been investigating two-dimensional (2D) antiscatter grids (2D ASGs) to reduce scatter fluence and improve image quality in cone beam computed tomography (CBCT). In this work, two different aspects of 2D ASGs, their scatter rejection and correction capability, were investigated in CBCT experiments. To correct residual scatter transmitted through the 2D ASG, it was used as a scatter measurement device with a novel method: grid-based scatter sampling.

METHODS

Three focused 2D ASG prototypes with grid ratios of 8, 12, and 16 were developed for linac-mounted offset detector CBCT geometry. In the first phase, 2D ASGs were used as a scatter rejection device, and the effect of grid ratio on CT number accuracy and contrast-to-noise ratio (CNR) evaluated in CBCT images. In the second phase, a grid-based scatter sampling method which exploits the signal modulation characteristics of the 2D ASG's septal shadows to measure and correct residual scatter transmitted through the grid was implemented. To evaluate CT number accuracy, the percent change in CT numbers was measured by changing the phantom from head to pelvis size and configuration.

RESULTS

When 2D ASG was used as a scatter rejection device, CT number accuracy increased and the CT number variation due to change in phantom dimensions was reduced from 23% to 2-6%. A grid ratio of 16 yielded the lowest CT number variation. All three 2D ASGs yielded improvement in CNR, up to a factor of two in pelvis-sized phantoms. When 2D ASG prototypes were used for both scatter rejection and correction, CT number variations were reduced further, to 1.3-2.6%. In comparisons with a clinical CBCT system and a high-performance radiographic ASG, 2D ASG provided higher CT number accuracy under the same imaging conditions.

CONCLUSIONS

When 2D ASG is used solely as a scatter rejection device, substantial improvement in CT number accuracy can be achieved by increasing the grid ratio. Two-dimensional ASGs also provided significant CNR improvement even at lower grid ratios. Two-dimensional ASGs used in conjunction with the grid-based scatter sampling method provided further improvement in CT number accuracy, irrespective of the grid ratio, while preserving 2D ASGs' capacity to improve CNR. The combined effect of scatter rejection and residual scatter correction by 2D ASG may accelerate implementation of new techniques in CBCT that require high quantitative accuracy, such as radiotherapy dose calculation and dual energy CBCT.

Quantitative analysis of metal artefacts of dental implant in CBCT image by correlation analysis to micro-CT: A microstructural study

OBJECTIVES

Quantification of dental implant metal artefacts in CBCT images using correlation analysis of trabecular microstructural parameters from CBCT and micro-CT, and analysis of the effect of varying the angular position of the subject.

METHODS

Polyurethane synthetic bone blocks were first scanned without implants by micro-CT and CBCT. Two dental implants were then placed parallel in the bone blocks and these specimens were scanned by CBCT with different alpha angles. Three volumes of interest (VOI) were set for further analysis. Six microstructural parameters were measured: trabecular thickness (Tb_Th), trabecular spacing (Th_Sp), bone volume per total volume (BV/TV), bone surface per total volume (BS/TV), connectivity density (CD) andfractal dimension (FD). Micro-CT measurements were used as a gold standard for CBCT. Spearman correlation coefficients for each microstructural parameter from CBCT and micro-CT were calculated and compared using Steiger's Z test.

RESULTS

Without the implants, in VOI₁, the Spearman correlation coefficients of Tb_Th, Tb_Sp, BV/TV, BS/TV, CD and FD were 0.599, 0.76, 0.552, 0.566, 0.664 and 0.607, respectively. With the implants, the correlation coefficients decreased sharply in VOI₁. As the alpha angle increased from zero to 90°, the correlation coefficients increased and became significant. Similar results appeared in VOI₂. In contrast, in VOI₃, the correlation coefficient decreased as the alpha angle increased.

CONCLUSIONS

Metal artefacts were successfully quantified using microstructural parameters in terms of the image quality of the CBCT. Changes in alpha angle affected the quality of the CBCT image.

Reducing metal artifacts between implants in cone-beam CT by adjusting angular position of the subject

OBJECTIVES

To reduce inter-implant metal artifacts in cone-beam CT (CBCT) imaging by adjusting angular position of the subject relative to the source-detector plane.

MATERIALS AND METHODS

Two dental implants were placed in a block made of homogeneous dental impression material. Using a custom-made apparatus, the specimen was scanned with a CBCT machine at seven different angles (0°, 15°, 30°, 45°, 60°, 75°, and 90°) along three different spatial axes, yielding 21 experimental groups. Thirteen volumes of interest (VOI) including inter- and peri-implant areas were selected from each axial reconstruction perpendicular to the implants. Gray values (GVs) of each pixel within these VOIs were measured. Mean differences in GV (ΔGV) between the VOIs and control area were calculated and expressed as a percentage. These ΔGVs from different spatial angle were compared and analyzed by Welch's analysis of variance (ANOVA) and linear regression using SPSS 25.0 software.

RESULTS

As alpha angle increased, the ΔGV of the inter-implant area increased from - 62.02% to near-zero while the standard deviation decreased. Welch's ANOVA and linear regression analysis revealed ΔGV increased significantly with alpha angle (p < 0.001, R² = 0.406). VOIs on the extension line of two implants showed similar results. After adjusting the beta and gamma angles, there was no significant change in ΔGV.

CONCLUSIONS

Increasing the alpha angle can reduce metal artifacts in the inter-implant area in CBCT images.

A new strategy for craniospinal axis localization and adaptive dosimetric evaluation using cone beam CT

BACKGROUND AND AIM

Computational complexities encountered in craniospinal irradiation (CSI) have been widely investigated with different planning strategies. However, localization of the entire craniospinal axis (CSA) and evaluation of adaptive treatment plans have traditionally been ignored in CSI treatment. In this study, a new strategy for CSI with comprehensive CSA localization and adaptive plan evaluation has been demonstrated using cone beam CT with extended longitudinal field-of-view (CBCTeLFOV).

MATERIALS AND METHODS

Multi-scan CBCT images were acquired with fixed longitudinal table translations (with 1 cm cone-beam overlap) and then fused into a single DICOM-set using the custom software coded in MatLab™. A novel approach for validation of CBCTeLFOV was demonstrated by combined geometry of Catphan-504 and Catphan-604 phantoms. To simulate actual treatment scenarios, at first, the end-to-end workflow of CSI with VMAT was investigated using an anthropomorphic phantom and then applied for two patients (based on random selection).

RESULTS

The fused CBCTeLFOV images were in excellent agreement with planning CT (pCT). The custom developed software effectively manages spatial misalignments arising out of the uncertainties in treatment/setup geometry. Although the structures mapped from pCT to CBCTeLFOV showed minimal variations, a maximum spatial displacement of up to 1.2 cm (and the mean of 0.8 ± 0.3 cm) was recorded in phantom study. Adaptive plan evaluation of patient paradigms showed the likelihood of under-dosing the craniospinal target.

CONCLUSION

Our protocol serves as a guide for precise localization of entire CSA and to ensure adequate dose to the large and complex targets. It can also be adapted for other complex treatment techniques such as total-marrow-irradiation and total-lymphoid-irradiation.

Correction of Cupping Artifacts in Axial Cone-Beam Computed Tomography Images by Using Image Processing Algorithms

Objectives:

One of the most important problems of cone-beam computed tomography (CBCT) imaging technique is the presence of dense objects, such as implants, amalgam fillings, and metal veneers, which result in beam-hardening artifacts. With an increase in the application of CBCT images and considering the problems in relation to cupping artifacts, some algorithms have been presented to reduce these artifacts. The aim was to present an algorithm to eliminate cupping artifacts from axial and other reconstructed CBCT images.

Materials and Methods:

We used CBCT images of NewTom VG imaging system (Verona, Italy, at Dentistry Faculty, Medical Sciences University, Tabriz, Iran) in which every image has a resolution of 366 × 320 in DICOM format. 50 images of patients with cupping artifacts were selected. Using Sobel edge detector and nonlinear gamma correction coefficient, the difference was calculated between the density of axial images in the main image and the image resulting from nonlinear gamma correction at the exact location of the radiopaque dental materials detected by Sobel. The points at which this density difference was out of a definite limit were treated as image artifacts and were eliminated from the main image by the inpainting method.

Results:

The resultant axial images, for producing reconstructed cross-sectional, panoramic images without cupping artifacts, were imported into NTT viewer V5.6 and utilized.

Conclusions:

With comparison, acquired images observed that the offering algorithm is practical and effective for reducing the cupping artifacts and preserving the quality of the reconstructed images. This algorithm does not need any additional equipment.

Shading artifact correction in breast CT using an interleaved deep learning segmentation and maximum-likelihood polynomial fitting approach

PURPOSE

The purpose of this work was twofold: (a) To provide a robust and accurate method for image segmentation of dedicated breast CT (bCT) volume data sets, and (b) to improve Hounsfield unit (HU) accuracy in bCT by means of a postprocessing method that uses the segmented images to correct for the low-frequency shading artifacts in reconstructed images.

METHODS

A sequential and iterative application of image segmentation and low-order polynomial fitting to bCT volume data sets was used in the interleaved correction (IC) method. Image segmentation was performed through a deep convolutional neural network (CNN) with a modified U-Net architecture. A total of 45 621 coronal bCT images from 111 patient volume data sets were segmented (using a previously published segmentation algorithm) and used for neural network training, validation, and testing. All patient data sets were selected from scans performed on four different prototype breast CT systems. The adipose voxels for each patient volume data set, segmented using the proposed CNN, were then fit to a three-dimensional low-order polynomial. The polynomial fit was subsequently used to correct for the shading artifacts introduced by scatter and beam hardening in a method termed "flat fielding." An interleaved utilization of image segmentation and flat fielding was repeated until a convergence criterion was satisfied. Mathematical and physical phantom studies were conducted to evaluate the dependence of the proposed algorithm on breast size and the distribution of fibroglandular tissue. In addition, a subset of patient scans (not used in the CNN training, testing or validation) were used to investigate the accuracy of the IC method across different scanner designs and beam qualities.

RESULTS

The IC method resulted in an accurate classification of different tissue types with an average Dice similarity coefficient > 95%, precision > 97%, recall > 95%, and F1-score > 96% across all tissue types. The flat fielding correction of bCT images resulted in a significant reduction in either cupping or capping artifacts in both mathematical and physical phantom studies as measured by the integral nonuniformity metric with an average reduction of 71% for cupping and 30% for capping across different phantom sizes, and the Uniformity Index with an average reduction of 53% for cupping and 34% for capping.

CONCLUSION

The validation studies demonstrated that the IC method improves Hounsfield Units (HU) accuracy and effectively corrects for shading artifacts caused by scatter contamination and beam hardening. The postprocessing approach described herein is relevant to the broad scope of bCT devices and does not require any modification in hardware or existing scan protocols. The trained CNN parameters and network architecture are available for interested users.

Correction of severe beam-hardening artifacts via a high-order linearization function using a prior-image-based parameter selection method

PURPOSE

Polychromatic x-rays are used in most computed tomography scanners. In this case, a beam-hardening effect occurs, which degrades the image quality and distorts the shapes of objects in the reconstructed images. When the beam-hardening artifact is not severe, conventional correction methods can reduce the artifact reasonably well. However, highly dense materials, such as iron and titanium, can produce more severe beam-hardening artifacts, which often cannot be corrected by conventional methods. Moreover, when the size of the metal is large, severe darks bands due to photon starvation as well as beam-hardening are generated. The purpose of our study was to develop a new method for correcting severe beam-hardening artifacts and severe dark bands using a high-order polynomial correction function and a prior-image-based linearization method.

METHODS

The initial estimate of an image free of beam-hardening (a prior image) was constructed from the initial reconstruction of the original projection data. Its corresponding beam-hardening-free projection data (a prior projection) were calculated by a projection operator onto the prior image. A new beam-hardening correction function G(praw ) with many high-order terms was effectively determined via a simple minimization process applied to the difference between the original projection data and the prior projection data. Using the determined correction function G(praw ), a corrected linearized sinogram pcorr can be obtained, which became effectively linear for the line integrals of the object. Final beam-hardening corrected images can be reconstructed from the linearized sinogram. The proposed method was evaluated in both simulation and real experimental studies.

RESULTS

All investigated cases in both simulations and real experiments showed that the proposed method effectively removed not only streaks for moderate beam-hardening artifacts but also dark bands for severe beam-hardening artifacts without causing structural and contrast distortion.

CONCLUSIONS

The prior-image-based linearization method exhibited better correction performance than conventional methods. Because the proposed method did not require time-consuming iterative reconstruction processes to obtain the optimal correction function, it can expedite the correction procedure and incorporate more high-order terms in the linearization correction function in comparison to the conventional methods.

Transmission characteristics of a two dimensional antiscatter grid prototype for CBCT

AIM

High fraction of scattered radiation in cone-beam CT (CBCT) imaging degrades CT number accuracy and visualization of low contrast objects. To suppress scatter in CBCT projections, we developed a focused, two-dimensional antiscatter grid (2DASG) prototype. In this work, we report on the primary and scatter transmission characteristics of the 2DASG prototype aimed for linac mounted, offset detector geometry CBCT systems in radiation therapy, and compared its performance to a conventional one-dimensional ASG (1DASG).

METHODS

The 2DASG is an array of through-holes separated by 0.1 mm septa that was fabricated from tungsten using additive manufacturing techniques. Through-holes' focusing geometry was designed for offset detector CBCT in Varian TrueBeam system. Two types of ASGs were evaluated: (a) a conventional 1DASG with a grid ratio of 10, (b) the 2DASG prototype with a grid ratio of 8.2. To assess the scatter suppression performance of both ASGs, Scatter-to-primary ratio (SPR) and scatter transmission fraction (Ts) were measured using the beam stop method. Scatter and primary intensities were modulated by varying the phantom thickness between 10 and 40 cm. Additionally, the effect of air gap and bow tie (BT) filter on SPR and Ts were evaluated. Average primary transmission fraction (T_P ) and pixel specific primary transmission were also measured for both ASGs. To assess the effect of transmission characteristics on projection image signal-to-noise ratio (SNR), SNR improvement factor was calculated. Improvement in contrast to noise ratio (CNR) was demonstrated using a low contrast object.

RESULTS

In comparison to 1DASG, 2DASG reduced SPRs by a factor of 3 to 6 across the range of phantom setups investigated. Ts values for 1D and 2DASGs were in the range of 21 to 29%, and 5 to 14% respectively. 2DASG continued to provide lower SPR and Ts at increased air gap and with BT filter. Tp of 1D and 2DASGs were 70.6% and 84.7% respectively. Due to the septal shadow of the 2DASG, its pixel specific primary transmission values varied between 32.5% and 99.1%. With respect to 1DASG, 2DASG provided up to factor of 1.7 more improvement in SNR across the SPR range investigated. Moreover, 2DASG provided improved visualization of low contrast objects with respect to 1DASG and NOASG setups.

CONCLUSIONS

When compared to a conventional 1DASG, 2DASG prototype provided noticeably lower SPR and Ts values, indicating its superior scatter suppression performance. 2DASG also provided 19% higher average primary transmission that was attributed to the absence of interseptal spacers and optimized grid geometry. Our results indicate that the combined effect of lower scatter and higher primary transmission provided by 2DASG may potentially translate into more accurate CT numbers and improved contrast resolution in CBCT images.

X-ray scatter correction for dedicated cone beam breast CT using a forward-projection model

PURPOSE

The quality of dedicated cone-beam breast CT (CBBCT) imaging is fundamentally limited by x-ray scatter contamination due to the large irradiation volume. In this paper, we propose a scatter correction method for CBBCT using a novel forward-projection model with high correction efficacy and reliability.

METHOD

We first coarsely segment the uncorrected, first-pass, reconstructed CBBCT images into binary-object maps and assign the segmented fibroglandular and adipose tissue with the correct attenuation coefficients based on the mean x-ray energy. The modified CBBCT are treated as the prior images toward scatter correction. Primary signals are first estimated via forward projection on the modified CBBCT. To avoid errors caused by inaccurate segmentation, only sparse samples of estimated primary are selected for scatter estimation. A Fourier-Transform based algorithm, herein referred to as local filtration hereafter, is developed to efficiently estimate the global scatter distribution on the detector. The scatter-corrected images are obtained by removing the estimated scatter distribution from measured projection data.

RESULTS

We evaluate the method performance on six patients with different breast sizes and shapes representing the general population. The results show that the proposed method effectively reduces the image spatial non-uniformity from 8.27 to 1.91% for coronal views and from 6.50 to 3.00% for sagittal views. The contrast-to-deviation ratio is improved by an average factor of 1.41. Comparisons on the image details reveal that the proposed scatter correction successfully preserves fine structures of fibroglandular tissues that are lost in the segmentation process.

CONCLUSION

We propose a highly practical and efficient scatter correction algorithm for CBBCT via a forward-projection model. The method is attractive in clinical CBBCT imaging as it is readily implementable on a clinical system without modifications in current imaging protocols or system hardware.

Population of 224 realistic human subject-based computational breast phantoms

PURPOSE

To create a database of highly realistic and anatomically variable 3D virtual breast phantoms based on dedicated breast computed tomography (bCT) data.

METHODS

A tissue classification and segmentation algorithm was used to create realistic and detailed 3D computational breast phantoms based on 230 + dedicated bCT datasets from normal human subjects. The breast volume was identified using a coarse three-class fuzzy C-means segmentation algorithm which accounted for and removed motion blur at the breast periphery. Noise in the bCT data was reduced through application of a postreconstruction 3D bilateral filter. A 3D adipose nonuniformity (bias field) correction was then applied followed by glandular segmentation using a 3D bias-corrected fuzzy C-means algorithm. Multiple tissue classes were defined including skin, adipose, and several fractional glandular densities. Following segmentation, a skin mask was produced which preserved the interdigitated skin, adipose, and glandular boundaries of the skin interior. Finally, surface modeling was used to produce digital phantoms with methods complementary to the XCAT suite of digital human phantoms.

RESULTS

After rejecting some datasets due to artifacts, 224 virtual breast phantoms were created which emulate the complex breast parenchyma of actual human subjects. The volume breast density (with skin) ranged from 5.5% to 66.3% with a mean value of 25.3% ± 13.2%. Breast volumes ranged from 25.0 to 2099.6 ml with a mean value of 716.3 ± 386.5 ml. Three breast phantoms were selected for imaging with digital compression (using finite element modeling) and simple ray-tracing, and the results show promise in their potential to produce realistic simulated mammograms.

CONCLUSIONS

This work provides a new population of 224 breast phantoms based on in vivo bCT data for imaging research. Compared to previous studies based on only a few prototype cases, this dataset provides a rich source of new cases spanning a wide range of breast types, volumes, densities, and parenchymal patterns.

Library based x-ray scatter correction for dedicated cone beam breast CT

PURPOSE

The image quality of dedicated cone beam breast CT (CBBCT) is limited by substantial scatter contamination, resulting in cupping artifacts and contrast-loss in reconstructed images. Such effects obscure the visibility of soft-tissue lesions and calcifications, which hinders breast cancer detection and diagnosis. In this work, we propose a library-based software approach to suppress scatter on CBBCT images with high efficiency, accuracy, and reliability.

METHODS

The authors precompute a scatter library on simplified breast models with different sizes using the geant4-based Monte Carlo (MC) toolkit. The breast is approximated as a semiellipsoid with homogeneous glandular/adipose tissue mixture. For scatter correction on real clinical data, the authors estimate the breast size from a first-pass breast CT reconstruction and then select the corresponding scatter distribution from the library. The selected scatter distribution from simplified breast models is spatially translated to match the projection data from the clinical scan and is subtracted from the measured projection for effective scatter correction. The method performance was evaluated using 15 sets of patient data, with a wide range of breast sizes representing about 95% of general population. Spatial nonuniformity (SNU) and contrast to signal deviation ratio (CDR) were used as metrics for evaluation.

RESULTS

Since the time-consuming MC simulation for library generation is precomputed, the authors' method efficiently corrects for scatter with minimal processing time. Furthermore, the authors find that a scatter library on a simple breast model with only one input parameter, i.e., the breast diameter, sufficiently guarantees improvements in SNU and CDR. For the 15 clinical datasets, the authors' method reduces the average SNU from 7.14% to 2.47% in coronal views and from 10.14% to 3.02% in sagittal views. On average, the CDR is improved by a factor of 1.49 in coronal views and 2.12 in sagittal views.

CONCLUSIONS

The library-based scatter correction does not require increase in radiation dose or hardware modifications, and it improves over the existing methods on implementation simplicity and computational efficiency. As demonstrated through patient studies, the authors' approach is effective and stable, and is therefore clinically attractive for CBBCT imaging.

A level set method for cupping artifact correction in cone-beam CT

PURPOSE

To reduce cupping artifacts and improve the contrast-to-noise ratio in cone-beam computed tomography (CBCT).

METHODS

A level set method is proposed to reduce cupping artifacts in the reconstructed image of CBCT. The authors derive a local intensity clustering property of the CBCT image and define a local clustering criterion function of the image intensities in a neighborhood of each point. This criterion function defines an energy in terms of the level set functions, which represent a segmentation result and the cupping artifacts. The cupping artifacts are estimated as a result of minimizing this energy.

RESULTS

The cupping artifacts in CBCT are reduced by an average of 90%. The results indicate that the level set-based algorithm is practical and effective for reducing the cupping artifacts and preserving the quality of the reconstructed image.

CONCLUSIONS

The proposed method focuses on the reconstructed image without requiring any additional physical equipment, is easily implemented, and provides cupping correction through a single-scan acquisition. The experimental results demonstrate that the proposed method successfully reduces the cupping artifacts.

Effects on image quality of a 2D antiscatter grid in x-ray digital breast tomosynthesis: Initial experience using the dual modality (x-ray and molecular) breast tomosynthesis scanner

PURPOSE

Radiation scattered from the breast in digital breast tomosynthesis (DBT) causes image degradation, including loss of contrast between cancerous and background tissue. Unlike in 2-dimensional (2D) mammography, an antiscatter grid cannot readily be used in DBT because changing alignment between the tube and detector during the scan would result in unacceptable loss of primary radiation. However, in the dual modality breast tomosynthesis (DMT) scanner, which combines DBT and molecular breast tomosynthesis, the tube and detector rotate around a common axis, thereby maintaining a fixed tube-detector alignment. This C-arm geometry raises the possibility of using a 2D (cellular) focused antiscatter grid. The purpose of this study is to assess change in image quality when using an antiscatter grid in the DBT portion of a DMT scan under conditions of fixed radiation dose.

METHODS

Two 2D focused prototype grids with 80 cm focal length were tested, one stack-laminated from copper (Cu) and one cast from a tungsten-polymer (W-poly). They were reciprocated using a motion scheme designed to maximize transmission of primary x-ray photons. Grid-in and grid-out scatter-to-primary ratios (SPRs) were measured for rectangular blocks of material simulating 30%, 50%, and 70% glandular tissue compositions. For assessment of changes in image quality through the addition of a grid, the Computerized Imaging Reference Systems, Inc., phantom Model 011A containing a set of 1 cm thick blocks simulating a range of glandular/adipose ratios from 0/100 to 100/0 was used. To simulate 6.5 and 8.5 cm thick compressed breasts, 1 cm thick slices of PMMA were added to the Model 011A phantom. DBT images were obtained with and without the grid, with exposure parameters fixed for a given compressed thickness. Signal-difference-to-noise ratios (SDNRs), contrast, and voxel value-based attenuation coefficients (μ) were measured for all blocks from reconstructed phantom images.

RESULTS

For 4, 6, and 8 cm tissue-equivalent block phantom thicknesses, the inclusion of the W-poly grid reduced the SPR by factors of 5, 6, and 5.8, respectively. For the same thicknesses, the copper grid reduced the SPR by factors of 3.9, 4.5, and 4.9. For the 011A phantom, the W-poly grid raised the SDNR of the 70/30 block from 0.8, -0.32, and -0.72 to 0.9, 0.76, and 0.062 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. It raised the SDNR of the 100/0 block from 3.78, 1.95, and 1.0 to 3.79, 3.67, and 3.25 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. Inclusion of the W-poly grid improved the accuracy of image-based μ values for all block compositions. However, smearing of attenuation across slices due to limited angular sampling decreases the sensitivity of voxel values to changing composition compared to theoretical μ values.

CONCLUSIONS

Under conditions of fixed radiation dose to the breast, use of a 2D focused grid increased contrast, SDNR, and accuracy of estimated attenuation for mass-simulating block compositions in all phantom thicknesses tested, with the degree of improvement depending upon material composition. A 2D antiscatter grid can be usefully incorporated in DBT systems that employ fully isocentric tube-detector rotation.

A general method for cupping artifact correction of cone-beam breast computed tomography images

PURPOSE

Cone-beam breast computed tomography (CBBCT), a promising breast cancer diagnostic technique, has been under investigation for the past decade. However, owing to scattered radiation and beam hardening, CT numbers are not uniform on CBBCT images. This is known as cupping artifact, and it presents an obstacle for threshold-based volume segmentation. In this study, we proposed a general post-reconstruction method for cupping artifact correction.

METHODS

There were four steps in the proposed method. First, three types of local region histogram peaks were calculated: adipose peaks with low CT numbers, glandular peaks with high CT numbers, and unidentified peaks. Second, a linear discriminant analysis classifier, which was trained by identified adipose and glandular peaks, was employed to identify the unidentified peaks as adipose or glandular peaks. Third, adipose background signal profile was fitted according to the adipose peaks using the least squares method. Finally, the adipose background signal profile was subtracted from original image to obtain cupping corrected image

RESULTS

In experimental study, standard deviation of adipose tissue CT numbers was obviously reduced and the CT numbers were more uniform after cupping correction by proposed method; in simulation study, root-mean-square errors were significantly reduced for both symmetric and asymmetric cupping artifacts, indicating that the proposed method was effective to both artifacts.

CONCLUSIONS

A general method without a circularly symmetric assumption was proposed to correct cupping artifacts in CBBCT images for breast. It may be properly applied to images of real patient breasts with natural pendent geometry.

A breast-specific, negligible-dose scatter correction technique for dedicated cone-beam breast CT: a physics-based approach to improve Hounsfield Unit accuracy

The purpose of this research was to develop a method to correct the cupping artifact caused from x-ray scattering and to achieve consistent Hounsfield Unit (HU) values of breast tissues for a dedicated breast CT (bCT) system. The use of a beam passing array (BPA) composed of parallel-holes has been previously proposed for scatter correction in various imaging applications. In this study, we first verified the efficacy and accuracy using BPA to measure the scatter signal on a cone-beam bCT system. A systematic scatter correction approach was then developed by modeling the scatter-to-primary ratio (SPR) in projection images acquired with and without BPA. To quantitatively evaluate the improved accuracy of HU values, different breast tissue-equivalent phantoms were scanned and radially averaged HU profiles through reconstructed planes were evaluated. The dependency of the correction method on object size and number of projections was studied. A simplified application of the proposed method on five clinical patient scans was performed to demonstrate efficacy. For the typical 10-18 cm breast diameters seen in the bCT application, the proposed method can effectively correct for the cupping artifact and reduce the variation of HU values of breast equivalent material from 150 to 40 HU. The measured HU values of 100% glandular tissue, 50/50 glandular/adipose tissue, and 100% adipose tissue were approximately 46, -35, and -94, respectively. It was found that only six BPA projections were necessary to accurately implement this method, and the additional dose requirement is less than 1% of the exam dose. The proposed method can effectively correct for the cupping artifact caused from x-ray scattering and retain consistent HU values of breast tissues.

Radiotherapy dose calculation on KV cone-beam CT image for lung tumor using the CIRS calibration

On-board kilovoltage (KV) cone-beam computed tomography (CBCT) images are used predominantly for the setup of patients' positioning. The image data can also potentially be used for dose calculation with the precise calibration of Hounsfield units (HU) to electron density (HU-density). CBCT calibration was analyzed in this study. A clinical treatment planning system was employed for CT and KV CBCT image to dose calculations and subsequent comparisons. Two HU-density tables were generated using the Computerized Imaging Reference Systems (CIRS) phantom. The results showed that a maximum ∼4% dose discrepancy was observed for inserts. The single field isodose curves were very close. The lung clinical patient study indicated that the volume of lung tumor that achieved the prescribed dose in CBCT was lower than in the CT plan. Our study showed that the dosimetric accuracy of CBCT-based dose calculation for lung tumor is acceptable only for the purpose of dosimetric checks with calibration applied. KV CBCT images cannot replace traditional CT images for dose calculation accuracy.

Personalized estimates of radiation dose from dedicated breast CT in a diagnostic population and comparison with diagnostic mammography

This study retrospectively analyzed the mean glandular dose (MGD) to 133 breasts from 132 subjects, all women, who participated in a clinical trial evaluating dedicated breast CT in a diagnostic population. The clinical trial was conducted in adherence to a protocol approved by institutional review boards and the study participants provided written informed consent. Individual estimates of MGD to each breast from dedicated breast CT was obtained by combining x-ray beam characteristics with estimates of breast dimensions and fibroglandular fraction from volumetric breast CT images, and using normalized glandular dose coefficients. For each study participant and for the breast corresponding to that imaged with breast CT, an estimate of the MGD from diagnostic mammography (including supplemental views) was obtained from the DICOM image headers for comparison. This estimate uses normalized glandular dose coefficients corresponding to a breast with 50% fibroglandular weight fraction. The median fibroglandular weight fraction for the study cohort determined from volumetric breast CT images was 15%. Hence, the MGD from diagnostic mammography was corrected to be representative of the study cohort. Individualized estimates of MGD from breast CT ranged from 5.7 to 27.8 mGy. Corresponding to the breasts imaged with breast CT, the MGD from diagnostic mammography ranged from 2.6 to 31.6 mGy. The mean (± inter-breast SD) and the median MGD (mGy) from dedicated breast CT exam were 13.9 ± 4.6 and 12.6, respectively. For the corresponding breasts, the mean (± inter-breast SD) and the median MGD (mGy) from diagnostic mammography were 12.4 ± 6.3 and 11.1, respectively. Statistical analysis indicated that at the 0.05 level, the distributions of MGD from dedicated breast CT and diagnostic mammography were significantly different (Wilcoxon signed ranks test, p = 0.007). While the interquartile range and the range (maximum-minimum) of MGD from dedicated breast CT was lower than diagnostic mammography, the median MGD from dedicated breast CT was approximately 13.5% higher than that from diagnostic mammography. The MGD for breast CT is based on a 1.45 mm skin layer and that for diagnostic mammography is based on a 4 mm skin layer; thus, favoring a lower estimate for MGD from diagnostic mammography. The median MGD from dedicated breast CT corresponds to the median MGD from four to five diagnostic mammography views. In comparison, for the same 133 breasts, the mean and the median number of views per breast during diagnostic mammography were 4.53 and 4, respectively. Paired analysis showed that there was approximately equal likelihood of receiving lower MGD from either breast CT or diagnostic mammography. Future work will investigate methods to reduce and optimize radiation dose from dedicated breast CT.

Scaling-law for the energy dependence of anatomic power spectrum in dedicated breast CT

PURPOSE

To determine the x-ray photon energy dependence of the anatomic power spectrum of the breast when imaged with dedicated breast computed tomography (CT).

METHODS

A theoretical framework for scaling the empirically determined anatomic power spectrum at one x-ray photon energy to that at any given x-ray photon energy when imaged with dedicated breast CT was developed. Theory predicted that when the anatomic power spectrum is fitted with a power curve of the form k f(-β), where k and β are fit coefficients and f is spatial frequency, the exponent β would be independent of x-ray photon energy (E), and the amplitude k scales with the square of the difference in energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues. Twenty mastectomy specimens based numerical phantoms that were previously imaged with a benchtop flat-panel cone-beam CT system were converted to 3D distribution of glandular weight fraction (f(g)) and were used to verify the theoretical findings. The 3D power spectrum was computed in terms of f(g) and after converting to linear attenuation coefficients at monoenergetic x-ray photon energies of 20-80 keV in 5 keV intervals. The 1D power spectra along the axes were extracted and fitted with a power curve of the form k f(-β). The energy dependence of k and β were analyzed.

RESULTS

For the 20 mastectomy specimen based numerical phantoms used in the study, the exponent β was found to be in the range of 2.34-2.42, depending on the axis of measurement. Numerical simulations agreed with the theoretical predictions that for a power-law anatomic spectrum of the form k f(-β), β was independent of E and k(E) = k(1)[μ(g)(E) - μ(a)(E)](2), where k(1) is a constant, and μ(g)(E) and μ(a)(E) represent the energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues, respectively.

CONCLUSIONS

Numerical simulations confirmed the theoretical predictions that in dedicated breast CT, the spatial frequency dependence of the anatomic power spectrum will be independent of x-ray photon energy, and the amplitude of the anatomic power spectrum scales by the square of difference in linear attenuation coefficients of fibroglandular and adipose tissues.

Generation of voxelized breast phantoms from surgical mastectomy specimens

PURPOSE

In the research and development of dedicated tomographic breast imaging systems, digital breast object models, also known as digital phantoms, are useful tools. While various digital breast phantoms do exist, the purpose of this study was to develop a realistic high-resolution model suitable for simulating three-dimensional (3D) breast imaging modalities. The primary goal was to design a model capable of producing simulations with realistic breast tissue structure.

METHODS

The methodology for generating an ensemble of digital breast phantoms was based on imaging surgical mastectomy specimens using a benchtop, cone-beam computed tomography system. This approach allowed low-noise, high-resolution projection views of the mastectomy specimens at each angular position. Reconstructions of these projection sets were processed using correction techniques and diffusion filtering prior to segmentation into breast tissue types in order to generate phantoms.

RESULTS

Eight compressed digital phantoms and 20 uncompressed phantoms from which an additional 96 pseudocompressed digital phantoms with voxel dimensions of 0.2 mm(3) were generated. Two distinct tissue classification models were used in forming breast phantoms. The binary model classified each tissue voxel as either adipose or fibroglandular. A multivalue scaled model classified each tissue voxel as percentage of adipose tissue (range 1%-99%). Power spectral analysis was performed to compare simulated reconstructions using the breast phantoms to the original breast specimen reconstruction, and fits were observed to be similar.

CONCLUSIONS

The digital breast phantoms developed herein provide a high-resolution anthropomorphic model of the 3D uncompressed and compressed breast that are suitable for use in evaluating and optimizing tomographic breast imaging modalities. The authors believe that other research groups might find the phantoms useful, and therefore they offer to make them available for wider use.

Cupping artifact correction and automated classification for high-resolution dedicated breast CT images

PURPOSE

To develop and test an automated algorithm to classify the different tissues present in dedicated breast CT images.

METHODS

The original CT images are first corrected to overcome cupping artifacts, and then a multiscale bilateral filter is used to reduce noise while keeping edge information on the images. As skin and glandular tissues have similar CT values on breast CT images, morphologic processing is used to identify the skin mask based on its position information. A modified fuzzy C-means (FCM) classification method is then used to classify breast tissue as fat and glandular tissue. By combining the results of the skin mask with the FCM, the breast tissue is classified as skin, fat, and glandular tissue. To evaluate the authors' classification method, the authors use Dice overlap ratios to compare the results of the automated classification to those obtained by manual segmentation on eight patient images.

RESULTS

The correction method was able to correct the cupping artifacts and improve the quality of the breast CT images. For glandular tissue, the overlap ratios between the authors' automatic classification and manual segmentation were 91.6% ± 2.0%.

CONCLUSIONS

A cupping artifact correction method and an automatic classification method were applied and evaluated for high-resolution dedicated breast CT images. Breast tissue classification can provide quantitative measurements regarding breast composition, density, and tissue distribution.

X-ray scatter correction method for dedicated breast computed tomography

PURPOSE

To improve image quality and accuracy in dedicated breast computed tomography (BCT) by removing the x-ray scatter signal included in the BCT projections.

METHODS

The previously characterized magnitude and distribution of x-ray scatter in BCT results in both cupping artifacts and reduction of contrast and accuracy in the reconstructions. In this study, an image processing method is proposed that estimates and subtracts the low-frequency x-ray scatter signal included in each BCT projection postacquisition and prereconstruction. The estimation of this signal is performed using simple additional hardware, one additional BCT projection acquisition with negligible radiation dose, and simple image processing software algorithms. The high frequency quantum noise due to the scatter signal is reduced using a noise filter postreconstruction. The dosimetric consequences and validity of the assumptions of this algorithm were determined using Monte Carlo simulations. The feasibility of this method was determined by imaging a breast phantom on a BCT clinical prototype and comparing the corrected reconstructions to the unprocessed reconstructions and to reconstructions obtained from fan-beam acquisitions as a reference standard. One-dimensional profiles of the reconstructions and objective image quality metrics were used to determine the impact of the algorithm.

RESULTS

The proposed additional acquisition results in negligible additional radiation dose to the imaged breast (∼0.4% of the standard BCT acquisition). The processed phantom reconstruction showed substantially reduced cupping artifacts, increased contrast between adipose and glandular tissue equivalents, higher voxel value accuracy, and no discernible blurring of high frequency features.

CONCLUSIONS

The proposed scatter correction method for dedicated breast CT is feasible and can result in highly improved image quality. Further optimization and testing, especially with patient images, is necessary to characterize its impact on clinical performance.

An analysis of the mechanical parameters used for finite element compression of a high-resolution 3D breast phantom

PURPOSE

The authors previously introduced a methodology to generate a realistic three-dimensional (3D), high-resolution, computer-simulated breast phantom based on empirical data. One of the key components of such a phantom is that it provides a means to produce a realistic simulation of clinical breast compression. In the current study, they have evaluated a finite element (FE) model of compression and have demonstrated the effect of a variety of mechanical properties on the model using a dense mesh generated from empirical breast data. While several groups have demonstrated an effective compression simulation with lower density finite element meshes, the presented study offers a mesh density that is able to model the morphology of the inner breast structures more realistically than lower density meshes. This approach may prove beneficial for multimodality breast imaging research, since it provides a high level of anatomical detail throughout the simulation study.

METHODS

In this paper, the authors describe methods to improve the high-resolution performance of a FE compression model. In order to create the compressible breast phantom, dedicated breast CT data was segmented and a mesh was generated with 4-noded tetrahedral elements. Using an explicit FE solver to simulate breast compression, several properties were analyzed to evaluate their effect on the compression model including: mesh density, element type, density, and stiffness of various tissue types, friction between the skin and the compression plates, and breast density. Following compression, a simulated projection was generated to demonstrate the ability of the compressible breast phantom to produce realistic simulated mammographic images.

RESULTS

Small alterations in the properties of the breast model can change the final distribution of the tissue under compression by more than 1 cm; which ultimately results in different representations of the breast model in the simulated images. The model properties that impact displacement the most are mesh density, friction between the skin and the plates, and the relative stiffness of the different tissue types.

CONCLUSIONS

The authors have developed a 3D, FE breast model that can yield high spatial resolution breast deformations under uniaxial compression for imaging research purposes and demonstrated that small changes in the mechanical properties can affect images generated using the phantom.

A general framework and review of scatter correction methods in x-ray cone-beam computerized tomography. Part 1: Scatter compensation approaches

Since scattered radiation in cone-beam volume CT implies severe degradation of CT images by quantification errors, artifacts, and noise increase, scatter suppression is one of the main issues related to image quality in CBCT imaging. The aim of this review is to structurize the variety of scatter suppression methods, to analyze the common structure, and to develop a general framework for scatter correction procedures. In general, scatter suppression combines hardware techniques of scatter rejection and software methods of scatter correction. The authors emphasize that scatter correction procedures consist of the main components scatter estimation (by measurement or mathematical modeling) and scatter compensation (deterministic or statistical methods). The framework comprises most scatter correction approaches and its validity also goes beyond transmission CT. Before the advent of cone-beam CT, a lot of papers on scatter correction approaches in x-ray radiography, mammography, emission tomography, and in Megavolt CT had been published. The opportunity to avail from research in those other fields of medical imaging has not yet been sufficiently exploited. Therefore additional references are included when ever it seems pertinent. Scatter estimation and scatter compensation are typically intertwined in iterative procedures. It makes sense to recognize iterative approaches in the light of the concept of self-consistency. The importance of incorporating scatter compensation approaches into a statistical framework for noise minimization has to be underscored. Signal and noise propagation analysis is presented. A main result is the preservation of differential-signal-to-noise-ratio (dSNR) in CT projection data by ideal scatter correction. The objective of scatter compensation methods is the restoration of quantitative accuracy and a balance between low-contrast restoration and noise reduction. In a synopsis section, the different deterministic and statistical methods are discussed with respect to their properties and applications. The current paper is focused on scatter compensation algorithms. The multitude of scatter estimation models will be dealt with in a separate paper.

Artefacts in CBCT: a review

Artefacts are common in today's cone beam CT (CBCT). They are induced by discrepancies between the mathematical modelling and the actual physical imaging process. Since artefacts may interfere with the diagnostic process performed on CBCT data sets, every user should be aware of their presence. This article aims to discuss the most prominent artefacts identified in the scientific literature and review the existing knowledge on these artefacts. We also briefly review the basic three-dimensional (3D) reconstruction concept applied by today's CBCT scanners, as all artefacts are more or less directly related to it.

Penalized maximum likelihood reconstruction for improved microcalcification detection in breast tomosynthesis

We examined the application of an iterative penalized maximum likelihood (PML) reconstruction method for improved detectability of microcalcifications (MCs) in digital breast tomosynthesis (DBT). Localized receiver operating characteristic (LROC) psychophysical studies with human observers and 2-D image slices were conducted to evaluate the performance of this reconstruction method and to compare its performance against the commonly used Feldkamp FBP algorithm. DBT projections were generated using rigorous computer simulations that included accurate modeling of the noise and detector blur. Acquisition dose levels of 0.7, 1.0, and 1.5 mGy in a 5-cm-thick compressed breast were tested. The defined task was to localize and detect MC clusters consisting of seven MCs. The individual MC diameter was 150 μm. Compressed-breast phantoms derived from CT images of actual mastectomy specimens provided realistic background structures for the detection task. Four observers each read 98 test images for each combination of reconstruction method and acquisition dose. All observers performed better with the PML images than with the FBP images. With the acquisition dose of 0.7 mGy, the average areas under the LROC curve (A(L)) for the PML and FBP algorithms were 0.69 and 0.43, respectively. For the 1.0-mGy dose, the values of A(L) were 0.93 (PML) and 0.7 (FBP), while the 1.5-mGy dose resulted in areas of 1.0 and 0.9, respectively, for the PML and FBP algorithms. A 2-D analysis of variance applied to the individual observer areas showed statistically significant differences (at a significance level of 0.05) between the reconstruction strategies at all three dose levels. There were no significant differences in observer performance for any of the dose levels.

An image quality comparison study between XVI and OBI CBCT systems

The purpose of this study is to evaluate and compare image quality characteristics for two commonly used and commercially available CBCT systems: the X‐ray Volumetric Imager and the On‐Board Imager. A commonly used CATPHAN image quality phantom was used to measure various image quality parameters, namely, pixel value stability and accuracy, noise, contrast to noise ratio (CNR), high‐contrast resolution, low contrast resolution and image uniformity. For the XVI unit, we evaluated the image quality for four manufacturer‐supplied protocols as a function of mAs. For the OBI unit, we did the same for the full‐fan and half‐fan scanning modes, which were respectively used with the full bow‐tie and half bow‐tie filters. For XVI, the mean pixel values of regions of interest were found to generally decrease with increasing mAs for all protocols, while they were relatively stable with mAs for OBI. Noise was slightly lower on XVI and was seen to decrease with increasing mAs, while CNR increased with mAs for both systems. For XVI and OBI, the high‐contrast resolution was approximately limited by the pixel resolution of the reconstructed image. On OBI images, up to 6 and 5 discs of 1% and 0.5% contrast, respectively, were visible for a high mAs setting using the full‐fan mode, while none of the discs were clearly visible on the XVI images for various mAs settings when the medium resolution reconstruction was used. In conclusion, image quality parameters for XVI and OBI have been quantified and compared for clinical protocols under various mAs settings. These results need to be viewed in the context of a recent study that reported the dose‐mAs relationship for the two systems and found that OBI generally delivered higher imaging doses than XVI.⁽¹⁾

PACS numbers: 85.57.C‐, 85.57.cj, 85.57.cm, 85.57.cf

Dosimetric characterization of a dedicated breast computed tomography clinical prototype

PURPOSE

To investigate the glandular dose magnitudes and characteristics resulting from image acquisition using a dedicated breast computed tomography (BCT) clinical prototype imaging system.

METHODS

The x-ray spectrum and output characteristics of a BCT clinical prototype (Koning Corporation, West Henrietta, NY) were determined using empirical measurements, breast phantoms, and an established spectrum model. The geometry of the BCT system was replicated in a Monte Carlo-based computer simulation using the GEANT4 toolkit and was validated by comparing the simulated results for exposure distribution in a standard 16 cm CT head phantom with those empirically determined using a 10 cm CT pencil ionization chamber and dosimeter. The computer simulation was further validated by replicating the results of a previous BCT dosimetry study. Upon validation, the computer simulation was modified to include breasts of varying sizes and homogeneous compositions spanning those encountered clinically, and the normalized mean glandular dose resulting from BCT was determined. Using the system's measured exposure output determined automatically for breasts of different size and density, the mean glandular dose for these breasts was computed and compared to the glandular dose resulting from mammography. Finally, additional Monte Carlo simulations were performed to study how the glandular dose values vary within the breast tissue during acquisition with both this BCT prototype and a typical craniocaudal (CC) mammographic acquisition.

RESULTS

This BCT prototype uses an x-ray spectrum with a first half-value layer of 1.39 mm Al and a mean x-ray energy of 30.3 keV. The normalized mean glandular dose for breasts of varying size and composition during BCT acquisition with this system ranges from 0.278 to 0.582 mGy/mGy air kerma with the reference air kerma measured in air at the center of rotation. Using the measured exposure outputs for the tube currents automatically selected by the system for the breasts of different sizes and densities, the mean glandular dose for a BCT acquisition with this prototype system varies from 5.6 to 17.5 mGy, with the value for a breast of mean size and composition being 17.06 mGy. The glandular dose throughout the breast tissue of this mean breast varies by up to +/- 50% of the mean value. During a typical CC view mammographic acquisition of an equivalent mean breast, which typically results in a mean glandular dose of 2.0-2.5 mGy, the glandular dose throughout the breast tissue varies from approximately 15% to approximately 400% of the mean value.

CONCLUSIONS

Acquisition of a BCT image with the automated tube output settings for a mean breast with the Koning Corp. clinical prototype results in mean glandular dose values approximately equivalent to three to five two-view mammographic examinations for a similar breast. For all breast sizes and compositions studied, this glandular dose ratio between acquisition with this BCT prototype and two-view mammography ranges from 1.4 to 7.2. In mammography, portions of the mean-sized breast receive a considerably higher dose than the mean value for the whole breast. However, only a small portion of a breast undergoing mammography would receive a glandular dose similar to that from BCT.

Reduction in x-ray scatter and radiation dose for volume-of-interest (VOI) cone-beam breast CT--a phantom study

With volume-of-interest (VOI) cone-beam computed tomography (CBCT) imaging, one set of projection images are acquired with the VOI collimator at a regular or high exposure level and the second set of projection images are acquired without the collimator at a reduced exposure level. The high exposure VOI scan data inside the VOI and the low exposure full-field scan data outside the VOI are then combined together to generate composite projection images for image reconstruction. To investigate and quantify scatter reduction, dose saving and image quality improvement in VOI CBCT imaging, a flat panel detector-based bench-top experimental CBCT system was built to measure the dose, the scatter-to-primary ratio (SPR), the image contrast, noise level, the contrast-to-noise ratio (CNR) and the figure of merit (FOM) in the CBCT reconstructed images for two polycarbonate cylinders simulating the small and the large phantoms. The results showed that, compared to the full field CBCT technique, radiation doses for the VOI CBCT technique were reduced by a factor of 1.20 and 1.36 for the small and the large phantoms at the phantom center, respectively, and from 2.7 to 3.0 on the edge of the phantom, respectively. Inside the VOI, the SPRs were substantially reduced by a factor of 6.6 and 10.3 for the small and the large phantoms, the contrast signals were improved by a factor of 1.35 and 1.8, and the noise levels were increased by a factor of 1.27 and 1.6, respectively. As a result, the CNRs were improved by a factor of 1.06 and 1.13 for the small and the large phantoms and the FOM improved by a factor of 1.4 and 1.7, respectively.

Methodology for generating a 3D computerized breast phantom from empirical data

The initial process for creating a flexible three-dimensional computer-generated breast phantom based on empirical data is described. Dedicated breast computed-tomography data were processed to suppress noise and scatter artifacts in the reconstructed image set. An automated algorithm was developed to classify the breast into its primary components. A preliminary phantom defined using subdivision surfaces was generated from the segmented data. To demonstrate potential applications of the phantom, simulated mammographic image data were acquired of the phantom using a simplistic compression model and an analytic projection algorithm directly on the surface model. The simulated image was generated using a model for a polyenergetic cone-beam projection of the compressed phantom. The methods used to create the breast phantom generate resulting images that have a high level of tissue structure detail available and appear similar to actual mammograms. Fractal dimension measurements of simulated images of the phantom are comparatively similar to measurements from images of real human subjects. A realistic and geometrically defined breast phantom that can accurately simulate imaging data may have many applications in breast imaging research.

Dose calculation on kV cone beam CT images: an investigation of the Hu-density conversion stability and dose accuracy using the site-specific calibration

Precise calibration of Hounsfield units (HU) to electron density (HU-density) is essential to dose calculation. On-board kV cone beam computed tomography (CBCT) imaging is used predominantly for patients' positioning, but will potentially be used for dose calculation. The impacts of varying 3 imaging parameters (mAs, source-imager distance [SID], and cone angle) and phantom size on the HU number accuracy and HU-density calibrations for CBCT imaging were studied. We proposed a site-specific calibration method to achieve higher accuracy in CBCT image-based dose calculation. Three configurations of the Computerized Imaging Reference Systems (CIRS) water equivalent electron density phantom were used to simulate sites including head, lungs, and lower body (abdomen/pelvis). The planning computed tomography (CT) scan was used as the baseline for comparisons. CBCT scans of these phantom configurations were performed using Varian Trilogy system in a precalibrated mode with fixed tube voltage (125 kVp), but varied mAs, SID, and cone angle. An HU-density curve was generated and evaluated for each set of scan parameters. Three HU-density tables generated using different phantom configurations with the same imaging parameter settings were selected for dose calculation on CBCT images for an accuracy comparison. Changing mAs or SID had small impact on HU numbers. For adipose tissue, the HU discrepancy from the baseline was 20 HU in a small phantom, but 5 times lager in a large phantom. Yet, reducing the cone angle significantly decreases the HU discrepancy. The HU-density table was also affected accordingly. By performing dose comparison between CT and CBCT image-based plans, results showed that using the site-specific HU-density tables to calibrate CBCT images of different sites improves the dose accuracy to approximately 2%. Our phantom study showed that CBCT imaging can be a feasible option for dose computation in adaptive radiotherapy approach if the site-specific calibration is applied.

Breast cancer imaging: a perspective for the next decade

Breast imaging is largely indicated for detection, diagnosis, and clinical management of breast cancer and for evaluation of the integrity of breast implants. In this work, a prospective view of techniques for breast cancer detection and diagnosis is provided based on an assessment of current trends. The potential role of emerging techniques that are under various stages of research and development is also addressed. It appears that the primary imaging tool for breast cancer screening in the next decade will be high-resolution, high-contrast, anatomical x-ray imaging with or without depth information. MRI and ultrasonography will have an increasingly important adjunctive role for imaging high-risk patients and women with dense breasts. Pilot studies with dedicated breast CT have demonstrated high-resolution three-dimensional imaging capabilities, but several technological barriers must be overcome before clinical adoption. Radionuclide based imaging techniques and x-ray imaging with intravenously injected contrast offer substantial potential as a diagnostic tools and for evaluation of suspicious lesions. Developing optical and electromagnetic imaging techniques hold significant potential for physiologic information and they are likely to be of most value when integrated with or adjunctively used with techniques that provide anatomic information. Experimental studies with breast specimens suggest that phase-sensitive x-ray imaging techniques can provide edge enhancement and contrast improvement but more research is needed to evaluate their potential role in clinical breast imaging. From the technological perspective, in addition to improvements within each modality, there is likely to be a trend towards multi-modality systems that combine anatomic with physiologic information. We are also likely to transition from a standardized screening, where all women undergo the same imaging exam (mammography), to selection of a screening modality or modalities based an individual-risk or other classification.