[OSGeoJapan-discuss] ラスタ作成時のスナップについて

[Mayumi Kubo]

赤木さま、岩崎さま

ラスタ上書きの情報ありがとうございます。
レイヤの重ね合わせ作業が多いので役立ちそうです。
上書きの方法でも試してみたいと思います。

取り急ぎ、御礼まで

久保まゆみ


2014年5月13日 23:41 Nobusuke Iwasaki <wata909 ＠ gmail.com>:

> 赤木さん、久保さん
>
> 岩崎です。
> 情報ありがとうございます。
>
> gdal_rasterlizeで既存のラスタ似合わせて出力するのが可能なのは知っていたので、おそらく他のコマンドでも可能ではないかと思って試してみたのですが、上手くいきませんでした。
> コメントを参考に、もう一度試してみようと思います。
>
> ありがとうございました。
>
>
> 2014年5月13日 15:15 Minoru Akagi <akaginch ＠ gmail.com>:
> > 久保さま、岩崎さま
> >
> > 型枠となるラスタと範囲・解像度が等しいラスタを作成して、それにgdalwarpで値を上書きする方法もあります。
> >
> > # セルの値がすべて0のGeoTIFFファイルを作成
> > gdal_translate -scale -9999 9999 0 0 -a_nodata none base.tif output.tif
> >
> > # 地理情報を持つ画像を用いる場合はバンド1だけに
> > # gdal_translate -ot Float32 -b 1 -scale 0 255 0 0 -a_nodata none
> > image.tif output.tif
> >
> > # 作成したGeoTIFFファイルに上書きされるようにgdalwarpを実行
> (不要な引数を指定すると新規作成モードになり実行されない場合がある。例えば-of)
> > gdalwarp -r binear D:\SRTM3\SRTM3.vrt output.tif
> >
> > ご参考まで
> >
> > 赤木 実
> >
> > 2014年5月13日 14:26 Mayumi Kubo <mayumitt ＠ gmail.com>:
> >> 岩崎さま
> >>
> >> どうもありがとうございます。
> >> -tapで原点を通るラスタと重なるように無事作成できました。
> >>
> >> -tapがあやしいと思い、何度も説明を読んだのですが読解力足りませんでした。
> >> -trの値に配列する＝原点を通るラスタと重なると言うことですね。
> >> -crop_to_cutlineも活用したいと思います。
> >> 大変助かりました。どうもありがとうございました。
> >>
> >> 久保まゆみ
> >>
> >>
> >>
> >> 2014年5月13日 10:41 Nobusuke IWASAKI <niwasaki ＠ affrc.go.jp>:
> >>>
> >>> 久保さん
> >>>
> >>> いわさきです。
> >>> gdalを使う場合、オプションに"-tap"というものがありました。
> >>> gdalwarpを使うときに、"-tap -tr 10 10"としてあげると、原点（0、0）を通っ
> >>> て、解像度10mのラスタができあがりました。
> >>> 範囲をあわせたい場合であれば、"-crop_to_cutline"オプションを使って、マッ
> >>> チさせる対象となるラスタファイルを指定すれば大丈夫のようです。
> >>>
> >>> たとえば、10mDEMを2451（平面直角９系）に変換する場合、
> >>>
> >>> gdalwarp -t_srs EPSG:2451 -tap  -tr 10 10 -crop_to_cutline
> >>> extent.tif（範囲を指定するtif） source_dem.tif dist_dem.tif
> >>>
> >>> といった感じでしょうか。
> >>> 試してみて下さい。
> >>>
> >>>
> >>> (2014/05/12 15:54), Mayumi Kubo wrote:
> >>> > お世話になります、久保です。
> >>> >
> >>> > 地理座標系（JGD2000）のDEMを平面直角座標系へ変換する際に
> >>> > 特定のラスタとセルがぴったり重なるようにすることは可能でしょうか。
> >>> > （平面直角座標系の原点を通る同じ解像度のラスタと重なるよう作成
> >>> > したいと思っています。）
> >>> >
> >>> > 今までは、gdalwarpで出力ファイルの範囲(-te)と解像度(-tr)を指定して
> >>> > セルが重なるようにしていたのですが、特定のラスタと重なるように
> >>> > 出力ファイルの範囲を確認する作業が大変なので簡単な方法はないか
> >>> > 探しております。
> >>> >
> >>> > ArcGISでは「スナップ」や「コントロールポイントの指定」を利用していま
> >>> > したが、GDALやGRASSでも簡単に指定する方法をご存知でしたらご教授下さい。
> >>> >
> >>> >
> >>> > 久保まゆみ
> >>> > 寒地土木研究所
> >>> >
> >>> 岩崎 亘典
> >>> (独)農業環境技術研究所
> >>> 生態系計測研究領域 主任研究員
> >>> 〒305-8604  茨城県つくば市観音台3-1-3
> >>>
> >>> 【歴史的農業環境閲覧システムを公開中です】
> >>>  http://habs.dc.affrc.go.jp/
> >>> _______________________________________________
> >>> OSGeoJapan-discuss mailing list
> >>> OSGeoJapan-discuss ＠ lists.osgeo.org
> >>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
> >>
> >>
> >>
> >> _______________________________________________
> >> OSGeoJapan-discuss mailing list
> >> OSGeoJapan-discuss ＠ lists.osgeo.org
> >> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
> > _______________________________________________
> > OSGeoJapan-discuss mailing list
> > OSGeoJapan-discuss ＠ lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>
>
>
> --
> 岩崎 亘典
>
-------------- next part --------------
HTMLの添付ファイルを保管しました...
URL: <http://lists.osgeo.org/pipermail/osgeojapan-discuss/attachments/20140514/c50d065b/attachment-0001.html>