使用Victor Velasquez编写的栅格降水数据提取代码时出现IndexError索引越界问题求助

嘿，我一眼就看出你遇到的问题大概率是坐标索引顺序搞反了，导致访问了完全超出栅格维度的索引值。咱们一步步拆解和解决：

错误根源分析 #

从报错信息IndexError: index 158116290 is out of bounds for axis 0 with size 198就能看出来：你用一个超大的数值（158116290）去访问只有198行的栅格数据（axis0对应栅格的行维度），这肯定越界啊。

问题就出在这两行代码：

X,Y = Data.index(self.x,self.y) #extraigo
Pre = Band[X,Y]

Rasterio的DatasetReader.index()方法返回的是**(行号, 列号)，对应栅格数组的Band[行, 列]。但这里有个很容易踩的坑：如果你的self.x是经度、self.y是纬度，那你把坐标顺序传反了——index()方法要求传入的是(x, y)，也就是(横坐标, 纵坐标)**，而很多地理数据场景里我们习惯先写纬度再写经度，或者你的坐标顺序和栅格的空间参考不匹配，直接导致返回的行号变成了一个离谱的大数。

具体解决方案 #

1. 修正坐标传入顺序 #

先搞清楚self.x和self.y的含义：

• 如果self.x是经度，self.y是纬度，那你需要把参数顺序反过来传入index()，同时把变量名改得更清晰（避免混淆）：

  # 交换x和y的传入顺序，同时用row/col代替X/Y，更直观
  row, col = Data.index(self.y, self.x)
  Pre = Band[row, col]
  

• 你也可以用Data.xy(row, col)方法反推验证：比如取栅格中心的行列，看看返回的坐标是否和你的点逻辑匹配，确认顺序是否正确。

2. 提前检查点是否在栅格范围内 #

就算顺序对了，如果你的点不在当前栅格的空间边界内，index()也会返回超出范围的索引值。建议在提取前加个范围校验：

with rasterio.open(self.Lis[i]) as Data:
    # 检查当前点是否在栅格的空间范围内
    bounds = Data.bounds
    if not (bounds.left <= self.x <= bounds.right and bounds.bottom <= self.y <= bounds.top):
        print(f"警告：点({self.x}, {self.y})不在文件{self.Lis[i]}的范围内")
        self.ListPre.append(None) # 用None标记缺失值，避免后续报错
        continue
    # 正确提取行和列（确认坐标顺序正确后使用）
    row, col = Data.index(self.x, self.y)
    Band = Data.read(1)
    Pre = Band[row, col]
    self.ListPre.append(Pre)

3. 验证栅格维度是否符合预期 #

你可以在循环里加一行打印，确认每个栅格的行/列数量，避免因为栅格本身的问题导致索引异常：

with rasterio.open(self.Lis[i]) as Data:
    print(f"文件{self.Lis[i]}的栅格形状（行, 列）：{Data.shape}")
    Band = Data.read(1)
    print(f"波段数组形状：{Band.shape}")

这样你就能直观看到，你要访问的row和col是否在合理的数值范围内。

额外小建议 #

• 把变量名改得更语义化，比如用row/col代替X/Y，能从根源上减少坐标和索引的混淆。
• 因为你处理的是逐日栅格，要确保所有.tif文件的空间范围、分辨率和投影是完全一致的，不然不同文件的行号列号可能对应不同的地理点，提取的数据会出错。

内容的提问来源于stack exchange，提问作者Gonzalo Lagas