Convolution of Sunpy maps

[afbattaglia]

Describe the feature

Hi everyone!

As discussed in the Element Sunpy group on 17-Oct-2023, I was wondering if there is a way to convolve Sunpy maps with externally provided functions. Specifically, I have an HMI map and for an optical design I'm working on, I need to convolve it with an artificial PSF of 10 arcsec FWHM. It seems that we cannot do this directly in Sunpy, but we can follow these steps: extract the array from the Sunpy map, apply scipy.ndimage.gaussian_filter, and then reconstruct the map.

This seems to be a common task, so it may be useful to have a direct method for doing it in Sunpy.

Proposed solution

A possible solution would be to create a method including the following steps: extract the array from the Sunpy map, apply scipy.ndimage.gaussian_filter, and then reconstruct the map.