我正在尝试完成一个过程,在这个过程中,我采用天文拟合文件,减去一个主平面文件,然后减去所得图像的对比度。
第一部分已成功完成但我的图像缺乏对比度。这是我的代码
from astropy.io.fits import getdata
import numpy
import numpy as np
import scipy
import Image
import PIL
import os
os.chdir("/localdir/")
from scipy import misc
import ImageEnhance
image = getdata('23484748.fts')
flat = getdata('Masterflat.fit')
normalized_flat = flat / numpy.mean(flat)
calibrated_image = image / normalized_flat
pix=numpy.fliplr(calibrated_image)
# the problem starts about here. How do I alter the contrast of pix?
from matplotlib import pyplot as plt
misc.imsave('saved image.gif', pix) # uses the Image module (PIL)
plt.imshow(pix, interpolation='nearest')
plt.show()
在你告诉我关于PIL功能和Matlib等的所有内容之前,我已经尝试过这些但没有成功。
我曾尝试使用image.fromarray将我的numpy数组转换为图像,但结果图像显示为纯白色。
如何拍摄我的numpy数组(pix)并改变其对比度?
出于测试目的,我已将两个示例文件放在http://members.optusnet.com.au/berrettp/
感谢你
彼得