我正在嘗試使用 astropy 2.0.11 和 python 2.7.15 通過對它應(yīng)用對數(shù)拉伸并更改對比度來編輯適合圖像,但我無法弄清楚。我一直在嘗試按照 astropy 網(wǎng)站上的教程打開和操作 fits 文件,但我想知道這些教程是否僅適用于最新版本的 astropy 和 python 3?抱歉我的代碼的組織。這是原型代碼,我只是想測試一些東西并讓它工作。import timeimport numpy as npimport matplotlibimport matplotlib.pyplot as pltfrom matplotlib.colors import LogNormimport astropy.visualizationfrom astropy.io import fitsfrom astropy.utils.data import download_filefrom astropy.visualization import astropy_mpl_styleplt.style.use(astropy_mpl_style)from astropy.utils.data import get_pkg_data_filenamedef main(): #My own fits file #fitsImage = get_pkg_data_filename("C:\\20180807T000456.fits") fitsImage = download_file('http://data.astropy.org/tutorials/FITS-images/HorseHead.fits', cache=True ) hdu_list = fits.open(fitsImage) hdu_list.info() #norm = ImageNormalize(stretch=LogStretch()) image_data = fits.getdata(fitsImage) print(type(image_data)) print(image_data.shape) hdu_list.close() plt.figure() plt.imshow(image_data, cmap='gray', norm=LogNorm()) plt.colorbar() # I chose the tick marks based on the histogram above cbar = plt.colorbar(ticks=[5.e3,1.e4,2.e4]) cbar.ax.set_yticklabels(['5,000','10,000','20,000']) time.sleep(10)我也無法讓圖像與 plt.imshow()任何見解都會有所幫助
1 回答

絕地?zé)o雙
TA貢獻(xiàn)1946條經(jīng)驗(yàn) 獲得超4個(gè)贊
你離得那么近!我在 Python 2.7 中運(yùn)行了你的代碼,你需要做的就是添加
plt.show()
之前time.sleep(10)
(你有什么理由包括這個(gè)?)然后你得到
另外,我認(rèn)為您不需要包含colorbar
and yticklabels
,plt.imshow
自動添加帶有 lognorm 比例的顏色條(我在獲得圖像時(shí)注釋掉了該部分)。
添加回答
舉報(bào)
0/150
提交
取消