ポリゴン フィーチャでクリップしたラスター レイヤーのセル最大値と最小値を取得

久しぶりに GIS らしい処理を ArcPy で実行したのでメモ。今回は、大きな範囲のラスター レイヤーの中にポリゴン フィーチャがいくつか存在し、各ポリゴン フィーチャの範囲内でセルの最大値・最小値を取得するというものです。

import arcpy

#以下を ArcGIS Pro の Python ウィンドウに貼り付けて実行

arcpy.env.addOutputsToMap = False #レイヤー追加を禁止

aprx = arcpy.mp.ArcGISProject("CURRENT") #現在のプロジェクトを取得
map1 = aprx.listMaps()[0] #最初のマップを取得

layer_target = map1.listLayers()[0]	#範囲ポリゴン レイヤー
layer_dem = map1.listLayers()[1]		#DEM レイヤー

cur = arcpy.da.SearchCursor(layer_target,["OID","SHAPE@"]

for row in cur:
  oid = row[0]
  geometry = row[1]
  extent = geometry.extent #範囲を取得
  in_raster = arcpy.Clip_management(layer_dem,str(extent),"in_memory\\" + str(oid).zfill(2),"#","#","NONE", "NOMAINTAIN_EXTENT") #ポリゴンの範囲でラスター レイヤーをクリップ

  #クリップしたラスター データセットのセル最小値と最大値を取得
  heightMax =  arcpy.GetRasterProperties_management(in_raster,"MAXIMUM")
  heightMin = arcpy.GetRasterProperties_management(in_raster,"MINIMUM")

  print("OID:{0} Max:{1} Min:{2}".format(oid, heightMax, heightMin))

これから ArcGIS Pro で処理できる場合は ArcPy も含め可能な限り ArcGIS Pro を使うようにしていきます。