Iterative Image Reconstruction Research Articles

Iterative image reconstruction algorithm analysis for optical CT radiochromic gel dosimetry

Background. Modern radiation therapy technologies aim to enhance radiation dose precision to the tumor and utilize hypofractionated treatment regimens. Verifying the dose distributions associated with these advanced radiation therapy treatments remains an active research area due to the complexity of delivery systems and the lack of suitable three-dimensional dosimetry tools. Gel dosimeters are a potential tool for measuring these complex dose distributions. A prototype tabletop solid-tank fan-beam optical CT scanner for readout of gel dosimeters was recently developed. This scanner does not have a straight raypath from source to detector, thus images cannot be reconstructed using filtered backprojection (FBP) and iterative techniques are required. Purpose. To compare a subset of the top performing algorithms in terms of image quality and quantitatively determine the optimal algorithm while accounting for refraction within the optical CT system. The following algorithms were compared: Landweber, superiorized Landweber with the fast gradient projection perturbation routine (S-LAND-FGP), the fast iterative shrinkage/thresholding algorithm with total variation penalty term (FISTA-TV), a monotone version of FISTA-TV (MFISTA-TV), superiorized conjugate gradient with the nonascending perturbation routine (S-CG-NA), superiorized conjugate gradient with the fast gradient projection perturbation routine (S-CG-FGP), superiorized conjugate gradient with with two iterations of CG performed on the current iterate and the nonascending perturbation routine (S-CG-2-NA). Methods. A ray tracing simulator was developed to track the path of light rays as they traverse the different mediums of the optical CT scanner. Two clinical phantoms and several synthetic phantoms were produced and used to evaluate the reconstruction techniques under known conditions. Reconstructed images were analyzed in terms of spatial resolution, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), signal non-uniformity (SNU), mean relative difference (MRD) and reconstruction time. We developed an image quality based method to find the optimal stopping iteration window for each algorithm. Imaging data from the prototype optical CT scanner was reconstructed and analysed to determine the optimal algorithm for this application. Results. The optimal algorithms found through the quantitative scoring metric were FISTA-TV and S-CG-2-NA. MFISTA-TV was found to behave almost identically to FISTA-TV however MFISTA-TV was unable to resolve some of the synthetic phantoms. S-CG-NA showed extreme fluctuations in the SNR and CNR values. S-CG-FGP had large fluctuations in the SNR and CNR values and the algorithm has less noise reduction than FISTA-TV and worse spatial resolution than S-CG-2-NA. S-LAND-FGP had many of the same characteristics as FISTA-TV; high noise reduction and stability from over iterating. However, S-LAND-FGP has worse SNR, CNR and SNU values as well as longer reconstruction time. S-CG-2-NA has superior spatial resolution to all algorithms while still maintaining good noise reduction and is uniquely stable from over iterating. Conclusions. Both optimal algorithms (FISTA-TV and S-CG-2-NA) are stable from over iterating and have excellent edge detection with ESF MTF 50% values of 1.266 mm−1 and 0.992 mm−1. FISTA-TV had the greatest noise reduction with SNR, CNR and SNU values of 424, 434 and 0.91 × 10−4, respectively. However, low spatial resolution makes FISTA-TV only viable for large field dosimetry. S-CG-2-NA has better spatial resolution than FISTA-TV with PSF and LSF MTF 50% values of 1.581 mm−1 and 0.738 mm−1, but less noise reduction. S-CG-2-NA still maintains good SNR, CNR, and SNU values of 168, 158 and 1.13 × 10−4, respectively. Thus, S-CG-2-NA is a well rounded reconstruction algorithm that would be the preferable choice for small field dosimetry.

Sparse-view cone-beam computed tomography iterative reconstruction based on new multi-gradient direction total variation.

