A deep learning approach for converting prompt gamma images to proton dose distributions: A Monte Carlo simulation study

PurposeIn proton therapy, imaging prompt gamma (PG) rays has the potential to verify proton dose (PD) distribution. Despite the fact that there is a strong correlation between the gamma-ray emission and PD, they are still different in terms of the distribution and the Bragg peak (BP) position. In this work, we investigated the feasibility of using a deep learning approach to convert PG images to PD distributions.MethodsWe designed the Monte Carlo simulations using 20 digital brain phantoms irradiated with a 100-MeV proton pencil beam. Each phantom was used to simulate 200 pairs of PG images and PD distributions. A convolutional neural network based on the U-net architecture was trained to predict PD distributions from PG images.ResultsOur simulation results show that the pseudo PD distributions derived from the corresponding PG images agree well with the simulated ground truths. The mean of the BP position errors from each phantom was less than 0.4 mm. We also found that 2000 pairs of PG images and dose distributions would be sufficient to train the U-net. Moreover, the trained network could be deployed on the unseen data (i.e. different beam sizes, proton energies and real patient CT data).ConclusionsOur simulation study has shown the feasibility of predicting PD distributions from PG images using a deep learning approach, but the reliable prediction of PD distributions requires high-quality PG images. Image-degrading factors such as low counts and limited spatial resolution need to be considered in order to obtain high-quality PG images.     

Therapy region monitoring based on PET using 478 keV single prompt gamma ray during BNCT: A Monte Carlo simulation study

We confirmed the feasibility of using our proposed system to extract two different kinds of functional images from a positron emission tomography (PET) module by using an insertable collimator during boron neutron capture therapy (BNCT). Coincidence events from a tumor region that included boron particles were identified by a PET scanner before BNCT; subsequently, the prompt gamma ray events from the same tumor region were collected after exposure to an external neutron beam through an insertable collimator on the PET detector. Five tumor regions that contained boron particles and were located in the water phantom and in the BNCT system with the PET module were simulated with Monte Carlo simulation code. The acquired images were quantitatively analyzed. Based on the receiver operating characteristic (ROC) curves in the five boron regions, A, B, C, D, and E, the PET and single-photon images were 10.2%, 11.7%, 8.2% (center region), 12.6%, and 10.5%, respectively. We were able to acquire simultaneously PET and single prompt photon images for tumor regions monitoring by using an insertable collimator without any additional isotopes.     

Comprehensive quality assurance phantom for the small animal radiation research platform (SARRP)

PurposeTo develop and test the suitability and performance of a comprehensive quality assurance (QA) phantom for the Small Animal Radiation Research Platform (SARRP).Methods and materialsA QA phantom was developed for carrying out daily, monthly and annual QA tasks including: imaging, dosimetry and treatment planning system (TPS) performance evaluation of the SARRP. The QA phantom consists of 15 (60 × 60 × 5 mm³) kV-energy tissue equivalent solid water slabs. The phantom can incorporate optically stimulated luminescence dosimeters (OSLD), Mosfet or film. One slab, with inserts and another slab with hole patterns are particularly designed for image QA.ResultsOutput constancy measurement results showed daily variations within 3%. Using the Mosfet in phantom as target, results showed that the difference between TPS calculations and measurements was within 5%. Annual QA results for the Percentage depth dose (PDD) curves, lateral beam profiles, beam flatness and beam profile symmetry were found consistent with results obtained at commissioning. PDD curves obtained using film and OSLDs showed good agreement. Image QA was performed monthly, with image-quality parameters assessed in terms of CBCT image geometric accuracy, CT number accuracy, image spatial resolution, noise and image uniformity.ConclusionsThe results show that the developed QA phantom can be employed as a tool for comprehensive performance evaluation of the SARRP. The study provides a useful reference for development of a comprehensive quality assurance program for the SARRP and other similar small animal irradiators, with proposed tolerances and frequency of required tests.     

Distal edge determination precision for a multi-slat prompt-gamma camera: A comprehensive simulation and optimization of the detection system

