傾斜地の面積の求め方

[primos colophon]

先週からQGISを使い始めた超初心者です。

以下の手順で、ある山の標高300m以上の面積を求めようとしていますができません。

初歩的な質問で申し訳ないのですが、ご教示いただけないでしょうか。

今は、3.のラスタデータをポリゴンに変換することができず、止まっています。

大まかな流れ（予定）

1.DEM作成

2.標高300ｍ以上の地域を抽出（ラスタデータ）

3.このラスタデータをポリゴン（ベクターデータ）に変換

4.DEMとこのポリゴンの重なる部分の面積を計算

具体的にやった作業は以下のとおりです。

1.

国土地理院の基盤地図情報ダウンロードページから数値標高モデル（DEM）をDL

今回は隠岐諸島のデータをDL、解凍。

これらの.xmiファイルを1つのフォルダに入れる。

 

エコリスの「基盤地図情報標高DEM変換ツール」（https://www.ecoris.co.jp/contents/demtool.html）を使用して、

すべての.xmlファイルを結合したGeoTIFF形式に変換する。

ファイル読み込みのときに、LatLog（緯度経度）0を選ばないようにする。今回は2を選んだ。

→merge.tif（投影法は緯度経度）とmerge_shade.tif（陰影起伏図）の2つのファイルが出力される。

 

これらのデータが入ったフォルダをC:\user\QGISに入れる。

ファイル名やフォルダ名は全角を使用しない。

2.

QGIS Desktop 3.10.1 with GRASS 7.8.1を起動。

merge.tifとmerge_shade.tifを読み込む。

 

300m以上の場所を抽出

以下を参照して、ラスタ計算機で300m以上を抽出

https://gis-oer.github.io/gitbook/book/materials/15/15.html#%E3%83%A9%E3%82%B9%E3%82%BF%E6%BC%94%E7%AE%97

→300ｍ.tifが出力。merge.tifと同じフォルダに入れて、QGISに読み込み。

 

3.

300ｍ以上のポリゴン（ベクターデータ）作成のため、

ラスタ―変換ーラスタのベクタ化

エラーになる

エラーメッセージは↓

メインウィンドウのログメッセージパネルのチェックの方法もわからないので、何が悪いのか検討が付きません。

[takanoko]

私がやった限りでは、ほぼそのとおりでベクタ化できましたね。

300ｍ.tifはベクタ化できるような形になってますか。

ビュー→パネル→メッセージパネル　（もしくはウインドウのステータスバー（右下）のメッセージアイコン）を見れば何かわかるかも

[水谷貴行]

こんにちは水谷です。

どうやらQGISの現時点での最新バージョン（3.4.14-1とか3.10.1-1）では「ポリゴン化」を実行するとエラーになるバグがあるようです。

（その他、pythonを使用するコマンドはダメっぽいです。）

詳細は以下のリンクをご覧ください。（英語ですがgoogle翻訳を通すとなんとなく分かります。）

https://github.com/qgis/QGIS/issues/33386

PYTHONPATHを設定するなどの対応方法も書いてありますが、前のバージョンでは問題ないとあるので、

とりあえず3.4.13とか3.10.0をインストールしてみるのが簡単ではないかと思います。

（次のマイナーアップデートでは、おそらく修正されると思いますが。）

[adachi]

エラーを見ると「gdal_polygonizeが見つからない」みたいなことを言っているようなので

オプション→システム

下のほうの「現在の環境変数」の「GDAL_DATA」のところで、ちゃんとgdalのパスが入っているか確認してみてはどうでしょうか。

なお求めたい「面積」というのは「水平面への投影面積」の意図で良いのでしょうか？

傾斜を考慮した表面面積もあるかと思いますが。

[primos colophon]

皆様

早速のご教示ありがとうございます。

1つずつ試してみたいと思います。

求めたいのは、後者です。前者も求めて比較してみたいとも思っています。

わかりにくくてすみません。

[primos colophon]

takanoko様

作業の流れとしては間違ってないことが確認できて安心しました。

水谷さんがご指摘くださったように、QGISのバージョンによってうまく行かないこともあるとのことでしたので、3.14.4をインストールしてやってみました。

ウインドウのステータスバーのPython警告（以下に全部をペースト）には以下のログが出ています。が、私にはなんのことかさっぱりわかりません。

最初の行に

>DeprecationWarning: invalid escape sequence \o

とありますので、このログが読み解ければ解決するかもとは思うのですが。

もしよろしければ解決方法をご教示いただけましたらありがたいです。

============

