无法复制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

• 第一个选项，获得GDAL 3.4，在这里解决了这个问题，GRIB在从GRIB转换为GeoTIFF时会自动从0-360转换为-180-180
• 第二个选项，使用NPM提供的geosub(npm -g install geosub(下载NOAA的GRIB。如果这是你正在使用的，它可以为你做到这一点
• 第三种选择，使用从早期就存在的gdalwarp --config CENTER_LONG 0

(免责声明：我是GDAL中GRIB 0-360翻译和geosub包的作者(