Multi-slat prompt-gamma camera is a promising tool for range monitoring during proton therapy. We report the results of a comprehensive simulation study analyzing the precision which is possible to reach with this camera in determination of the position of the distal edge of the Bragg peak. For the first time we include simulation of optical photons. The proton beam (single pencil beam, 130 MeV, 10 ns bunch period, total of 1·10⁸ protons) is interacting with a polymethyl methacrylate (PMMA) phantom, which is a cylinder of 200 mm in diameter and length. The prompt gamma rays generated in the phantom are collimated with a multi-slat collimator and detected using a combination of yttrium aluminum perovskite (YAP) scintillators, installed in the collimator apertures, and light sensors. Two scintillator packing schemes, with one and with two scintillator plates per aperture, are considered. The collimator configuration (the septal thickness, aperture and height), resulting in the best precision, is determined using two methods of detector optimization. Precision of 2.1 mm (full width at half maximum) in the edge position determination is demonstrated.     

Monte Carlo characterizations mapping of the (γ,n) and (n,γ) photonuclear reactions in the high energy X-ray radiation therapy

Aim

The aim of this work was to map the characteristics of (n,γ) and (γ,n) reactions in a high energy photon radiation therapy.

Background

Photoneutrons produced in the high energy X-Ray radiation therapy may damage patients and staff. It is due to high RBE of the produced neutrons according to their energy and isotropic emission. Characterization of the photoneutrons can help us in appropriate shielding.

Materials and methods

This study focused on the photoneutron and capture gamma ray phenomena. Characteristics such as dose value, fluence and spectra of both the neutrons and the by produced prompt gamma ray were described.

Results and discussion

Neutron and prompt gamma spectra in different points showed the neutrons to be thermalized when increasing the distance from the linac. Energy of the neutrons changed from about 0.6 MeV at the isocentre to around 10⁻⁰⁸ MeV at the outer door position. Although the neutrons were found as fast neutrons, their spectra showed they were thermal neutrons at the outer door position. Additionally, it was seen that the energy of the gamma rays is higher than the scattered X-ray energy. The energy of gamma rays was seen to be up to 10 MeV while the linac photons had energy lower than 1 MeV. Neutron source strength obtained in this work was in good agreement with the published data, which may be a confirmation of our simulation accuracy.

Conclusion

The study showed that the Monte Carlo simulation can be applied in the radiotherapy and industrial radiation works as a useful and precise estimator. We also concluded that the dose from the prompt gamma ray at the outer door location is higher than the scattered radiation from the linac and should be considered in the shielding.     

Prompt-gamma emission in GEANT4 revisited and confronted with experiment

PurposeThe field of online monitoring of the beam range is one of the most researched topics in proton therapy over the last decade. The development of detectors that can be used for beam range verification under clinical conditions is a challenging task. One promising possible solution are modalities that record prompt-gamma radiation produced by the interactions of the proton beam with the target tissue. A good understanding of the energy spectra of the prompt gammas and the yields in certain energy regions is crucial for a successful design of a prompt-gamma detector. Monte-Carlo simulations are an important tool in development and testing of detector concepts, thus the proper modelling of the prompt-gamma emission in those simulations are of vital importance. In this paper, we confront a number of GEANT4 simulations of prompt-gamma emission, performed with different versions of the package and different physics lists, with experimental data obtained from a phantom irradiation with proton beams of four different energies in the range 70–230 MeV.MethodsThe comparison is made on different levels: features of the prompt-gamma energy spectrum, gamma emission depth profiles for discrete transitions and the width of the distal fall-off in those profiles.ResultsThe best agreement between the measurements and the simulations is found for the GEANT4 version 10.4.2 and the reference physics list QGSP_BIC_HP.ConclusionsModifications to prompt-gamma emission modelling in higher versions of the software increase the discrepancy between the simulation results and the experimental data.     

A new analytical formula for neutron capture gamma dose calculations in double-bend mazes in radiation therapy

BackgroundPhotoneutrons are produced in radiation therapy with high energy photons. Also, capture gamma rays are the byproduct of neutrons interactions with wall material of radiotherapy rooms.AimIn the current study an analytical formula was proposed for capture gamma dose calculations in double bend mazes in radiation therapy rooms.Materials and methodsA total of 40 different layouts with double-bend mazes and a 18 MeV photon beam of Varian 2100 Clinac were simulated using MCNPX Monte Carlo (MC) code. Neutron capture gamma ray dose equivalent was calculated by the MC method along the maze and at the maze entrance door of all the simulated rooms. Then, all MC resulted data were fitted to an empirical formula for capture gamma dose calculations. Wu–McGinley analytical formula for capture gamma dose equivalent at the maze entrance door in single-bend mazes was also used for comparison purposes.ResultsFor capture gamma dose equivalents at the maze entrance door, the difference of 2–11% was seen between MC and the derived equation, while the difference of 36–87% was found between MC and the Wu–McGinley methods.ConclusionOur results showed that the derived formula results were consistent with the MC results for all of 40 different geometries. However, as a new formula, further evaluations are required to validate its use in practical situations. Finally, its application is recommend for capture gamma dose calculations in double-bend mazes to improve shielding calculations.     

