Importance of point-by-point back projection correction for isocentric motion in digital breast tomosynthesis: relevance to morphology of structures such as microcalcifications

Digital breast tomosynthesis is a three-dimensional imaging technique that provides an arbitrary set of reconstruction planes in the breast from a limited-angle series of projection images acquired while the x-ray tube moves. Traditional shift-and-add (SAA) tomosynthesis reconstruction is a common mathematical method to line up each projection image based on its shifting amount to generate reconstruction slices. With parallel-path geometry of tube motion, the path of the tube lies in a plane parallel to the plane of the detector. The traditional SAA algorithm gives shift amounts for each projection image calculated only along the direction of x-ray tube movement. However, with the partial isocentric motion of the x-ray tube in breast tomosynthesis, small objects such as microcalcifications appear blurred (for instance, about 1-4 pixels in blur for a microcalcification in a human breast) in traditional SAA images in the direction perpendicular to the direction of tube motion. Some digital breast tomosynthesis algorithms reported in the literature utilize a traditional one-dimensional SAA method that is not wholly suitable for isocentric motion. In this paper, a point-by-point back projection (BP) method is described and compared with traditional SAA for the important clinical task of evaluating morphology of small objects such as microcalcifications. Impulse responses at different three-dimensional locations with five different combinations of imaging acquisition parameters were investigated. Reconstruction images of microcalcifications in a human subject were also evaluated. Results showed that with traditional SAA and 45 degrees view angle of tube movement with respect to the detector, at the same height above the detector, the in-plane blur artifacts were obvious for objects farther away from x-ray source. In a human subject, the appearance of calcifications was blurred in the direction orthogonal to the tube motion with traditional SAA. With point-by-point BP, the appearance of calcifications was sharper. The point-by-point BP method demonstrated improved rendition of microcalcifications in the direction perpendicular to the tube motion direction. With wide angles or for imaging of larger breasts, this point-by-point BP rather than the traditional SAA should also be considered as the basis of further deblurring algorithms that work in conjunction with the BP method.     

A computer simulation platform for the optimization of a breast tomosynthesis system

In breast tomosynthesis there is a compromise between resolution, noise, and acquisition speed for a given glandular dose. The purpose of the present work is to develop a simulation platform to investigate the potential imaging performance for the many possible tomosynthesis system configurations. The simulation platform was used to investigate the dependence of image blur and signal difference to noise ratio (SDNR) for several different tomosynthesis acquisition configurations. Simulated projections of a slanted thin tungsten wire placed in different object planes were modified according to the detector's modulation transfer function (MTF), with or without pixel binning. In addition, the focal spot blur (FSB), which depends on the location of the wire, the system geometry, the source-detector movement speed, and the exposure time, was also incorporated into the projections. Both expectation maximization (EM) and filtered back projection (FBP) based algorithms were used for 3D image reconstruction. The in-plane MTF was calculated from the reconstructed image of the wire. To evaluate the noise performance, simulated noiseless projections of calcification and tumor in uniform breast tissue were modified with the noise power spectrum (NPS) calculated from a cascaded linear system model for the detector for a given x-ray dose. The SDNR of the reconstructed images was calculated with different tomosynthesis configurations, e.g., pixel binning, view number, and angular range. Our results showed that for a source-to-imager distance (SID) of 66 cm, pixel binning (2 x 2) caused more degradation to the in-plane MTF than the blur caused by the moving focal spot and reconstruction. The in-depth resolution can be improved by increasing the angular range.     

Thick Slices from Tomosynthesis Data Sets: Phantom Study for the Evaluation of Different Algorithms

PURPOSE  

