栅格计算 python

Raster calculations python

我在 python 中使用 Gtif 文件进行栅格计算。我想应用公式:

干旱指数=降水量/(气温+10).

当我将它用于 Gtiff 时,出现以下错误: “类型错误:/ 不支持的操作数类型:'DatasetReader' 和'DatasetReader”。

我是 python 的新手。谢谢

这是我的代码:

tmp_tif= rasterio.open('temp.tif')
pcp_tif = rasterio.open('pcp.tif')

AI_DM = pcp_tif/tmp_tif + 10

正如@alex 指出的(见评论),您应该先阅读 DatasetReaders 的内容(见 Documentation

因此,你要做的是:

tmp_tif= rasterio.open('temp.tif').read()
pcp_tif = rasterio.open('pcp.tif').read()

AI_DM = pcp_tif/tmp_tif + 10