Deconvolution algorithms have been used successfully for optimization/restoration/deblurring of astronomical images in a variety of wavelengths, especially in the optical band (e.g., for HST). We present here an iterative Richardson-Lucy type method
designed for treatment of images obtained with the H.E.S.S. array of ground-based gamma-ray telescopes. Its application to shell-type supernova remnant images yields refined details relevant for the study of correlations with other wavelengths, and hence for interpretation in terms either of hadronic or leptonic origin of the observed VHE gamma-ray emission.
