Introduction In this study, the optimal input parameters point spread function (PSF) and the number of iterations of the Richardson–Lucy algorithm were experimentally determined to restore Tc-99 m methyl diphosphonate (MDP) whole-body bone scan images. Materials and methods The experiment was performed on 60 anonymized Tc-99 m MDP whole-body bone scan images. Ten images were used for estimating the optimum value of PSF and the number of iterations to restore scintigraphic images. The remaining 50 images were used for validation of estimated parameters. The image quality of observed and restored images was assessed objectively using blind/referenceless image spatial quality evaluator (BRISQUE), mean brightness (MB), discrete entropy (DE), and edge-based contrast measure (EBCM) image quality metrics. Image quality was subjectively assessed by two nuclear medicine physicians (NMPs) by comparing the restored image quality with observed image quality and assigning a score to each image on the scale of 0–5. Results Based on BRISQUE, MB, DE, and EBCM scores, the restored images were significantly sharper, less bright, had more detailed information, and had less contrast around edges compared to the input images. The restored images had improved resolution based on visual assessment as well; NMPs assigned an average image quality score of 4.00 to restored images. Maximum resolution enhancement was noticed at PSF (size: 11 pixels, sigma: 1.75 pixels) and the number of iterations = 10. With the increase in the number of iterations, noise also gets amplified along with resolution enhancement and affects the detectability of small lesions; in the case of relatively low noisy input images, the number of iterations = 5 gave better results. Conclusion Tc-99 m MDP bone scan images were restored to improve image quality using the Richardson–Lucy algorithm. The optimum value of the PSF parameter was found to be of size = 11 pixels and sigma = 1.75 pixels.

Introduction A DnCNN for image denoising trained with natural images is available in MATLAB. For Tc-99m DMSA images, any loss of clinical details during the denoising process will have serious consequences since denoised image is to be used for diagnosis. The objective of the study was to find whether this pre-trained DnCNN can be used for denoising Tc-99m DMSA images and compare its performance with block matching 3D (BM3D) filter. Materials and methods Two hundred forty-two Tc-99m DMSA images were denoised using BM3D filter (at sigma = 5, 10, 15, 20, and 25) and DnCNN. The original and denoised images were reviewed by two nuclear medicine physicians and also assessed objectively using the image quality metrics: SSIM, FSIM, MultiSSIM, PIQE, Blur, GCF, and Brightness. Wilcoxon signed-rank test was applied to find the statistically significant difference between the value of image quality metrics of the denoised images and the corresponding original images. Results Nuclear medicine physicians observed no loss of clinical information in DnCNN denoised image and superior image quality compared to its original and BM3D denoised images. Edges/boundaries of the scar were found to be well preserved, and doubtful scar became obvious in the denoised image. Objective assessment also showed that the quality of DnCNN denoised images was significantly better than that of original images at P-value <0.0001. Conclusion The pre-trained DnCNN available with MATLAB Deep Learning Toolbox can be used for denoising Tc-99m DMSA images, and the performance of DnCNN was found to be superior in comparison with BM3D filter.

Introduction: In this pilot study, we have proposed and evaluated pipelined application of the dynamic stochastic resonance (DSR) algorithm and block-matching 3D (BM3D) filter for the enhancement of nuclear medicine images. The enhanced images out of the pipeline were compared with the corresponding enhanced images obtained using individual applications of DSR and BM3D algorithm. Materials and Methods: Twenty 99m-Tc MDP bone scan images acquired on SymbiaT6 SPECT/CT gamma camera system fitted with low-energy high-resolution collimators were exported in DICOM format to a personal computer and converted into PNG format. These PNG images were processed using the proposed algorithm in MATLAB . Two nuclear medicine physicians visually compared each input and its corresponding three enhanced images to select the best-enhanced image. The image quality metrics ( Brightness , Global Contrast Factor (GCF) , Contrast per pixel (CPP), and Blur ) were used to assess the image quality objectively. The Wilcoxon signed test was applied to find a statistically significant difference in Brightness , GCF, CPP, and Blur of enhanced and its input images at a level of significance. Results: Images enhanced using the pipelined application of SR and BM3D were selected as the best images by both nuclear medicine physicians. Based on Brightness , Global Contrast Factor (GCF), CPP, and Blur , the image quality of our proposed pipeline was significantly better than enhanced images obtained using individual applications of DSR and BM3D algorithm. The proposed method was found to be very successful in enhancing details in the low count region of input images. The enhanced images were bright, smooth, and had better target-to-background ratio compared to input images. Conclusion: The pipelined application of DSR and BM3D algorithm produced enhancement in nuclear medicine images having following characteristics: bright, smooth, better target-to-background ratio, and improved visibility of details in the low count regions of the input image, as compared to individual enhancements by application of DSR or BM3D algorithm.

