傾斜区分図が作成できません

[ti131...@gmail.com]

　初めて質問させていただきます。

田舎の百姓のおじさんにもわかるように教えてください。

近年、施設園芸では環境制御技術として日射量が重視されていることからQGISで確認しようと試みたのですが…。

さっそく、表題のように傾斜区分図で詰まってしまいました。

【範囲】

二次メッシュ＿503326区画

【手順】

1).数値標高モデル(5mメッシュ)をDLし、適当なフォルダに解凍

2).「標高DEMデータ変換ツール(㈱エコリス)」によりGeoTiff化

　投影法＝2(平面直角座標系)　,系番号＝4(四国)

　→「merge.tif」を作成。「Plane_5m.tif」に名前を変更。

3).QGISv3.6を起動し、レイヤ->ラスタレイヤの追加->でラスタデータセットに「Plane_5m.tif」を指定。

　→グレーの画像が表示される。

4).再投影(warp)処理[図]

　変換元：JGD2000

　変換先：JGD2000 / Japan Plane …　IV

　再投影の保存先：指定

5).傾斜区分処理

　先の4)の手順でできたレイヤー「再投影したラスタファイル」を選択し、傾斜(slope)処理を実行[図]

　　↓

　結果、出力された「傾斜(slope)」というレイヤーは0,90(度)の2値のみの白黒画像になってしまいます。

どうも手順4)の再投影がおかしい気がしていますが、対処方法がわかりません。

　そもそもですが投影法の使い分けがまだわかっておりません。

傾斜区分のときに直角座標系の投影法をつかう理由は何でしょう？

UTMとかだと曲面のためDEMがm×nの関係になってないから？でしょうか。

　日射量の計算までまだ先が長そうですが、ともかく傾斜区分図の作り方を教えていただきたいです。

よろしくお願いします。

▲

[wata909]

こんにちは

手順を見ていると，

2).で平面直角座標系として作成したDEMを，4).で，JGD2000（地理座標系）→JGD2000 / Japan Plane …　IV（平面直角４）へと，変換しているようにみえます．

エコリスさんのツールを使っている段階で，平面直角座標系になっているので，４）の手順は必要ないのでは？

２）でできたデータで，傾斜を計算してみてはいかがでしょうか？

試しに，FG-GML-5033-26-DEM5Aでやってみた例が以下になります．

参考までに

[ti131...@gmail.com]

※！．回答メールに直接返信してしまいました。

　失礼しました。ただ、そちらにも別の情報をくださったので以下に載せておきます。

> また、質問にあった

>> 傾斜区分のときに直角座標系の投影法をつかう理由は何でしょう？

>については、実はUTMでも問題ありません。地理座標（単位が緯度経度）ではなく、投影座標（単位がm）であれば大丈夫です。

>ただ、平面直角座標系の方が国内であれば歪みが少ない（と思う）ので、そちらを使う場合が多いのかもしれません。

>参考までに

---*---*---*---*---*---*---*---*---*---

解決しました。

「そういうことか！」と納得しました。

日照時間のほうに取り掛かりたいと思います。

ありがとうございました。

[ti131...@gmail.com]

(案の定)日射量の計算で詰まっています

QGISv3.6から、プロセッシングツールボックス中にあるGRASSのr.sunの各種パラメータを設定して実行させるとエラーが表示され、出力されません。

ログをたどると、

警告: データム <Japanese_Geodetic_Datum_2011> はGRASSに認識されていません そしてパラメータが見つかりません

という記載が見えます。

計算に用いる各種レイヤー(標高,傾斜,傾斜方位)のCRSを確認すると、

「標高DEMデータ変換ツール(㈱エコリス)」によりGeoTiff化した大本の標高レイヤーが、

EPSG:6672 - JGD2011 / Japan Plane Rectangular CS IV - 投影された

になっていることが分かったので、これを再投影すると、

EPSG:4612 - JGD2000 - 地理学

に変換できました。

これを元に、傾斜と傾斜方位のレイヤーを解析するわけですが、

傾斜のほうが例によって失敗します(白黒二値で水部と陸地になってしまいます)。

エラーへの対処方法が間違っているのかもしれません。

対処方法をご存知の方がいらっしゃいましたら、教えていただきたいです。

[福岡]

ご覧の皆様へ

検証しようと思って「r.sun」を実行しようとしたところ、私の環境ではうまく動作しませんでした。

Windows10 Pro 64bit

QGIS3.4.8と3.6.3で検証中

まず、「r.sun」を立ち上げようとすると「このアルゴリズムを実行できません　The specified GRASS 7 folder "C:\PROGRA~1\QGIS3~1.4\bin\bin" does not contain a valid set of GRASS 7 modules. Please, go to the Processing settings dialog, and check that the GRASS 7 folder is correctly configured」というメッセージが出ます。

そこで情報を探してみると「grass76」というフォルダ名を「grass-7.6.0」に変更すると良いという情報をみつけました。

インストールされているGRASSバージョンをみると7.6.1なので、「grass-7.6.1」に変更してみました。

これで「r.sun」を呼び出すことはできました。

GRASS not working with QGIS 3.6

https://gis.stackexchange.com/questions/313724/grass-not-working-with-qgis-3-6

次に実行してみると「次のレイヤーは正しく生成されませんでした・・・」というメッセージが出て結果が出力されません。

以前あった「incidence angle raster map」の出力の問題がまだ残っている可能性もあるかとスキップしてみても同様の結果です。

プロセッシングツールのGRASSで何らかの不具合が発生しているのでしょうか？

それとも私の環境だけの問題でしょうか？