Search

Intravoxel incoherent motion (IVIM) is a type of diffusion-weighted imaging (DWI), and IVIM model parameters (water molecule diffusion rate Dt, pseudo diffusion coefficient Dp, and tissue perfusion fraction Fp) have been widely used in the diagnosis and characterization of malignant lesions.This study proposes a deep-learning model with synthetic-to-real domain adaptation to fit the IVIM model parameters of DWI.Ninety-eight consecutive patients diagnosed with hepatocellular carcinoma between January 2017 and September 2020 were included in the study, and routine IVIM-DWI serial examinations were performed using a 3.0 T magnetic resonance imaging system in preoperative MR imaging. The proposed method is mainly composed of two modules: a convolutional neural network-based IVIM model fitting network to map b-value images to the IVIM parameter maps and a domain discriminator to improve the accuracy of the IVIM parameter maps in the real data. The proposed method was compared with previously reported fitting methods, including the nonlinear least squares (NLS), IVIM-NETThe DWI reconstruction performance demonstrates that the proposed method has better reconstruction accuracy for DWI with a low signal-to-noise ratio, which implies that the proposed method improves the fitting accuracy of the IVIM parameters. Noise-corrupt experiments show that the proposed method is more robust against noise-corrupted signals. With the proposed method, no outliers were found in Dt, and outliers were reduced for Fp in the abnormal regions (proposed method: 1.85%; NLS: 5.90%; IVIM-NETIVIM parameters can be estimated using a synthetic-to-real domain-adaptation framework with deep learning, and the proposed method outperforms previously reported methods. This article is protected by copyright. All rights reserved.

Purpose: In transmitted X-ray tomography imaging, projections are sometimes incomplete due to a variety of reasons, such as geometry inaccuracy, defective detector cells, etc. To address this issue, we have derived a direct consistency condition based on John's Equation, and proposed a method to effectively restore incomplete projections based on this consistency condition. Methods: Through parameter substitutions, we have derived a direct consistency condition equation from John's equation, in which the left side is only projection derivative of view and the right side is projection derivative of other geometrical parameters. Based on this consistency condition, a projection restoration method is proposed, which includes five steps: 1) Forward projecting reconstructed image and using linear interpolation to estimate the incomplete projections as the initial result; 2) Performing Fourier transform on the projections; 3) Restoring the incomplete frequency data using the consistency condition equation; 4) Performing inverse Fourier transform; 5) Repeat step 2)∼4) until our criteria is met to terminate the iteration. Results: A beam-blocking-based scatter correction case and a bad-pixel correction case were used to demonstrate the efficacy and robustness of our restoration method. The mean absolute error (MAE), signal noise ratio (SNR) and mean square error (MSE) were employed as our evaluation metrics of the reconstructed images. For the scatter correction case, the MAE is reduced from 63.3% to 71.7% with 4 iterations. Compared with the existing Patch's method, the MAE of our method is further reduced by 8.72%. For the bad-pixel case, the SNR of the reconstructed image by our method is increased from 13.49% to 21.48%, with the MSE being decreased by 45.95%, compared with linear interpolation method. Conclusion: Our studies have demonstrated that our restoration method based on the new consistency condition could effectively restore the incomplete projections, especially for their high frequency component.

Purpose: Sparse-view computed tomography (CT) reconstruction is an effective strategy to reduce the radiation dose delivered to patients. Due to its insufficiency of measurements, traditional non-local means (NLM) based reconstruction methods often lead to over-smoothness in image edges. To address this problem, an adaptive NLM reconstruction method based on rotational invariance (RIANLM) is proposed. Methods: The method consists of four steps: 1) Initializing parameters; 2) Algebraic reconstruction technique (ART) reconstruction using raw projection data; 3) Positivity constraint of the image reconstructed by ART; 4) Update reconstructed image by using RIANLM filtering. In RIANLM, a novel similarity metric that is rotational invariance is proposed and used to calculate the distance between two patches. In this way, any patch with similar structure but different orientation to the reference patch would win a relatively large weight to avoid over-smoothed image. Moreover, the parameter h in RIANLM which controls the decay of the weights is adaptive to avoid over-smoothness, while it in NLM is not adaptive during the whole reconstruction process. The proposed method is named as ART-RIANLM and validated on Shepp-Logan phantom and clinical projection data. Results: In our experiments, the searching neighborhood size is set to 15 by 15 and the similarity window is set to 3 by 3. For the simulated case with a resolution of 256 by 256 Shepp-Logan phantom, the ART-RIANLM produces higher SNR (35.38dB

Purpose: For traditional geometric calibration, the calibration accuracy relies on both accuracy of geometry phantom manufacture and accuracy of ball bearings (BB) location estimation. In this work, we have developed a method to perform phantom calibration and geometry calibration iteratively and accurately in a whole procedure. Methods: We have designed and manufactured a geometry phantom with BB and an evaluation phantom of a crystal ball contained in a cubic gel box. Our calibration method consists of five steps: 1) Estimate BB locations using spiral CT image, which are then used to initialize the particles in Particle Swarm Optimization (PSO) algorithm; 2) Perform geometric calibration; 3) Reconstruct the images of the evaluation phantom based on the current geometry calibration; 4) Evaluate the reconstructed images using a geometry evaluate index; 5) Update BB locations in PSO algorithm. Repeat step2)-5) until our stopping criteria is met. The edge of the crystal ball in the calibration phantom on CBCT images is detected by Hough transform to define two circular rings outside and inside the ball. The evaluation index used in step 4) is defined as the difference of the averaged image intensities of these two circular rings. Results: We have demonstrated the feasibility and performance of our method on a benchtop CBCT system. It is observed that inaccurate BB locations lead to severe image distortion and relative small evaluation index. With our method, streak artifacts are reduced and the structure becomes sharper and clearer. The evaluation index increases fast within 10 iterations, and then becomes stable gradually. Conclusion: Our method can perform accurate phantom calibration and geometry calibration together in a whole procedure. It helps to mitigate the impact of the geometry phantom manufacture errors on the calibration, which could hence save the cost of the geometry phantom manufacture.

