Focused beam-stop array for the measurement of scatter in megavoltage portal and cone beam CT imaging

We describe a focused beam-stop array (BSA) for the measurement of object scatter in imaging systems that utilize x-ray beams in the megavoltage (MV) energy range. The BSA consists of 64 doubly truncated tungsten cone elements of 0.5 cm maximum diameter that are arranged in a regular array on an acrylic slab. The BSA is placed in the accessory tray of a medical linear accelerator at a distance of approximately 50 cm from the focal spot. We derive an expression that allows us to estimate the scatter in an image taken without the array present, given image values in a second image with the array in place. The presence of the array reduces fluence incident on the imaged object. This leads to an object-dependent underestimation bias in the scatter measurements. We apply corrections in order to address this issue. We compare estimates of the flat panel detector response to scatter obtained using the BSA to those derived from Monte Carlo simulations. We find that the two estimates agree to within 10% in terms of RMS error for 30 cm x 30 cm water slabs in the thickness range of 10-30 cm. Larger errors in the scatter estimates are encountered for thinner objects, probably owing to extrafocal radiation sources. However, RMS errors in the estimates of primary images are no more than 5% for water slab thicknesses in the range of 1-30 cm. The BSA scatter estimates are also used to correct cone beam tomographic projections. Maximum deviations of central profiles of uniform water phantoms are reduced from 193 to 19 HU after application of corrections for scatter, beam hardening, and lateral truncation that are based on the BSA-derived scatter estimate. The same corrections remove the typical cupping artifact from both phantom and patient images. The BSA proves to be a useful tool for quantifying and removing image scatter, as well as for validating models of MV imaging systems.     

A general framework and review of scatter correction methods in cone beam CT. Part 2: scatter estimation approaches

The main components of scatter correction procedures are scatter estimation and a scatter compensation algorithm. This paper completes a previous paper where a general framework for scatter compensation was presented under the prerequisite that a scatter estimation method is already available. In the current paper, the authors give a systematic review of the variety of scatter estimation approaches. Scatter estimation methods are based on measurements, mathematical-physical models, or combinations of both. For completeness they present an overview of measurement-based methods, but the main topic is the theoretically more demanding models, as analytical, Monte-Carlo, and hybrid models. Further classifications are 3D image-based and 2D projection-based approaches. The authors present a system-theoretic framework, which allows to proceed top-down from a general 3D formulation, by successive approximations, to efficient 2D approaches. A widely useful method is the beam-scatter-kernel superposition approach. Together with the review of standard methods, the authors discuss their limitations and how to take into account the issues of object dependency, spatial variance, deformation of scatter kernels, external and internal absorbers. Open questions for further investigations are indicated. Finally, the authors refer on some special issues and applications, such as bow-tie filter, offset detector, truncated data, and dual-source CT.     

An iterative approach to the beam hardening correction in cone beam CT

In computed tomography (CT), the beam hardening effect has been known to be one of the major sources of deterministic error that leads to inaccuracy and artifact in the reconstructed images. Because of the polychromatic nature of the x-ray source used in CT and the energy-dependent attenuation of most materials, Beer's law no longer holds. As a result, errors are present in the acquired line integrals or measurements of the attenuation coefficients of the scanned object. In the past, many studies have been conducted to combat image artifacts induced by beam hardening. In this paper, we present an iterative beam hardening correction approach for cone beam CT. An algorithm that utilizes a tilted parallel beam geometry is developed and subsequently employed to estimate the projection error and obtain an error estimation image, which is then subtracted from the initial reconstruction. A theoretical analysis is performed to investigate the accuracy of our methods. Phantom and animal experiments are conducted to demonstrate the effectiveness of our approach.     

Tilted plane Feldkamp type reconstruction algorithm for spiral cone beam CT

An approximate image reconstruction method for spiral cone beam computed tomography (CT), called tilted plane Feldkamp type reconstruction algorithm (TPFR), is presented in this paper, which extends Feldkamp cone beam reconstruction algorithm to deal with its inaccuracy and artifact problems caused by large cone angle. This is done by tilting the reconstructing planes to minimize the cone angle and optimally fit the spiral segment of the source. The tilted plane image reconstruction requires reforming the three-dimensional projection data set for the tilted plane and application of Feldkamp algorithm to the reformed data set. Analytical and computational results can show that the image reconstruction performance of the proposed TPFR algorithm is superior to that of the Feldkamp reconstruction algorithm in the image quality, volume coverage speed, maximum achievable pitch value, and slice sensitivity profiles. Moreover, it provides more accurate image reconstruction than the existing two-dimensional reconstruction algorithms.     

