TINやラスターから断面データを取得する

EZ ProfilerというArcScript（フリーソフトみたいなもの）を紹介いたします。

EZ Profiler for ArcMap 9.x
http://arcscripts.esri.com/details.asp?dbid=13688

このプログラムはArcMapの機能を拡充させるためのプラグインで、Spatial Analystや3D Analyst等のエクステンションが無くても、polyline、grid、TINなどからデータを取得できるようになります。

TINやラスターデータから線分に沿った値を取得し、グラフ描画やエクセルへのエクスポートの出来るツールで、複数のレイヤからもまとめて取得することも可能。

【インストール】

メニューより、
ツール > カスタマイズ > コマンド
とすすみ、『ファイルから追加』をクリックし、『profile91.dll』を選択。

そうすると、『Easy Profiler 9.1』ツールバーがArcGISに追加されます。

【アンインストール】

アンインストールするときは、付属の『EZPuninstall.bat』を実行してください。

【使い方

PSボタン
細かい設定を行うためのボタン

Choose one method and input the value
> pointを発生させる個数の設定

　Specify the number of points along the profile
　> 指定した線などを等分割してポイントを発生させたい場合に選択

　Specity a set interval distance along the profile
　> 指定した距離でポイントを発生させたい場合に選択

Graph Properties
> グラフの設定

Make Profile From
> 解析の対象とするレイヤが１つかどうか
　既存polylineに沿ってラスタの値を取得したい場合などは、Multiple Layersを選択

Dボタン
ラスタ又はTINの上に、自由に線を書いて、ラスタから取得した属性値を持ったpointシェイプファイルを作成します。
使い方は押せば分かります。

mボタン
既存のポリラインに沿ってラスタから取得した属性値を持ったpointシェイプファイルを作成します。
使い方としては、既存ポリラインを選択してからボタンを押します。そうすると対象とするレイヤのダイアログが表示されますので、解析の対象とするレイヤを選択します。

Eボタン
選択したレイヤのテーブルをExcelに出力します。ArcGISの標準機能でexport等を行うよりも楽チンです。

その他細かいことは、zipファイルに入っているマニュアルに記載されています（英語ですが）。

大量のポイントデータを元にメッシュデータを作成する

ArcGISでは、大量のデータを含むシェイプファイルを扱うと、動きが緩慢になったり、下手すると勝手にソフトが終了してしまったりします。

そこで、大量のデータを扱う場合は、ラスタに変換すると良いでしょう。

まず準備として、『Spatial Analyst』でpointファイルをラスタに変換します。
また、適当なサイズのメッシュを『タイルポリゴン作成ツール』により作成します。さらに、各メッシュに適当なゾーンフィールド（ID番号のコピーで良い）を作ります。

ArcToolBoxより、
Spatial Analystツール > ゾーン > ゾーン統計
を起動

入力ラスタ、またはフィーチャーデータ欄
>メッシュファイル

ゾーンフィールド
>最初に作ったフィールド

入力値ラスタ
>最初に作ったラスタ

統計情報の種類（オプション）
>自分が必要とする情報を選ぶ。
　MEAN = 平均値
　MAJORITY = 最頻値
　MAXIMUM = 最大値
　MEDIAN = 中央値
　MINIMUM = 最小値
　MINORITY = 最も頻度が低い値
　RANGE = 最大値と最小値の幅
　STD = 標準偏差
　SUM = 合計
　VERIETY = ユニークな値の数

として、OKボタンを押す。

そうすると、メッシュの大きさが1つのセルになったラスタが作成されます。

ArcToolBoxより、
Spatial Analystツール > ゾーン > ゾーン統計をテーブルに出力
を起動

入力ラスタ、またはフィーチャーデータ欄
>メッシュファイル

ゾーンフィールド
>最初に作ったフィールド

入力値ラスタ
>ゾーン統計で作ったラスタ

として、OKボタンを押す。

そうすると、ラスタのセル毎に平均値などが算出されたテーブルが作成されますので、そのファイルをマップに追加してください。

ここまで来たら、あとは『最初に作ったメッシュ』に『ゾーン統計で作ったテーブル』をテーブル結合すると、膨大な数のデータをメッシュ毎に集計した結果を知ることができます。

シェープファイルのサイズを小さくする

ArcVIEWでは、１つのファイルとして扱うシェイプファイルは、エクスプローラーで見てみると、dbf,prj,sbn,sbx,shxなどのファイル群から構成されています。

ArcMap上で、レイヤを右クリックし、テーブルを表示するとズラズラズラズラと不要な属性値が入っていませんか？

CADのdxfファイル等をshpファイルに変換した場合、使わない属性値が沢山付属してくるため、拡張子がdbfのファイルサイズが抜けて大きくなります。

そこで、その属性値（フィールド）を一括で削除する方法を紹介いたします。

ArcToolBoxより、
データ管理ツール > フィールド > フィールドの削除

このコマンド（Delete Field）を使うと、テーブル上では1つづつしか削除できないフィールドを、一気にまとめて削除することができます。

カシミールの情報をArcGISで読み込む

カシミールで作ったルートなどの情報を、ArcGISで読み込むには、「gpx2shp」というフリーソフトを用います。
http://gpx2shp.sourceforge.jp/index.html.ja

gpx2shp は GPXファイル（GPS Exchange Format: *.gpx）を ESRIシェープファイルに変換するコンバーターで、下記HPよりダウンロードできます。数字が大きいほど最新版なので、それをダウンロードします。

http://sourceforge.jp/projects/gpx2shp/files/

使い方は、下記の通り超簡単です。

1. ダウンロードしたファイルを、適当なフォルダに解凍
2. カシミールで作ったルートやウェイポイントなどを、gpx形式で保存
3. 解凍したフォルダにあるgpx2shpのアイコンにgpxファイルをドラッグ
4. gpxファイルがあったフォルダに、シェイプファイル（dbf,shp,shx）が作成されます。

昔の地形と今の地形を比較する

昔と今で地形がどの程度違うのかを比較する場合、写真と地図のどちらを使っても比較することが出来ますが、写真の場合は位置あわせが面倒なので、私の場合は古地図を使います。

昭和初期の古い地形図は、国土地理院より購入することができます。
国土地理院（旧版地図）
http://www.gsi.go.jp/MAP/HISTORY/5-25/index5-25.html

古い1/25000地形図は、日本測地系で作成されているため、標準地域メッシュ・ポリゴン作成ユーティリティ（http://www.esrij.com/support/arcgis_desktop/downloads/ArcGIS_Standard.html）を使って、日本測地系の２次メッシュを作成し、四隅を合わせます。

そのままだと日本測地系で使いにくいので、ArcToolboxを使って世界測地系に変換しておくとよいでしょう。その際、ラスタ名が日本語だとエラーが出るようです。

データ管理ツール＞投影変換と座標変換＞ラスタ＞ラスタの投影変換

そうして作った古地図と最新の地図を比較すれば、地形がどの程度変化しているのか一目瞭然です。