2015-10-13 45 views
1

我有5个python天文图像,每个图像的波长不同,因此它们具有不同的角度分辨率和网格大小,并且为了比较它们以便我可以创建我需要的温度图他们是相同的角度分辨率和网格大小。在python中重新调整拟合图像

我已经成功地将每幅图像卷成与最差图像相同的角度分辨率,但是我很难找到一种方法来重新格网化每个图像,并想知道是否有人知道如何去做这件事?

我希望将图像重新网格化为与最差质量图像相同的网格大小,因此如果需要,我可以将它用作参考图像。谢谢

+0

你使用什么软件包?你用什么来形象化? matplotlib? – wgwz

+0

我已经安装了所有基本软件包,numpy,scipy等以及其他如astropy和APLpy。是的,我正在使用matplotlib – taylor993

+0

嗯,我已经尝试从scipy imresize功能和从numpy调整大小,但这些不适合我适合的文件,我还没有能够找到另一种方法到目前为止。 – taylor993

回答

2

如果图像头有正确的世界坐标系统的数据,你可以使用重新投影包重新采样图像: http://reproject.readthedocs.org/en/stable/

+0

'reproject'是大多数情况下使用的正确工具。如果你正在处理多维数据集,至少现在FITS_tools(http://fits-tools.readthedocs.org/en/latest/)有一些更方便的包装。 – keflavich

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)以获取更多信息。