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

Question

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

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

Answer 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']