Use astropy.nddata for convolution. We need to re-implement smoothing with an asymmetric kernel, and we should also be using FFT convolution for large arrays/kernels.
Removed unused file