Generalize reproject_and_coadd for N-dimensional data, and add option to specify blank pixel value and progress bar

Generalize reproject_and_coadd for N-dimensional data, and add option to specify blank pixel value and progress bar #865