Tomosynthesis is a 3-dimensional mammography technique that generates thin slices separated one to the other by typically 1 mm from source data sets. The relatively high image noise in these thin slices raises the value of 1-cm thick slices computed from the set of reconstructed slices for image interpretation. In an initial evaluation, we investigated the potential of different algorithms for generating thick slices from tomosynthesis source data (maximum intensity projection—MIP; average algorithm—AV, and image generation by means of a new algorithm, so-called softMip). The three postprocessing techniques were evaluated using a homogeneous phantom with one textured slab with a total thickness of about 5 cm in which two 0.5-cm-thick slabs contained objects to simulate microcalcifications, spiculated masses, and round masses. The phantom was examined by tomosynthesis (GE Healthcare). Microcalcifications were simulated by inclusion of calcium particles of four different sizes. The slabs containing the inclusions were examined in two different configurations: adjacent to each other and close to the detector and with the two slabs separated by two 1-cm thick breast equivalent material slabs. The reconstructed tomosynthesis slices were postprocessed using MIP, AV, and softMip to generate 1-cm thick slices with a lower noise level. The three postprocessing algorithms were assessed by calculating the resulting contrast versus background for the simulated microcalcifications and contrast-to-noise ratios (CNR) for the other objects. The CNRs of the simulated round and spiculated masses were most favorable for the thick slices generated with the average algorithm, followed by softMip and MIP. Contrast of the simulated microcalcifications was best for MIP, followed by softMip and average projections. Our results suggest that the additional generation of thick slices may improve the visualization of objects in tomosynthesis. This improvement differs from the different algorithms for microcalcifications, speculated objects, and round masses. SoftMip is a new approach combining features of MIP and average showing image properties in between MIP and AV.     

Noise and signal decoupling in maximum-likelihood reconstructions and Metz filters for PET brain images

Images reconstructed with the maximum-likelihood-by-expectation-maximization (ML) algorithm have lower noise in some regions, particularly low count areas, compared with images reconstructed with filtered backprojection (FBP). The use of statistically correct noise model coupled with the positivity constraint in the ML algorithm provides this noise improvement, but whether this model confers a general advantage for ML over FBP with no noise model and any reconstruction filter, is unclear. We have studied the quantitative impact of the correct noise model in the ML algorithm applied to simulated and real PET fluoro-deoxyglucose (FDG) brain images, given a simplified but accurate reconstruction model with spatially invariant resolution. For FBP reconstruction, several Metz filters were chosen and images with different resolution were obtained depending on the order (1-400) of the Metz filters. Comparisons were made based on the mean Fourier spectra of the projection amplitudes, the noise-power spectra, and the mean region-of-interest signal and noise behaviour in the images. For images with resolution recovery beyond the intrinsic detector resolution, the noise increased significantly for FBP compared with ML. This indicates that in the process of signal recovery using ML, the noise is decoupled from the signal. Such noise decoupling is not possible for FBP. However, for image resolution equivalent to or less than the intrinsic detector resolution, FBP with Metz filters of various orders can achieve a performance similar to ML. The significance of the noise decoupling advantage in ML is dependent on the reconstructed image resolution required for specific imaging tasks.     

Contrast Detail Phantom Comparison on a Commercially Available Unit. Digital Breast Tomosynthesis (DBT) versus Full-Field Digital Mammography (FFDM)

The performance of a commercial digital mammographic system working in 2D planar versus tomosynthesis mode was evaluated in terms of the image signal difference to noise ratio (SDNR). A contrast detail phantom was obtained embedding 1 cm Plexiglas, including 49 holes of different diameter and depth, between two layers containing a breast-simulating material. The phantom was exposed with the details plane perpendicular to the X-ray beam using the manufacturer’s standard clinical breast acquisition parameters. SDNR in the digital breast tomosynthesis (DBT) images was higher than that of the full-field digital mammography (FFDM) for 38 out of 49 details in complex background conditions. These differences (p < 0.05) are statistically significant for 19 details out of 38. The relative SDNR results for DBT and FFDM images showed a dependence on the diameter of the details considered. This paper proposes an initial framework for a global image quality evaluation for commercial systems that can operate with different image acquisition modality using the same detector.     

Evaluation of Back Projection Methods for Breast Tomosynthesis Image Reconstruction

