Geometric Calibration Using Line Fiducials for Cone-Beam CT with General, Non-Circular Source-Detector Trajectories

Line-based iterative geometric calibration method for a tomosynthesis system

BACKGROUND

A next generation tomosynthesis (NGT) system, capable of two-dimensional source motion, detector motion in the perpendicular direction, and magnification tomosynthesis, was constructed to investigate different acquisition geometries. Existing position-based geometric calibration methods proved ineffective when applied to the NGT geometries.

PURPOSE

A line-based iterative calibration method is developed to perform accurate geometric calibration for the NGT system.

METHODS

The proposed method calculates the system geometry through virtual line segments created by pairs of fiducials within a calibration phantom, by minimizing the error between the line equations computed from the true and estimated fiducial projection pairs. It further attempts to correct the 3D fiducial locations based on the initial geometric calibration. The method's performance was assessed via simulation and experimental setups with four distinct NGT geometries: X, T, XZ, and TZ. The X geometry resembles a conventional DBT acquisition along the chest wall. The T geometry forms a "T"-shaped source path in mediolateral (ML) and posteroanterior (PA) directions. A descending detector motion is added to both X and T geometries to form the XZ and TZ geometries, respectively. Simulation studies were conducted to assess the robustness of the method to geometric perturbations and inaccuracies in fiducial locations. Experimental studies were performed to assess the impact of phantom magnification and the performance of the proposed method for various geometries, compared to the traditional position-based method. Star patterns were evaluated for both qualitative and quantitative analyses; the Fourier spectral distortions (FSDs) graphs and the contrast transfer function (CTF) were extracted. The limit of spatial resolution (LSR) was measured at 5% modulation of the CTF.

RESULTS

The proposed method presented is highly robust to geometric perturbation and fiducial inaccuracies. After the line-based iterative method, the mean distance between the true and estimated fiducial projections was [X, T, XZ, TZ]: [0.01, 0.01, 0.02, 0.01] mm. The impact of phantom magnification was observed; a contact-mode acquisition of a calibration phantom successfully provided an accurate geometry for 1.85× magnification images of a star pattern, with the X geometry. The FSD graphs for the contact-mode T geometry acquisition presented evidence of super-resolution, with the LSR of [0°-quadrant: 8.57, 90°-quadrant: 8.47] lp/mm. Finally, a contact-mode XZ geometry acquisition and a 1.50× magnification TZ geometry acquisition were reconstructed with three calibration methods-position-based, line-based, and iterative line-based. As more advanced methods are applied, the CTF becomes more isotropic, the FSD graphs demonstrate less spectral leakage as super-resolution is achieved, and the degree of blurring artifacts reduces significantly.

CONCLUSIONS

This study introduces a robust calibration method tailored to the unique requirements of advanced tomosynthesis systems. By employing virtual line segments and iterative techniques, we ensure accurate geometric calibration while mitigating the limitations posed by the complex acquisition geometries of the NGT system. Our method's ability to handle various NGT configurations and its tolerance to fiducial misalignment make it a superior choice compared to traditional calibration techniques.

A nonconvex model-based combined geometric calibration scheme for micro cone-beam CT with irregular trajectories

BACKGROUND

Many dedicated cone-beam CT (CBCT) systems have irregular scanning trajectories. Compared with the standard CBCT calibration, accurate calibration for CBCT systems with irregular trajectories is a more complex task, since the geometric parameters for each scanning view are variable. Most of the existing calibration methods assume that the intrinsic geometric relationship of the fiducials in the phantom is precisely known, and rarely delve deeper into the issue of whether the phantom accuracy is adapted to the calibration model.

PURPOSE

A high-precision phantom and a highly robust calibration model are interdependent and mutually supportive, and they are both important for calibration accuracy, especially for the high-resolution CBCT. Therefore, we propose a calibration scheme that considers both accurate phantom measurement and robust geometric calibration.

METHODS

Our proposed scheme consists of two parts: (1) introducing a measurement model to acquire the accurate intrinsic geometric relationship of the fiducials in the phantom; (2) developing a highly noise-robust nonconvex model-based calibration method. The measurement model in the first part is achieved by extending our previous high-precision geometric calibration model suitable for CBCT with circular trajectories. In the second part, a novel iterative method with optimization constraints based on a back-projection model is developed to solve the geometric parameters of each view.

RESULTS

The simulations and real experiments show that the measurement errors of the fiducial ball bearings (BBs) are within the subpixel level. With the help of the geometric relationship of the BBs obtained by our measurement method, the classic calibration method can achieve good calibration based on far fewer BBs. All metrics obtained in simulated experiments as well as in real experiments on Micro CT systems with resolutions of 9 and 4.5 μm show that the proposed calibration method has higher calibration accuracy than the competing classic method. It is particularly worth noting that although our measurement model proves to be very accurate, the classic calibration method based on this measurement model can only achieve good calibration results when the resolution of the measurement system is close to that of the system to be calibrated, but our calibration scheme enables high-accuracy calibration even when the resolution of the system to be calibrated is twice that of the measurement system.

CONCLUSIONS

The proposed combined geometrical calibration scheme does not rely on a phantom with an intricate pattern of fiducials, so it is applicable in Micro CT with high resolution. The two parts of the scheme, the "measurement model" and the "calibration model," prove to be of high accuracy. The combination of these two models can effectively improve the calibration accuracy, especially in some extreme cases.

Correction of geometric artifact in cone-beam computed tomography through a deep neural network

Cone-beam computed tomography is a noninvasive detection system that can obtain the three-dimensional structure of objects in a way that does not damage the object. It is widely applied in precision instruments, medical detection, and other fields. However, in the actual process, if a geometric artifact appears in the results, it will affect the quality of reconstructed images, including detail loss and decreased spatial resolution, which leads to inaccurate distinction of defects in detection. We propose a method for correcting a geometric artifact by means of data-driven projection and neural networks. The network designed is a deep neural network with six convolutional layers and six deconvolutional layers that can correct a geometric artifact by training a large number of labeled data and unlabeled data. Compared with other networks that require prior information for reconstructed images, the proposed method uses a projection data-driven approach that can avoid the requirement for prior information. The simulation data have been tested under varying degrees of noise, and satisfactory geometric artifact correction results have been obtained. Meanwhile, we use the actual data of line pairs and ball grid array solder joints to conduct experiments. The results obtained by our method are compared with two other phantom-based method and the U-net method, respectively. The results of similarity and spatial resolution show that the proposed method can achieve the comparable results as the two types of methods. At the same time, we apply a projection data-driven approach to avoid the requirement for prior information, which is more conducive to the correction of the geometric artifact in practical situations where prior information is lacking.