Feasibility study on the use of 230 MeV proton cyclotron in proton therapy centers as a spallation neutron source for BNCT

AimThe feasibility of using 230 MeV proton cyclotrons in proton therapy centers as a spallation neutron source for Boron Neutron Capture Therapy (BNCT) was investigated.BackgroundBNCT is based on the neutron irradiation of a ¹⁰B-containing compound located selectively in tumor cells. Among various types of neutron generators, the spallation neutron source is a unique way to generate high-energy and high-flux neutrons.Materials and MethodsNeutron beam was generated by a proton accelerator via spallation reactions and then the produced neutron beam was shaped to be appropriate for BNCT. The proposed Beam Shaping Assembly (BSA) consists of different moderators, a reflector, a collimator, as well as thermal and gamma filters. In addition, the simulated Snyder head phantom was utilized to evaluate the dose distribution in tumor and normal tissue due to the irradiation by the designed beam. MCNPX2.6 Monte Carlo code was used to optimize BSA as well as evaluate dose evaluation.ResultsA BSA was designed. With the BSA configuration and a beam current of 104 nA, epithermal neutron flux of 3.94 × 10⁶ [n/cm²] can be achieved, which is very low. Provided that we use the beam current of 5.75 μA, epithermal neutron flux of 2.18 × 10⁸ [n/cm²] can be obtained and the maximum dose of 38.2 Gy-eq can be delivered to tumor tissue at 1.4 cm from the phantom surface.ConclusionsResults for 230 MeV protons show that with proposed BSA, proton beam current about 5.75 μA is required for this purpose.     

Characterization of prompt gamma-ray emission with respect to the Bragg peak for proton beam range verification: A Monte Carlo study

In this paper we report a Geant4 simulation study to investigate the characteristic prompt gamma (PG) emission in a water phantom for real-time monitoring of the Bragg peak (BP) during proton beam irradiation. The PG production, emission spatial correlation with the BP, and position preference for detection with respect to the BP have been quantified in different PG energy windows as a function of proton pencil-beam energy from 100 to 200 MeV. The PG response to small BP shifts was evaluated using a 2 cm-thick slab with different human body materials embedded in a water phantom. Our results show that the prominent characteristic PG emissions of 4.44, 5.21 and 6.13 MeV exhibit distinctive correlation with the dose deposition curve. The accuracy in BP position identification using these characteristic PG rays is highly consistent as the beam energy increases from 100 to 200 MeV. There exists a position preference for PG detection with respect to the BP position, which has a strong dependence on the proton beam energy and PG energies. It was also observed that a submillimeter shift of the BP position can be realized by using PG signals. These results indicate that the characteristic PG signal is sensitive and reliable for BP tracking. Although the maximization of the PG measurement associated with the BP is difficult, it can be optimized with energy and detection position preferences.     

Use of uniform shots for robust planning of mask-based treatment in Gamma Knife Icon

PurposeTo verify whether Icon automatic correction is robust in preserving plan quality.Materials/methodsAn end-to-end phantom was used to verify Icon’s correction accuracy qualitatively. For quantitative assessment, two plans, a composite- and a uniform-shot-only, were created for an elliptical- (E) and a sausage-shaped (S) lesion inside a PseudoPatient head phantom with a film insert. The phantom was irradiated in the planned and three other positions under each plan: 14° pitch (B); 14° rotation + 8° pitch (C); 95° rotation + 4-cm shift (D).ResultsIcon accurately corrects the locations of the shots. For the uniform-shot plans: all gamma index passing rates were >97%, and the differences between the planned and the delivery doses (minimum, maximum, and mean) were all ≤0.1 Gy. For the composite-shot plans, however, the dose differences increased as the phantom was shifted through positions B-D, with a gamma index passing rate of 61% for lesion-E in position D, and 92%, 79%, and 45% for lesion-S in positions B, C, and D, respectively.ConclusionsPlans using only uniform shots are more robust to deviations in treatment position. The tolerance for such deviations may be lower for plans using composite shots.     