The accurate reconstruction of cone-beam computed tomography (CBCT) from sparse projections is one of the most important areas for study. The compressed sensing theory has been widely employed in the sparse reconstruction of CBCT. However, the total variation (TV) approach solely uses information from the i-coordinate, j-coordinate, and k-coordinate gradients to reconstruct the CBCT image. It is well recognized that the CBCT image can be reconstructed more accurately with more gradient information from different directions. Thus, this study introduces a novel approach, named the new multi-gradient direction total variation minimization method. The method uses gradient information from the ij-coordinate, ik-coordinate, and jk-coordinate directions to reconstruct CBCT images, which incorporates nine different types of gradient information from nine directions. This study assessed the efficacy of the proposed methodology using under-sampled projections from four different experiments, including two digital phantoms, one patient's head dataset, and one physical phantom dataset. The results indicated that the proposed method achieved the lowest RMSE index and the highest SSIM index. Meanwhile, we compared the voxel intensity curves of the reconstructed images to assess the edge structure preservation. Among the various methods compared, the curves generated by the proposed method exhibited the highest level of consistency with the gold standard image curves. In summary, the proposed method showed significant potential in enhancing the quality and accuracy of CBCT image reconstruction.

Improved Complex Convolutional Neural Network Based on SPIRiT and Dense Connection for Parallel MRI Reconstruction

To accelerate the data acquisition speed of magnetic resonance imaging (MRI) and improve the reconstructed MR images’ quality, we propose a parallel MRI reconstruction model (SPIRiT-Net), which combines the iterative self-consistent parallel imaging reconstruction model (SPIRiT) with the cascaded complex convolutional neural networks (CCNNs). More specifically, this model adopts the SPIRiT model for reconstruction in the k-space domain and the cascaded CCNNs with dense connection for reconstruction in the image domain. Meanwhile, this model introduces the data consistency layers for better reconstruction in both the image domain and the k-space domain. The experimental results on two clinical knee datasets as well as the fastMRI brain dataset under different undersampling patterns show that the SPIRiT-Net model achieves better reconstruction performance in terms of visual effect, peak signal-to-noise ratio, and structural similarity over SPIRiT, Deepcomplex, and DONet. It will be beneficial to the diagnosis of clinical medicine.

Development and validation of the effective CNR analysis method for evaluating the contrast resolution of CT images.

Contrast resolution is an important index for evaluating the signal detectability of computed tomographic (CT) images. Recently, various noise reduction algorithms, such as iterative reconstruction (IR) and deep learning reconstruction (DLR), have been proposed to reduce the image noise in CT images. However, these algorithms cause changes in the image noise texture and blurred image signals in CT images. Furthermore, the contrast-to-noise ratio (CNR) cannot be accurately evaluated in CT images reconstructed using noise reduction methods. Therefore, in this study, we devised a new method, namely, "effective CNR analysis," for evaluating the contrast resolution of CT images. We verified whether the proposed algorithm could evaluate the effective contrast resolution based on the signal detectability of CT images. The findings showed that the effective CNR values obtained using the proposed method correlated well with the subjective visual impressions of CT images. To investigate whether signal detectability was appropriately evaluated using effective CNR analysis, the conventional CNR analysis method was compared with the proposed method. The CNRs of the IR and DLR images calculated using conventional CNR analysis were 13.2 and 10.7, respectively. By contrast, those calculated using the effective CNR analysis were estimated to be 0.7 and 1.1, respectively. Considering that the signal visibility of DLR images was superior to that of IR images, our proposed effective CNR analysis was shown to be appropriate for evaluating the contrast resolution of CT images.

Image reconstruction based on nonlinear diffusion model for limited-angle computed tomography

The problem of limited-angle computed tomography (CT) imaging reconstruction has a wide range of practical applications. Due to various factors such as high x-ray absorption, structural characteristics of the scanned object, and equipment limitations, it is often impractical to obtain a complete angular scan, resulting in limited-angle scan data. In this paper, we propose an iterative image reconstruction algorithm for limited-angle CT. The algorithm carries out a traditional CT reconstruction and a nonlinear diffusion process alternatively. Specifically, a subtle partial differential equation is constructed to guide the nonlinear diffusion process to eliminate limited-angle artifacts in the reconstructed image. Numerical experiments on both analytic data and real data validate the efficacy of the proposed nonlinear diffusion reconstruction algorithm. Furthermore, a linear diffusion reconstruction algorithm which combines a traditional CT reconstruction algorithm and a linear diffusion process is also presented in the paper.

