Python:GDAL 错误 - 访问窗口超出 RasterIO() 范围

Question

我的 GDAL/遇到此错误问题

INPUT_raster = "E:\myraster.asc"
ds = gdal.Open(INPUT_raster, gdal.GA_ReadOnly)
band = ds.GetRasterBand(1)
data = band.ReadAsArray(9658, 11599, 1, 1)
ERROR 5: E:\myraster.asc, band 1: Access window out of range in RasterIO().  Requested
(9658,11599) of size 1x1 on raster of 9658x16934.

是否可以使用try:和expect:来避免这个错误？

非常感谢您的意见/建议。

问候，詹尼

Answer 1

尝试/期望不要避免错误，它们会处理它们。

您正尝试从栅格中不存在的位置读取数据。 x 维度为 9658 个元素，如果您想要最外面的元素，您应该使用 9657，因为索引从零开始。

所以使用:

data = band.ReadAsArray(9657, 11599, 1, 1)

或最后一个 x 位置:

data = band.ReadAsArray(ds.RasterXSize-1, 11599, 1, 1)