Breast cancer is the most common cancer among women in the USA. Compared to mammography, digital breast tomosynthesis is a new imaging technique that may improve the diagnostic accuracy by removing the ambiguities of overlapped tissues and providing 3D information of the breast. Tomosynthesis reconstruction algorithms generate 3D reconstructed slices from a few limited angle projection images. Among different reconstruction algorithms, back projection (BP) is considered an important foundation of quite a few reconstruction techniques with deblurring algorithms such as filtered back projection. In this paper, two BP variants, including α-trimmed BP and principal component analysis-based BP, were proposed to improve the image quality against that of traditional BP. Computer simulations and phantom studies demonstrated that the α-trimmed BP may improve signal response performance and suppress noise in breast tomosynthesis image reconstruction.     

2D linear and iterative reconstruction algorithms for a PET-insert scanner

We are developing novel insert devices for existing whole body PET scanners to achieve better resolution in selected regions of interest such as the head, neck, breast or abdomen. The insert considered here is a full ring of high resolution detectors, which can be placed around the object of interest. Adding the insert inside the scanner leads to three different types of coincidences: insert-insert, insert-scanner and scanner-scanner. The insert-insert and scanner-scanner coincidences are similar to the coincidences obtained in a traditional PET system. The insert-scanner coincidences have an inherent fan-beam geometry for which a spatially variant system matrix is proposed. The system matrix is computed using the intersection of a fan beam with a pixel. A filtered back-projection (FBP) algorithm for the insert-scanner geometry is presented. This FBP algorithm yields images with significantly reduced artifacts compared to FBP reconstructions on insert-scanner data rebinned into parallel beams. This is demonstrated using simulated point source data acquired using SimSET. It is proposed to use a penalized ML-EM (PML-EM) algorithm using a log-cosh roughness penalty function to reconstruct a single activity distribution from all three data sets. This is demonstrated qualitatively using simulated point source data. A quantitative comparison of PML-EM and FBP was performed on data acquired from insert-scanner coincidences using a phantom with hot and cold tumors imaged in an experimental setup. The quantitative studies demonstrate that the resolution/noise tradeoff of PML-EM is improved relative to FBP.     

A comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis

Digital tomosynthesis mammography (DTM) is a promising new modality for breast cancer detection. In DTM, projection-view images are acquired at a limited number of angles over a limited angular range and the imaged volume is reconstructed from the two-dimensional projections, thus providing three-dimensional structural information of the breast tissue. In this work, we investigated three representative reconstruction methods for this limited-angle cone-beam tomographic problem, including the backprojection (BP) method, the simultaneous algebraic reconstruction technique (SART) and the maximum likelihood method with the convex algorithm (ML-convex). The SART and ML-convex methods were both initialized with BP results to achieve efficient reconstruction. A second generation GE prototype tomosynthesis mammography system with a stationary digital detector was used for image acquisition. Projection-view images were acquired from 21 angles in 3 degrees increments over a +/- 30 degrees angular range. We used an American College of Radiology phantom and designed three additional phantoms to evaluate the image quality and reconstruction artifacts. In addition to visual comparison of the reconstructed images of different phantom sets, we employed the contrast-to-noise ratio (CNR), a line profile of features, an artifact spread function (ASF), a relative noise power spectrum (NPS), and a line object spread function (LOSF) to quantitatively evaluate the reconstruction results. It was found that for the phantoms with homogeneous background, the BP method resulted in less noisy tomosynthesized images and higher CNR values for masses than the SART and ML-convex methods. However, the two iterative methods provided greater contrast enhancement for both masses and calcification, sharper LOSF, and reduced interplane blurring and artifacts with better ASF behaviors for masses. For a contrast-detail phantom with heterogeneous tissue-mimicking background, the BP method had strong blurring artifacts along the x-ray source motion direction that obscured the contrast-detail objects, while the other two methods can remove the superimposed breast structures and significantly improve object conspicuity. With a properly selected relaxation parameter, the SART method with one iteration can provide tomosynthesized images comparable to those obtained from the ML-convex method with seven iterations, when BP results were used as initialization for both methods.     

