[OSGeoJapan-discuss] 【教えて下さい】Qgisのラスタ範囲と計算機について

2014年 3月 10日 (月) 05:29:54 PDT

朝日様、岩崎様

島田です。お世話になります。
本日、外におりましたので、お返事が遅くなりましたことお詫びします。

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

アドバイスを参考に、もう一度こちらで作業してみます。
それでも解決されない場合は、データをお送りしたいと思います。

まずは頑張ってみます。


> また、不等号がmo ＠ 1>01 となっていますが、mo ＠ 1>1 とした場合はいかがでしょうか

すみません、書き間違いでした。いずれも0より大きい、としました。

mo ＠ 1>0 AND higata ＠ 1>0 
(mo ＠ 1>0) AND (higata ＠ 1>0)

>計算させようとしているラスタの形式ですが、何になりますでしょうか。grasでしょうか？

Qgisでベクトル→ラスタ変換したものですから、デフォルトでgeotiff・・・だと思うので
すが・・・これも確認します。

島田＠調布市
----
Nobusuke Iwasaki <wata909 ＠ gmail.com>さん：
> 島田さん、朝日さん
> 
> 岩崎です。
> 範囲などは朝日さんが答えてくれた通りだと思います。
> #つか、ソースは見てもわからないので、参考になりました。
> 
> 計算が上手くいかない点ですが、計算させようとしているラスタの形式ですが、何にな
りますでしょうか。grasでしょうか？
> 
> また、不等号がmo ＠ 1>01 となっていますが、mo ＠ 1>1 とした場合はいかがでしょうか
> 
> 2014/03/10 13:49 "ASAHI Kosuke" <asahi ＠ hcc.co.jp>:
> 
> 島田様
> 
> 朝日です。
> 
> ぱっとソース見てみてQGISの挙動だけ。
> 
> > ●ベクタ→ラスタ変換をする時、範囲設定はどうなっているんでしょうか？
> これは変換するベクタレイヤの範囲となるようです。
> 
> >> ●ラスタ計算機について
> ・出力ラスタの範囲、ピクセル数
>  ダイアログ表示時にデフォルトで表示される、範囲、ピクセル数は、最初にプ
> ロジェクトに読み込まれているラスタレイヤのものが表示されるようです。
>  変更しない場合、この範囲、ピクセル数が出力ラスタに使用されます。
> 
> ・出力ラスタと、入力ラスタの範囲、ピクセル数が違う場合
> 　計算前に入力ラスタを出力ラスタに合わせる処理が入っています。
> 　範囲外の部分にはnodataが入りますし、出力ラスタに合わせたピクセルの値は
> NearestNeighbourで求めるようです。
> 
> 島田様の処理がなぜうまく行ってないのかはわかりませんね。。。
> 差し支えなければ、データを別途送付頂ければ見てみますが。
> 
> 以上です。
> 
> (2014/03/09 13:38), SGW00412 ＠ nifty.com wrote:
> > 皆様
> >
> > Qgis初心者の島田です。
> >
> > Qgisのラスタ計算機について、二点ほど分からないことがあります。
> >
> > ●ベクタ→ラスタ変換をする時、範囲設定はどうなっているんでしょうか？
> > いくつかのレイヤをラスタ群に変換して、ラスタ計算機を使っていろいろやりたいよ
うな
> > 場合、範囲を統一すると思います。
> >
> > grass　だと　g.region vect/rast=hoge で設定するのですが、Qgisnのベクタ→ラス
タ変
> > 換には、それがないように思います。。。。
> >
> > gdal解説をみると、「指定しない限り、ベクトルレイヤ群(layer[s])の範囲だ」とい
うニ
> > ュアンスで書かれています・・・つまり、変換したいレイヤ群をすべて表示しておけ
ば、
> > その範囲で統一される、ということなのでしょうか？
> >
> > あるいは、単に変換しようとするベクタレイヤの「範囲」なのでしょうか？ただ、こ
れだ
> > とレイヤ毎に範囲が異なってしまうので、ラスタ計算できないような気がします。。。
> >
> > -te xmin ymin xmax ymax :
> > (GDAL >= 1.8.0) set georeferenced extents. The values must be expressed in
> > georeferenced units. If not specified, the extent of the output file will be 
the
> > extent of the vector layers.
> >
> > ●ラスタ計算機について
> > 前の質問で用いた「藻場」と「干潟」で、約500m（0.005度）のラスタに変換し、両者
が共
> > 通して出現するラスタを求めてみました(不在0,存在1のラスタ）。
> >
> > 両者の範囲は異なっていました。
> >
> > ncols        145
> > nrows        132
> > xllcorner    129.497745118728
> > yllcorner    33.009945147708
> > cellsize     0.005000000000
> > ・・・・以下、略
> >
> > ncols        185
> > nrows        146
> > xllcorner    129.497751548845
> > yllcorner    32.838408821121
> > cellsize     0.005000000000
> > ・・・・以下、略
> >
> > ただ、範囲の異なる二つのラスタを用いて「ラスタ計算機」で処理しても、計算して
くれ
> > るようです。しかし・・・
> >
> > mo ＠ 1>01 AND higata ＠ 1>0 計算しますが、不正解（これは何故なのでしょう？）
> > (mo ＠ 1>01) AND (higata ＠ 1>0) 計算しますが、不正解（これも何故なのでしょう？）
> >
> > そこで・・・
> > mo ＠ 1 + higata ＠ 1
> > としてから
> > sum_mo_higata ＠ 1=2
> > とすると正解が得られました。
> >
> > なかなか難しいですね。。。
> >
> > ---
> > 島田＠調布市
> > _______________________________________________
> > OSGeoJapan-discuss mailing list
> > OSGeoJapan-discuss ＠ lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
> >
> 
> 
> --
> 北海道地図株式会社 研究開発課
> 朝日　孝輔
> mailto:asahi ＠ hcc.co.jp
> Tel : 0166-62-4711
> 
> _______________________________________________
> OSGeoJapan-discuss mailing list
> OSGeoJapan-discuss ＠ lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss