astropy

    0熱度

    1回答

    我正在爲一個項目做一些目錄匹配,發現使用astropy和其他座標匹配軟件包(例如Topcat)相比,我得到的結果是不同的。特別是,我不明白爲什麼下面的代碼會產生不對稱的結果 import numpy as np from astropy.coordinates import match_coordinates_sky, SkyCoord # Get some random coordinat

    1熱度

    1回答

    我已經嘗試調試我的代碼,並且我已經意識到,當我嘗試將我的AltAz座標保存到.csv文件時,它最終會崩潰,因爲它不是一個numpy數組,它是一個SkyCoord對象。有人可能會提出一個簡單的方法來將大型赤道座標錶轉換爲AltAz,或者我可以如何將我的代碼保存到文件中。 # Get time now time = astropy.time.Time.now() time.delta_ut1_ut

    -1熱度

    1回答

    我想我搞砸了的東西從擬合文件選擇數據時...基本上我所做的是: 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

    1熱度

    1回答

    我想在Python中讀取一個ascii表格中的Numpy/Pandas/Astropy數組/數據框/表格。表中的每一行看起來是這樣的: 329444.6949 0.0124 -6.0124 3 97.9459 15 32507 303 7 3 4 8 2 7 HDC-13-O 的問題是,有列之間沒有明顯的分隔符/分隔符,所以對於一些行有兩列之間沒有空格,就像這樣: 332174.928

    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

    0熱度

    1回答

    將np.ndarray(將numpy導入aspopy.coordinates.Angle類)最快/最有效的方法是什麼?我無法將其保存爲np.ndarray,因爲.wrap_at()操作不起作用。

    0熱度

    1回答

    我已經有一個寫好的文件(有人發給我),我想添加校驗和和datasum頭。我使用astropy.io.fits添加校驗和的唯一示例涉及在將每個部分添加到HDU時構建新的擬合HDU and verifying it。我可以做的,但似乎它會有更多的開銷,然後是需要的。 有沒有辦法將校驗和和datasum添加到現有的HDU?

    1熱度

    1回答

    我正在使用astropy.table.write(filename,path = run_dir)將一個astropy表寫入名爲dat.h5的文件。但是我收到文件存在的錯誤,我在下面用pdb跟蹤顯示它沒有。到底是怎麼回事? (Pdb) run_dir '/Users/ms/run0' (Pdb) os.system("ls " + run_dir) param.txt temp_in.da

    6熱度

    1回答

    我有兩個來源與赤道座標(ra, dec)和(ra_0, dec_0)位於距離r和r_0,我需要計算它們之間的3D距離。 我使用兩種方法,應該給我相同的結果,據我瞭解,但不。 第一種方法是應用astropy的separation_3d函數。第二種方法是使用給出了兩個源之間的距離與球面座標的表達式: 如圖here。 在下面的MCVE,返回的值是: 91.3427173002 pc 93.847049

    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