[OSGeoJapan-discuss] 電子地形図25000をQGISで利用する方法(その2)

村石眞澄 ataualpa-requi @ ta2.so-net.ne.jp
2013年 9月 21日 (土) 18:03:31 PDT


 いわさきさん、Tetsuya Obaさんにアドバイを頂いたおかげで、課題を解決す 
ることができました。本当にありがとうございました。

 不具合の第一の原因は、Geotiff画像を編集加工する場合に、レイヤのCRSとプ 
ロジェクトのCRSを一致させていなかったためです。
レイヤの「座標参照システム」を「JGD2000/UTMzone54N」に指定すると画像が傾 
いてしまい、矩形で切り出せなくなるので、一 旦「JDG2000」にワープしてか 
ら、クリップしました。一気に「JDG2000/Japan Plane Rectangular CSVIII」に 
ワープしたかたのですが、これだとクリップする値が座標値となるため一度 
「JDG2000」にワープする方法を採りました。二度ワープす るので、望ましくな 
いとは思いましたが、画質は変化がなかったので良しとしました。

 以下に、正しくない部分もあるかもしれませんが、一連の作業の過程をまとめ 
ました。お気づきの点がありましたら、ご教示いただければ幸いで す。

 
国土地理院の電子地形図25000をQuantumGISで利用する

1. オンラインで定形図郭版(TIFF形式)を購入
・ 定形図郭版 tiff形式、300dpi(通常印刷用) @170円
・ 甲府市周辺

2. QGISに「電子地形図25000」を読み込む
a. メニュー:レイヤ>「ラスタレイヤの追加」

b. メニュー:レイヤ」>「レイヤのCRSの指定」
・ レイヤの「座標参照システム」を指定
・ 地理座標系>「JGD2000/UTMzone54N」を指定

3. 地図投影法を変換
a. メニュー:ラスタ>Projection>ワープ(再投影)
・ 出力ファイル名の指定
・ ソースSRS::JDG2000/UTMzone54N
・ ターゲットSRS:JDG2000
・ サンプリングメソッド:キュービック

4. 地図の余白を切り落とす
a. 新規プロジェクトからワープしたファイルを開く
c. メニュー:レイヤ>「レイヤのCRSをプロジェクトに設定する」
d. メニュー:ラスタ>Extraction >クリップ
・ 出力ファイル名の指定
・ 図郭の左上と右下の緯度経度を指定

5. 地図投影法を再度変換
a. 新規プロジェクトからクリップしたファイルを開く
b. メニュー:ラスタ>Projection>ワープ(再投影)
・ 出力ファイル名の指定
・ ソースSRS:: JDG2000
・ ターゲットSRS:JDG2000/Japan Plane Rectangular CSVIII
・ サンプリングメソッド:キュービック

6. 国土地理院基盤地図25000「行政界」と照合確認















