將北極立體二進制重投影到 Lon/Lat 網格（python 中的光柵或 gdal）

Question

現在，我正在研究第四階段數據集（CONUS 中的降水數據集）。
https://data.eol.ucar.edu/dataset/21.093

這是形成 (1121 x 881) 陣列和局部 4 公里極地立體網格的二進制數據。
header 文件包括： “極地立體聲：Lat1 23.117000 Long1 -119.023000 Orient -105.000000”。

我想將此數據轉換為 Lon/Lat 網格，以便在空間上與其他數據集匹配。
所以，我需要先將此二進制文件轉換為 geoTif。
通常，我使用 rasterio 將二進制（數組）轉換為 tif，如下所示；

def convert_np2tif(data, out_file,
                   width_of_pixel, height_of_pixel, lon_upper_left, lat_upper_left,
                   crs):
    affine = Affine(width_of_pixel, 0, lon_upper_left, 0, -1*height_of_pixel, lat_upper_left)
    with rasterio.open(out_file, 'w',
                        driver='GTiff',
                        height=data.shape[0],
                        width=data.shape[1],
                        count=1,
                        dtype=data.dtype,
                        crs=crs,
                        transform=affine
                        ) as dst:
        dst.write(data.reshape(1,data.shape[0],data.shape[1]))

但是，這種方式需要指定左上 lon/lat。 我不確定極地立體網格中的這個值。

有誰知道將此二進制文件轉換為 GeoTif 並轉換為 Lon/Lat 網格的好方法。
不僅 rasterio，而且 gdal 的方式都是受歡迎的。 預先感謝您的建議。

Answer 1

這個問題包含您需要的一切： https://gis.stackexchange.com/questions/278627/using-gdal-warp-and-gdal-warpoptions-of-gdal-python-api/341693

在WarpOptions dstSRS你不需要其他任何東西