2020-01-16T18:59:40     WARNING    warning:C:\PROGRA~1\QGIS3~1.4\apps\Python37\lib\site-packages\requests\status_codes.py:18: DeprecationWarning: invalid escape sequence \o
              """
            
             traceback: File "", line 1, in
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 335, in startPlugin
              plugins[packageName] = package.classFactory(iface)
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python/plugins\MetaSearch\__init__.py", line 29, in classFactory
              from MetaSearch.plugin import MetaSearchPlugin
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 685, in _import
              mod = _builtin_import(name, globals, locals, fromlist, level)
              File "", line 983, in _find_and_load
              File "", line 967, in _find_and_load_unlocked
              File "", line 677, in _load_unlocked
              File "", line 728, in exec_module
              File "", line 219, in _call_with_frames_removed
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python/plugins\MetaSearch\plugin.py", line 34, in
              from MetaSearch.dialogs.maindialog import MetaSearchDialog
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 685, in _import
              mod = _builtin_import(name, globals, locals, fromlist, level)
              File "", line 983, in _find_and_load
              File "", line 967, in _find_and_load_unlocked
              File "", line 677, in _load_unlocked
              File "", line 728, in exec_module
              File "", line 219, in _call_with_frames_removed
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python/plugins\MetaSearch\dialogs\maindialog.py", line 50, in
              from owslib.csw import CatalogueServiceWeb # spellok
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 685, in _import
              mod = _builtin_import(name, globals, locals, fromlist, level)
              File "", line 983, in _find_and_load
              File "", line 967, in _find_and_load_unlocked
              File "", line 677, in _load_unlocked
              File "", line 728, in exec_module
              File "", line 219, in _call_with_frames_removed
              File "C:\PROGRA~1\QGIS3~1.4\apps\Python37\lib\site-packages\owslib\csw.py", line 27, in
              from owslib.util import OrderedDict
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 685, in _import
              mod = _builtin_import(name, globals, locals, fromlist, level)
              File "", line 983, in _find_and_load
              File "", line 967, in _find_and_load_unlocked
              File "", line 677, in _load_unlocked
              File "", line 728, in exec_module
              File "", line 219, in _call_with_frames_removed
              File "C:\PROGRA~1\QGIS3~1.4\apps\Python37\lib\site-packages\owslib\util.py", line 32, in
              import requests
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 685, in _import
              mod = _builtin_import(name, globals, locals, fromlist, level)
              File "", line 983, in _find_and_load
              File "", line 967, in _find_and_load_unlocked
              File "", line 677, in _load_unlocked
              File "", line 728, in exec_module
              File "", line 219, in _call_with_frames_removed
              File "C:\PROGRA~1\QGIS3~1.4\apps\Python37\lib\site-packages\requests\__init__.py", line 114, in
              from .models import Request, Response, PreparedRequest
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 685, in _import
              mod = _builtin_import(name, globals, locals, fromlist, level)
              File "", line 983, in _find_and_load
              File "", line 967, in _find_and_load_unlocked
              File "", line 677, in _load_unlocked
              File "", line 728, in exec_module
              File "", line 219, in _call_with_frames_removed
              File "C:\PROGRA~1\QGIS3~1.4\apps\Python37\lib\site-packages\requests\models.py", line 43, in
              from .status_codes import codes
              File "C:/PROGRA~1/QGIS3~1.4/apps/qgis-ltr/./python\qgis\utils.py", line 685, in _import
              mod = _builtin_import(name, globals, locals, fromlist, level)
              File "", line 983, in _find_and_load
              File "", line 967, in _find_and_load_unlocked
              File "", line 677, in _load_unlocked
              File "", line 724, in exec_module
              File "", line 860, in get_code
              File "", line 791, in source_to_code
              File "", line 219, in _call_with_frames_removed
            

=============

[primos colophon]

水谷様

ご教示ありがとうございます。

バージョンでうまく行かなかったりするとはまったく考えてもいませんでした。

ですが、長期安定版の3.14.4をインストールしてやり直しましたが、同じ結果でした。

[primos colophon]

adachi様

パスを調べてみました。

写真の通り、

C:\PROGRA~1\QGIS3~1.4\apps\gdal2\share\gdal

と入力されており、パスは通っているようです。

[水谷貴行]

こんにちは

マイナーバージョン（例えば3.4.xx のxx部分）の違いでバグの有無が異なります。

今回のバグは、最新版(3.10)および長期リリース版(3.4)ともに、現在公開中のバージョン(3.10.1-1および3.4.14-1)にバグがあるようなので、

それ以外を以下からダウンロードしてインストールしてみてください。(おそらくwindows版限定のバグだと思います。)

https://qgis.org/downloads/

私が試したところ、3.4.11-1では問題ありませんでした。

64ビットパソコンであれば以下のファイル名がそれに該当します。

QGIS-OSGeo4W-3.4.11-1-Setup-x86_64.exe

（公開中のものが早く修正されると良いのですが...）

[水谷貴行]

すみません。本題ではないですが「マイナーバージョン」と言う呼び方が間違ってたので訂正です。

リビジョンとかビルドと呼ぶのが正しいですね。

参考：

https://www.change-your-tomorrow.com/software-version-general-rules/

[primos colophon]

水谷様

詳しく教えていただき、誠にありがとうございます。

教えていただいたマイナーバージョンのものをインストールしたら、かんたんにできてしまいました。

まさか長期リリース版でもこういう不具合があるとは思っていませんでしたので、とても助かりました。以後はいくつかのマイナーバージョンでも同じ操作をやってみます。

本題の面積の求め方についても不具合？がありますので、別に質問を立てようと思います。