tiffとtfwの投影変換について
Up to GISフォーラム
初歩的質問かと思いますが、ご教授下さい。
Qgisで地図画像の投影変換を行いたいと思っております。
kasimirでとってきた地図画像と、その座標(左上と右上)から
http://www.yaskey.cside.tv/mapserver/note/geotiff.html
を用いて、tfwファイルを作成しました(緯度経度座標)。
地図画像とtfwファイルの名前を同一にしてQgisに読み込むと
もちろん緯度経度座標で表示されるのですが、これをUTMに投影変換したいのですが
何か方法はありますでしょうか?
(ベクトルデータなら、CRSを変更しつつ、保存すれば次に開くときに変わっている。。。ように思います)
たとえば、gdalでtfwを含めて投影変換できるのでしょうか。。。
こんにちは、質問有難う御座います。
まずは、GDALのtfw対応についていなのですが、GDALのホームページを読んだところ、tifファイルに空間参照系に関するのヘッダーがついていない場合はGDALはtfwファイルを自動的に検索するようになっています。
If no georeferencing information is available in the TIFF file itself, GDAL will also check for, and use an ESRI world file with the extention .tfw, .tifw/.tiffw or .wld, as well as a MapInfo .tab file (only control points used, Coordsys ignored). (http://www.gdal.org/frmt_gtiff.html より引用)
ということでGDALで読み込んで、gdalwarpを使えば投影系の変換はできることになります(自分でこのようなケースにぶつかったことが無いのであくまで理論的にですが)。
もう一つの方法として、tifとtfwファイルから、GeoTiffファイルを作ってしまうという方法もあります。作り方としては、
- libtiffライブラリをダウンロード (http://www.remotesensing.org/libtiff/) ウィンドウズだとバイナリのインストーラがあります。
- libtiffライブラリーをインストール
- ライブラリー内のgeotifcpを使って、tifファイルにヘッダーとしてtwfの内容を入れ込む。
という順番になります。geotifcpの文法は、以下のように、-eというオプションを付け、その後スペースで区切りながら、入力tfwファイル、入力tifファイル、出力GeoTiffファイル、と言う順番になります。(http://www.omg.unb.ca/~jonnyb/processing/geotiff_tifw_format.html を参照)
geotifcp -e test.tfw test.tif test_geotiff.tif
公式の文法は、
geotifcp [options] [-e esri_worldfile] [-g metadata] input output
となっています(http://www.remotesensing.org/geotiff/geotifcp.html)。
いったんGeoTiffファイルを作ってしまえばgdalwarpで投影系の変換ができます。
まず最初の方法を試してみて、うまく行かなかったら2番目の方法を試してみてはどうでしょうか?うまく行かなかったらまた連絡下さい。
今木
今木様
>..tifとtfwファイルから、GeoTiffファイルを作ってしまう
ありがとうございます。なるほど、そうか。勉強になります。
さっそく試してみたいと思います。どうもありがとうございました。
