[OSGeoJapan-discuss] ＧＲＡＳＳについて

2022年 1月 17日 (月) 17:12:00 PST

いわささん、みなさん

いわさきです。

> ロケーションとマップセットについての考え方について教えて頂ければ助かります。
>
>　ロケーション：shikoku(CRSは6672)
>　マップセット：komi(現場名)
>
> 自分の場合、作業対象は四国内＝CRS６６７２なのでロケーションは１つとしマップセットを現場ごとに増やしていくような
感じでやっていくのがいいのでしょうか？

実は、その辺のやり方は人それぞれだったりします。
かなり古い資料ですが、私の作った資料が以下になるのですが、私の場合は以下のように使い分けています。
https://www.slideshare.net/wata909/20gisgrass/11

・LOCATION：座標系、分析対象範囲、ラスタのデフォルト解像度を定義。異なる座標系は混在できない
・MAPSET：実際の分析や処理の対象や、データの型によってmapsetを作り分ける。たとえば、DEMの処理結果、ベクタデータ等々
　ただし、共通して使うものは、PARMANENTに格納する

ということで、もし自分がやればですが、以下のようになるかと思います。
・LOCATION：komi(現場名)
・MAPSET：PERMANENT（もとになるオルソ画像、DSM）
　　　　　　mapset1（例えば、DMSを処理したデータ）
　　　　　　mapset2（例えば、ベクタデータ）
r.fillnullsをするだけであれば、個別のマップセットは作らず、PERMANENTだけでもいいと思います。

また作業に詰まっている部分の「３．ＧＲＡＳＳだと計算対象範囲指定方法」ですが、locationを作るときに、処理対象とするDSMを読みこんでつくると、自動的にそのDSMの範囲、解像度が処理対象として設定されます。
以下のスライドや、GRASSビギナーズマニュアルの「ロケーションとマップセットの作成」でいうと、「新しいロケーションの作成方法を選択」のところで、「地理参照ファイルを選択」にチェックを入れて、処理対象となるDSMを選択すると、DSMにあわせてロケーションが作成されると思います。
https://www.slideshare.net/wata909/20gisgrass/25
https://gis-oer.github.io/gitbook/book/materials/GRASS/GRASS.html
https://gis-oer.github.io/gitbook/book/materials/GRASS/pic/grasspic_5.png

ここで、GRASSのメニューの「設定」→「Computatinal region」→「Show current region
[g.region -p]」を選んでいただくと、計算範囲や解像度の数値が表示されます。
以上の手順、ざっくりキャプチャしたものを共有します。
https://www.dropbox.com/s/ohbdtiehy4fpedk/GRASS_Mapset%E4%BD%9C%E6%88%90.docx?dl=0

QGISから扱うのは、ちょっと時間がないので、また機会があれば…

