我在python中有5个天文图像,每个图像都有不同的波长,因此它们具有不同的角度分辨率和网格尺寸,为了比较它们,我可以创建温度图,我需要它们具有相同的角度分辨率和网格大小。
我已经设法将Gaussian卷积成每个图像的角度分辨率与最差的图像相同,但是我无法找到一种方法来重新格式化python中的每个图像,并想知道是否有人知道如何去做这个?
我希望将图像重新网格化为与最差质量图像相同的网格尺寸,因此如果需要,我可以将其用作参考图像。谢谢
答案 0 :(得分:2)
如果图像标题具有正确的世界坐标系数据,则可以使用reproject包重新采样图像: http://reproject.readthedocs.org/en/stable/
答案 1 :(得分:1)
您可以使用FITS_tools(https://pypi.python.org/pypi/FITS_tools,它也可以通过$ pip install FITS_TOOLS
安装在python的anaconda发行版中)。
两个图像必须在标题信息中具有wcs。
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)
它返回:两个图像投射到同一空间,并且可选择用于投影它们的标题。安装后,您可以提供帮助(FITS_tools.match_fits)以获取更多信息。