ArcGIS 栅格计算器

一、概述

ArcGIS 栅格计算器是一个广泛使用的基于栅格数据的计算工具，它可以帮助用户进行各种复杂的栅格数据处理和分析。使用栅格计算器，可以对栅格数据进行加、减、乘、除、逻辑运算等计算，同时也可以进行复杂的数学和统计计算。

二、基本用法

在 ArcGIS 软件中打开需要进行栅格计算的地图，然后打开栅格计算器窗口。在窗口中，用户可以输入需要计算的表达式，表达式中可以使用栅格图层的名称、BMP 单元、常量等元素进行运算。下面是一个简单的栅格计算器表达式示例：

("Elevation" - 1000) * 0.1

这个表达式将栅格图层 "Elevation" 中的所有单元格值都减去 1000，然后乘以 0.1。

三、常用函数

ArcGIS 栅格计算器中提供了众多的函数，可以方便地进行各种栅格计算和统计分析。下面是一些常用函数示例：

• Sin()：计算给定角度的正弦值。

  Sin("Aspect" * 3.14159265 / 180)
  

• Log10()：计算给定数字的以 10 为底的对数。

  Log10("Slope" * 100 + 1)
  

• Exp()：计算给定数字的以 e 为底的指数。

  Exp("Slope" / 100)
  

• Raster()：将栅格数据类型转换成其他类型。

  Raster("Elevation" > 1000)
  

四、条件函数

ArcGIS 栅格计算器中提供了多个条件函数，可以根据单元格值进行判断和计算。下面是一些条件函数示例：

• Con()：根据条件返回结果，如果条件为真，则返回真部分的值，否则返回假部分的值。

  Con(("Elevation" >= 500) & ("Elevation" < 1000), "Elevation" * 2, "Elevation" / 2)
  

• SetNull()：根据条件将指定单元格的值设为 NoData。

  SetNull("Elevation" < 0, "Elevation")
  

• FocalStatistics()：对指定领域内的值进行统计分析。

  FocalStatistics("Elevation", NbrCircle(3, "CELL"), "MEAN")
  

五、用户自定义函数

如果 ArcGIS 栅格计算器中提供的函数不能满足用户的需求，用户还可以通过编写 Python 脚本来自定义函数。用户自定义函数可以在 ArcGIS 环境中被调用，并且可以像普通函数一样进行计算。下面是一个用户自定义函数的示例：

import arcpy

def myfunc(in_raster):
  mean_value = arcpy.GetRasterProperties_management(in_raster, "MEAN")
  return mean_value

my_raster = "Elevation"
output_raster = arcpy.sa.ExtractByMask(my_raster, "Mask")
result = myfunc(output_raster)

六、总结

ArcGIS 栅格计算器是一个功能强大的栅格数据处理和分析工具。它提供了基本运算、函数计算、条件计算等多种功能，可以满足用户各种不同的需求。同时，用户还可以通过自定义 Python 函数来扩展计算器的功能。