Python利用ArcPy与NumPy高效提取栅格图像边界像元

键盘旁的猫

发布时间：2026-10-12 | 3人浏览过

原创

本文介绍在GIS中利用Python提取栅格图像边界的方法。通过ArcPy的RasterToNumpyArray将栅格转为数组，结合设定的nodata值（如0），遍历行列识别首尾非背景像元，从而快速生成边界掩膜。

在地理信息系统（GIS）应用中，栅格数据包含丰富的元信息，其中背景值（nodata value）对于准确识别图像边界像元至关重要。我们的核心目标是通过定位每一行和每一列中首次出现的非nodata像元，以及最后一次出现nodata像元的前一个像元，来精确提取边界。

针对栅格数据，可以利用ArcPy库中的RasterToNumpyArray函数，将栅格数据转换为NumPy数组。转换后，即可按照预设逻辑读取每一行和每一列中位于首尾位置的像元数据。

以下是提取边界像元的核心算法代码片段，其实现思路较为直观（此处假设0代表nodata value）。

a=[[0 for col in range(Raster.width)]for row in range(Raster.height)]
 
for i in range(0,Raster.width):
... for j in range(0,Raster.height):
...  if(myRaster[j][i]!=0 and myRaster[j-1][i]==0):
...    a[j][i]=myRaster[j][i]
...  if(myRaster[j][i]==0 and myRaster[j-1][i]!=0):
...    a[j-1][i]=myRaster[j-1][i]
...    
 
 
for i in range(0,myRaster.height):
... for j in range(0,myRaster.width):
...  if(arr[i][j]!=0 and arr[i][j-1]==0):
...    a[i][j]=arr[i][j]
...  if(arr[i][j]==0 and arr[i][j-1]!=0):
...    a[i][j-1]=arr[i][j-1]

以上关于Python简单栅格图像边界提取的方法分享完毕，希望能为读者提供有价值的参考。欢迎大家支持本站，获取更多技术内容。