UMC-PET: a fast and flexible Monte Carlo PET simulator.

Objective.The GPU-based Ultra-fast Monte Carlo positron emission tomography simulator (UMC-PET) incorporates the physics of the emission, transport and detection of radiation in PET scanners. It includes positron range, non-colinearity, scatter and attenuation, as well as detector response. The objective of this work is to present and validate UMC-PET as a a multi-purpose, accurate, fast and flexible PET simulator.Approach.We compared UMC-PET against PeneloPET, a well-validated MC PET simulator, both in preclinical and clinical scenarios. Different phantoms for scatter fraction (SF) assessment following NEMA protocols were simulated in a 6R-SuperArgus and a Biograph mMR scanner, comparing energy histograms, NEMA SF, and sensitivity for different energy windows. A comparison with real data reported in the literature on the Biograph scanner is also shown.Main results.NEMA SF and sensitivity estimated by UMC-PET where within few percent of PeneloPET predictions. The discrepancies can be attributed to small differences in the physics modeling. Running in a 11 GB GeForce RTX 2080 Ti GPU, UMC-PET is ∼1500 to ∼2000 times faster than PeneloPET executing in a single core Intel(R) Xeon(R) CPU W-2155 @ 3.30 GHz.Significance.UMC-PET employs a voxelized scheme for the scanner, patient adjacent objects (such as shieldings or the patient bed), and the activity distribution. This makes UMC-PET extremely flexible. Its high simulation speed allows applications such as MC scatter correction, faster SRM estimation for complex scanners, or even MC iterative image reconstruction.

Optimization of the "Perth CT" Protocol for Preoperative Planning and Postoperative Evaluation in Total Knee Arthroplasty.

Background and Objectives: Total knee arthroplasty (TKA) has become the treatment of choice for advanced osteoarthritis. The aim of this paper was to show the possibilities of optimizing the Perth CT protocol, which is highly effective for preoperative planning and postoperative assessment of alignment. Materials and Methods: The cross-sectional study comprised 16 patients for preoperative planning or postoperative evaluation of TKA. All patients were examined with the standard and optimized Perth CT protocol using advance techniques, including automatic exposure control (AEC), iterative image reconstruction (IR), as well as a single-energy projection-based metal artifact reduction algorithm for eliminating prosthesis artifacts. The effective radiation dose (E) was determined based on the dose report. Imaging quality is determined according to subjective and objective (values of signal to noise ratio (SdNR) and figure of merit (FOM)) criteria. Results: The effective radiation dose with the optimized protocol was significantly lower compared to the standard protocol (p < 0.001), while in patients with the knee prosthesis, E increased significantly less with the optimized protocol compared to the standard protocol. No significant difference was observed in the subjective evaluation of image quality between protocols (p > 0.05). Analyzing the objective criteria for image quality optimized protocols resulted in lower SdNR values and higher FOM values. No significant difference of image quality was determined using the SdNR and FOM as per the specified protocols and parts of extremities, and for the presence of prothesis. Conclusions: Retrospecting the ALARA ('As Low As Reasonably Achievable') principles, it is possible to optimize the Perth CT protocol by reducing the kV and mAs values and by changing the collimation and increasing the pitch factor. Advanced IR techniques were used in both protocols, and AEC was used in the optimized protocol. The effective dose of radiation can be reduced five times, and the image quality will be satisfactory.

Five-view three-dimensional reconstruction for ultrafast dynamic imaging of pulsed radiation sources

