An attenuation field network for dedicated cone beam breast CT with short scan and offset detector geometry

The feasibility of full-scan, offset-detector geometry cone-beam CT has been demonstrated for several clinical applications. For full-scan acquisition with offset-detector geometry, data redundancy from complementary views can be exploited during image reconstruction. Envisioning an upright breast CT system, we propose to acquire short-scan data in conjunction with offset-detector geometry. To tackle the resulting incomplete data, we have developed a self-supervised attenuation field network (AFN). AFN leverages the inherent redundancy of cone-beam CT data through coordinate-based representation and known imaging physics. A trained AFN can query attenuation coefficients using their respective coordinates or synthesize projection data including the missing projections. The AFN was evaluated using clinical cone-beam breast CT datasets (n = 50). While conventional analytical and iterative reconstruction methods failed to reconstruct the incomplete data, AFN reconstruction was not statistically different from the reference reconstruction obtained using full-scan, full-detector data in terms of image noise, image contrast, and the full width at half maximum of calcifications. This study indicates the feasibility of a simultaneous short-scan and offset-detector geometry for dedicated breast CT imaging. The proposed AFN technique can potentially be expanded to other cone-beam CT applications.

Technical Note: Intrinsic raw data-based CT misalignment correction without redundant data

PURPOSE

CT image reconstruction requires accurate knowledge of the used geometry or image quality might be degraded by misalignment artifacts. To overcome this issue, an intrinsic method, that is, a method not requiring a dedicated calibration phantom, to perform a raw data-based misalignment correction for CT is proposed herein that does not require redundant data and hence is applicable to measurements with less than 180 ∘ plus fan-angle of data.

METHODS

The forward projection of a volume reconstructed from a misaligned geometry resembles the acquired raw data if no redundant data are used, that is, if less than 180 ∘ plus fan-angle are used for image reconstruction. Hence, geometric parameters cannot be deduced from such data by an optimization of the geometry-dependent raw data fidelity. We propose to use a nonlinear transform applied to the reconstructed volume to introduce inconsistencies in the raw data that can be employed to estimate geometric parameters using less than 180 ∘ plus fan-angle of data. The proposed method is evaluated using simulations of the FORBILD head phantom and using actual measurements of a contrast-enhanced scan of a mouse acquired using a micro-CT.

RESULTS

Noisy simulations and actual measurements demonstrate that the proposed method is capable of correcting for artifacts arising from a misaligned geometry without redundant data while ensuring raw data fidelity.

CONCLUSIONS

The proposed method extends intrinsic raw data-based misalignment correction methods to an angular range of 180 ∘ or less and is thus applicable to systems with a limited scan range.

A novel concept for CT with fixed anodes (FACT): Medical imaging based on the feasibility of thermal load capacity

Focussing primarily on thermal load capacity, we describe the performance of a novel fixed anode CT (FACT) compared with a 100 kW reference CT. Being a fixed system, FACT has no focal spot blurring of the X-ray source during projection. Monte Carlo and finite element methods were used to determine the fluence proportional to thermal capacity. Studies of repeated short-time exposures showed that FACT could operate in pulsed mode for an unlimited period. A virtual model for FACT was constructed to analyse various temporal sequences for the X-ray source ring, representing a circular array of 1160 fixed anodes in the gantry. Assuming similar detector properties at a very small integration time, image quality was investigated using an image reconstruction library. Our model showed that approximately 60 gantry rounds per second, i.e. 60 sequential targetings of the 1160 anodes per second, were required to achieve a performance level equivalent to that of the reference CT (relative performance, RP = 1) at equivalent image quality. The optimal projection duration in each direction was about 10 μs. With a beam pause of 1 μs between projections, 78.4 gantry rounds per second with consecutive source activity were thermally possible at a given thermal focal spot. The settings allowed for a 1.3-fold (RP = 1.3) shorter scan time than conventional CT while maintaining radiation exposure and image quality. Based on the high number of rounds, FACT supports a high image frame rate at low doses, which would be beneficial in a wide range of diagnostic and technical applications.

Artifact-resistant motion estimation with a patient-specific artifact model for motion-compensated cone-beam CT

PURPOSE

