如何使用 R 中的雙線性插值降低 netCDF 的分辨率（Regrid）？

Question

我從這里下載了 netCDF 文件。 它們的分辨率為 0.5*0.5。 我想以更粗的 1*1 分辨率重新網格這些文件。 我找到了一些鏈接。 第一個鏈接討論 R 中的重新網格，但不使用雙線性插值。 第二個鏈接處理雙線性插值，但使用氣候數據運算符（我不太熟悉）。 然后我遇到了一個 R 包HiClimR 。 在這個包中，一個命令coarseR降低了數據的分辨率。 我將 netCDF 文件轉換為 excel 文件並使用coarseR 但是得到結果后我發現這個命令實際上以某種方式跳過了經緯度並將分辨率降低到1 * 1。 簡而言之，我的問題是

(1) 使用coarseR來降低分辨率是否正確？ (2) 如何將雙線性變換用於我在 R 中的特定問題？

提前謝謝了。

Answer 1

這個解決方案不在 R 中，只是為了指出 CDO 內置了一些非常好的重新網格化功能。例如，按照您的要求使用雙線性插值重新網格化為 1x1 規則網格，您可以簡單地執行以下操作：

cdo remapbil,r720x360 in.nc out.nc

但是，如果您要轉換為較粗的網格，則建議使用保守的重新映射技術，否則您可能會在重新映射期間遺漏點。 這對於高度不均勻的領域（例如降水）尤其重要。 在這種情況下，CDO 提供一階和二階保守重映射技術。 使用一階技術

cdo remapcon,r720x360 in.nc out.nc

（請注意，有時您可能會發現 CDO 由於轉換過程中的精度損失而導致抖動，在這種情況下，它會建議您使用選項“-b32”或“-b64”。除非被迫，否則您不想這樣做以增加文件大小。）

要考慮的常見重新網格選項有：

• remapbil：雙線性插值
• remapnn：最近鄰插值（即從最近的單元格取值）
• remapcon：一階保守重映射
• remapcon2：二階保守重映射