2022年1月17日(月) 15:18 <y-iwasa ＠ limecon.co.jp>:
>
> 岩崎様
>
> やっと時間ができてGRASSを今日試したのですが午前中一杯かけてどうにかできるように
> なったのですがロケーションとマップセットについての考え方について教えて頂ければ助かります。
>
> ロケーション：shikoku(CRSは6672)
> マップセット：komi(現場名)
>
> 自分の場合、作業対象は四国内＝CRS６６７２なのでロケーションは１つとしマップセットを現場ごとに増やしていくような
> 感じでやっていくのがいいのでしょうか？
>
> なおＧＲＡＳＳでの処理計算ですが以下のようか通りでした（QGISは３．２２です）
> ①プロセッシングツールで"r.fillnulls"を実行：エラー
> ②プラグインのＧＲＡＳＳで実行：入力ファイル（処理対象のファイル）が候補に出ないも一度GRASSの単独プログラム側でラスタファイルを対象マップセットで
> 表示させるとその後はQGIS側で変換対象で表示されます、また対象範囲はＱＧＩＳのプラグイン側で設定しないとできませんでした。
> ③ＧＲＡＳＳ単独：範囲指定方法がわからず、そのまま実行するとエラーで終了、ただし事前にＱＧＩＳのプラグインで範囲が指定されいると処理計算は問題なく完了
>
> なので現在の状況では②のＱＧＩＳプラグインだとどうにかＧＲＡＳＳでの計算ができた状況です。ただし正しい手順での作業かは自信がありません。
> 何となく何も設定をしていない状態で①（プロセッシングツール）だとうまく動かないかと、結局②か③になるのですが②プラグインだと入力対象ファイルがＱＧＩＳ上では登録できない
> ③ＧＲＡＳＳだと計算対象範囲指定方法がわからず結局ＱＧＩＳとＧＲＡＳＳを行ったり来たりという状況です。
>
> 何かやり方が間違っていただアドバイスを頂ければ助かります。
>
> いわさ＠松山
>
> -----Original Message-----
> From: IWASAKI Nobusuke <niwasaki ＠ naro.affrc.go.jp>
> Sent: Thursday, January 6, 2022 8:12 PM
> To: y-iwasa ＠ limecon.co.jp; osgeojapan-discuss ＠ lists.osgeo.org
> Subject: Re: [OSGeoJapan-discuss] ＤＥＭデータ加工（水面部の処理修正）
>
> いわささま
>
> GRASSは単独でも使う事が出来ますし、QGISのprocessingから使うこともできます。
> 個人的にはGRASS単独で使う方が安定しているからなと思うのですが、数枚のデータを処理するだけであれば、QGISのprocesingから使った方が簡単でいいと思います。
>
> 以下、参考のサイトです
> GRASSビギナーズマニュアル
> https://gis-oer.github.io/gitbook/book/materials/GRASS/GRASS.html
>
> チュートリアル: GRASSツールを用いて河川と集水域の境界線を求める
> （使っている処理は違いますが、QGISからやる場合、同じ手順でできると思います）
> https://courses.gisopencourseware.org/mod/book/tool/print/index.php?id=213
>
>
> On 2022/01/06 19:57, y-iwasa ＠ limecon.co.jp wrote:
> > 補足です。
> >
> > GRASSは別のプログラムなんですね。
> > まずはGRASSの設定をしないといけないようなのでまずはGRASSの設定＆使い方を
> > 調べてから週末にこの機能について確認させて頂きます。
> >
> > いわさ＠松山
> >
> > -----Original Message-----
> > From: y-iwasa ＠ limecon.co.jp <y-iwasa ＠ limecon.co.jp>
> > Sent: Thursday, January 6, 2022 7:52 PM
> > To: 'IWASAKI Nobusuke' <niwasaki ＠ naro.affrc.go.jp>;
> > 'osgeojapan-discuss ＠ lists.osgeo.org'
> > <osgeojapan-discuss ＠ lists.osgeo.org>
> > Subject: RE: [OSGeoJapan-discuss] ＤＥＭデータ加工（水面部の処理修正）
> >
> > 岩崎様
> >
> > 追加コメントありがとうございます。
> > 早速確認させて頂きます。
> >
> > ただしGRASSのコマンド自体使うのが初めてなのでちょっと不安です。
> >
> > いわさ＠松山
> >
> > -----Original Message-----
> > From: OSGeoJapan-discuss <osgeojapan-discuss-bounces ＠ lists.osgeo.org>
> > On Behalf Of IWASAKI Nobusuke
> > Sent: Thursday, January 6, 2022 11:11 AM
> > To: osgeojapan-discuss ＠ lists.osgeo.org
> > Subject: Re: [OSGeoJapan-discuss] ＤＥＭデータ加工（水面部の処理修正）
> >
> > いわささま、田鎖さま
> >
> > 岩崎です。
> > GRASSのコマンドで、r.fillnullsというのがあります。
> > https://grass.osgeo.org/grass78/manuals/r.fillnulls.html
> >
> > みんなの自動翻訳で翻訳させると、
> > r.fillnullsは、入力ラスターマップのNULLピクセル(データ領域なし)を塗りつぶし、塗りつぶされたデータを新しい出力ラスターマップに格納します。塗りつぶし領域は、v.surf.rst正則化スプライン補間(方法=rst)またはr.resamp.bsplineキュービックまたは線形スプライン補間(Tykhonov正則化)を使用して、非データ領域境界バッファから補間されます。
> >
> > という感じです。
> > GRASSの機能なので走らせるまでがちょっと手間かもしれませんが、参考までに。
> >
> >
> > On 2022/01/05 6:31, y-iwasa ＠ limecon.co.jp wrote:
> >> 田鎖様
> >>
> >> 実はコメントのリンクは自分も今一度調べた時に出てきて確認させて頂きました。
> >>
> >> 前半部分の補完処理はうまくいったりいかなかったりで結局挫折してしまいました。
> >>
> >> もしかすると対象が比較的川幅が細かったかもしれません。
> >>
> >> 読んでみて意外と１０ｍＤＥＭならば粗いですが欠測がないのでドローンでの地形把握ならば
> >>
> >> これが一番簡単なのかもしれません。
> >>
> >> ５ｍと１０ｍの合成も試したのですが境目の部分がまたー９９９９もしくは０のままでした。
> >>
> >> いわさ＠松山
> >>
> >> *From:*OSGeoJapan-discuss
> >> <osgeojapan-discuss-bounces ＠ lists.osgeo.org>
> >> *On Behalf Of *Junta TAGUSARI
> >> *Sent:* Tuesday, January 4, 2022 3:02 PM
> >> *To:* osgeojapan-discuss ＠ lists.osgeo.org
> >> *Subject:* Re: [OSGeoJapan-discuss] ＤＥＭデータ加工（水面部の処理修正）
> >>
> >> いわさ様
> >>
> >> はじめまして，田鎖と申します。
> >> 地理情報が専門というわけではないのですが，業務で必要に迫られてGISを触っている者です。
> >>
> >> 水面データはー９９９９で入っていますが，これは「欠測」として扱い，周囲の値を使って補間するのが良いのではないかと思います。
> >> 過去によく似たやりとりが行われているのを見つけましたのでご参照ください。
> >>
> >> https://groups.google.com/g/qgisshitumon01/c/mTnnGJGk0-U?pli=1
> >> <https://groups.google.com/g/qgisshitumon01/c/mTnnGJGk0-U?pli=1>
> >>
> >> ちなみに，このやり取りの中でも触れられていますが，海面の高さまで勝手に補間されるので注意が必要です。
> >> たとえば，海が陸地で挟まれている湾のような地形で，かつ海に接している陸地の標高がそこそこ大きいと，湾内の標高が数mとかで補間されてしまいます。
> >> このような場合，目的によっては，水面を0mで設定した方が良いかもしれません。
> >>
> >> ご参考になれば幸いです。
> >>
> >> 〒060-8628　札幌市北区北13条西8丁目
> >>
> >> 北海道大学大学院工学研究院　環境工学部門
> >>
> >> 地域環境研究室
> >>
> >> 田鎖順太　Junta Tagusari
> >>
> >> Tel&Fax: 011-706-6872
> >>
> >> Email: j.tagusari ＠ eng.hokudai.ac.jp
> >> <mailto:j.tagusari ＠ eng.hokudai.ac.jp>
> >>
> >> On 2022/01/03 9:00, y-iwasa ＠ limecon.co.jp <mailto:y-iwasa ＠ limecon.co.jp> wrote:
> >>
> >>      いわさ＠松山です。
> >>
> >>      本来の本年最初の投稿をさせて頂きます。
> >>
> >>      仕事はじめは１／５なのでいきなり現場観測なので．．．
> >>
> >>      実は最近ドローンでの撮影する頻度が増えてきて事前に可能ならばＤＥＭデータを地理院からダウンロードして
> >>
> >>      確認しているのですが川などの水面部分をできれば川辺の高さに表示できるとイメージがわきやすいのですが。
> >>
> >>      いまだと単純に水面部が“-9999”にしているので渓谷に見えてしまいます．．．
> >>
> >>      基盤地図情報の変換ソフトでは水面は０かー９９９９に指定できるのですが例えば後者（－９９９９）に指定すれば
> >>
> >>      ラスタ計算機でできるのかなと思っていますがその方法がわからず。
> >>
> >>      自分で思いついたのはその場所の標高あたりに適当に設定してｉｆ文で“－９９９９”の場合“標高試算値”に置き換えですが
> >>
> >>      これだと対象画像が同じ高さに修正されています。できれば水面と陸地の境の標高値にスムーズにかつ自動的に合わせてくれると
> >>
> >>      いいのですがこのような機能はありませんでしょうか？
> >>
> >>      もしかしたらＦＵＳＩＯＮなどであるかもしれないのですがご存知の方がいらっしゃいましたら教えて頂けると助かります。
> >>
> >>      いわさ＠松山
> >>
> >>
> >>
> >>      _______________________________________________
> >>
> >>      OSGeoJapan-discuss mailing list
> >>
> >>      OSGeoJapan-discuss ＠ lists.osgeo.org
> >> <mailto:OSGeoJapan-discuss ＠ lists.osgeo.org>
> >>
> >>      https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
> >> <https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss>
> >>
> >>
> >> _______________________________________________
> >> OSGeoJapan-discuss mailing list
> >> OSGeoJapan-discuss ＠ lists.osgeo.org
> >> https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
> >
> > --
> > ☆通常勤務中です☆
> > #staypositive
> >
> > ---
> > Nobusuke IWASAKI
> > Representative, OSGeo Japan Chapter
> >
> > Institute for Agro-Environmental Sciences, NARO 3-1-1, Kannondai,
> > Tsukuba-shi, Ibaraki-ken 305-8604, JAPAN Tel / Fax +81298388227
> > _______________________________________________
> > OSGeoJapan-discuss mailing list
> > OSGeoJapan-discuss ＠ lists.osgeo.org
> > https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
> --
> ☆通常勤務中です☆
> #staypositive
>
> ---
> Nobusuke IWASAKI
> Representative, OSGeo Japan Chapter
>
> Institute for Agro-Environmental Sciences, NARO 3-1-1, Kannondai, Tsukuba-shi, Ibaraki-ken 305-8604, JAPAN Tel / Fax +81298388227
> _______________________________________________
> OSGeoJapan-discuss mailing list
> OSGeoJapan-discuss ＠ lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss


-- 
岩崎 亘典