2014年01月12日

林班の外周と中心をカシミール3Dで表示1

TNTmipsで,カシミール3D用の林班外周線(.e00)と林班中心点(.gdbのウェイポイント)を作成する手順をメモしておく。

例示する林班データは,森林GISからエクスポートしたShapeファイルをTNTmipsでインポートし,加工したものである。

まず初めにShapeファイルインポート時の注意点を書く。

Shapeファイルのインポート窓↓で,[次]をクリックし,
aインポート1.gif

開く窓で,「ポリゴナル」を選び,標準属性テーブルの生成にチェックを入れること↓。
aインポート2.gif
標準属性テーブルがないと,林班ポリゴンの中心座標をカシミール3Dへ渡せないので。

TNTmipsで林班ポリゴンデータを開いたところ↓。
b_mips_01.gif
左下が林班名が格納された属性データベース,右下(POLYSTATS)が,Shapeインポート時にチェックした,標準属性テーブルで,ポリゴンの面積や中心x座標(CentX),中心y座標(CentY)が格納されている。

このデータはこの時点ではJGD2000測地系の19座標系第VII系の座標を持っている↓。
b_mips_02.gif

これを,リプロジェクトモジュールで,WGS84測地系の経緯度座標系(ジオグラフィック)へ変換する↓。
b_mips_03.gif
カシミール3Dは.e00ファイルの座標値が,WGS84測地系の経緯度座標系であると解釈して描画するので,この変換が必要となる。

変換後のオブジェクトの範囲↓。
b_mips_04.gif

今度は右下のPOLYSTATSの中心座標値が,経度(CentX)と緯度(CentY)になった↓。
b_mips_05.gif

林班名が格納されているデータベースで,テーブル/プロパティをクリック↓。
b_mips_06.gif

テーブルプロパティ窓が開く↓。
b_mips_07.gif

左上のフィールドの追加アイコンをクリック↓
b_mips_08.gif

フィールド名を,lat(緯度)とし,フィールドタイプを数値処理,幅10,小数点以下5桁まで表示と設定し,スクリプトの編集をクリック↓。
b_mips_09.gif

開いたスクリプト編集窓で,フィールドの挿入アイコンをクリック↓。
b_mips_10.gif

こんな窓が開く↓。
b_mips_11.gif

Internal右側の下向き矢印をクリックし,POLYSTATSを選ぶ↓。
b_mips_12.gif

Areaの右側矢印をクリックし,CentY(緯度なのでY座標)を選ぶ↓。
b_mips_13.gif

選び終わったところで挿入をクリック↓。
b_mips_14.gif

スクリプト編集窓がこんな風になる↓。
b_mips_15.gif
これは,データベース名とフィールド名をピリオドで連結したもの。

latフィールドに,中心y座標を入れる準備ができたので,okを押す↓。
b_mips_16.gif

林班データベースにlat列が追加され,中心の緯度値が転写されたところ↓。
b_mips_17.gif
latの文字が青色で表示されている。

latをマウスで右クリックし,計算フィールドを固定化を選ぶ↓。
b_mips_18.gif

開く窓で「はい」をクリック↓
b_mips_19.gif

latの文字が黒色で表示され,中心緯度値が固定化された↓。
b_mips_20.gif

再びテーブルプロパティでフィールドの追加を選び↓,
b_mips_21.gif

今度はlon(経度)と名前を付け,スクリプト編集をクリックし↓,
b_mips_22.gif

CentX(経度なのでX座標)を指定したところ↓。
b_mips_23.gif

lonフィールドが計算されたところ↓。
b_mips_24.gif

lonフィールドを固定したところ↓。
b_mips_25.gif

テーブル/別名保存をクリック↓。
b_mips_26.gif

フォーマットはcsv(カシミール3DのGPSファイルツールプラグイン用),見出し左にチェックしてok↓。
b_mips_27.gif

ここでは,「魚津_林班中心」と命名してok↓。
b_mips_28.gif

エクスポートモジュールで,ARC-E00-VECTOR(e00)を選んで次へ↓。
b_mips_29.gif

WGS84の経緯度に変換した林班ポリゴンを選んでOK↓。
b_mips_30.gif

エクスポートパラメータ窓で,ラインにチェックを入れて,出力↓。
b_mips_31.gif
カシミール3Dはe00のラインにのみ対応している。

「魚津_林班外周」と命名してOK↓。
b_mips_32.gif

エクスポート中↓
b_mips_33.gif

エクスポート終了↓
b_mips_34.gif

ここまでで,
・魚津_林班中心.csv(さらに加工が必要)
・魚津_林班外周.e00
の2つのファイルができた。

第1部おわり

2014年01月11日

森林計画図をカシミール3Dで表示

講演会のスライド用に,キャプチャ画像をいくつか作ったので,こっちにも載せておきます。

TNTmipsで,森林計画画像をカシミール3D用画像に変換する手順の抜粋です。

スキャンした地図画像の四隅に位置情報を与えたところ↓。
立山計画図_01.gif

左下隅を拡大したこところ↓。
立山計画図_02.gif
平面直角座標系,第VII系の,
東西座標=16(km),南北座標=66(km),第VII系の原点は岐阜県内にある。

座標値はメートル(m)単位で入力する↓。
立山計画図_03.gif

地図画像本体(図郭内)を切り出したところ↓。
立山計画図_04.gif
実際は,切り出す前に一度リサンプルを掛けて,スキャン時の地図の方向ずれ(紙が少し斜めにスキャナに吸い込まれてしまった場合など)を修正している。

