<reponame>kapteyn-astro/kapteyn
from kapteyn import maputils
from matplotlib import pyplot as plt
import numpy

# Read first image as base 
Basefits = maputils.FITSimage(promptfie=maputils.prompt_fitsfile)
print(type(Basefits), isinstance(Basefits, maputils.FITSimage))

# Get data from a second image. This is the data that 
# should be reprojected to fit the header of Basefits.
Secondfits = maputils.FITSimage(promptfie=maputils.prompt_fitsfile)
#Secondfits.set_imageaxes(promptfie=maputils.prompt_imageaxes)
#Secondfits.set_limits(promptfie=maputils.prompt_box)

# Now we want to overlay the data of this Base fits object onto
# the wcs of the second fits object. This is done with the 
# reproject_to() method of
# the first FITSimage object (the data object) with the second
# FITSimage object as parameter. This results in a new fits file

#Reprojfits = Basefits.reproject_to(Secondfits.hdr, plimlo=(2,1), plimhi=(2,1))
#Reprojfits = Basefits.reproject_to(Secondfits.hdr, pxlim=(100,400),  pylim=(100,400))
Reprojfits = Basefits.reproject_to(Secondfits.hdr)
Reprojfits.writetofits("reproj.fits", clobber=True)

from kapteyn import maputils
from matplotlib import pyplot as plt
import numpy

# Read first image as base 
Basefits = maputils.FITSimage(promptfie=maputils.prompt_fitsfile)
print(type(Basefits), isinstance(Basefits, maputils.FITSimage))

# Get data from a second image. This is the data that 
# should be reprojected to fit the header of Basefits.
Secondfits = maputils.FITSimage(promptfie=maputils.prompt_fitsfile)
#Secondfits.set_imageaxes(promptfie=maputils.prompt_imageaxes)
#Secondfits.set_limits(promptfie=maputils.prompt_box)

# Now we want to overlay the data of this Base fits object onto
# the wcs of the second fits object. This is done with the 
# reproject_to() method of
# the first FITSimage object (the data object) with the second
# FITSimage object as parameter. This results in a new fits file

#Reprojfits = Basefits.reproject_to(Secondfits.hdr, plimlo=(2,1), plimhi=(2,1))
#Reprojfits = Basefits.reproject_to(Secondfits.hdr, pxlim=(100,400),  pylim=(100,400))
Reprojfits = Basefits.reproject_to(Secondfits.hdr)
Reprojfits.writetofits("reproj.fits", clobber=True)
# Read first image as base # Get data from a second image. This is the data that # should be reprojected to fit the header of Basefits. #Secondfits.set_imageaxes(promptfie=maputils.prompt_imageaxes) #Secondfits.set_limits(promptfie=maputils.prompt_box) # Now we want to overlay the data of this Base fits object onto # the wcs of the second fits object. This is done with the # reproject_to() method of # the first FITSimage object (the data object) with the second # FITSimage object as parameter. This results in a new fits file #Reprojfits = Basefits.reproject_to(Secondfits.hdr, plimlo=(2,1), plimhi=(2,1)) #Reprojfits = Basefits.reproject_to(Secondfits.hdr, pxlim=(100,400),  pylim=(100,400))