Python调用gdal_calc与命令行结果不一致，顶部栅格空白

我有一个Z值均为正数的插值文件，希望使用gdal_calc将其转为负值。在GDAL命令行执行以下命令可正常得到结果：

gdal_calc -a interpolated.tif --outfile gdal_negative.tif --calc="-abs(-a)"

我尝试通过Python的subprocess模块集成该逻辑，代码如下（插值文件已存入interpolated_file变量）：

print('Pre calc')
negative_file = output_file + '_negative.tif'
calculation = '--calc=-abs(-a)'
subprocess.call(['gdal_calc.py', '-a', interpolated_file, '--outfile', negative_file, '--overwrite', calculation])
print('Post calc')

脚本运行无报错，但生成的栅格仅约3/4区域符合预期，顶部区域为空白，而命令行执行则无此问题，请问原因是什么？

核心原因是subprocess参数拆分错误：

• 命令行中--calc="-abs(-a)"的引号是给系统shell用的，作用是避免shell把计算式拆分成多个独立参数；
• 但用subprocess的列表形式传参时，不需要额外加引号，且必须把--calc和它的取值-abs(-a)拆成两个独立的列表元素。

你的代码把--calc=-abs(-a)作为单个参数传递，gdal_calc.py解析参数时出现异常，导致计算范围处理错误，最终生成的栅格顶部出现空白。

修正后的代码：

print('Pre calc')
negative_file = output_file + '_negative.tif'
# 将--calc与计算式拆分为两个独立参数
subprocess.call(['gdal_calc.py', '-a', interpolated_file, '--outfile', negative_file, '--overwrite', '--calc', '-abs(-a)'])
print('Post calc')

如果一定要用字符串形式传参（不推荐，存在命令注入风险），可以开启shell=True，此时要保留命令行的引号格式：

print('Pre calc')
negative_file = output_file + '_negative.tif'
cmd = f'gdal_calc.py -a {interpolated_file} --outfile {negative_file} --overwrite --calc="-abs(-a)"'
subprocess.call(cmd, shell=True)
print('Post calc')

内容的提问来源于stack exchange，提问作者user09873920134