Noise simulation in cone beam CT imaging with parallel computing

We developed a computer noise simulation model for cone beam computed tomography imaging using a general purpose PC cluster. This model uses a mono-energetic x-ray approximation and allows us to investigate three primary performance components, specifically quantum noise, detector blurring and additive system noise. A parallel random number generator based on the Weyl sequence was implemented in the noise simulation and a visualization technique was accordingly developed to validate the quality of the parallel random number generator. In our computer simulation model, three-dimensional (3D) phantoms were mathematically modelled and used to create 450 analytical projections, which were then sampled into digital image data. Quantum noise was simulated and added to the analytical projection image data, which were then filtered to incorporate flat panel detector blurring. Additive system noise was generated and added to form the final projection images. The Feldkamp algorithm was implemented and used to reconstruct the 3D images of the phantoms. A 24 dual-Xeon PC cluster was used to compute the projections and reconstructed images in parallel with each CPU processing 10 projection views for a total of 450 views. Based on this computer simulation system, simulated cone beam CT images were generated for various phantoms and technique settings. Noise power spectra for the flat panel x-ray detector and reconstructed images were then computed to characterize the noise properties. As an example among the potential applications of our noise simulation model, we showed that images of low contrast objects can be produced and used for image quality evaluation.     

High-quality 3D correction of ring and radiant artifacts in flat panel detector-based cone beam volume CT imaging

The use of an x-ray flat panel detector is increasingly becoming popular in 3D cone beam volume CT machines. Due to the deficient semiconductor array manufacturing process, the cone beam projection data are often corrupted by different types of abnormalities, which cause severe ring and radiant artifacts in a cone beam reconstruction image, and as a result, the diagnostic image quality is degraded. In this paper, a novel technique is presented for the correction of error in the 2D cone beam projections due to abnormalities often observed in 2D x-ray flat panel detectors. Template images are derived from the responses of the detector pixels using their statistical properties and then an effective non-causal derivative-based detection algorithm in 2D space is presented for the detection of defective and mis-calibrated detector elements separately. An image inpainting-based 3D correction scheme is proposed for the estimation of responses of defective detector elements, and the responses of the mis-calibrated detector elements are corrected using the normalization technique. For real-time implementation, a simplification of the proposed off-line method is also suggested. Finally, the proposed algorithms are tested using different real cone beam volume CT images and the experimental results demonstrate that the proposed methods can effectively remove ring and radiant artifacts from cone beam volume CT images compared to other reported techniques in the literature.     

基于Contourlet变换的CT和锥形束CT图像配准算法

目的提出一种基于Contourlet变换,用于放射治疗定位的CT与锥形束CT（cone beam CT,CBCT）图像配准的方法。方法利用Contourlet变换多尺度多方向的分辨特性,将待配准图像进行Contourlet变换分解,分解后的高频方向子带合成梯度图像,采用归一化互信息作为相似性测度,把梯度图像与低频方向子带以加权函数结合,进行临床医学图像的刚性配准,有效弥补了互信息配准中缺少空间信息的不足。结果通过已知空间变换参数图像的配准结果验证了算法的准确性。配准后lO幅图像变换参数的误差极小,且均方根误差接近于0。结论该图像配准算法精确度高,并具有很好的鲁棒性,有助于提高图像引导放射治疗（image guid edradiation therapy,IGRT）中解剖组织结构和靶区的定位精度。     

Image registration algorithm in CT and cone beam CT based on Contourlet transform

目的提出一种基于Contourlet变换,用于放射治疗定位的CT与锥形束CT(cone beam CT,CBCT)图像配准的方法.方法 利用Contourlet变换多尺度多方向的分辨特性,将待配准图像进行Contourlet变换分解,分解后的高频方向子带合成梯度图像,采用归一化互信息作为相似性测度,把梯度图像与低频方向子带以加权函数结合,进行临床医学图像的刚性配准,有效弥补了互信息配准中缺少空间信息的不足.结果 通过已知空间变换参数图像的配准结果验证了算法的准确性.配准后10幅图像变换参数的误差极小,且均方根误差接近于0.结论 该图像配准算法精确度高,并具有很好的鲁棒性,有助于提高图像引导放射治疗(image guided radiation therapy,IGRT)中解剖组织结构和靶区的定位精度.     

Respiratory correlated cone beam CT

