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

[ASAHI Kosuke]

島田様

朝日です。

ぱっとソース見てみて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