ArcGIS大师之路500技---031栅格计算器-育师

文章目录

前言
一、基本概念
二、常用函数
- 2.1 数学函数
- 2.2 统计函数
- 2.3 条件函数
- 2.4 逻辑函数
- 2.5 重分类函数
- 2.6 表面分析函数
- 2.7 距离函数
- 2.8 水文分析函数
- 2.9 插值函数
三、最佳实践建议
总结

前言

ArcGIS栅格计算器（Raster Calculator）是空间分析工具箱中的强大工具，用于对栅格数据进行数学运算和逻辑分析。

一、基本概念

作用
执行栅格数据的数学运算（加、减、乘、除等）
进行条件判断和逻辑运算
多个栅格数据的叠加分析
地图代数表达式计算
访问方式
ArcMap: 空间分析工具 → 地图代数 → 栅格计算器
: 通过arcpy.sa.RasterCalculator()调用
栅格引用方式
直接双击图层列表中的栅格名称
手动输入栅格名称（需加引号）
使用完整路径：“C:/data/dem.tif”

二、常用函数

2.1 数学函数

基本数学运算

# 绝对值Abs("temperature_change")#平方根Sqrt("area")#幂运算Power("radius",2)# 半径的平方Square("distance")# 平方#对数Log("population")# 自然对数Log10("light_value")# 以10为底的对数Exp("growth_rate")# 指数函数

2 三角函数

#角度转弧度("angle_degree"*3.1415926535)/180#三角函数计算Sin("slope_radians")# 正弦Cos("aspect_radians")# 余弦Tan("angle_radians")# 正切#反三角函数ASin("value")# 反正弦ACos("value")# 反余弦ATan("value")# 反正切实例：计算日照强度 假设slope为坡度，aspect为坡向，latitude为纬度 slope_rad=("slope"*3.14159)/180aspect_rad=("aspect"*3.14159)/180lat_rad=(30*3.14159)/180# 北纬30度#简化的日照强度计算solar_radiation=Cos(slope_rad)*Cos(lat_rad)+Sin(slope_rad)*Sin(lat_rad)*Cos(aspect_rad)

2.2 统计函数

局部统计

#邻域统计FocalStatistics("elevation",NbrRectangle(3,3),"MEAN")# 3x3窗口均值FocalStatistics("rainfall",NbrCircle(5,"CELL"),"MAX")# 5像元半径最大值#块统计BlockStatistics("ndvi",NbrRectangle(10,10),"MEAN")# 10x10块统计

区域统计

#分区统计ZonalStatistics("watershed","VALUE","slope","MEAN")# 流域平均坡度ZonalStatistics("landuse","VALUE","temperature","RANGE")# 土地利用类型内温度范围

2.3 条件函数

Con函数（最常用）

#基本语法Con(condition,true_value,false_value)Con(condition,true_value,false_value,where_clause)#简单重分类Con("elevation">1000,1,0)# 高于1000米为1，否则为0#多条件重分类Con("landuse"==1,"耕地",Con("landuse"==2,"林地",Con("landuse"==3,"水域","其他")))

Pick函数

#根据位置索引选择值Pick("index_raster",["raster1","raster2","raster3"])#实例：根据月份选择温度数据month="month_index"# 值从1到12temperature=Pick(month,["jan_temp","feb_temp","mar_temp","apr_temp","may_temp","jun_temp","jul_temp","aug_temp","sep_temp","oct_temp","nov_temp","dec_temp"])

2.4 逻辑函数

比较运算

#等于"landuse"==1# 不等于"slope"!=0# 范围判断("temperature">=10)&("temperature"<=30)# 多个条件组合("soil_type"==3)&("slope"<15)&("aspect">90)&("aspect"<270)

布尔运算

# 与运算("rainfall">1000)&("temperature">15)# 或运算("elevation"<50)|("distance_sea"<10)# 非运算~("protected_area"==1)# 非保护区# 异或运算("urban"==1)^("industrial"==1)# 城市或工业区，但不是两者都是

2.5 重分类函数

Reclassify函数

# 直接重分类（在ArcGIS Pro中）Reclassify("slope",RemapRange([[0,5,1],[5,15,2],[15,25,3],[25,90,4]]))# 使用栅格计算器实现类似功能slope_class=Con("slope"<=5,1,Con("slope"<=15,2,Con("slope"<=25,3,4)))

自定义重分类

# 连续值离散化temperature_class=Int(("temperature"+10)/5)# 每5度一个等级# 归一化到0-1范围normalized=("value"-FocalStatistics("value",NbrRectangle(100,100),"MINIMUM"))/(FocalStatistics("value",NbrRectangle(100,100),"MAXIMUM")-FocalStatistics("value",NbrRectangle(100,100),"MINIMUM"))

2.6 表面分析函数

地形分析

# 坡度计算（需要Spatial Analyst扩展）Slope("dem","DEGREE")# 度为单位Slope("dem","PERCENT_RISE")# 百分比坡度# 坡向计算Aspect("dem")# 曲率计算Curvature("dem")

2.7 距离函数

欧氏距离

# 到最近道路的距离EucDistance("roads")# 成本距离CostDistance("source","cost_raster")# 成本路径CostPath("destination","cost_distance","backlink")

2.8 水文分析函数

基本水文分析

# 流向计算flow_dir=FlowDirection("dem","NORMAL")# 汇流累积量flow_acc=FlowAccumulation(flow_dir)# 盆地（流域）划分basin=Basin(flow_dir)# 水流长度flow_length=FlowLength(flow_dir,"DOWNSTREAM")

2.9 插值函数

空间插值

# IDW反距离权重插值Idw("rain_gauges","RAINFALL",2000,2)# 克里金插值Kriging("temperature_stations","TEMP","Spherical")# 样条函数插值Spline("elevation_points","ELEVATION","REGULARIZED")