Set-up error validation with EPID images: Measurements vs Egs_cbct simulation

AimIn this study, the egs_cbct code’s ability to replicate an electronic portal imaging device (EPID) is explored.BackgroundWe have investigated head and neck (H&N) setup verification on an Elekta Precise linear accelerator. It is equipped with an electronic portal imaging device (EPID) that can capture a set of projection images over different gantry angles.Methods and materialsCone-beam computed tomography (CBCT) images were reconstructed from projection images of two different setup scenarios. Projections of an Anthropomorphic Rando head phantom were also simulated by using the egs_cbct Monte Carlo code for comparison with the measured projections.Afterwards, CBCT images were reconstructed from this data. Image quality was evaluated against a metric defined as the image acquisition interval (IAI). It determines the number of projection images to be used for CBCT image reconstruction.ResultsFrom this results it was established that phantom shifts could be determined within 2 mm and rotations within one degree accuracy using only 20 projection images (IAI = 10 degrees). Similar results were obtained with the simulated data.ConclusionIn this study it is demonstrated that a head and neck setup can be verified using substantially fewer projection images. Bony landmarks and air cavities could still be observed in the reconstructed Rando head phantom. The egs_cbct code can be used as a tool to investigate setup errors without tedious measurements with an EPID system.     

3D-printed patient-specific pelvis phantom for dosimetry measurements for prostate stereotactic radiotherapy with dominant intraprostatic lesion boost

AimDeveloping and assessing the feasibility of using a three-dimensional (3D) printed patient-specific anthropomorphic pelvis phantom for dose calculation and verification for stereotactic ablative radiation therapy (SABR) with dose escalation to the dominant intraprostatic lesions.Material and methodsA 3D-printed pelvis phantom, including bone-mimicking material, was fabricated based on the computed tomography (CT) images of a prostate cancer patient. To compare the extent to which patient and phantom body and bones overlapped, the similarity Dice coefficient was calculated. Modular cylindrical inserts were created to encapsulate radiochromic films and ionization chamber for absolute dosimetry measurements at the location of prostate and at the boost region. Gamma analysis evaluation with 2%/2mm criteria was performed to compare treatment planning system calculations and measured dose when delivering a 10 flattening filter free (FFF) SABR plan and a 10FFF boost SABR plan.ResultsDice coefficients of 0.98 and 0.91 were measured for body and bones, respectively, demonstrating agreement between patient and phantom outlines. For the boost plans the gamma analysis yielded 97.0% of pixels passing 2%/2mm criteria and these results were supported by the chamber average dose difference of 0.47 ± 0.03%. These results were further improved when overriding the bone relative electron density: 97.3% for the 2%/2mm gamma analysis, and 0.05 ± 0.03% for the ionization chamber average dose difference.ConclusionsThe modular patient-specific 3D-printed pelvis phantom has proven to be a highly attractive and versatile tool to validate prostate SABR boost plans using multiple detectors.     

Evaluation of the effects of implementing a diode transmission device into the clinical workflow

PurposeThis work investigated effects of implementing the Delta⁴ Discover diode transmission detector into the clinical workflow.MethodsPDD and profile scans were completed with and without the Discover for a number of photon beam energies. Transmission factors were determined for all beam energies and included in Eclipse TPS to account for the attenuation of the Discover. A variety of IMRT plans were delivered to a Delta⁴ Phantom+ with and without the Discover to evaluate the Discover’s effects on IMRT QA. An imaging QA phantom was used to assess the detector’s effects on MV image quality. OSLDs placed on the Phantom+ were used to determine the detector’s effects on superficial dose.ResultsThe largest effect on PDDs after d_max was 0.5%. The largest change in beam profile symmetry and flatness was 0.2% and 0.1%, respectively. An average difference in gamma passing rates (2%/2 mm) of 0.2% was observed between plans that did not include the Discover in the measurement and calculation to plans that did include the Discover in the measurement and calculation. The Discover did not significantly change the MV image quality, and the largest observed increase in the relative superficial dose when the Discover was present was 1%.ConclusionsThe effects the Discover has on the linac beam were found to be minimal. The device can be implemented into the clinic without the need to alter the TPS beam modeling, other than accounting for the device’s attenuation. However, a careful workflow review to implement the Discover should be completed.     

