[OSGeoJapan-discuss] 基盤地図情報数値標高（GML）の変換結果について

[Nobusuke Iwasaki]

赤木様

岩崎です。
まず、分類結果がおかしくなる部分ですが、私の勘違いで、10m DEMでなく、50m DEMで行ったときに見られる結果でした。10m
DEMでは問題ありませんでした。
なお、参考までに対象ですが、千葉県を対象にしたときに、平野部でbinary positiveとbinary negativeの結果が傾斜があるとなっていました。

なお、同様の地域を10mDEMで分析したところ、50mDEMで見られたような現象はありませんでした。

また、誤差についてもありがとうございます。
確かに、上記の通り10m DEMでは起きませんでしたので、おっしゃるとおりこの誤差が問題ではありませんでした。

以上、確認とお礼までに。


2015年4月29日 11:39 Minoru Akagi <akaginch ＠ gmail.com>:
> 岩崎さま
>
> 山地地域の地形データを使ってr.geomorphonモジュールを試してみましたがラスタデータを構成する下位のメッシュの境界部で地形分類結果がおかしくなっているといったことは見られませんでした。図郭部が地形分類結果に斜面として表れたというのは平野の地形データでしょうか?
>
> ちなみにご指摘の誤差についてはデータ変換ソフトで標高値を4バイト(32ビット)の浮動小数点数型の変数に格納したり、このデータ型のGeoTIFFファイルに保存したりするために生じる誤差です。Pythonでは次のようにすればこの誤差が生じることを確認することができます。
>
>>>> import struct
>>>> a = 439.40
>
> # 数値を32ビット浮動小数点数バイト列へ変換
>>>> b = struct.pack("f", a)
>>>> b
> '3\xb3\xdbC'
>
> # 数値へ戻す
>>>> struct.unpack("f", b)
> (439.3999938964844,)
>
> この誤差はごく小さいので、境界部の分類結果がおかしい原因がこれだとは言えないのではないかと個人的には思います。
>
> 赤木実
>
> 2015年4月24日 17:53 Nobusuke Iwasaki <wata909 ＠ gmail.com>:
>> リストの皆様
>>
>> いわさき＠OSGeo.JPです。
>> 基盤地図情報数値標高（GML）の変換について、質問があります。
>>
>> このデータを変換するためのツールとしては、
>> ・基盤地図情報数値標高（GML）対応テスト版GDAL
>> http://lists.osgeo.org/pipermail/osgeojapan-discuss/2014-January/002044.html
>>
>> ・基盤地図情報 標高DEMデータ変換ツール
>> http://www.ecoris.co.jp/contents/demtool.html
>>
>> ・fgddemImporter (QGISプラグイン)
>> http://space.geocities.jp/bischofia_vb/qgis-plugins/fgddemImporter/
>>
>> などなどあります。
>> 先日、このツールを使って10m DEMを変換し、GRASS7
>> のr.geomorphonで分析したところ、図郭の境の部分に斜面があるという結果が出てしまいました。
>>
>> 気になって確認したところ、国土地理院から公開されている基盤地図情報変換ツールを使って変換したshpファイルの標高と、DEMの標高が本当に微妙に（0.0001mmのオーダーで）ずれていることが分かりました。
>> 以下、参考のスクリーンショットです。
>> https://www.dropbox.com/s/77g1ega1add4r9l/img.png?dl=0
>>
>> おそらく、変換に使っているgdalが原因だと思うのですが、どなたか同じような問題にあたった方はいらっしゃいますでしょうか？実質的な解決としては、ラスタを100倍して整数値にすればいいと思いますが、少々気になりました。
>> 何か情報がありましたら、提供頂ければ幸いです。
>>
>> よろしくお願いします。
>>
>> --
>> 岩崎 亘典
>> _______________________________________________
>> OSGeoJapan-discuss mailing list
>> OSGeoJapan-discuss ＠ lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss



-- 
岩崎 亘典