我定義了以下函式來插入網格資料:
def interpolate_quadrant(file, r_out):
f = open(file, 'r')
x = []
y = []
z = []
for line in f:
line = line.strip()
columns = line.split()
x.append(float(columns[0]))
y.append(float(columns[1]))
z.append(float(columns[2]))
f.close()
xold = np.asarray(x)
yold = np.asarray(y)
zold = np.asarray(z)
xx = np.linspace(np.min(xold), np.max(xold), num= r_out)
yy = np.linspace(np.min(yold), np.max(yold), num= r_out)
xx, yy = np.meshgrid(xx, yy)
valsq1 = interpolate.griddata((xold, yold), zold, (xx, yy), method = 'cubic')
valsq1 = valsq1[1:,:]
valsq2 = valsq1[::-1,:]
valsq3 = valsq1[:,::-1]
valsq4 = valsq1[::-1,::-1]
hdu = fits.PrimaryHDU(data = r)
hdu.writeto('FITS_file_(1).fits')
當我在 Python 單元格中 Shift Enter 時,它會作業并通過。但是,當我通過以下方式在新單元格中呼叫該函式時:
interpolate_quadrant(/home/ben/Documents/Planetary_Astrophysics/PE_MHD_Project/faceon1cm_mhd_5_30_ml-8.6_VE.out, 200)
我收到附加影像中顯示的語法錯誤。Jupyter 實驗室中的語法錯誤。
我不確定我用作函式引數的檔案鏈接有什么問題。要輸入的檔案型別是“.out”。
uj5u.com熱心網友回復:
似乎解決方案是用引號對路徑進行編碼。google collab 就是這種情況,不過我從未使用過 jupyter notebook。
interpolate_quadrant('/home/ben/Documents/Planetary_Astrophysics/PE_MHD_Project/faceon1cm_mhd_5_30_ml-8.6_VE.out', 200)
轉載請註明出處,本文鏈接:https://www.uj5u.com/caozuo/348415.html
上一篇:在使用函式時如何結束我的代碼?
下一篇:如何將此jquery方法-.text()、.height()和.css()轉換為vanillajavascript?
