Simple Image Reprojection

reproject.py

from reproject import reproject_interp
from astropy.io import fits
from glob import glob
import numpy as np

# get the header for the file we want to reproject to
hdr = fits.getheader("NEP5281_r_0_new.fits")
# get the filenames as a list for each filter
rfiles = glob("*r*_new.fits")
ifiles = glob("*i*_new.fits")
zfiles = glob("*z*_new.fits")

# here is an example projection and reduction for just the r frames. You'll need to do the rest
reprojected_r = [reproject_interp(fname, hdr, return_footprint=False) for fname in rfiles]
combined_r = np.median(reprojected_r, axis=0)
fits.writeto("NEP5281_r_master.fits", combined_r, header=hdr)