Multiaxial neutron/x-ray imaging and three-dimensional (3D) reconstruction techniques play a crucial role in gaining valuable insights into the generation and evolution mechanisms of pulsed radiation sources. Owing to the short emission time (∼200 ns) and drastic changes of the pulsed radiation source, it is necessary to acquire projection data within a few nanoseconds in order to achieve clear computed tomography 3D imaging. As a consequence, projection data that can be used for computed tomography image reconstruction at a certain moment are often available for only a few angles. Traditional algorithms employed in the process of reconstructing 3D images with extremely incomplete data may introduce significant distortions and artifacts into the final image. In this paper, we propose an iterative image reconstruction method using cylindrical harmonic decomposition and a self-supervised denoising network algorithm based on the deep image prior method. We augment the prior information with a 2D total variation prior and a 3D deep image prior. Single-wire Z-pinch imaging experiments have been carried out at Qin-1 facility in five views and four frames, with a time resolution of 3 ns for each frame and a time interval of 40 ns between adjacent frames. Both numerical simulations and experiments verify that our proposed algorithm can achieve high-quality reconstruction results and obtain the 3D intensity distribution and evolution of extreme ultraviolet and soft x-ray emission from plasma.

Iterative Image Reconstruction with Under-Sampled Data Assisted by a Neural Network.

Image reconstruction with under-sampled data is usually achieved by an iterative algorithm, which minimizes an objective function. The objective function commonly contains a data fidelity term and one or more Bayesian terms. A popular Bayesian term is the total variation (TV) norm of the image. This paper suggests an addition Bayesian term that is generated by a neural network. This neural network is essentially a classifier. This classifier can recognize the artifacts caused by the incomplete data. This classifier is trained by patient images reconstructed by complete and incomplete data sets. This newly introduced Bayesian term is referred to as the CNN score, which is a real number in (-∞, ∞). Patient studies show the good correlation between the CNN score and the severeness of the artifacts due to the incomplete measurements. A neural network can extract features from images that are suffering from incomplete measurements and convert the features into a CNN score. An iterative image reconstruction algorithm can be developed to minimize this CNN score to suppress the artifacts in the reconstructed image.

Innovative advances in pediatric radiology: computed tomography reconstruction techniques, photon-counting detector computed tomography, and beyond.

In pediatric radiology, balancing diagnostic accuracy with reduced radiation exposure is paramount due to the heightened vulnerability of younger patients to radiation. Technological advancements in computed tomography (CT) reconstruction techniques, especially model-based iterative reconstruction and deep learning image reconstruction, have enabled significant reductions in radiation doses without compromising image quality. Deep learning image reconstruction, powered by deep learning algorithms, has demonstrated superiority over traditional techniques like filtered back projection, providing enhanced image quality, especially in pediatric head and cardiac CT scans. Photon-counting detector CT has emerged as another groundbreaking technology, allowing for high-resolution images while substantially reducing radiation doses, proving highly beneficial for pediatric patients requiring frequent imaging. Furthermore, cloud-based dose tracking software focuses on monitoring radiation exposure, ensuring adherence to safety standards. However, the deployment of these technologies presents challenges, including the need for large datasets, computational demands, and potential data privacy issues. This article provides a comprehensive exploration of these technological advancements, their clinical implications, and the ongoing efforts to enhance pediatric radiology's safety and effectiveness.

Application of Elastic Wave CT in Geotechnical Investigation of Suning Square

Elastic wave CT technology, also known as elastic wave tomography technology. It is on the premise of not destroying the internal structure of the rock and soil medium in the test area, according to the basic principle that elastic waves propagate at different speeds in different rock and soil media, the travel time and wave velocity of the elastic wave in the rock and soil between two points are obtained through the experiment, and through a special inversion method, a large amount of elastic wave velocity information is used to obtain the distribution law and the elastic wave velocity chromatogram of the medium in the test area, in this way, the lithology, structural development and velocity parameters of the rock-soil medium in the test area can be inferred. This article briefly introduces the basic theory of elastic wave CT, including Radon transform and inverse transform, forward ray tracing method and inversion iterative image reconstruction technology. Then combined with the engineering survey example of the Suning Plaza project, it briefly introduced the field work method, data parameter collection, geological interpretation and result analysis. The analysis shows that the results of the three-dimensional velocity map obtained by the elastic wave CT method are intuitive, the test accuracy is high, and the development state, size and location of the karst in the survey area can be clearly seen. It can describe the spatial distribution characteristics of karst in detail, has outstanding advantages in karst exploration, and has good application prospects.

