我正在使用aplypy从三个FITS文件创建Eagle星云的RGB图像,表示图像的红色,绿色和蓝色分量。 FITS文件可用here,红色为673nm,绿色为656nm,蓝色为502nm。
import aplpy
aplpy.make_rgb_cube(['673nmos.fits','656nmos.fits','502nmos.fits'], 'nmod_cube.fits')
aplpy.make_rgb_image('nmod_cube.fits','nmod.png')
f = aplpy.FITSFigure('../data/nmod_cube_2d.fits')
f.show_rgb('../data/output/nmod.png')
图像应该看起来像this,但看起来像下面的图像。蓝色组件显然与红色和绿色不对齐。
根据make_rgb_cube
,docs应根据每个FITS标题中的WCS信息将三个图像重新对齐到相同的投影中,但在这种情况下它似乎不起作用。
我是否需要使用星形匹配库来获得对齐精度?
Blue FITS文件的标题中的WCS坐标出错。我已经联系了spacetelescope.org以提醒他们。
来自M17的the site的另一个例子是使用问题中显示的aplpy脚本正确创建的。按照注释中的建议,使用DS9闪烁每个过滤器,确认正确对齐。
使用STScI软件TweakReg对齐图像并在组合之前获得良好的对齐。更多信息和工作示例可以在DrizzlePac website上找到。
TweakReg:
使用星形图的组合图像要求每个输入图像的标题中的WCS信息在子像素精度内对齐。 tweakreg任务允许用户将图像集彼此对齐和/或与外部天体测量参考帧或图像对齐。
它针对HST仪器进行了优化,但可以适用于其他天文台/仪器。