Apply deconvolution filter on the image to correct blurring

Hello Dear Developers and Users

I have images acquired from a radiation field using Electronic Portal Imaging Device (EPID) in radiotherapy.
Due to scattered beams in EPID, there was a blurring effect. As expected on the edges, the image quality is reduced.
It seems that deconvolution filters should be used to improve image quality.
For it, I have calculated the EPID response to a pencil beam through the GATE Monte Carlo code. In other words, I have determined the EPID response to a pencil beam (PSF system).
The following images are related to the image acquired with the EPID and the logarithm of the EPID response to a pencil beam.

Specifications of EPID image:
Image Dimensions: 1024 1024 1
Image spacing: 0.25mm 0.25mm 1.00mm
Scalar Type: double
Scalar Range: 530 to 304446
EPID

Specification of scatter kernel:
Image Dimensions: 250 250 1
Image spacing: 0.40mm 0.40mm 1.00mm
Scalar Type: double
Scalar Range: -27 to 0
scatterKernel

I’ve noticed that there are several deconvolution filters in the “itk Simple Filters” module such as:
InverseDeconvolutionImageFilter
LandweberDeconvolutionImageFilter
ProjectedLandweberDeconvolutionImageFilter
RichardsonLucyDeconvolutionImageFilter
TikhonovDeconvolutionImageFilter
WienerDeconvolutionImageFilter

When I use these filters with these image as input volume (with default settings), the output image is badly damaged.
With this image and kenel, what solution you recommend to me?
Please guide me.
Thanks a lot.
Shahrokh

It is expected that the device manufacturer performs all possible image quality improvement filtering that optimizes overall signal-to-noise ratio in the image. This includes deconvolution filtering, if applicable. If you have doubts that the image quality is optimal then report this error to the manufacturer. If you work with images or features that the imaging device was not optimized for then you might be able to achieve further improvements by enhancing those specific features that you are interested in (at the cost of decreasing overall image quality).

Unsharp masking only improves image quality with optimal parameter settings but if the output image is completely off then you may need to convert the pixel type of the input image to float (using Cast Scalar Volume module). Some filters may not work on single-slice volumes, so you may need to resample it to have multiple slices (using Resample Scalar Volume, adjust Spacing parameters to have at least 10 slices). Make sure to also crop any black borders in the input image.

I was looking at some papers about improving the accuracy of segmented models from CT scans. One paper suggested that the slightly blurred edge between bone and soft tissue could be sharpened using a Richardson Lucy Deconvolution filter with some knowledge of the point spread function.

It would be good to experiment with this filter to see if edges could be un-blurred to potentially make more accurate bone surface models for modelling implants.

I was looking at the Richardson Lucy Filter in Simple Filters and was wondering if anyone knows how to use this filter? I couldn’t find much info on the ITK website either. Why are there two volume inputs? Also I was under the impression that this filter needed an estimate of the PSF? Where would this be input?

I haven’t had any luck using this filter. The output is always a greyed out volume with a black dot in the middle. I have tried playing with all the settings and different iterations. I have also tried converting the volume to float as well as double.

Thanks,

RichardsonLucy

We worked on orbital bone reconstruction from CT images where segmentation of thin bones was very challenging. Deconvolution helped but simple unsharp masking worked about just as well in practice.

See these posts for some more details:

I would recommend not to worry about achieve optimal results by preprocessing. If you oversample the volume (so that you can represent thinner walls) and apply some edge enhancement filter then it should help. You can take care of the remaining issues in further segmentation steps (e.g., you can create a continuous bone surface from segmented contours by using WrapSolidify segment editor effect).