Radial magnetic resonance image reconstruction with a deep unrolled projected fast iterative soft-thresholding network

Radially sampling of magnetic resonance imaging (MRI) is an effective way to accelerate the imaging. How to preserve the image details in reconstruction is always challenging. In this work, a deep unrolled neural network is designed to emulate the iterative sparse image reconstruction process of a projected fast soft-threshold algorithm (pFISTA). The proposed method, an unrolled pFISTA network for Deep Radial MRI (pFISTA-DR), include the preprocessing module to refine coil sensitivity maps and initial reconstructed image, the learnable convolution filters to extract image feature maps, and adaptive threshold to robustly remove image artifacts. Experimental results show that, among the compared methods, pFISTA-DR provides the best reconstruction and achieved the highest PSNR, the highest SSIM and the lowest reconstruction errors.

MCGPU-PET: An open-source real-time Monte Carlo PET simulator

Monte Carlo (MC) simulations are commonly used to model the emission, transmission, and/or detection of radiation in Positron Emission Tomography (PET). In this work, we introduce a new open-source MC software for PET simulation, MCGPU-PET, which has been designed to fully exploit the computing capabilities of modern GPUs to simulate the acquisition of more than 100 million coincidences per second from voxelized sources and material distributions. The new simulator is an extension of the PENELOPE-based MCGPU code previously used in cone-beam CT and mammography applications. We validated the accuracy of the accelerated code by comparing it to GATE and PeneloPET simulations achieving an agreement within 10 percent approximately. As an example application of the code for fast estimation of PET coincidences, a scan of the NEMA IQ phantom was simulated. A fully 3D sinogram with 6382 million true coincidences and 731 million scatter coincidences was generated in 54 s in one GPU. MCGPU-PET provides an estimation of true and scatter coincidences and spurious background (for positron-gamma emitters such as 124I) at a rate 3 orders of magnitude faster than CPU-based MC simulators. This significant speed-up enables the use of the code for accurate scatter and prompt-gamma background estimations within an iterative image reconstruction process. Program SummaryProgram Title: MCGPU-PETCPC Library link to program files:https://doi.org/10.17632/k5x2bsf27m.2Licensing provisions: CC0 1.0Programming language: C (with NVIDIA CUDA extensions)Journal reference of the previous version: Badal A, Sharma D, Graff C G, Zeng R and Badano A, Mammography and breast tomosynthesis simulator for virtual clinical trials, Computer Physics Communications 261, 107,779 (2021); https://doi.org/10.1016/j.cpc.2020.107779Reasons for the new version: The open source x-ray imaging simulator MCGPU was modified to simulate Positron Emission Tomography (PET), extending the medical imaging applications that can be modeled with the system. The main radiation transport algorithm has not been changed, but new radiation source and detector models were implemented.Nature of problem: Due to their stochastic nature, MC methods require a large number of computations to obtain low-noise estimations, which makes them usually very time-consuming. This has limited their application in PET imaging to offline tasks where the simulation time is not a critical factor, leaving out other relevant cases that would benefit from accurate calculations.Solution method: The development of computers with multiple CPUs and GPUs for scientific computing has enabled tackling real-time MC simulations, i.e. performing simulations fast enough to obtain results at the same time as other processes of the imaging pipeline. We introduce a new open source software, MCGPU-PET, for simulating PET acquisitions with GPUs that simulates more than 100 million detected coincidences per second, which makes it even faster than many real PET acquisitions. The software samples gamma rays emitted from inside a voxelized anatomical model following a user input radiotracer distribution. The gamma ray coincidences arriving at the ideal ring detector are tallied in either a sinogram or list mode, and can be reconstructed using standard PET reconstruction algorithms.

99mTc/123I Dual-Radionuclide Correction for Self-Scatter, Down-Scatter, and Tailing Effect for a CZT SPECT with Varying Tracer Distributions.