A cone beam computed tomography (CBCT) scanner integrated with a linear accelerator is a powerful tool for image guided radiotherapy. Respiratory motion, however, induces artifacts in CBCT, while the respiratory correlated procedures, developed to reduce motion artifacts in axial and helical CT are not suitable for such CBCT scanners. We have developed an alternative respiratory correlated procedure for CBCT and evaluated its performance. This respiratory correlated CBCT procedure consists of retrospective sorting in projection space, yielding subsets of projections that each corresponds to a certain breathing phase. Subsequently, these subsets are reconstructed into a four-dimensional (4D) CBCT dataset. The breathing signal, required for respiratory correlation, was directly extracted from the 2D projection data, removing the need for an additional respiratory monitor system. Due to the reduced number of projections per phase, the contrast-to-noise ratio in a 4D scan reduced by a factor 2.6-3.7 compared to a 3D scan based on all projections. Projection data of a spherical phantom moving with a 3 and 5 s period with and without simulated breathing irregularities were acquired and reconstructed into 3D and 4D CBCT datasets. The positional deviations of the phantoms center of gravity between 4D CBCT and fluoroscopy were small: 0.13 +/- 0.09 mm for the regular motion and 0.39 +/- 0.24 mm for the irregular motion. Motion artifacts, clearly present in the 3D CBCT datasets, were substantially reduced in the 4D datasets, even in the presence of breathing irregularities, such that the shape of the moving structures could be identified more accurately. Moreover, the 4D CBCT dataset provided information on the 3D trajectory of the moving structures, absent in the 3D data. Considerable breathing irregularities, however, substantially reduces the image quality. Data presented for three different lung cancer patients were in line with the results obtained from the phantom study. In conclusion, we have successfully implemented a respiratory correlated CBCT procedure yielding a 4D dataset. With respiratory correlated CBCT on a linear accelerator, the mean position, trajectory, and shape of a moving tumor can be verified just prior to treatment. Such verification reduces respiration induced geometrical uncertainties, enabling safe delivery of 4D radiotherapy such as gated radiotherapy with small margins.     

Projection correlation based view interpolation for cone beam CT: primary fluence restoration in scatter measurement with a moving beam stop array

Scatter correction is an open problem in x-ray cone beam (CB) CT. The measurement of scatter intensity with a moving beam stop array (BSA) is a promising technique that offers a low patient dose and accurate scatter measurement. However, when restoring the blocked primary fluence behind the BSA, spatial interpolation cannot well restore the high-frequency part, causing streaks in the reconstructed image. To address this problem, we deduce a projection correlation (PC) to utilize the redundancy (over-determined information) in neighbouring CB views. PC indicates that the main high-frequency information is contained in neighbouring angular projections, instead of the current projection itself, which provides a guiding principle that applies to high-frequency information restoration. On this basis, we present the projection correlation based view interpolation (PC-VI) algorithm; that it outperforms the use of only spatial interpolation is validated. The PC-VI based moving BSA method is developed. In this method, PC-VI is employed instead of spatial interpolation, and new moving modes are designed, which greatly improve the performance of the moving BSA method in terms of reliability and practicability. Evaluation is made on a high-resolution voxel-based human phantom realistically including the entire procedure of scatter measurement with a moving BSA, which is simulated by analytical ray-tracing plus Monte Carlo simulation with EGSnrc. With the proposed method, we get visually artefact-free images approaching the ideal correction. Compared with the spatial interpolation based method, the relative mean square error is reduced by a factor of 6.05-15.94 for different slices. PC-VI does well in CB redundancy mining; therefore, it has further potential in CBCT studies.     

A geometric calibration method for cone beam CT systems

Cone beam CT systems are being deployed in large numbers for small animal imaging, dental imaging, and other specialty applications. A new high-precision method for cone beam CT system calibration is presented in this paper. It uses multiple projection images acquired from rotating point-like objects (metal ball bearings) and the angle information generated from the rotating gantry system is also used. It is assumed that the whole system has a mechanically stable rotation center and that the detector does not have severe out-of-plane rotation (<2 degrees). Simple geometrical relationships between the orbital paths of individual BBs and five system parameters were derived. Computer simulations were employed to validate the accuracy of this method in the presence of noise. Equal or higher accuracy was achieved compared with previous methods. This method was implemented for the geometrical calibration of both a micro CT scanner and a breast CT scanner. The reconstructed tomographic images demonstrated that the proposed method is robust and easy to implement with high precision.     

Effect of electron contamination on scatter correction factors for photon beam dosimetry.

