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.