SPECT systems distinguish radionuclides by using multiple energy windows. For CZT detectors, the energy spectrum has a low energy tail leading to additional crosstalk between the radionuclides. Previous work developed models to correct the scatter and crosstalk for CZT-based dedicated cardiac systems with similar 99mTc/123I tracer distributions. These models estimate the primary and scatter components by solving a set of equations employing the MLEM approach. A penalty term is applied to ensure convergence. The present work estimates the penalty term for any 99mTc/123I activity level. An iterative approach incorporating Monte Carlo into the iterative image reconstruction loops was developed to estimate the penalty terms. We used SIMIND and XCAT phantoms in this study. Distribution of tracers in the myocardial tissue and blood pool were varied to simulate a dynamic acquisition. Evaluations of the estimated and the real penalty terms were performed using simulations and large animal data. The myocardium to blood pool ratio was calculated using ROIs in the myocardial tissue and the blood pool for quantitative analysis. All corrected images yielded a good agreement with the gold standard images. In conclusion, we developed a CZT crosstalk correction method for quantitative imaging of 99mTc/123I activity levels by dynamically estimating the penalty terms.

Multi-shaping sparse-continuous reconstruction for an optical coherence tomography sidelobe suppression.

Optical coherence tomography (OCT) images are commonly affected by sidelobe artifacts due to spectral non-uniformity and spectral leakage. Conventional frequency domain spectral shaping methods widen the mainlobe and compromise axial resolution. While image-domain deconvolution techniques can address the trade-off between axial resolution and artifact suppression, their reconstruction quality relies on accurate measurement or estimation of system point spread function (PSF). Inaccurate PSF estimation leads to loss of details in the reconstructed images. In this Letter, we introduce multi-shaping sparse-continuous reconstruction (MSSCR) for an OCT image, a novel, to the best of our knowledge, framework that combines spectral multi-shaping and iterative image reconstruction with sparse-continuous priors. The MSSCR achieves sidelobe suppression without requiring any PSF measurement or estimation and effectively preserving the axial resolution. The experimental results demonstrate that the MSSCR achieves sidelobe suppression of more than 8 dB. We believe that the MSSCR holds potential for addressing sidelobe artifacts in OCT.

Coil sketching for computationally efficient MR iterative reconstruction.

Parallel imaging and compressed sensing reconstructions of large MRI datasets often have a prohibitive computational cost that bottlenecks clinical deployment, especially for three-dimensional (3D) non-Cartesian acquisitions. One common approach is to reduce the number of coil channels actively used during reconstruction as in coil compression. While effective for Cartesian imaging, coil compression inherently loses signal energy, producing shading artifacts that compromise image quality for 3D non-Cartesian imaging. We propose coil sketching, a general and versatile method for computationally-efficient iterative MR image reconstruction. We based our method on randomized sketching algorithms, a type of large-scale optimization algorithms well established in the fields of machine learning and big data analysis. We adapt the sketching theory to the MRI reconstruction problem via a structured sketching matrix that, similar to coil compression, considers high-energy virtual coils obtained from principal component analysis. But, unlike coil compression, it also considers random linear combinations of the remaining low-energy coils, effectively leveraging information from all coils. First, we performed ablation experiments to validate the sketching matrix design on both Cartesian and non-Cartesian datasets. The resulting design yielded both improved computatioanal efficiency and preserved signal-to-noise ratio (SNR) as measured by the inverse g-factor. Then, we verified the efficacy of our approach on high-dimensional non-Cartesian 3D cones datasets, where coil sketching yielded up to three-fold faster reconstructions with equivalent image quality. Coil sketching is a general and versatile reconstruction framework for computationally fast and memory-efficient reconstruction.

The effects of various penalty parameter values in Q.Clear algorithm for rectal cancer detection on 18F-FDG images using a BGO-based PET/CT scanner: a phantom and clinical study