Automated detection of microcalcification clusters for digital breast tomosynthesis using projection data only: a preliminary study

Digital breast tomosynthesis (DBT) is a promising modality for breast imaging in which an anisotropic volume image of the breast is obtained. We present an algorithm for computerized detection of microcalcification clusters (MCCs) for DBT. This algorithm operates on the projection views only. Therefore it does not depend on reconstruction, and is computationally efficient. The algorithm was developed using a database of 30 image sets with microcalcifications, and a control group of 30 image sets without visible findings. The patient data were acquired on the first DBT prototype at Massachusetts General Hospital. Algorithm sensitivity was estimated to be 0.86 at 1.3 false positive clusters, which is below that of current MCC detection algorithms for full-field digital mammography. Because of the small number of patient cases, algorithm parameters were not optimized and one linear classifier was used. An actual limitation of our approach may be that the signal-to-noise ratio in the projection images is too low for microcalcification detection. Furthermore, the database consisted of predominantly small MCC. This may be related to the image quality obtained with this first prototype.     

Eigenvector decomposition of full-spectrum x-ray computed tomography

Energy-discriminated x-ray computed tomography (CT) data were projected onto a set of basis functions to suppress the noise in filtered back-projection (FBP) reconstructions. The x-ray CT data were acquired using a novel x-ray system which incorporated a single-pixel photon-counting x-ray detector to measure the x-ray spectrum for each projection ray. A matrix of the spectral response of different materials was decomposed using eigenvalue decomposition to form the basis functions. Projection of FBP onto basis functions created a de facto image segmentation of multiple contrast agents. Final reconstructions showed significant noise suppression while preserving important energy-axis data. The noise suppression was demonstrated by a marked improvement in the signal-to-noise ratio (SNR) along the energy axis for multiple regions of interest in the reconstructed images. Basis functions used on a more coarsely sampled energy axis still showed an improved SNR. We conclude that the noise-resolution trade off along the energy axis was significantly improved using the eigenvalue decomposition basis functions.     

Scatter radiation in digital tomosynthesis of the breast

Digital tomosynthesis of the breast is being investigated as one possible solution to the problem of tissue superposition present in planar mammography. This imaging technique presents various advantages that would make it a feasible replacement for planar mammography, among them similar, if not lower, radiation glandular dose to the breast; implementation on conventional digital mammography technology via relatively simple modifications; and fast acquisition time. One significant problem that tomosynthesis of the breast must overcome, however, is the reduction of x-ray scatter inclusion in the projection images. In tomosynthesis, due to the projection geometry and radiation dose considerations, the use of an antiscatter grid presents several challenges. Therefore, the use of postacquisition software-based scatter reduction algorithms seems well justified, requiring a comprehensive evaluation of x-ray scatter content in the tomosynthesis projections. This study aims to gain insight into the behavior of x-ray scatter in tomosynthesis by characterizing the scatter point spread functions (PSFs) and the scatter to primary ratio (SPR) maps found in tomosynthesis of the breast. This characterization was performed using Monte Carlo simulations, based on the Geant4 toolkit, that simulate the conditions present in a digital tomosynthesis system, including the simulation of the compressed breast in both the cranio-caudal (CC) and the medio-lateral oblique (MLO) views. The variation of the scatter PSF with varying tomosynthesis projection angle, as well as the effects of varying breast glandular fraction and x-ray spectrum, was analyzed. The behavior of the SPR for different projection angle, breast size, thickness, glandular fraction, and x-ray spectrum was also analyzed, and computer fit equations for the magnitude of the SPR at the center of mass for both the CC and the MLO views were found. Within mammographic energies, the x-ray spectrum was found to have no appreciable effect on the scatter PSF and on the SPR. Glandular fraction and compressed breast size were found to have a small effect, while compressed breast thickness and projection angle, as expected, introduced large variations in both the scatter PSF and SPR. The presence of the breast support plate and the detector cover plate in the simulations introduced important effects on the SPR, which are also relevant to the scatter content in planar mammography.     

