Use astropy.nddata for convolution

I need to re-implement smoothing with an asymmetric kernel, and we should also be using FFT convolution for large arrays/kernels. However, I had issues with the FFT convolution when testing, so I will investigate some more.