Out-of-field doses from secondary radiation produced in proton therapy and the associated risk of radiation-induced cancer from a brain tumor treatment

PurposeTo determine out-of-field doses produced in proton pencil beam scanning (PBS) therapy using Monte Carlo simulations and to estimate the associated risk of radiation-induced second cancer from a brain tumor treatment.MethodsSimulations of out-of-field absorbed doses were performed with MCNP6 and benchmarked against measurements with tissue-equivalent proportional counters (TEPC) for three irradiation setups: two irradiations of a water phantom using proton energies of 78–147 MeV and 177–223 MeV, and one brain tumor irradiation of a whole-body phantom. Out-of-field absorbed and equivalent doses to organs in a whole-body phantom following a brain tumor treatment were subsequently simulated and used to estimate the risk of radiation-induced cancer. Additionally, the contribution of absorbed dose originating from radiation produced in the nozzle was calculated from simulations.ResultsOut-of-field absorbed doses to the TEPC ranged from 0.4 to 135 µGy/Gy. The average deviation between simulations and measurements of the water phantom irradiations was about 17%. The absorbed dose contribution from radiation produced in the nozzle ranged between 0 and 70% of the total dose; the contribution was however small in absolute terms. The absorbed and equivalent doses to the organs ranged between 0.2 and 60 µGy/Gy and 0.5–151 µSv/Gy. The estimated lifetime risk of radiation-induced second cancer was approximately 0.01%.ConclusionsThe agreement of out-of-field absorbed doses between measurements and simulations was good given the sources of uncertainties. Calculations of out-of-field organ doses following a brain tumor treatment indicated that proton PBS therapy of brain tumors is associated with a low risk of radiation-induced cancer.     

Chromosome aberration frequencies produced by a 70-MeV proton beam

The effectiveness of a 70-MeV proton beam in the induction of chromosome aberrations was studied. We employed peripheral lymphocytes and analyzed the frequencies of dicentrics and rings after irradiation at doses ranging from 0.1 to 8.0 Gy at various depths within a Lucite phantom. The frequency of chromosome aberrations after irradiation with an unmodulated proton beam at 5 mm showed a dose-response relationship similar to that of 60Co gamma rays. However, irradiation at greater depths with the spread-out Bragg peak induced higher aberration frequencies at doses lower than those with gamma rays. Furthermore, the distribution curve of chromosome aberration frequencies as a function of depth was found to be slightly different from the physically measured depth-dose curve. With the spread-out Bragg peak the biological effects were more marked at greater depths, resulting in a distribution of relative biological effectiveness values. The results obtained from chromosome aberration analysis may not be related directly to those for the relationship between dose and cell killing. Slight differences in values for relative biological effectiveness due to the change of dose and site of proton beam irradiation may not be important for practical proton beam therapy, but may be important in the prevention of late radiation injuries.     

Exit fluence analysis using portal dosimetry in volumetric modulated arc therapy

AimIn measuring exit fluences, there are several sources of deviations which include the changes in the entrance fluence, changes in the detector response and patient orientation or geometry. The purpose of this work is to quantify these sources of errors.BackgroundThe use of the volumetric modulated arc therapy treatment with the help of image guidance in radiotherapy results in high accuracy of delivering complex dose distributions while sparing critical organs. The transit dosimetry has the potential of Verifying dose delivery by the linac, Multileaf collimator positional accuracy and the calculation of dose to a patient or phantom.Materials and methodsThe quantification of errors caused by a machine delivery is done by comparing static and arc picket fence test for 30 days. A RapidArc plan, created for the pelvis site was delivered without and with Rando phantom and exit portal images were acquired. The day to day dose variation were analysed by comparing the daily exit dose images during the course of treatment. The gamma criterion used for analysis is 3% dose difference and 3 mm distance to agreement with a threshold of 10% of maximum dose.ResultsThe maximum standard deviation for the static and arc picket fence test fields were 0.19 CU and 1.3 CU, respectively. The delivery of the RapidArc plans without a phantom shows the maximum standard deviation of 1.85 CU and the maximum gamma value of 0.59. The maximum gamma value for the RapidArc plan delivered with the phantom was found to be 1.2. The largest observed fluence deviation during the delivery to patient was 5.7% and the maximum standard deviation was 4.1 CU.ConclusionIt is found from this study that the variation due to patient anatomy and interfraction organ motion is significant.     

