CODEWITHSUNDEEP
pythonimage-resizingfitsastropy
taylor993
taylor993

Reputation: 49

re-sizing a fits image in python

I have 5 astronomy images in python, each for a different wavelength, therefore they are of different angular resolutions and grid sizes and in order to compare them so that i can create temperature maps i need them to be the same angular resolution and grid size.

I have managed to Gaussian convolve each image to the same angular resolution as the worst one, however i am having trouble finding a method to re-grid each image in python and wondered if anyone knew how to go about doing this?

I wish to re-grid the images to the same grid size as the worst quality image and so i can use that as a reference image if required. Thank you

Upvotes: 1

Views: 5778

Answers (2)

user3835290
user3835290

Reputation: 196

You can use FITS_tools (https://pypi.python.org/pypi/FITS_tools, it can also be installed via $ pip install FITS_TOOLS in the anaconda distribution of python). Both images must have wcs in the header information. 

import FITS_tools
to_be_projected = 'fits_file_to_be_projected.fits'
reference_fits  = 'fits_file_serving_as_reference.fits'
im1,im2 = FITS_tools.match_fits(to_be_projected,reference_fits)

It returns: Two images projected into the same space, and optionally the header used to project them. Once installed, you can do help(FITS_tools.match_fits) for more information.

Upvotes: 1

user3148185
user3148185

Reputation: 524

If the image headers have the correct World Coordinate System data, you can use the reproject package to resample the images: http://reproject.readthedocs.org/en/stable/

Upvotes: 3

Related Questions