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

      2018/01/12

久しぶりに GIS らしい処理を ArcPy で実行したのでメモ。これから ArcGIS Pro で処理できる場合は ArcPy も含め可能な限り ArcGIS Pro を使うようにしていきます。

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

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_manageme  nt(in_raster,"MINIMUM")

  print("OID:{0} Max:{1} Min:{2}".format(oid, heightMax, heightMin))<span data-mce-type="bookmark" style="display: inline-block; width: 0px; overflow: hidden; line-height: 0;" class="mce_SELRES_start"></span>

 - プログラミング, ArcGIS ,