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 函數來擴展計算器的功能。