基于微波热声层析成像定量重建生物组织电导率的改进方法

目的 生物电磁学参数中的电导率与组织的功能性信息直接相关,精准重建生物组织电导率在医学成像技术和医学诊断领域中有着重要意义。本文改进定量微波热声层析成像（microwave-induced thermoacoustic tomography,MTAT）算法,使组织电导率的重建精度提高。方法 本文在利用有限元离散法求解热声波动方程和亥姆霍兹方程的基础之上,提出了一种基于正则化牛顿迭代法（regularized Newton iteration method,RNIM）定量重建组织电导率的改进方法。结果 通过数值模拟实验和含不同浓度NaCl溶液的仿体实验,验证了算法改进的有效性。组织仿体实验结果表明,目标在不同位置、不同大小、不同对比度情况下,相比于定量微波热声层析成像采用拟合（fitting）的方法,采用正则化牛顿法定量重建的仿体电导率相对误差明显降低,重建目标精度提高。在仿体实验中采用RNIM方法重建相同浓度的单目标在不同位置的电导率变化幅度更小,以及重建多目标电导率的相对比值与实际更接近,实验结果验证了改进方法的稳定性。结论 研究结果表明优化算法能更加准确地定量重建组织仿体的电导率,...     

Calibration of MTT assay in proton beams using radiochromic films

PurposeThis study provides methodology of calibrating as well as controlling the output for an MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) colorimetric assay irradiated in a low energy proton beam using EBT3-model GAFCHROMIC^TM film, without correcting for quenching effect.MethodsA calibrated Markus ionization chamber was used to measure the depth dose and beam output for 26.5 MeV protons produced by a CS30 cyclotron. A time-controlled aluminum cylinder was added in front of the horizontal beam-exit serving as a radiation shutter. Following the TRS-398 reference dosimetry protocol for proton beams, the output was calibrated in water at a reference depth of 3 mm. EBT3 film was calibrated for doses up to 8 Gy at the same depth. To verify the dose distribution for each 96-well MTT assay plate, EBT3 film was placed at the reference depth during irradiation and cell doses were scaled by measured percent depth dose (PDD) data.ResultsThe radiochromic film dosimetry system in this study provides dose measurements with an uncertainty better than 3.3% for doses higher than 1 Gy. From a single exposure and utilizing the Gaussian shape of the beam, multiple dose points can be obtained within different wells of the same plate ranging from 6.9 Gy (sigma ∼4%) in the central well, and 2 Gy (sigma ∼8%) for wells positioned closer to the periphery.ConclusionsWe described a methodology for radiochromic film-based dose monitoring system, using low-energy protons, which can be used for the MTT assay in any proton beam, except within Bragg peak region.     

An iterative deconvolution algorithm for image recovery in clinical CT: A phantom study

PurposeTo study the feasibility of using an iterative reconstruction algorithm to improve previously reconstructed CT images which are judged to be non-diagnostic on clinical review. A novel rapidly converging, iterative algorithm (RSEMD) to reduce noise as compared with standard filtered back-projection algorithm has been developed.Materials and methodsThe RSEMD method was tested on in-silico, Catphan^®500, and anthropomorphic 4D XCAT phantoms. The method was applied to noisy CT images previously reconstructed with FBP to determine improvements in SNR and CNR. To test the potential improvement in clinically relevant CT images, 4D XCAT phantom images were used to simulate a small, low contrast lesion placed in the liver.ResultsIn all of the phantom studies the images proved to have higher resolution and lower noise as compared with images reconstructed by conventional FBP. In general, the values of SNR and CNR reached a plateau at around 20 iterations with an improvement factor of about 1.5 for in noisy CT images. Improvements in lesion conspicuity after the application of RSEMD have also been demonstrated. The results obtained with the RSEMD method are in agreement with other iterative algorithms employed either in image space or with hybrid reconstruction algorithms.ConclusionsIn this proof of concept work, a rapidly converging, iterative deconvolution algorithm with a novel resolution subsets-based approach that operates on DICOM CT images has been demonstrated. The RSEMD method can be applied to sub-optimal routine-dose clinical CT images to improve image quality to potentially diagnostically acceptable levels.