python - 如何在 Python 中使用 GDAL 打开 geotiff 图像？

Question

我正在尝试运行以下代码:

from osgeo import gdal
import sys

# This allows GDAL to throw Python exceptions
src_ds = gdal.Open("fused.tif")
src_ds.show()

但是我收到以下错误:

Traceback (most recent call last):
    File ".../gdalopen1.py", line 5, in module src_ds.show()
AttributeError: 'Dataset' object has no attribute 'show'

为什么会这样？

Answer 1

以下代码打开一个光栅文件并将光栅的波段读入 NumPy数组。

from osgeo import gdal
ds = gdal.Open('input.tif', gdal.GA_ReadOnly)
rb = ds.GetRasterBand(1)
img_array = rb.ReadAsArray()