我正在爲一個項目做一些目錄匹配,發現使用astropy和其他座標匹配軟件包(例如Topcat)相比,我得到的結果是不同的。特別是,我不明白爲什麼下面的代碼會產生不對稱的結果 import numpy as np
from astropy.coordinates import match_coordinates_sky, SkyCoord
# Get some random coordinat
我已經嘗試調試我的代碼,並且我已經意識到,當我嘗試將我的AltAz座標保存到.csv文件時,它最終會崩潰,因爲它不是一個numpy數組,它是一個SkyCoord對象。有人可能會提出一個簡單的方法來將大型赤道座標錶轉換爲AltAz,或者我可以如何將我的代碼保存到文件中。 # Get time now
time = astropy.time.Time.now()
time.delta_ut1_ut
我想我搞砸了的東西從擬合文件選擇數據時...基本上我所做的是: File='/SomePath/xxx.fits'
hdulist=fits.open(File)
tbdata=hdulist[1].data
和應用選擇規則: for i in range (0, len(tbdata)):
if tbdata[i]['z']<0.2:
A.append(tbdata
我嘗試讀取世界座標系(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
我已經有一個寫好的文件(有人發給我),我想添加校驗和和datasum頭。我使用astropy.io.fits添加校驗和的唯一示例涉及在將每個部分添加到HDU時構建新的擬合HDU and verifying it。我可以做的,但似乎它會有更多的開銷,然後是需要的。 有沒有辦法將校驗和和datasum添加到現有的HDU?