乳腺断层摄影合成技术的临床价值

目的:探讨数字乳腺断层合成X线成像(DBT)结合合成2D图像(SM)对乳腺微钙化的检出和诊断效能。方法:回顾性分析228例乳腺影像及病理资料。3名影像医师独立阅读DBT结合全视野数字化乳腺摄影(FFDM)、DBT结合SM、FFDM、SM 4种模式下影像资料,记录微钙化有无,根据BI-RADS 2013版对微钙化进行分类,分析不同密度乳腺类型中良、恶性微钙化的检出率及诊断效能。结果:不管在致密型乳腺或所有腺体类型乳腺中,4种阅片模式对微钙化检出敏感度的差异无统计学意义(P>0.05),特异度均为100%。DBT结合SM与DBT结合FFDM对微钙化诊断敏感度、特异度及ROC曲线下面积的差异无统计学意义(P>0.05);FFDM的敏感度高于SM,特异度低于SM,ROC曲线下面积高于SM,差异均具有统计学意义(P<0.05)。结论:DBT结合SM与DBT结合FFDM对乳腺微钙化的检出、诊断效能相似。     

Megavoltage CT on a tomotherapy system.

A megavoltage computed tomography (MVCT) system was developed on the University of Wisconsin tomotherapy benchtop. This system can operate either axially or helically, and collect transmission data without any bounds on delivered dose. Scan times as low as 12 s per slice are possible, and scans were run with linac output rates of 100 MU min(-1), although the system can be tuned to deliver arbitrarily low dose rates. Images were reconstructed with clinically reasonable doses ranging from 8 to 12 cGy. These images delineate contrasts below 2% and resolutions of 3.0 mm. Thus, the MVCT image quality of this system should be sufficient for verifying the patient's position and anatomy prior to radiotherapy. Additionally, synthetic data were used to test the potential for improved MVCT contrast using maximum-likelihood (ML) reconstruction. Specifically, the maximum-likelihood expectation-maximization (ML-EM) algorithm and a transmission ML algorithm were compared with filtered backprojection (FBP). It was found that for expected clinical MVCT doses enough imaging photons are used such that little benefit is conferred by the improved noise model of ML algorithms. For significantly lower doses, some quantitative improvement is achieved through ML reconstruction. Nonetheless, the image quality at those lower doses is not satisfactory for radiotherapy verification.     

Imaging performance of an amorphous selenium digital mammography detector in a breast tomosynthesis system

In breast tomosynthesis a rapid sequence of N images is acquired when the x-ray tube sweeps through different angular views with respect to the breast. Since the total dose to the breast is kept the same as that in regular mammography, the exposure used for each image of tomosynthesis is 1/N. The low dose and high frame rate pose a tremendous challenge to the imaging performance of digital mammography detectors. The purpose of the present work is to investigate the detector performance in different operational modes designed for tomosynthesis acquisition, e.g., binning or full resolution readout, the range of view angles, and the number of views N. A prototype breast tomosynthesis system with a nominal angular range of +/-25 degrees was used in our investigation. The system was equipped with an amorphous selenium (a-Se) full field digital mammography detector with pixel size of 85 microm. The detector can be read out in full resolution or 2 x 1 binning (binning in the tube travel direction). The focal spot blur due to continuous tube travel was measured for different acquisition geometries, and it was found that pixel binning, instead of focal spot blur, dominates the detector modulation transfer function (MTF). The noise power spectrum (NPS) and detective quantum efficiency (DQE) of the detector were measured with the exposure range of 0.4-6 mR, which is relevant to the low dose used in tomosynthesis. It was found that DQE at 0.4 mR is only 20% less than that at highest exposure for both detector readout modes. The detector temporal performance was categorized as lag and ghosting, both of which were measured as a function of x-ray exposure. The first frame lags were 8% and 4%, respectively, for binning and full resolution mode. Ghosting is negligible and independent of the frame rate. The results showed that the detector performance is x-ray quantum noise limited at the low exposures used in each view of tomosynthesis, and the temporal performance at high frame rate (up to 2 frames per second) is adequate for tomosynthesis.     

