Arcpy入门教程02：计算两景影像或DEM的差值并输出差值大于给定阈值的像素位置

目录

1 测试数据制作

2 代码编写

2.1 用到的接口

2.2 接口拼接实现功能

2.3 完整代码

3 报错及处理

---

1 测试数据制作

随便找到一幅影像数据，使用Arcmap的Toolbox工具生成随机栅格：

数据管理工具 --  栅格 -- 栅格数据集 -- 创建随机栅格

创建两次，我们即可得到两幅范围相同的栅格数据，我们将这个数据作为我们接下来编码的测试数据。

2 代码编写

2.1 用到的接口

编写这个代码我们要用到一下几个红盒子工具：

• 空间分析工具 -- 数学分析 -- 减 (Minus)
• 空间分析工具 -- 数学分析 -- 绝对值(Abs)
• 空间分析工具 -- 提取分析 -- 按属性提取(ExtractByAttributes)
• 转换工具 -- 由栅格转出 -- 栅格转点 (RasterToPoint_conversion)

这个四个工具的接口定义如下：​​​​​​​

• Abs (in_raster_or_constant)

• Minus (in_raster_or_constant1, in_raster_or_constant2)

• ExtractByAttributes (in_raster, where_clause)

  <