In case you want to convolve the image later with a given point spread function, make sure to use a larger image size. After convolution, the profiles become larger and a profile that is normally completely outside of the image might fall within it.
On one axis, if you want your final (convolved) image to be \(m\) pixels and your PSF is \(2n+1\) pixels wide, then when calling MakeProfiles, set the axis size to \(m+2n\), not \(m\). You also have to shift all the pixel positions of the profile centers on the that axis by \(n\) pixels to the positive.
After convolution, you can crop the outer \(n\) pixels with the section crop box specification of Crop: --section=n+1:*-n,n+1:*-n (according to the FITS standard, counting is from 1 so we use n+1
) assuming your PSF is a square, see Crop section syntax.
This will also remove all discrete Fourier transform artifacts (blurred sides) from the final image.
To facilitate this shift, MakeProfiles has the options --xshift, --yshift and --prepforconv, see Invoking MakeProfiles.
JavaScript license information
GNU Astronomy Utilities 0.23 manual, July 2024.