Physical quantities for use in megavoltage photon beam dose calculations which are defined at the depth of maximum absorbed dose are sensitive to electron contamination and are difficult to measure and to calculate. Recently, formalisms have therefore been presented to assess the dose using collimator and phantom scatter correction factors, Sc and Sp, defined at a reference depth of 10 cm. The data can be obtained from measurements at that depth in a miniphantom and in a full scatter phantom. Equations are presented that show the relation between these quantities and corresponding quantities obtained from measurements at the depth of the dose maximum. It is shown that conversion of Sc and Sp determined at a 10 cm depth to quantities defined at the dose maximum such as (normalized) peak scatter factor, (normalized) tissue-air ratio, and vice versa is not possible without quantitative knowledge of the electron contamination. The difference in Sc at dmax resulting from this electron contamination compared with Sc values obtained at a depth of 10 cm in a miniphantom has been determined as a multiplication factor, Scel, for a number of photon beams of different accelerator types. It is shown that Scel may vary up to 5%. Because in the new formalisms output factors are defined at a reference depth of 10 cm, they do not require Scel data. The use of Sc and Sp values, defined at a 10 cm depth, combined with relative depth-dose data or tissue-phantom ratios is therefore recommended. For a transition period the use of the equations provided in this article and Scel data might be required, for instance, if treatment planning systems apply Sc data normalized at d(max).     

Small angle scatter imaging from wide beam diffraction patterns

In this paper we report on the extension of the technique of mapping small angle x-ray scatter (SAXS) across a soft material specimen several millimetres square. In the conventional SAXS mapping technique a pencil beam of x-rays is raster scanned over the specimen with the scatter pattern recorded from each point in the raster. In our technique a wide, parallel beam is used, speeding up the data collection time considerably. An image processing algorithm is used to separate the scatter pattern features from individual points along the line of the beam. To test the efficacy of the technique a phantom was constructed using gelatin and rat tail tendon collagen. Collagen fibres in the phantom were arranged in quarters horizontally, diagonally and vertically leaving one quarter with just gelatin. The phantom was used to collect both raster scanned sets of SAXS patterns spaced at 0.25 mm horizontally and vertically and also a wide beam data set. The width of the beam in this case was approximately 7 mm. Using the third-order diffraction of rat tail tendon intensity data were gathered from each SAXS pattern and used to construct a map. Data from the raster scan image and that from the wide beam are compared. Finally using a phantom made from dehydrated rat tail tendon and paraffin wax a tomographic slice constructed using data from SAXS patterns is shown.     

Helical cone beam CT with an asymmetrical detector

If a multislice or other area detector is shifted to one side to cover a larger field of view, then the data are truncated on one side. We propose a method to restore the missing data in helical cone-beam acquisitions that uses measured data on the longer side of the asymmetric detector array. The method is based on the idea of complementary rays, which is well known in fan beam geometry; in this paper we extend this concept to the cone-beam case. Different cases of complementary data coverage and dependence on the helical pitch are considered. The proposed method is used in our prototype 16-row CT scanner with an asymmetric detector and a 700 mm field of view. For evaluation we used scanned body phantom data and computer-simulated data. To simulate asymmetric truncation, the full, symmetric datasets were truncated by dropping either 22.5% or 45% from one side of the detector. Reconstructed images from the prototype scanner with the asymmetrical detector show excellent image quality in the extended field of view. The proposed method allows flexible helical pitch selection and can be used with overscan, short-scan, and super-short-scan reconstructions.     

A megavoltage scatter correction technique for cone-beam CT images acquired during VMAT delivery

Kilovoltage cone-beam CT (kV CBCT) can be acquired during the delivery of volumetric modulated arc therapy (VMAT), in order to obtain an image of the patient during treatment. However, the quality of such CBCTs is degraded by megavoltage (MV) scatter from the treatment beam onto the imaging panel. The objective of this paper is to introduce a novel MV scatter correction method for simultaneous CBCT during VMAT, and to investigate its effectiveness when compared to other techniques. The correction requires the acquisition of a separate set of images taken during VMAT delivery, while the kV beam is off. These images--which contain only the MV scatter contribution on the imaging panel--are then used to correct the corresponding kV/MV projections. To test this method, CBCTs were taken of an image quality phantom during VMAT delivery and measurements of contrast to noise ratio were made. Additionally, the correction was applied to the datasets of three VMAT prostate patients, who also received simultaneous CBCTs. The clinical image quality was assessed using a validated scoring system, comparing standard CBCTs to the uncorrected simultaneous CBCTs and a variety of correction methods. Results show that the correction is able to recover some of the low and high-contrast signal to noise ratio lost due to MV scatter. From the patient study, the corrected CBCT scored significantly higher than the uncorrected images in terms of the ability to identify the boundary between the prostate and surrounding soft tissue. In summary, a simple MV scatter correction method has been developed and, using both phantom and patient data, is shown to improve the image quality of simultaneous CBCTs taken during VMAT delivery.