使用numpy和scipy將wrf數據的網格點插值到python中的lat lon

Question

我是python的新手，我正在嘗試編寫一個函數，它將從一個帶有尺寸[時間，高度，經度，緯度]的netcdf文件中獲取一個numpy數組，並將該函數內插到指定的lat和lon。 我已經研究過scipy.interpolate但我仍然不確定從那里去哪里。 任何幫助？

Answer 1

如果您希望（a）將整個數據集內插到一組新的緯度/經度坐標或（b）在單個緯度/經度位置獲取該值，我無法從問題中看出。

無論哪種方式，如果你沒有在python中這樣做，那么在Climate Data Operators中使用remapgrid會更加容易（確保你使用netCDF支持安裝）。 這確實需要netCDF文件遵循CF約定（請參閱此文章，以使您的WRF文件符合CF標准）。

有關remapgrid的文檔可在此處找到。 有關網格選項的信息，請參閱第1.3節（第9-12頁）;有關插值技術的信息，請參見第2.10節（第106頁開始）。 如何重映射有多種選擇，您使用的命令取決於您的應用程序。

（a）這是使用雙線性插值重新映射到高斯（128x64）網格的示例

cdo remapbil,n32 wrfout_d01_1999-01-01-01_00:00:00.nc out.nc

（b）這是使用最近鄰居映射到一個緯度經度點的示例重映射：

cdo remapnn,lon=-107.0_lat=34.0 wrfout_d01_1999-01-01-01_00:00:00.nc out.nc

Answer 2

您可以使用XESMF重新編譯它。

它有GEOSChem的教程 ，但對於WRF來說幾乎是一樣的。

我寫了一個關於如何重新調整WRF數據的筆記本 。