无法复制gdal输出

Question

我有一组被翻转的GRIB文件(经度跨度从0到365)，我使用gdal首先将数据转换为GeoTIFF，然后将网格数据扭曲为标准的WGS84经度(-180到180)。到目前为止，我一直在使用命令行中的gdal_translate和gdalwarp的组合，并使用parallel快速地遍历所有文件。这些是bash脚本中的函数：

gdal_multiband_transform(){
    FILEPATH=$1
    SAVEPATH=$2

    NUM_BANDS=$(gdalinfo $FILEPATH | grep 'Band' | wc -l)

    if [[ $NUM_BANDS -eq 1 ]]
    then
        echo "Extracting 1 band from $FILEPATH"
        gdal_translate -of GTiff -b 1 $FILEPATH $SAVEPATH
    else
        echo "Extracting 2 bands from $FILEPATH"
        gdal_translate -of GTiff -b 1 -b 2 $FILEPATH $SAVEPATH
    fi
}

warp_raster(){

    echo "Rewarp all rasters in $PATH_TO_GRIB"
    find $PATH_TO_GRIB -type f -name '*.tif' | parallel -j 5 -- gdalwarp -t_srs WGS84 {} {.}_warp.tif \
        -wo SOURCE_EXTRA=1000 --config CENTER_LONG 0 -overwrite
}

warp_raster

现在，我想使用osgeo库在Python中复制同样的行为。我跳过了翻译部分，因为我意识到osgeo.gdal可以直接扭曲GRIB文件，而不必转换/转换为GeoTIFF格式。为此，我使用了以下Python代码：

from osgeo import gdal

OPTS = gdal.WarpOptions(dstSRS='WGS84',
                        warpOptions=['SOURCE_EXTRA=1000'],
                        options=['CENTER_LONG 0'])

       
try:
   ds = gdal.Open(filename)
except RuntimeError:
   ds = gdal.Open(str(filename))

if os.path.getsize(filename):
   ds_transform = gdal.Warp(file_temp_path,
                             ds,
                             options=OPTS)
   # is this a hack?
   ds_transform = None
else:
   print(f'{filename} is an empty file. No GDAL transform')

其中，我使用gdal.WarpOptions从bash脚本中定义了相同的选项。结果在视觉上是一样的；代码实现了主要目标:在-180到180之间扭曲经度。但是，当我拿本地统计数据时，差别是很大的。只是整个网格数据的平均值相差4摄氏度(即表面温度数据)。我在osgeo中缺少的任何GDAL选项都会产生这些差异？我不想使用bash脚本，因为我只想要一个-Python实现。

Answer 1

• First选项，在这个问题解决的地方获得GDAL3.4，当从GRIB转换为GeoTIFF
• 第二个选项时，GRIBs会自动从0-360转换到-180-180，使用NPM (npm -g install geosub)提供的geosub下载NOAA的GRIBs，如果这是您正在使用的，它可以为您提供
• 第三选项，使用自早期就一直存在的gdalwarp --config CENTER_LONG 0

f 210。

(免责声明:我是GDAL和geosub软件包中Grib0-360翻译的作者)