普段の生活で”モザイク処理”というと、”公然と見せてはいけないものをぼかす”ような処理を想像してしまいますが、GISでは複数の画像をくっつける処理のことを言います。

ArcToolboxより、モザイクと検索すると幾つかのコマンドが出てきますので、お好みの方法で処理してください。
関連しているかもしれない記事：

• 表示画面を回転させる
• ラスタデータを複雑な形で抽出する方法
• AutoCADで作成した図面を画像に変換する
• 配置した画像のレクティファイ
• ラスタの色を薄くする

Create Raster Index Shapefile

インストールは、インストーラーがついていますので、setup.exeをダブルクリックし、インストール先フォルダを選択するだけです。

ボタンを表示するには、
ツール＞カスタマイズ＞コマンド＞Query Tools＞Create Raster Index Shapefile
とすすみ、適当な所にボタンを配置します。

設定は以下の通り、誰でも分かる簡単な内容です。

• Input Folder：ラスタが入っているフォルダを選択
• Output Shapefile：保存先及び名前を指定
• Process Subfolders：サブフォルダも検索してくれます
• Project polygons WGS84：座標系をWGS84に変換してくれます

あとは、Create Indexボタンを押せば、サブフォルダにあるラスタデータも全て統合した、使いやすいインデックスポリゴンを作成してくれます。進捗率だけではなく、どのファイルを処理中かも表示されますので、本当に作業中なのか不安になることがありません。

で、出来上がるシェイプファイルも扱いやすい。ジオメトリが壊れているのでそのままでは使えなかったり、パス名が長いと途中でパス名が切れてしまうような、困った仕様ではありません。

属性テーブルも扱いやすいです。

• Shape：当然Polygon
• NAME：ファイル名
• PATH：ファイルの保存先（ファイルパス）＆ファイル名
• FORMAT：ファイル形式（TIFFとかJPEG）
• COMP：圧縮タイプ（無圧縮ならNone）
• PYRAMID：ピラミッドの有無（作成済みだとTrue）
• BANDS：バンド数
• PIXEL_TYPE：画素タイプ（8ビットとか）
• SP_REF：そのラスタが設定されている空間参照
• UNIT：単位
• SIZE_X：ピクセルサイズX（Ｘ方向の画素数）
• SIZE_Y：ピクセルサイズY（Ｙ方向の画素数）
• COLS：何かを基準にした列数
• ROWS：何かを基準にした行数
• X_MIN：X座標の最小値
• X_MAX：X座標の最大値
• Y_MIN：Y座標の最小値
• Y_MAN：Y座標の最大値

これだけの情報がリスト化されていれば、属性で表示すべきラスタをフィルタリングしたり、ハイパーリンクを付けたりが容易ですね。コレは、誰にでもお勧めできて、ぜひ常用したいArcScriptです。
関連しているかもしれない記事：

• 複数のラスタの配置図を作る
• フォルダ内のテキストファイルを全て結合する方法
• ラスタデータを複雑な形で抽出する方法
• ラインを指定間隔に分割する
• 数値地図50mメッシュ標高を使う

そこで、３次元の点の集まりであるレーザープロファイラデータの処理に特化したフリーソフトの紹介です。また、同サイトでは、xyz座標が記載されたcsvテキストファイルをASCII Gridに変換するフリーソフトも公開されています。ArcGISのプラグインではなく、単体で動くソフトです。

このソフトが無料とは信じられません。これはすごいです。

以下に紹介するどちらのソフトも、日本語化されていませんが、ごくごく簡単な英語なので、誰でも使えます。また、ライセンス取得等も必要がないので、ソフトをダウンロードしたらすぐに使えます。インストール自体は、特にオプションはないので迷わず行えます。

■超高速3次元ポイントビューアー
LViz – LiDAR / ALSM 3D Visualization tool
http://lidar.asu.edu/LViz.html (LViz-1.4.1.1.exe)

●LVizの使い方

1. x,y,z座標が記載されたテキストを用意します。
2. LVizを起動します。
3. メニューの『Import』から、『Import Multi-Column Point Cloud File』を選択
   Filename　＞　座標を記載したファイルを指定
   File Delimeter　＞　csvファイルの区切りを指定
   Get Headersをクリック
   X,Y,Z　＞　表示したい列を選択する。csvファイルにidとかが記載されている場合は、適宜指定する。
4. Importが終わると、ArcGISやAutoCADでは考えられないような描画スピードで、膨大なポイントデータを表示してくれます。

■テキストからラスタ（ASCII Grid）に高速変換するソフト
GEON Points2Grid Utility
http://lidar.asu.edu/points2grid.html (GEON-points2grid-Utility-Setup-v1_3.exe)

●GEON Points2Gridの使い方

1. x,y,z座標が記載されたテキストを用意します。
2. GEON points2Gridを起動します。
3. いきなり入力ダイアログが表示されるので、適宜入力する。
   Input File　＞　上記テキストを指定
   Grid Type　＞　Min,Max,Mean,IDW,Point Countのうち好きなのを選択
   Algorithm Parameters　＞　Grid Resolution　例えば10mメッシュが欲しければ10と入力
   Search Radius Value　特に変更不要
   Output Directory　＞　出力先フォルダを指定
   Output File ＞　ファイル名を指定
4. あとは、OKボタンを押すだけで、あっという間にラスタ変換が終わります。

これでフィルタリング処理もできれば最高ですが、フリーソフトでここまで出来れば文句ありませんね。
関連しているかもしれない記事：

• FLASHから画像を抽出する
• ラスタデータを複雑な形で抽出する方法
• 配置した画像のレクティファイ
• 国土交通省がAjax地図ビューアを利用した「航空写真画像情報所在検索・案内システム」
• ArcGISでアスタリスクを使ってデータを検索する

この作業には、ArcGISとSpatialAnalystというエクステンションを使用します。

○必要なデータ
ベースとなるラスタ（数値地図・航空写真・衛星写真・DEMなど）と抽出したい形のポリゴンシェープファイル（どんなに複雑な形状でもOK）を用意する。

○作業手順
ArctoolBoxより、Spatial Analystツール＞抽出＞マスクで抽出（Extract by Mask）を実行。

【入力ラスタ】欄に、先ほどのラスタを指定
【入力ラスタ、またはフィーチャマスクデータ】欄に、先ほどのshapeファイルを指定
【出力ラスタ】欄は、出力先なのでお好きなフォルダ＆ファイル名を指定

これだけの手順で、ポリゴンと重なった部分のみで作られたラスタデータを作成することができます。
関連しているかもしれない記事：

• PDFに使われているイラストを抽出する
• 標高データから、河川および流域界を抽出する
• FLASHから画像を抽出する
• ラスタ変換の方法は目的にあわせて変える
• ラスタの色を薄くする

昭和初期の古い地形図は、国土地理院より購入することができます。
国土地理院（旧版地図）

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

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

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

• 国土基本図や森林基本図を背景にしたい
• 表示画面を回転させる
• 数値地図を読み込んで、簡単につなげる
• 既存のマップに数値地図画像を追加配置するとずれる
• CADデータや数値地図のデータの位置関係が合わない