<\/HEAD>
これを動かすことはできましたが、元のpython構文のどこに問題があったのかまだわかりません。<\/P>
<\/P>
ある地域でいくつか提案された変更に対して表流水解析を実行し、洪水をどれだけ防げるかを確認したかったので、2015年のプロジェクトからpythonスクリプトを入手し、それを修正して別のワークスペースで実行し、複数の異なる表面でこれらのステップを実行できるようにしました。 <\/P>
<SPAN class="comment token"># -*- coding: utf-8 -*-<\/SPAN>
<SPAN class="comment token"># ---------------------------------------------------------------------------<\/SPAN>
<SPAN class="comment token"># Fill2Stream.py<\/SPAN>
<SPAN class="comment token"># 作成日: 2015-03-06 10:28:41.00000<\/SPAN>
<SPAN class="comment token"># (ArcGIS/ModelBuilderによって生成)<\/SPAN>
<SPAN class="comment token"># 説明: 標高ラスターから流向および蓄積グリッドを生成するため。 <\/SPAN>
<SPAN class="comment token"># 2017-11-28にPHによってMud Lakeスピルウェイ解析用にこれらのステップを実行するよう修正<\/SPAN>
<SPAN class="comment token"># ---------------------------------------------------------------------------<\/SPAN>
<SPAN class="comment token"># arcpyモジュールのインポート<\/SPAN>
<SPAN class="keyword token">import<\/SPAN> arcpy
<SPAN class="keyword token">from<\/SPAN> arcpy<SPAN class="punctuation token">.<\/SPAN>sa <SPAN class="keyword token">import<\/SPAN> <SPAN class="operator token">*<\/SPAN>
<SPAN class="comment token"># 必要なライセンスのチェックアウト<\/SPAN>
arcpy<SPAN class="punctuation token">.<\/SPAN>CheckOutExtension<SPAN class="punctuation token">(<\/SPAN><SPAN class="string token">"spatial"<\/SPAN><SPAN class="punctuation token">)<\/SPAN>
arcpy<SPAN class="punctuation token">.<\/SPAN>env<SPAN class="punctuation token">.<\/SPAN>workspace <SPAN class="operator token">=<\/SPAN> <SPAN class="string token">"G:\/elevation\/lidar_2005\/MudLake\/MudLakeTerrians.gdb"<\/SPAN>
arcpy<SPAN class="punctuation token">.<\/span>env<SPAN class="punctuation token">.<\/span>overwrite <span class="operator token">=<\/span> <span class="string token">"True"<\/span>
<span class="comment token"># ローカル変数:<\/span>
<span class="comment token">#DEMrs = "G:\\elevation\\lidar_2005\\MudLake\\MudLakeTerrians\\DEMFloodway"<\/span>
DEMrs <span class="operator token">=<\/span> <span class="string token">"G:\elevation\lidar_2005\MudLake\MudLakeTerrians.gdb\DEMFloodway"<\/span>
<span class="comment token">#DEMrs = "DEMFloodway"<\/span>
Fill <span class="operator token">=<\/span> <span class="string token">"DEMFWfill"<\/span>
flowdir <span class="operator token">=<\/span> <span class="string token">"DEMFWdir"<\/span>
drop_raster <span class="operator token">=<\/span> <span class="string token">"DEMFWDrop"<\/span>
Accum <span class="operator token">=<\/span> <span class="string token">"DEMFWaccum"<\/span>
<span class="keyword token">print<\/span> DEMrs
<span class="comment token">#処理: Fill<\/span>
<span class="comment token">#outFill = Fill(DEMrs)<\/span>
<span class="comment token">#outFill.save(Fill)<\/span>
arcpy.gp.Fill_sa(DEMrs, Fill, "")
<span class="comment token"># 処理: Flow Direction<\/span>
arcpy.gp.FlowDirection_sa(Fill, flowdir, "NORMAL", drop_raster)
<span class="comment token"># 処理: Flow Accumulation<\/span>
arcpy.gp.FlowAccumulation_sa(flowdir, Accum, "", "FLOAT")
<span class="comment token"># ストリームラインラスターの処理<\/span>
strm = Con(Raster(Accum) > 200000, 1, 0)
strm10 = Con(Raster(Accum) > 10000, 1, 0)
strm.save("FWstrm")
strm10.save("FWstrm10")
<span class="line-numbers-rows"><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span><\/span><span></ span >< span >< < EMOJI_19 > >< \/ span >< span >< < EMOJI_20 > >< \/ span >< span >< < EMOJI_21 > >< \/ span >< span >< < EMOJI_22 > >< \/ span >< span >< < EMOJI_23 > >< \/ span >< span >< < EMOJI_24 > >< \/ span >< span >< < EMOJI_25 > >< \/ span >< span >< < EMOJI_26 > >< \/ span >< span >< < EMOJI_27 > >< \/ span >< span >< < EMOJI_28 > >< \/ span >< span >< < EMOJI_29 > >< \/ span >< span >< < EMOJI_30 > >< \/ span >< span >< < EMOJI_31 > >< \/ span >< span >< < EMOJI_32 > >< \/ span >< span >< < EMOJI_33 > >< \/ span >< span >< < EMOJI_34 > >< \/ span >< span >< < EMOJI_35 > >< \/ span >< span >< < EMOJI_36 > >< \/ span >< span >< < EMOJI_37 > >< \/ span >< span >< < EMOJI_38 > >< \/ span >< span >< < EMOJI_39 > >< \/ span >< span >< < EMOJI_40 > >< \/ span >< span >< < EMOJI_41 > >< \/ span >< span >< < EMOJI_42 > >( \ / )