Introduction: Wavelet transforms of an image result in set of wavelet coefficients. Thresholding eliminates insignificant coefficients while retaining the significant ones (resulting in matrix having few nonzero elements that need to be stored). The compressed image is reconstructed by applying inverse wavelet transform. The quality of compressed image deteriorates with increase in compression. Hence, finding optimum value of scale and threshold is a challenging task. The objective of the study was to find the optimum value of scale and threshold for compressing 99mTc-methylene diphosphonate (99 mTc-MDP) bone scan images using Haar wavelet transform. Materials and Methods: Haar wavelet transform at scale 1–8 was applied on 106 99 mTc-MDP whole-body bone scan images, and wavelet coefficients were threshold at 90, 95, 97, and 99 percentiles, followed by inverse wavelet transform to get 3392 compressed images. Nuclear medicine physician (NMP) compared compressed image with its corresponding input to label it as acceptable or unacceptable. The values of scale and threshold that resulted in majority of acceptable images were considered to be optimum. The quality of compressed image was also evaluated using perception image quality evaluator (PIQE) image quality metrics. Compression ratio was calculated by dividing the number of nonzero elements after thresholding wavelet coefficients by the number of nonzero elements in Haar decomposed matrix. Results: NMP found quality of compressed images (obtained at scale 2 and 90 percentile threshold) identical to the quality of the corresponding input images. As per PIQE score, quality of compressed images was perceptually better than that of the corresponding input images. Conclusions: The optimum values of scale and threshold were determined to be 2 and 90 percentiles, respectively.

Objective: The objective of the study was to develop a Personal Computer (PC) based tool to estimate the center of rotation (COR) offsets from COR projection datasets using methods mentioned in IAEA-TECDOC-602. Materials and Methods: Twenty-four COR studies were acquired on Discovery NM 630 Dual head gamma camera fitted with parallel hole collimator, and COR offsets were estimated with the software available at the terminal for processing a COR study. These COR projection images were exported in DICOM. A MATLAB script (software program) was written to estimate COR offset using Method A (using opposite pair of projections) and Method B (using curve fitting method) as mentioned in IAEA-TECDOC-602. Our program read the COR study (in DICOM) and estimated COR offsets using Method A and Method B. The accuracy of the program was verified using simulated projection dataset of a point source object acquired at 6° interval in the range of 0°–360° angle. Bland Altman plot was used for analyzing the agreement between the COR offsets estimated using (1) Method A and Method B mentioned in IAEA-TECDOC-602, and (2) Our program and vendor program available at Discovery NM 630 acquisition terminal. Results: On simulated data, offset from center of gravity (COG) in X direction (COGX) and COG in Y direction (COGY) estimated using Method A was constant (same) at each pair of angles while using Method B, it was found to be in the range (−2 × 10 −10 , 1 × 10 −10) which is negligible. Most of the differences (23 out of 24) between the result of Method A and Method B, and between the results of our program and vendor program was found to be within 95% confidence interval (mean ± 1.96 standard deviation). Conclusions: Our PC-based tool to estimate COR offsets from COR projection datasets using methods mentioned in IAEA-TECDOC-602 was found to be accurate and provides results in agreement with vendor's program. It can be used as an independent tool to estimate COR offset for standardization and calibration purposes.