(2013/09/20 10:56), Nobusuke Iwasaki wrote:
> 村石様
>
> いわさきです。
> クリップを実施する際の座標値はどうなっているでしょうか?
>
> 2013年9月20日 7:40 村石眞澄 <ataualpa-requi @ ta2.so-net.ne.jp>:
>> いわさき様
>>
>>  ご教示ありがとうございます。
>>
>>  地図画像のレイヤのみで、クリップを試してみましたが、以前と同じ以下の メッセージが出ました。
>>
>>
>>  Computed -srcwin falls outside raster side of 12096×9094
>>
>>  図中の凡例のスケールの距離の計測を試して、確認をしています。
>>
>>   取り急ぎ、お礼と報告まで
>>
>>
>>
>>> 村石さま
>>>
>>> いわさきです。
>>> ちょっと想像の部分が入るのですが、WGS84緯度経度で表示されている状態でクリップを行っていますでしょうか?
>>> もしそうだとすると、GUIからクリップを行うとに入力した座標値(緯度経度)と、元のデータの座標値(UTM54)が異なっている可能性があります。
>>>
>>> まずは、地図画像だけ表示してクリップを行ってはいかがでしょうか。
>>> よろしくお願いします。
>>>
>>>
>>> 2013年9月19日 2:27 村石眞澄 <ataualpa-requi @ ta2.so-net.ne.jp>:
>>>> いわさき様
>>>>
>>>>  早速のご教示ありがとうございます。
>>>>
>>>> 1.メニュー > ラスタ > ワープ を実行してみました。
>>>>
>>>>  ソースSRS::JDG2000/UTMzone54N
>>>>  ターゲット:SRS:WGS84/UTMzone54N
>>>>
>>>>  に変換したところ、処理速度が格段に向上しました。 満足です。
>>>>
>>>>
>>>> 2.しかし、クリッパーは上手く動きません。
>>>>
>>>>  メッセージは、Computed -srcwin falls outside raster side of 12096×9094
>>>>
>>>>
>>>>  一歩前進しましたが、もう一歩進みたいところです。
>>>>
>>>>  とりあえず、現状の報告です!
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> (2013/09/18 9:32), Nobusuke Iwasaki wrote:
>>>>> 村石さま
>>>>>
>>>>> いわさきです。
>>>>> クリップの方は上手く行ったということでしょうか?
>>>>>
>>>>> さて、オーバービューの方ですが、実は使っているプログラムが上の二つで違います。
>>>>>
>>>>>
>>>>> 1.のやり方はgdalというライブラリを使っています。タイムスタンプが変わっていないとの事ですが、同じフォルダにovrという拡張子のファイルが無いでしょうか?
>>>>> これが、オーバービューデータになります。
>>>>>
>>>>> ちなみに、表示が早くならないのは別の原因だと思います。
>>>>> まずは、ラスター→ワープで、座標系UTM54のGeoTiff形式に変換してからオーバービューを作成してみてください。
>>>>>
>>>>>
>>>>>
>>>>> 2.の方はQGISの独自機能を使ったオーバービューになります。私の方でも確認してみたのですが、上記と同様にワープをかけた後の画像であればオーバービューの作成が可能でした。
>>>>>
>>>>> とりいそぎ。
>>>>>
>>>>>
>>>>>
>>>>> 2013年9月18日 7:32 村石眞澄 <ataualpa-requi @ ta2.so-net.ne.jp>:
>>>>>> Tetsuya Oba 様
>>>>>>
>>>>>>  ご教示ありがとうございます。
>>>>>>
>>>>>>  ファイル名は、注意していましたが、フォルダ名は忘れていました。早速に修正しました。
>>>>>>
>>>>>>  その結果、少し変化がありましたが、問題解決はまだです。
>>>>>>
>>>>>> 1.メニュー > ラスタ > Miscellaneous > オーバービューの作成 の実行はできるようになりました。
>>>>>>
>>>>>>  処理ができたとのメッセージは出ましたが、あまり軽くなったように見えますが、ファイルのタイムスタンプは変わっていません。
>>>>>>
>>>>>> 2.メニュー >  レイヤ > プロパティ  からの オーバービューの作成 はエラーメッセージが出て、実行できません。
>>>>>>
>>>>>>  現状の報告です。もう少し試してみますが、どなたかヒントがあれば、ご指導ください。
>>>>>>
>>>>>>  以上
>>>>>>
>>>>>>
>>>>>>
>>>>>> (2013/09/16 17:02), Tetsuya Oba wrote:
>>>>>>
>>>>>> 村石 様
>>>>>>
>>>>>> もしかして、日本語のfolder名が悪さしてるんじゃないでしょうか?
>>>>>>
>>>>>> 英数字で空白無しのfolder名に変えてみてはどうでしょうか?
>>>>>>
>>>>>> 以 上
>>>>>> (2013/09/16 16:40), 村石眞澄 wrote:
>>>>>>
>>>>>> GIS初心者の村石です。
>>>>>>
>>>>>>  QGISで利用する方法でQGISに取り込むことが出来た「電子地形図25000」のgeotiffの画像加工で困ってい ま す。
>>>>>>
>>>>>> 1.「電子地形図25000」の余白を切り取りたいので、GDALTools のクリッパーを使ってみましたが、以下のよう なエ ラー
>>>>>> メッセージが出て処
>>>>>> 理できません。
>>>>>>
>>>>>>   ※オーバービューの作成も同様です。
>>>>>>
>>>>>> ERROR 4: `I:\OPEN\MapData\“dŽq’nŒ` }25000\–³ˆ³ k\533844b.tif' does not
>>>>>> exist
>>>>>> in the file system,
>>>>>> and is not recognised as a supported dataset name.
>>>>>>
>>>>>>
>>>>>> 2.tiff形式に問題があるのかと思い、無圧縮のもので試してみましたが結果は同じでした。
>>>>>>
>>>>>>
>>>>>> 3.ファイル名もエラーメッセージにあるように英数字のみです。
>>>>>>
>>>>>>  手詰まり状態です、どなたか、Geotiffの編集加工について、ご教示ください。
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> OSGeoJapan-discuss mailing list
>>>>>> OSGeoJapan-discuss @ lists.osgeo.org
>>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>>>
>>>>>>
>>>>>>
>>>>>> --
>>>>>> ==================================
>>>>>> 合同会社 TYS
>>>>>> 技術開発部
>>>>>> 大庭哲哉
>>>>>> Tetsuya Oba
>>>>>>
>>>>>> 〒252-0318
>>>>>> 相模原市 南区 上鶴間本町5-32-25
>>>>>>
>>>>>> TEL : 042-812-5920
>>>>>> FAX : 042-812-5920
>>>>>> E-mail : t.oba @ tys-yokohama.co.jp
>>>>>> URL    : http://www.godo-tys.jp/
>>>>>>           : http://blog.godo-tys.jp/
>>>>>>           : http://www.disaster-software.net/
>>>>>>           : http://www.tys-yokohama.co.jp/
>>>>>>
>>>>>> ==================================
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> OSGeoJapan-discuss mailing list
>>>>>> OSGeoJapan-discuss @ lists.osgeo.org
>>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>>>
>>>>>>
>>>>>>
>>>>>> --
>>>>>> ****************************
>>>>>> 甲府市音羽町 自宅PC
>>>>>>          村石眞澄
>>>>>> ****************************
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> OSGeoJapan-discuss mailing list
>>>>>> OSGeoJapan-discuss @ lists.osgeo.org
>>>>>> http://lists.osgeo.org/mailman/listinfo/osgeojapan-discuss
>>>>>
>>>>>
>>>> --
>>>> ****************************
>>>> 甲府市音羽町 自宅PC
>>>>          村石眞澄
>>>> ****************************
>>>>
>>>
>>
>> --
>> ****************************
>> 甲府市音羽町 自宅PC
>>          村石眞澄
>> ****************************
>>
>
>


-- 
****************************
甲府市音羽町 自宅PC
         村石眞澄
****************************



More information about the OSGeoJapan-discuss mailing list