In image-guided radiation therapy (IGRT) valuable information for patient positioning, dose verification, and adaptive treatment planning is provided by an additional kV imaging unit. However, due to the limited gantry rotation speed during treatment the typical acquisition time is quite long. Tomographic images of the thorax suffer from motion blurring or, if a gated 4D reconstruction is performed, from significant streak artifacts. Our purpose is to provide a method that reliably estimates respiratory motion in presence of severe artifacts. The estimated motion vector fields are then used for motion-compensated image reconstruction to provide high quality respiratory-correlated 4D volumes for on-board cone-beam CT (CBCT) scans.

METHODS

The proposed motion estimation method consists of a model that explicitly addresses image artifacts because in presence of severe artifacts state-of-the-art registration methods tend to register artifacts rather than anatomy. Our artifact model, e.g., generates streak artifacts very similar to those included in the gated 4D CBCT images, but it does not include respiratory motion. In combination with a registration strategy, the model gives an error estimate that is used to compensate the corresponding errors of the motion vector fields that are estimated from the gated 4D CBCT images. The algorithm is tested in combination with a cyclic registration approach using temporal constraints and with a standard 3D-3D registration approach. A qualitative and quantitative evaluation of the motion-compensated results was performed using simulated rawdata created on basis of clinical CT data. Further evaluation includes patient data which were scanned with an on-board CBCT system.

RESULTS

The model-based motion estimation method is nearly insensitive to image artifacts of gated 4D reconstructions as they are caused by angular undersampling. The motion is accurately estimated and our motion-compensated image reconstruction algorithm can correct for it. Motion artifacts of 3D standard reconstruction are significantly reduced, while almost no new artifacts are introduced.

CONCLUSIONS

Using the artifact model allows to accurately estimate and compensate for patient motion, even if the initial reconstructions are of very low image quality. Using our approach together with a cyclic registration algorithm yields a combination which shows almost no sensitivity to sparse-view artifacts and thus ensures both high spatial and high temporal resolution.

Panoramic cone beam computed tomography

PURPOSE

Cone-beam computed tomography (CBCT) is the main imaging tool for image-guided radiotherapy but its functionality is limited by a small imaging volume and restricted image position (imaged at the central instead of the treatment position for peripheral lesions to avoid collisions). In this paper, the authors present the concept of "panoramic CBCT," which can image patients at the treatment position with an imaging volume as large as practically needed.

METHODS

In this novel panoramic CBCT technique, the target is scanned sequentially from multiple view angles. For each view angle, a half scan (180° + θ(cone) where θ(cone) is the cone angle) is performed with the imaging panel positioned in any location along the beam path. The panoramic projection images of all views for the same gantry angle are then stitched together with the direct image stitching method (i.e., according to the reported imaging position) and full-fan, half-scan CBCT reconstruction is performed using the stitched projection images. To validate this imaging technique, the authors simulated cone-beam projection images of the Mathematical Cardiac Torso (MCAT) thorax phantom for three panoramic views. Gaps, repeated/missing columns, and different exposure levels were introduced between adjacent views to simulate imperfect image stitching due to uncertainties in imaging position or output fluctuation. A modified simultaneous algebraic reconstruction technique (modified SART) was developed to reconstruct CBCT images directly from the stitched projection images. As a gold standard, full-fan, full-scan (360° gantry rotation) CBCT reconstructions were also performed using projection images of one imaging panel large enough to encompass the target. Contrast-to-noise ratio (CNR) and geometric distortion were evaluated to quantify the quality of reconstructed images. Monte Carlo simulations were performed to evaluate the effect of scattering on the image quality and imaging dose for both standard and panoramic CBCT.

RESULTS

Truncated images with artifacts were observed for the CBCT reconstruction using projection images of the central view only. When the image stitching was perfect, complete reconstruction was obtained for the panoramic CBCT using the modified SART with the image quality similar to the gold standard (full-scan, full-fan CBCT using one large imaging panel). Imperfect image stitching, on the other hand, lead to (streak, line, or ring) reconstruction artifacts, reduced CNR, and/or distorted geometry. Results from Monte Carlo simulations showed that, for identical imaging quality, the imaging dose was lower for the panoramic CBCT than that acquired with one large imaging panel. For the same imaging dose, the CNR of the three-view panoramic CBCT was 50% higher than that of the regular CBCT using one big panel.

