[OSGeoJapan-discuss] QGIS1.4.0のトラブルについて（解決しました！）

[Hayato IIJIMA]

嘉山様、今木様、皆様

今木様に質問をさせていただきました、山梨県森林総合研究所の飯島と申します。まず、非常に曖昧な質問をしてしまったことに対してお詫びを申し上げます。そして、昨日QGISを再インストールしてみたところ、問題としていた点が解決してしまいました。お騒がせして申し訳ありませんでした。ですが、何かの参考になればと思い、以下長文になりますが、昨日まで解決していなかった問題の状況を述べさせていただきます。

【QGISで実行したいこと】
・異なる由来の地図をQGIS上で重ね合わせる
・重ね合わせた地図について、片方の地図（ベクタ形式の点データ）の点の周囲にバッファを発生させ、バッファに含まれる別の地図のデータ（ベクタ形式のポリゴンデータ）の面積を集計したい

このうち、QGIS上で異なる由来の地図を重ねて表示させることは、On the
Fly投影をOnにすることによって、再インストールする前のQGIS ver.
1.4.0でも実行できました。2番目の、片方のレイヤの点から発生させたバッファに含まれる他のレイヤの値を集計することが上手くいっていませんでした。この要因は、それぞれの地図の測地系が統一されていないことが原因であると考えました。そこで、まず両方のレイヤの測地系を平面直角座標で統一しようとしたのですが、再インストールするまでは測地系の定義が上手くいっていないようでした。

【実行環境】
・OS: Windows XP Service Pack 3
・筐体：Panasonic Let’s note Light CF-W2
・QGISのバージョン：1.4.0 Enceladus
・QGISの入手方法：QGISの公式ホームページ（http://qgis.org/）から、オンラインインストーラ（Windows -
OSGeo4W）をダウンロードし、インストール

【使用した地図データ】
・国土数値情報（http://nlftp.mlit.go.jp/ksj/jpgis/jpgis_datalist.html）で、JPGIS準拠のデータから取得した、平成18年度世界測地系版土地利用細分メッシュの山梨県部分（L03-b-06_5338-jgd.zip）。地図の測地系は、JGD2000で緯度経度。
・GPSで落とした調査地（14地点）の位置データ
これらを、以下に述べる方法でQGISに読み込んだ。

【地図データのQGISへの取り込み方法】
・国土数値情報のデータ：同サイトで配布されている、ksjtoolを用い、shpファイルに変換
・平面直角座標・緯度経度オンライン変換（http://www.n-survey.com/online/xybl.htm）サイトで、GPSの緯度経度を平面直角座標に変換し、それらをMicrosoft
Excelに打ち込んで、csvファイルとした。これを、「デリミティッドテキストレイヤの読み込み」プラグインでQGISに読み込んだ。

【作業手順】
最終的には、GPSの点データの周辺に「ftoolsプラグイン」に含まれる「空間演算ツール」の「バッファ」でバッファを発生させ、発生させたバッファに含まれる国土数値情報のポリゴンの面積を計算させたかったため、XY座標で投影されている平面直角座標（JGD2000）に変換する方がいいと考えたため、緯度経度形式で配布されている国土数値情報のデータを平面直角座標に変換しようと考えた（GPSデータは上記の方法ですでに平面直角座標になっている）。

国土数値情報のデータは、ダウンロードしたサイトによればJGD2000で緯度経度で表現されている。そのため、座標参照系はEPSG:
4612であるJGD2000の測地系であると判断。

1. QGISを立ち上げ、On the Fly投影をOnにし、座標参照系をEPSG: 4612とした。
2. 国土数値情報のshpファイルを読み込んだ。
3. 読み込んだ国土数値情報のレイヤを右クリックし、プロパティからCRSの指定で現在の座標参照系をEPSG: 4612とした。
4. ftoolsのデータマネージメントツールの「新しい投影法へエクスポートする」を選ぶ

ここで、前の座標参照系の定義が上手くいっていれば、「Input spatial reference
system」の所に、3.で定義した座標参照系（+proj=longlat +ellp=GRS80
+towgs84=0,0,0,0,0,0,0 +no_def）←つまりはEPSG:
4612であるJGD2000が入力されているはずなのですが、QGISを再インストールするまではここが空欄になっていました。

そのままでも一応新しい投影法へのエクスポートは成功するのですが、緯度経度形式から例えばJGD2000/UTM zone 54N（EPSG:
3100）形式に変換すると、緯度経度の値がそのままメートル（QGISの下部に表示される座標が緯度経度の値のまま）になってしまっていました。

このとき考えたのは、「新しい投影法へエクスポートする」のときに、Input spatial reference
systemが空欄になっていたことから、座標参照系の定義が上手くいっていないのではないかと言うことです。

ですが、昨日QGIS1.4.0をオンラインインストーラーで再インストールしたところ、同じ動作を行ったときに、Input spatial
reference systemのところにプロパティで指定した座標参照系が入っていました。ということで、問題が解決しました。

今回の顛末は以上です。お騒がせして申し訳ありませんでした。なぜ再インストールで上手くいったのかわかりませんが、もしも同じような現象が発生している方がいらっしゃいましたら、参考になれば幸いです。今後ともよろしくお願い致します。


-- 
**********************************************************
山梨県森林総合研究所森林研究部森林保護科

研究員　飯島勇人（いいじまはやと）

〒400-0502
山梨県南巨摩郡富士川町最勝寺2290-1
TEL:0556-22-8005（直）8001（代）
E-mail: iijima-akks ＠ pref.yamanashi.lg.jp
**********************************************************