A comparison of reconstruction algorithms for C-arm mammography tomosynthesis

Digital tomosynthesis is an imaging technique to produce a tomographic image from a series of angular digital images in a manner similar to conventional focal plane tomography. Unlike film focal plane tomography, the acquisition of the data in a C-arm geometry causes the image receptor to be positioned at various angles to the reconstruction tomogram. The digital nature of the data allows for input images to be combined into the desired plane with the flexibility of generating tomograms of many separate planes from a single set of input data. Angular datasets were obtained of a low contrast detectability (LCD) phantom and cadaver breast utilizing a Lorad stereotactic biopsy unit with a coupled source and digital detector in a C-arm configuration. Datasets of 9 and 41 low-dose projections were collected over a 30 degrees angular range. Tomographic images were reconstructed using a Backprojection (BP) algorithm, an Iterative Subtraction (IS) algorithm that allows the partial subtraction of out-of-focus planes, and an Algebraic Reconstruction (AR) algorithm. These were compared with single view digital radiographs. The methods' effectiveness at enhancing visibility of an obscured LCD phantom was quantified in terms of the Signal to Noise Ratio (SNR), and Signal to Background Ratio (SBR), all normalized to the metric value for the single projection image. The methods' effectiveness at removing ghosting artifacts in a cadaver breast was quantified in terms of the Artifact Spread Function (ASF). The technology proved effective at partially removing out of focus structures and enhancing SNR and SBR. The normalized SNR was highest at 4.85 for the obscured LCD phantom, using nine projections and IS algorithm. The normalized SBR was highest at 23.2 for the obscured LCD phantom, using 41 projections and an AR algorithm. The highest normalized metric values occurred with the obscured phantom. This supports the assertion that the greatest value of tomosynthesis is in imaging fibroglandular breasts. The ASF performance was best with the AR technique and nine projections.     

Dual-energy contrast-enhanced breast tomosynthesis: optimization of beam quality for dose and image quality

Dual-energy contrast-enhanced breast tomosynthesis is a promising technique to obtain three-dimensional functional information from the breast with high resolution and speed. To optimize this new method, this study searched for the beam quality that maximized image quality in terms of mass detection performance. A digital tomosynthesis system was modeled using a fast ray-tracing algorithm, which created simulated projection images by tracking photons through a voxelized anatomical breast phantom containing iodinated lesions. The single-energy images were combined into dual-energy images through a weighted log subtraction process. The weighting factor was optimized to minimize anatomical noise, while the dose distribution was chosen to minimize quantum noise. The dual-energy images were analyzed for the signal difference to noise ratio (SdNR) of iodinated masses. The fast ray-tracing explored 523?776 dual-energy combinations to identify which yields optimum mass SdNR. The ray-tracing results were verified using a Monte Carlo model for a breast tomosynthesis system with a selenium-based flat-panel detector. The projection images from our voxelized breast phantom were obtained at a constant total glandular dose. The projections were combined using weighted log subtraction and reconstructed using commercial reconstruction software. The lesion SdNR was measured in the central reconstructed slice. The SdNR performance varied markedly across the kVp and filtration space. Ray-tracing results indicated that the mass SdNR was maximized with a high-energy tungsten beam at 49 kVp with 92.5 μm of copper filtration and a low-energy tungsten beam at 49 kVp with 95 μm of tin filtration. This result was consistent with Monte Carlo findings. This mammographic technique led to a mass SdNR of 0.92 ± 0.03 in the projections and 3.68 ± 0.19 in the reconstructed slices. These values were markedly higher than those for non-optimized techniques. Our findings indicate that dual-energy breast tomosynthesis can be performed optimally at 49 kVp with alternative copper and tin filters, with reconstruction following weighted subtraction. The optimum technique provides best visibility of iodine against structured breast background in dual-energy contrast-enhanced breast tomosynthesis.     