BackgroundThe Q.Clear algorithm is a fully convergent iterative image reconstruction technique. We hypothesize that different PET/CT scanners with distinct crystal properties will require different optimal settings for the Q.Clear algorithm. Many studies have investigated the improvement of the Q.Clear reconstruction algorithm on PET/CT scanner with LYSO crystals and SiPM detectors. We propose an optimum penalization factor (β) for the detection of rectal cancer and its metastases using a BGO-based detector PET/CT system which obtained via accurate and comprehensive phantom and clinical studies.Methods18F-FDG PET-CT scans were acquired from NEMA phantom with lesion-to-background ratio (LBR) of 2:1, 4:1, 8:1, and 15 patients with rectal cancer. Clinical lesions were classified into two size groups. OSEM and Q.Clear (β value of 100–500) reconstruction was applied. In Q.Clear, background variability (BV), contrast recovery (CR), signal-to-noise ratio (SNR), SUVmax, and signal-to-background ratio (SBR) were evaluated and compared to OSEM.ResultsOSEM had 11.5–18.6% higher BV than Q.Clear using β value of 500. Conversely, RC from OSEM to Q.Clear using β value of 500 decreased by 3.3–7.7% for a sphere with a diameter of 10 mm and 2.5–5.1% for a sphere with a diameter of 37 mm. Furthermore, the increment of contrast using a β value of 500 was 5.2–8.1% in the smallest spheres compared to OSEM. When the β value was increased from 100 to 500, the SNR increased by 49.1% and 30.8% in the smallest and largest spheres at LBR 2:1, respectively. At LBR of 8:1, the relative difference of SNR between β value of 100 and 500 was 43.7% and 44.0% in the smallest and largest spheres, respectively. In the clinical study, as β increased from 100 to 500, the SUVmax decreased by 47.7% in small and 31.1% in large lesions. OSEM demonstrated the least SUVmax, SBR, and contrast. The decrement of SBR and contrast using OSEM were 13.6% and 12.9% in small and 4.2% and 3.4%, respectively, in large lesions.ConclusionsImplementing Q.Clear enhances quantitative accuracies through a fully convergent voxel-based image approach, employing a penalization factor. In the BGO-based scanner, the optimal β value for small lesions ranges from 200 for LBR 2:1 to 300 for LBR 8:1. For large lesions, the optimal β value is between 400 for LBR 2:1 and 500 for LBR 8:1. We recommended β value of 300 for small lesions and β value of 500 for large lesions in clinical study.

Taking into account a priori information in the iterative reconstruction of images of foundry products

Methods of restoring images and properties of non-destructive testing objects based on solving inverse problems (problems of restoring distribution functions of unknown characteristics of an object based on the results of indirect measurements) are considered. Management methods are based on solving inverse problems and allow you to get the most complete information about the distributed properties of an object. The need to attract additional information imposes serious restrictions on the development of universal applied algorithms for solving incorrectly set tasks. As a rule, individual additional information is available for each specific non-destructive testing task. An effective numerical algorithm for solving an incorrectly posed problem should be focused on taking this information into account at each stage of the solution search. When solving an applied problem, it is also necessary that the algorithm corresponds to both the measuring capabilities and the capabilities of available computing tools. The problem of low-projection X-ray tomography is always associated with a lack of initial data and can only be solved using a priori information. To introduce the necessary additional information into the numerical algorithm, the methods of iterative reconstruction of tomographic images are identified as the most suitable. One of the approaches to the presentation of this kind of information is described. A practical solution to this problem will expand the scope of the X-ray tomography method.

Deep learning imaging reconstruction of reduced-dose 40 keV virtual monoenergetic imaging for early detection of colorectal cancer liver metastases

