pyfits

    -1熱度

    1回答

    從適合文件中讀取數據時遇到問題。通常我可以從適合的文件中讀取數據,但現在我有一些數據文件,這給我一些問題。 當我試圖讀取文件,其中應包含三列和一個標題,我得到的是一個數字表格: [[64 195 45 ...,0 0 0]] 現在,這是一個整數,而數據文件確實應該包括含有雙三列和三列這樣的花車: [9819.3801,0.00000,0.00000] [9820.0025,5.50764e-16,

    -1熱度

    1回答

    我在PyCharm(社區版本4.0.2)從.fits文件中檢索數據。 z = data.field [SDSS_Z]結果爲 TypeError:'instancemethod'對象沒有屬性'getitem'。 同.fits其他外地文件我沒有問題 import pyfits import pylab as plt import numpy as np plt.close('all') hdu

    2熱度

    1回答

    我在Python中打開FITS文件時遇到問題。我收到以下錯誤消息: File "G:\Anaconda\lib\site-packages\pyfits\file.py", line 416, in _open_filelike % self.mode) IOError: File-like object does not have a 'write' method, required for

    5熱度

    1回答

    我看過並尋找解決這個問題的方法,並且什麼也沒有。 我通過matplotlib生成矩形FITS圖像,然後使用AstroPy(或PyFITS)向它們應用WCS座標。我的圖像在銀河的經度和緯度,所以適合我的地圖的標題關鍵字應該是GLON-CAR和GLAT-CAR(用於笛卡爾投影)。我看過other maps,它在SAO DS9中使用這個相同的地圖投影,並且座標工作得很好......網格完全正交,因爲它應

    0熱度

    1回答

    我想用兩個較舊的使用PyFITS創建一個新的FITS文件。 import pyfits from sys import stdout from sys import argv import time file1 = argv[1] file2 = argv[2] hdu1 = pyfits.open(file1) hdu2 = pyfits.open(file2) new0 =

    1熱度

    4回答

    我有一個包含許多列的FITS文件。該文件的一個簡化的例子中的一部分是這樣的: A B C 100 1 90 100 2 90 100 3 90 100 4 90 211 40 70 211 41 70 211 42 70 211 43 70 211 44 70 如果

    1熱度

    1回答

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

    2熱度

    1回答

    我需要製作輪廓圖並覆蓋圖像上的輪廓。我用aplpy庫來覆蓋天文圖像上的輪廓。 我已經下載在APlpy網站(https://github.com/aplpy/aplpy-examples/tree/master/data)的2MASS數據和寫了下面的代碼: import numpy as np import aplpy import atpy from pyavm import AVM im

    0熱度

    1回答

    我有一個包含二進制表的60GB FITS文件。我想一次讀取(並處理)這個表一行/入口/行/塊*。 (*我不能確定正確的命名) 我使用pyfits什麼,我想這樣做可以歸結爲簡單: import pyfits hdulist = = pyfits.open("file.fits") # the binary table has to be in the 2nd extension # henc

    1熱度

    3回答

    我正在嘗試更改FITS文件標題中關鍵字的值。 相當簡單,這是代碼: import pyfits hdulist = pyfits.open('test.fits') # open a FITS file prihdr = hdulist[1].header print prihdr['AREASCAL'] effarea = prihdr['AREASCAL']/5. print