Evaluation of an Automated Information Extraction Tool for Imaging Data Elements to Populate a Breast Cancer Screening Registry

Breast cancer screening is central to early breast cancer detection. Identifying and monitoring process measures for screening is a focus of the National Cancer Institute’s Population-based Research Optimizing Screening through Personalized Regimens (PROSPR) initiative, which requires participating centers to report structured data across the cancer screening continuum. We evaluate the accuracy of automated information extraction of imaging findings from radiology reports, which are available as unstructured text. We present prevalence estimates of imaging findings for breast imaging received by women who obtained care in a primary care network participating in PROSPR (n = 139,953 radiology reports) and compared automatically extracted data elements to a “gold standard” based on manual review for a validation sample of 941 randomly selected radiology reports, including mammograms, digital breast tomosynthesis, ultrasound, and magnetic resonance imaging (MRI). The prevalence of imaging findings vary by data element and modality (e.g., suspicious calcification noted in 2.6 % of screening mammograms, 12.1 % of diagnostic mammograms, and 9.4 % of tomosynthesis exams). In the validation sample, the accuracy of identifying imaging findings, including suspicious calcifications, masses, and architectural distortion (on mammogram and tomosynthesis); masses, cysts, non-mass enhancement, and enhancing foci (on MRI); and masses and cysts (on ultrasound), range from 0.8 to1.0 for recall, precision, and F-measure. Information extraction tools can be used for accurate documentation of imaging findings as structured data elements from text reports for a variety of breast imaging modalities. These data can be used to populate screening registries to help elucidate more effective breast cancer screening processes.     

三维坐标下任意方向的投影变换及反投影重建方法

提出一种三维坐标下任意方向的投影变换及反投影重建的方法,并且给出一个基于多个二维投影变换及反投影重建的简单实现方式,计算机仿真实验的结果证明了该方法的可行性.本文提出的三维投影变换及反投影重建,不仅可以应用在放疗计划的设计过程中,也可用在多层CT的成像过程、数字减影等多种三维图像重建中.     

Application of boundary detection information in breast tomosynthesis reconstruction

Digital tomosynthesis mammography (DTM) is one of the most promising techniques that can potentially improve early detection of breast cancers. DTM can provide three-dimensional (3D) structural information by reconstructing the whole imaged volume from a sequence of projection-view (PV) mammograms that are acquired at a small number of projection angles over a limited angular range. Our previous study showed that simultaneous algebraic reconstruction technique (SART) can produce satisfactory tomosynthesized image quality compared to maximum likelihood-type algorithms. To improve the efficiency of DTM reconstruction and address the problem of boundary artifacts, we have developed methods to incorporate both two-dimensional (2D) and 3D breast boundary information within the SART reconstruction algorithm in this study. A second generation GE prototype tomosynthesis mammography system with a stationary digital detector was used for PV image acquisition from 21 angles in 3 degrees increments over a +/- 30 degrees angular range. The 2D breast boundary curves on all PV images were obtained by automated segmentation and were used to restrict the SART reconstruction to be performed only within the breast volume. The computation time of SART reconstruction was reduced by 76.3% and 69.9% for cranio-caudal and mediolateral oblique views, respectively, for the chosen example. In addition, a 3D conical trimming method was developed in which the 2D breast boundary curves from all PVs were back projected to generate the 3D breast surface. This 3D surface was then used to eliminate the multiple breast shadows outside the breast volume due to reconstruction by setting these voxels to a constant background value. Our study demonstrates that, by using the 2D and 3D breast boundary information, all breast boundary and most detector boundary artifacts can be effectively removed on all tomosynthesized slices.