オルソ手順書_はじめに

簡易オルソ作成手順書のうち，「第1章はじめに」を書きました↓。


Wordで目次をちゃんと作っておくと，pdfにしてもしおり？になってくれて便利ですよ。

地理院DEMをdcm形式へ(2)

(1)の最後に示した，接合後のDEMの左上の経緯度（JGD2000測地系）を，Tokyo測地系の経緯度へ変換しておきます。変換にはTNTmipsのMap Calculatorを使います↓。

東経136度45分00.000秒（JGD2000）
→136度45分10.746秒（Tokyo）

北緯37度00分00.000秒（JGD2000）
→36度59分49.196秒（Tokyo）
と，概ね11秒ずつずれます。

（注意）ここで作成しているDEMは，後にカシミール3Dから”USGS30秒/GTOTP30”形式で読み込みます。
USGSは米国地質調査所のことなので，左上の経緯度はWGS84（≒JGD2000）の経緯度を指定すればよい，つまり変換せずに地理院DEMの元々の経緯度を使えばよいと思い込み，変換なしのそのままの経緯度を使っていたら，何度やっても山の頂上が400〜500mずれ，途方に暮れていました。
で，左上の経緯度の値を，上述のように，Tokyo測地系（旧日本測地系）の値に変えてやったらうまく山頂がDEMの極大地点と一致しました。どうやらカシミール3Dは，”USGS30秒/GTOTP30”形式のDEMの経緯度を，Tokyo測地系のものに決め打ちして読み込むようです。これが最大の注意点です。

接合したDEMは，SIMPLE-ARRAY形式でエクスポートします↓。

Byte Orderは，High-Lowの順番を指定します。

バイト並びの根拠は赤い下線の通りです↓。

カシミール3DのQ&Aより

フォルダを作成し（ここではexp_simple_array），ファイル名はデフォルトの”MOSAIC”のままOKを押します↓。


エクスポート後のファイルのエクスプローラ表示です↓。

拡張子はなく，容量は約130MBです。

ここでExcelを使います↓。

赤枠内は，1画素当たりの経度の増分（XDIM，10進数の度単位）を計算しています（前掲Q&A図参照）。
同様に，1画素当たりの緯度の増分も計算しておきます。

また，Map Calculatorで算出した，左上隅のTokyo測地系での経緯度を，10進数度単位へ変換しておきます↓。


これまで得た情報を元に，以下のヘッダファイルをテキストエディタで作成します↓。

項目の説明は前掲Q&Aを参照してください。

DEM本体のMOSAICに，拡張子".bil"を追加します↓。


(3)へつづく。

地理院DEMをdcm形式へ(1)

国土地理院の10mメッシュDEMは，基盤地図情報（標高）プラグインでカシミール3D用に変換できますが，ファイル数が多くなる（2次メッシュ単位（＝1/25,000地形図図幅単位））ので，管理が面倒だと感じるときがあります。

そこで今回は，この地理院の10mDEMの富山県全体を含む範囲を，1つの.dcmファイル（圧縮された標高データ）に変換する方法を載せておきます。

今回紹介する変換には，TNTmipsのプロ版（約80万円）が必要ですが，自分用の備忘録も兼ねているので，持っていない方，ごめんなさい。

まずはTNTmipsの日本代理店のホームページからDEMをダウンロードします↓。


富山県全体を含む2次メッシュの矩形範囲（私はこれを「小富山」，「ST」と呼ぶ）のDEMファイル群です↓。77個あります。


TNTmipsで全てを表示したところです↓。


ひとつの画像の範囲情報です↓。地理院オリジナルと同じ，JGD2000測地系の経緯度座標系です。


モザイク（接合）用コントロール窓(Input)です↓。


モザイク用View窓です↓。


モザイク用コントロール窓(Output)です↓。

元データは浮動小数点ですが，カシミール3Dへは16bit整数で渡すので，ここで16bit符号なし整数を選びます。

新しいファイルを作り（地理院dem10.rvc），デフォルトのMOSAICというオブジェクト名を指定します↓。

[OK]を押します。

はい，接合できました↓。


接合後の画像範囲と行列数を表示して↓，ハードコピーを取ります。後で必要です。


(2)へつづく