CONCLUSIONS

The authors have developed a panoramic CBCT technique and demonstrated with simulation data that it can image tumors of any location for patients of any size at the treatment position with comparable or less imaging dose and time. However, the image quality of this CBCT technique is sensitive to the reconstruction artifacts caused by imperfect image stitching. Better algorithms are therefore needed to improve the accuracy of image stitching for panoramic CBCT.

Low-dose cardio-respiratory phase-correlated cone-beam micro-CT of small animals

PURPOSE

Micro-CT imaging of animal hearts typically requires a double gating procedure because scans during a breath-hold are not possible due to the long scan times and the high respiratory rates, Simultaneous respiratory and cardiac gating can either be done prospectively or retrospectively. True five-dimensional information can be either retrieved with retrospective gating or with prospective gating if several prospective gates are acquired. In any case, the amount of information available to reconstruct one volume for a given respiratory and cardiac phase is orders of magnitud lower than the total amount of information acquired. For example, the reconstruction of a volume from a 10% wide respiratory and a 20% wide cardiac window uses only 2% of the data acquired. Achieving a similar image quality as a nongated scan would therefore require to increase the amount of data and thereby the dose to the animal by up to a factor of 50.

METHODS

To achieve the goal of low-dose phase-correlated (LDPC) imaging, the authors propose to use a highly efficient combination of slightly modified existing algorithms. In particular, the authors developed a variant of the McKinnon-Bates image reconstruction algorithm and combined it with bilateral filtering in up to five dimensions to significantly reduce image noise without impairing spatial or temporal resolution.

RESULTS

The preliminary results indicate that the proposed LDPC reconstruction method typically reduces image noise by a factor of up to 6 (e.g., from 170 to 30 HU), while the dose values lie in a range from 60 to 500 mGy. Compared to other publications that apply 250-1800 mGy for the same task [C. T. Badea et al., "4D micro-CT of the mouse heart," Mol. Imaging 4(2), 110-116 (2005); M. Drangova et al., "Fast retrospectively gated quantitative four-dimensional (4D) cardiac micro computed tomography imaging of free-breathing mice," Invest. Radiol. 42(2), 85-94 (2007); S. H. Bartling et al., "Retrospective motion gating in small animal CT of mice and rats," Invest. Radiol. 42(10), 704-714 (2007)], the authors' LDPC approach therefore achieves a more than tenfold dose usage improvement.

CONCLUSIONS

The LDPC reconstruction method improves phase-correlated imaging from highly undersampled data. Artifacts caused by sparse angular sampling are removed and the image noise is decreased, while spatial and temporal resolution are preserved. Thus, the administered dose per animal can be decreased allowing for long-term studies with reduced metabolic inference.

The influence of the heel effect in cone-beam computed tomography: artifacts in standard and novel geometries and their correction

For decades, the heel effect has been known to cause an angular dependence of the emitted spectrum of an x-ray tube. In radiography, artifacts were observed and attributed to the heel effect. However, no problems due to the heel effect were discerned in multi-slice computed tomography (MSCT) so far. With flat-detector CT (FDCT), involving larger cone angles and different system geometries, the heel effect might cause new artifacts. These artifacts were analyzed in this paper for system geometries different from the ones widely used nowadays. Simulations and measurements were performed. Simulations included symmetric as well as asymmetric detector layouts and different x-ray tube orientations with respect to the detector plane. The measurements were performed on a micro-CT system in an asymmetric detector layout. Furthermore, an analytical correction scheme is proposed to overcome heel effect artifacts. It was shown that the type of artifact greatly depends on the orientation of the x-ray tube and also on the type of detector alignment (i.e. symmetric or different types of asymmetric alignment). Certain combinations exhibited almost no significant artifact while others greatly influenced the quality of the reconstructed images. The proposed correction scheme showed good results that were further improved when also applying a scatter correction. When designing CT systems, care should be taken when placing the tube and the detector. Orientation of the x-ray tube like in most MSCT systems seems advisable in asymmetric detector layouts. However, a different type of tube orientation can be overcome with suitable correction schemes.