Purpose: Image artifacts are usually evaluated qualitatively via visual observation of the reconstructed images, which is susceptible to subjective factors due to the lack of an objective evaluation criterion. In this work, we propose a Helgason-Ludwig consistency condition (HLCC) based evaluation method to quantify the severity level of different image artifacts in dental CBCT. Methods: Our evaluation method consists of four step: 1) Acquire Cone beam CT(CBCT) projection; 2) Convert 3D CBCT projection to fan-beam projection by extracting its central plane projection; 3) Convert fan-beam projection to parallel-beam projection utilizing sinogram-based rebinning algorithm or detail-based rebinning algorithm; 4) Obtain HLCC profile by integrating parallel-beam projection per view and calculate wave percentage and variance of the HLCC profile, which can be used to describe the severity level of image artifacts. Results: Several sets of dental CBCT projections containing only one type of artifact (i.e. geometry, scatter, beam hardening, lag and noise artifact), were simulated using gDRR, a GPU tool developed for efficient, accurate, and realistic simulation of CBCT Projections. These simulated CBCT projections were used to test our proposed method. HLCC profile wave percentage and variance induced by geometry distortion are about 3∼21 times and 16∼393 times as large as that of the artifact-free projection, respectively. The increase factor of wave percentage and variance are 6 and133 times for beam hardening, 19 and 1184 times for scatter, and 4 and16 times for lag artifacts, respectively. In contrast, for noisy projection the wave percentage, variance and inconsistency level are almost the same with those of the noise-free one. Conclusion: We have proposed a quantitative evaluation method of image artifact based on HLCC theory. According to our simulation results, the severity of different artifact types is found to be in a following order: Scatter>Geometry>Beam hardening>Lag>Noise>Artifact-free in dental CBCT

Purpose: In computed tomography (CT) system, CT images with ring artifacts will be reconstructed when some adjacent bins of detector don’t work. The ring artifacts severely degrade CT image quality. We present a useful CT ring artifacts reduction based on projection data correction, aiming at estimating the missing data of projection data accurately, thus removing the ring artifacts of CT images. Methods: The method consists of ten steps: 1) Identification of abnormal pixel line in projection sinogram; 2) Linear interpolation within the pixel line of projection sinogram; 3) FBP reconstruction using interpolated projection data; 4) Filtering FBP image using mean filter; 5) Forwarding projection of filtered FBP image; 6) Subtraction forwarded projection from original projection; 7) Linear interpolation of abnormal pixel line area in the subtraction projection; 8) Adding the interpolated subtraction projection on the forwarded projection; 9) FBP reconstruction using corrected projection data; 10) Return to step 4 until the pre-set iteration number is reached. The method is validated on simulated and real data to restore missing projection data and reconstruct ring artifact-free CT images. Results: We have studied impact of amount of dead bins of CT detector on the accuracy of missing data estimation in projection sinogram. For the simulated case with a resolution of 256 by 256 Shepp-Logan phantom, three iterations are sufficient to restore projection data and reconstruct ring artifact-free images when the dead bins rating is under 30%. The dead-bin-induced artifacts are substantially reduced. More iteration number is needed to reconstruct satisfactory images while the rating of dead bins increases. Similar results were found for a real head phantom case. Conclusion: A practical CT image ring artifact correction scheme based on projection data is developed. This method can produce ring artifact-free CT images feasibly and effectively.