fits

    1熱度

    2回答

    我有5個python天文圖像,每個圖像的波長不同,因此它們具有不同的角度分辨率和網格大小,並且爲了比較它們以便我可以創建我需要的溫度圖他們是相同的角度分辨率和網格大小。 我已經成功地將每幅圖像捲成與最差圖像相同的角度分辨率,但是我很難找到一種方法來重新格網化每個圖像,並想知道是否有人知道如何去做這件事? 我希望將圖像重新網格化爲與最差質量圖像相同的網格大小,因此如果需要,我可以將它用作參考圖像。謝

    1熱度

    1回答

    我嘗試讀取世界座標系(WCS)從使用satrapy這個代碼FITS文件: from astropy.wcs import WCS from astropy.io import fits data = 'file.fits' hdu = fits.open(data) w = WCS(hdu[0].header) 我得到的錯誤: WARNING: FITSFixedWarning: R

    1熱度

    1回答

    我想用PyFITS在FITS文件中編寫一個布爾結構數組。 我有一些問題。這是一個簡單的例子。 我創建測試字典並將其轉換爲結構化數組。 In [241]: test = {'p':np.array([True]*10+[False]*10,dtype='b')} In [242]: test = np.core.records.fromarrays(list(test.values()), nam

    2熱度

    2回答

    我正在製作一個程序,用於在計算機上搜索.fts和.fits文件,在該文件中打開文件並檢索與標題中特定關鍵字相對應的信息,並將該文件重命名爲該關鍵字。 我有一個問題,我不斷收到一個KeyError,因爲我在文件中找不到頭標關鍵字。有沒有解決的辦法?我希望能夠搜索各種關鍵字並執行一些操作,即使該文件中不存在該關鍵字。 這裏是代碼: from astropy.io import fits import

    4熱度

    2回答

    在previous posts我看到改變dtype的recarray可以使用astype執行。然而,我無法設法用recarray這個在其中一列中有一個數組的方法。 我recarray來自一個FITS文件記錄: > f = fits.open('myfile.fits') > tbdata = f[1].data > tbdata # FITS_rec([ (0.275916799999999

    0熱度

    1回答

    我想基於從頭部提供的緯度拼接擬合數組。但是,我似乎無法用我對python的知識和文檔中的astropy來做到這一點。我的代碼是這樣的: from astropy.io import fits import numpy as np Wise1 = fits.open('Image1.fits') im1 = Wise1[0].data im1 = np.where(im1 > *lati

    0熱度

    2回答

    我遇到了一個問題,機智適合文件操作在astropy包中,我需要一些幫助。 我基本上想拍攝一張我適合文件格式的圖像,並創建一個新文件,我需要開始輸入校正因子,然後可以使用校正因子和原始圖像生成新圖像以生成修正圖像。每個這些將具有相同的尺寸。 與此開始: from astropy.io import fits # Compute the size of the images (you can al

    1熱度

    3回答

    從FITS-Table讀取數據,並獲取包含該表的Struct,其中每個標記表示一列。 有沒有辦法將數組的結構重新格式化爲結構數組? 這樣,數組中的一個Struct代表一個Row? 通過@mgalloy做了一般性的解決方案(見下文): function SoA2AoS, table if (table eq !NULL) then return, !NULL tnames =

    1熱度

    1回答

    我有一個合適的圖像與許多天文物體。我正在嘗試從感興趣的物體上創建小型4x4「郵票」(物體周圍的部分/圖像)。我已經計算了原始適配文件中對象的像素座標,並創建了一個包含座標的文檔。我知道imshow()命令可能是最好的選擇,但我很難理解如何使用像素座標來完成任務。 from pylab import * import numpy as np import pyfits import matpl

    3熱度

    3回答

    我試圖加載到內存中的幾個2 000 FITS利用astropy.io.fits: def readfits(filename): with fits.open(filename) as ft: # the fits contain a single HDU data = ft[0].data return data data_sci = [] fo