python - 无法在 python 中读取地理空间数据（img、rrd、ige 文件）

Question

我有一个矢量图保存在一个文件夹中，其中包括三个文件（.img、.rrd、.ige）。我想将向量值（整数）保存在 NumPy 数组中。我已经通过以下代码将这张地图读入python，它没有报告错误。

from osgeo import gdal
from osgeo.gdalconst import *
import numpy as np

driver = gdal.GetDriverByName('HFA')
driver.Register()
img = gdal.Open('./5m DEM/tw_dem.img', GA_ReadOnly)
geoinfo = img.GetGeoTransform()
ltX = geoinfo[0]
pixX = geoinfo[1]
ltY = geoinfo[3]
pixY = geoinfo[5]

但在我添加以下代码后，Spyder 报告错误（显示在代码下方）。有人能告诉我如何阅读这张地图中的矩阵吗？非常感激。

x,y = 333960.43196014385,2776767.60157236

xOffset = int((x-ltX) / pixX)
yOffset = int((y-ltY) / pixY)
band = img.GetRasterBand(1)
data = band.ReadAsArray(xOffset, yOffset, 1,1 )

错误 4：无法打开外部数据文件：./5m DEM\tw_dem.ige 错误 1：./5m DEM/tw_dem.img，波段 1：IReadBlock 在 X 偏移量 0、Y 偏移量 0 处失败：无法打开外部数据文件: ./5m DEM\tw_dem.ige 错误 4: 无法打开外部数据文件: ./5m DEM\tw_dem.ige 错误 1: ./5m DEM/tw_dem.img, band 1: IReadBlock 在 X 偏移量 0, Y 处失败偏移量 0：无法打开外部数据文件：./5m DEM\tw_dem.ige 错误 5：./5m DEM/tw_dem.img: GDALDataset::GetRasterBand(0) - Illegal band # 错误 4：无法打开外部数据文件：./5m DEM\tw_dem.ige 错误 1：./5m DEM/tw_dem.img，波段 1：IReadBlock 在 X 偏移 569，Y 偏移 74 处失败：无法打开外部数据文件：./5m DEM\tw_dem.ige错误 4：无法打开外部数据文件：./5m DEM\tw_dem.ige 错误 1：./5m DEM/tw_dem.img，波段 1：IReadBlock 在 X 偏移 569，Y 偏移 74 处失败：无法打开外部数据文件：./5m DEM\tw_dem.ige

标签： pythonnumpyreadfile