立山町全体のモザイク(接合)前↓
立山計画図_05.gif

モザイク後↓
立山計画図_06.gif

この状態では,計画図オリジナルの,Tokyo測地系,19座標系,第VII系になっている↓。
立山計画図_07.gif

これを,カシミール3Dが唯一サポートする,測地座標系(geographic)に変換する。
リサンプルモジュールでモザイク画像を選び↓
立山計画図_08.gif

出力参照系として,Tokyo測地系のgeographic(ジオグラフィック)を選ぶ↓。
立山計画図_09.gif
リサンプル処理は,全ての画素を並べ替えるので,時間が掛かるが,必須の処理である。

リサンプル後のモザイク画像の範囲↓。
立山計画図_10.gif
青四角で囲んだ,左上と右下の緯度,経度を後でカシミール3Dで手入力する必要がある。

ビットマップ形式(.bmp)でエクスポートするところ↓。
立山計画図_11.gif
実際は,エクスポート前に,白黒1bit(白=0,黒=1)を,グレースケール8bit(白=0,黒=255)に変換している。これは,カシミール3Dが1bitのビットマップ画像を読めないからである。

カシミール3Dから,ビットマップの計画図を初めて開くところ↓。
立山計画図_12.gif

先ほどの緯度,経度値を赤四角内に手入力する↓。
立山計画図_13.gif
地図測地系はTokyoを選ぶ。

立山町の森林計画図が開いたところ↓
立山計画図_14.gif
ここで,市町村界のe00ファイルを重ねて表示し,境界付近を眺め回すと,先ほどの緯度,経度の入力間違いがないかどうかチェックできる。
この後,bmpをcmpに圧縮してやれば,ファイルサイズが小さくなる。

以上,だいぶ端折ったが,おわりです。

2014年01月06日

道路と河川のe00

地理院の基盤地図情報のサイトから道路や河川のベクトルデータをダウンロードして,地理院のコンバートソフト(FGDV)でxml→shapeに変換し,TNTmipsでshape→e00にすればカシミール3DのGISツールプラグインで表示できるのですが,道路が2重線(道路縁),河川も2重線(水涯線)で表示されるのがうっとうしいです。

道路も河川も1本線(道路中心線,河川中心線)のデータで,有料で,しかも多少古くても構わない場合には,数値地図25000(空間情報基盤)が便利です。

大多数の都道府県版が売り切れですが,日本地図センターから復刻版が購入できます↓。
http://www.jmc.or.jp/data/25000spa/juchu.html

これらのCDからカシミール3Dで表示するまでの手順をメモしておきます。

まずは元データのCDの写真から↓
a_CDの写真.gif
CDの日付は平成14年6月1日です。(6月以降か以前かで少し違うらしい)

(1)地理院のデコーダ:生データ→xml
(2)ESRIジャパンの変換ソフト:xml→shape
(3)TNTmips:shape→rvc→e00
(4)カシミール3D:GISツールプラグインでe00を表示

(1)CDの中のDECODE.exeを起動する↓。
bデコード_01.gif

解凍ソフトの初期画面で次へ↓
bデコード_02.gif

参照ボタンを押して,CD内のCSVファイルを選ぶと,DTDとデータフォルダが自動で設定され,出力フォルダを指定する↓。
bデコード_03.gif

ファイルセットから全てを選び,右向き矢印で右へ送って次へ↓。
bデコード_04.gif

地物選択の画面でそのまま次へ↓
bデコード_05.gif

出力ファイル画面で,地物別ファイルとシフトJISをチェックして次へ↓
bデコード_06.gif
ESRIジャパンの変換ソフトのマニュアルに,「シフトJIS」を選べと書いてある。

変換が始まる↓。
bデコード_07.gif

変換中のプログレスバー↓。
bデコード_08.gif

しばらく(20分くらい?)待つと終わる↓。
bデコード_09.gif

出力されたxmlファイル群↓。
bデコード_10.gif

(2)ESRIジャパンの変換ソフトを起動し,数値地図25000(空間データ基盤)を選ぶ↓。
cコンバータ_1.gif

xmlファイルが入っている,親フォルダを指定してok↓。
cコンバータ_2.gif

開くウィンドウで全て追加をクリック↓。
cコンバータ_3.gif

出力フォルダを指定し,6月以降にチェックがあるのを確認して変換をクリック↓
cコンバータ_4.gif
CDの日付が平成14年6月1日なので,デフォルトの6月以降でOK

xml→shapeの変換中(席を外して別の用事をした方が良い)↓。
cコンバータ_5.gif

変換終わり↓
cコンバータ_6.gif

出力されたshapeファイル群↓。
cコンバータ_7.gif

(3)TNTmipsでshapeをインポート↓
d_mips_01.gif

道路(道路区間)を表示したところ↓
d_mips_02.gif

河川(河川区間)を表示したところ↓
d_mips_03.gif

データは,JGD2000測地系の経緯度座標系(測地座標,地理座標)です↓。
d_mips_04.gif

これを念のため,WGS84測地系の経緯度座標系に変換し,e00でエクスポートする。
エクスポートはこの記事参照↓。
http://gpsrsgis.seesaa.net/article/383205920.html

(4)カシミール3DのGISツールプラグインで道路を開いて色を変えたところ↓
e_k3d_道路e00.gif

河川の色と線の太さを変えたところ↓
e_k3d_河川e00.gif

以上で説明終わり。