使用gdalinfo从栅格中提取栅格属性表

时间:2017-02-04 05:26:09

标签: python gdal

_rat = subprocess.check_output('gdalinfo -json ' + dataset_uri, shell=True)

我想提取.tif文件的栅格属性表。在上面的命令中,我能够将信息转换为_rat但不确定如何从_rat中提取rat部分。有什么建议?

1 个答案:

答案 0 :(得分:1)

在您的代码中,_rat是一个有效JSON的字符串。您可以将该JSON转换为python dict,以便您轻松访问元素。我记得gdalinfo没有提供完整的栅格属性表,但仍然存在有意义的带级统计数据。

import json
import subprocess

dataset_uri = 'input.tif'
_rat = subprocess.check_output('gdalinfo -json ' + dataset_uri, shell=True)
data = json.loads(_rat) # load json string into dictionary
print data

# to get band-level data
bands = data['bands']