在这种情况下,我使用astropy可视化制作M66的彩色图像。 在做任何事情之前,我了解到我必须使用numpy.float(()
forCasting = np.float_()
### READING
b = fits.open("data/"+"M66-Blue.fts")[0].data
r = fits.open("data/"+"M66-Red.fts")[0].data
g = fits.open("data/"+"M66-Green.fts")[0].data
### CASTING
r = np.array(r,forCasting)
g = np.array(g,forCasting)
b = np.array(b,forCasting)
这样我就可以继续做伸展运动,比如:
stretch = SqrtStretch() + ZScaleInterval()
r = stretch(b)
g = stretch(r)
b = stretch(g)
plt.imshow(r, origin='lower')
plt.show()
plt.imshow(g, origin='lower')
plt.show()
plt.imshow(b, origin='lower')
plt.show()
然后我使用astropy.visualizaion中的方法制作rgb,如下所示,但我有一张超暗图像,我无法分辨任何东西。有人知道我为什么在这里有一个黑暗的最终图像吗?你有什么建议吗
### SAVING
# rgb_default = make_lupton_rgb(r, g, b, minimum=1000, stretch=900, Q=100, filename="provafinale.png")
rgb_default = make_lupton_rgb(r, g, b, filename="provafinale.png")
plt.imshow(rgb_default, origin='lower')
plt.show()
谢谢
看起来您必须设置
make_lupton_rgb
的stretch
和Q
参数。默认值是
stretch=5
和Q=8
,这会产生暗结果我没有
astropy
或天文学方面的经验。我只是玩了一下参数,然后用
stretch=1
和Q=0
得到了明亮的图像我尝试使用
np.percentile
计算minimum
和stretch
,以线性拉伸输出我使用来自index_fits的
m8_050507_9i9m
图像测试了代码以下是我用于测试的代码:
结果:
相关问题 更多 >
编程相关推荐