ObjectiveTo explore whether reduced-dose (RD) gemstone spectral imaging (GSI) and deep learning image reconstruction (DLIR) of 40 keV virtual monoenergetic image (VMI) enhanced the early detection and diagnosis of colorectal cancer liver metastases (CRLM). MethodsThirty-five participants with pathologically confirmed colorectal cancer were prospectively enrolled from March to August 2022 after routine care abdominal computed tomography (CT). GSI mode was used for contrast-enhanced CT, and two portal venous phase CT images were obtained [standard-dose (SD) CT dose index (CTDIvol) = 15.51 mGy, RD CTDIvol = 7.95 mGy]. The 40 keV-VMI were reconstructed via filtered back projection (FBP) and iterative reconstruction (ASIR-V 60 %, AV60) of both SD and RD images. RD medium-strength deep learning image reconstruction (DLIR-M) and RD high-strength deep learning image reconstruction (DLIR-H) were used to reconstruct the 40 keV-VMI. The contrast-to-noise ratio (CNR) and signal-to-noise ratio (SNR) of the liver and the lesions were objectively evaluated. The overall image quality, lesion conspicuity, and diagnostic confidence were subjectively evaluated, to compare the differences in evaluation results among the different images. ResultsAll 35 participants (mean age: 59.51 ± 11.01 years; 14 females) underwent SD and RD GSI portal venous-phase CT scans. The dose-length product of the RD GSI scan was reduced by 49–53 % lower than that of the SD GSI scan (420.22 ± 31.95) vs (817.58 ± 60.56). A total of 219 lesions were identified, including 55 benign lesions and 164 metastases, with an average size of 7.37 ± 4.14 mm. SD-FBP detected 207 lesions, SD-AV60 detected 201 lesions, and DLIR-M and DLIR-H detected 199 and 190 lesions, respectively. For lesions ≤ 5 mm, there was no statistical difference between SD-FBP vs DLIR-M (χ2McNemar = 1.00, P = 0.32) and SD-AV60 vs DLIR-M (χ2McNemar = 0.33, P = 0.56) in the detection rate. The CNR, SNR, and noise of DLIR-M and DLIR-H 40 keV-VMI images were better than those of SD-FBP images (P < 0.01) but did not differ significantly from those of SD-AV60 images (P > 0.05). When the lesions ≤ 5 mm, there were statistical differences in the overall diagnostic sensitivity of lesions compared with SD-FBP, SD-AV60, DLIR-M and DLIR-H (P＜0.01). There were no statistical differences in the sensitivity of lesions diagnosis between SD-FBP, SD-AV60 and DLIR-M (both P＞0.05). However, the DLIR-M subjective image quality and lesion diagnostic confidence were higher for SD-FBP (both P < 0.01). ConclusionReduced dose DLIR-M of 40 keV-VMI can be used for routine follow-up care of colorectal cancer patients, to optimize evaluations and ensure CT image quality. Meanwhile, the detection rate and diagnostic sensitivity and specificity of small lesions, early liver metastases is not obviously reduced.

A quantitative CBCT pipeline based on 2D antiscatter grid and grid-based scatter sampling for image-guided radiation therapy.

Quantitative accuracy is critical for expanding the role of cone beam CT (CBCT) imaging from target localization to quantitative treatment monitoring and plan adaptations in radiation therapy. Despite advances in CBCT image quality improvement methods, quantitative accuracy gap between CBCT and multi-detector CT (MDCT) remains. In this work, a physics-driven approach was investigated that combined robust scatter rejection, raw data correction and iterative image reconstruction to further improve CBCT image quality and quantitative accuracy, referred to as quantitative CBCT (qCBCT). QCBCT approach includes tungsten 2D antiscatter grid hardware, residual scatter correction with grid-based scatter sampling, image lag, and beam hardening correction for offset detector geometry linac-mounted CBCT. Images were reconstructed with iterative image reconstruction to reduce image noise. qCBCT was evaluated using a variety of phantoms to investigate the effect of object size and its composition on image quality, and image quality was benchmarked against clinical CBCT and gold standard MDCT images used for treatment planning. QCBCT provided statistically significant improvement in CT number accuracy and reduced image artifacts when compared to clinical CBCT images. When compared to gold standard MDCT, mean HU errors in qCBCT and clinical CBCT were 17±9 and 38±29 HU, respectively. Magnitude of phantom size dependent HU variations were comparable between MDCT and qCBCT images. With iterative reconstruction, contrast-to-noise ratio improved by 25% when compared to clinical CBCT protocols. Combination of novel scatter suppression techniques and other data correction methods in qCBCT provided CT number accuracy comparable to gold standard MDCT used for treatment planning. This approach may potentially improve CBCT's promise in fulfilling the tasks that demand high quantitative accuracy, such as online dose calculations and treatment response assessment, in image guided radiation therapy.