Python のバックアップ(No.2)

Arcpyはじめのいっぽ

Arcpyはじめのいっぽ †

こんなときにはArcpyを使って処理するとべんりです。

バッチ処理をしたいとき
GUIで処理しようとしたら意味不明なエラーが出てきたとき
マシンのスペックが低くてGUIを起動するだけで精いっぱいなとき
結構簡単なのでお試しください。
python用のエディタを用意するとよいです(Pythonインストール時に付いてくるIDLEでもよいです)
テキストエディタでもよいですがpython用のエディタはインデント(pythonでは重要です)や関数などをわかりやすく整理・表示することができます。

シンプルな処理をしてみる †

# arcpyをインポートする
import arcpy
# 必要なエクステンションがあるかどうかをチェックする
# spatial analyst = "spatial"
# 3D analyst = "3D"
arcpy.CheckOutExtension("spatial")
# 入出力する変数を定義する
# パスとファイル名はダブルコーテーションで囲う
# パス区切りは\\で行う
hexagons_shp = "D:\\test\\hexagons.shp"
dem50gsi_j2ku54_img = "D:\\test\\dem50gsi_j2ku54.img"
output_table = "D:\\temp\\zst_slope_02.dbf"
# 処理コマンド中に上記変数を入力する
# 上記の変数以外の設定値はダブルコーテーションで囲って指定する
# どのような変数が必須かはHELPを参照する
arcpy.gp.ZonalStatisticsAsTable_sa(hexagons_shp, "Unique_ID", dem50gsi_j2ku54_img, output_table, "NODATA", "ALL")
# 処理が終了したらfinishという文字を表示させる
print "finish"

↑

たくさんのファイルをあつめてマージしてみる †

ワイルドカード(*)を用いて複数ファイルを指定して実行する手順

# globモジュールをインポートする
import glob
# arcpyをインポートする
import arcpy
# 出力ファイルのパスと名前を指定
outputdata = "C:\\temp\\ps_merge_all.shp"
# ワイルドカードを使用して入力データリストを指定
files = glob.glob('C:\\temp\\*.shp')
# Process: マージ（Merge）
arcpy.Merge_management(files, outputdata)
# 処理終了時にfinishと出力する
print = "finish"

↑

Arcpyはじめのいっぽ †

シンプルな処理をしてみる †

たくさんのファイルをあつめてマージしてみる †

バッチ処理をしてみる †