我一直在尝试在ipython中使用光栅计算来获取我上传的tif文件，但是我无法找到该函数的整个代码。我一直在寻找下面的例子，但我不确定如何使用它。
gdal_calc.py -A input.tif --outfile=result.tif --calc="A*(A>0)" --NoDataValue=0
然后我通过分配部分尝试了另一个过程，但是这仍然不起作用（下面的代码）
a = '/iPythonData/cstone/prec_7.tif'
outfile = '/iPythonData/cstone/prec_result.tif'
expr = 'A<125'
gdal_calc.py -A=a --outfile=outfile --calc='expr' --NoDataValue=0
不断提出无法分配给运营商。有人可以帮忙整个代码。
答案 0 :(得分：0)
查看gdal_calc.py的源代码，该文件只有大约300行。这是该文件的链接。
https://raw.githubusercontent.com/OSGeo/gdal/trunk/gdal/swig/python/scripts/gdal_calc.py
他们只是在main中创建一个OptionParser对象并将其传递给doit()方法（第63行）。您可以根据通过命令行传递给它的相同参数生成相同的OptionParser实例，并直接调用它们的doit方法。
也就是说，系统调用对@thomas-k完全有效。只有当你真的想要留在Python环境中时才会这样。