森林土木memoある自治体職員が、仕事で使うGIS、GPS、CADなどについてメモします。基盤地図対応GDAL/OGRの使い方FaceBookで質問を受けたけど、自分でも忘れていたのでメモ。OsGeo.jpのWebサイトで、基盤地図情報を扱えるGDAL/OGRのバイナリが公開されています。http://www.osgeo.jp/foss4g-mext/このシステムを使用すると、基盤地図情報のJPGIS形式のXMLファイルを一括でShapefileやGeotiffに変換したり、QGISで直接基盤地図情報が読み込めたりします。使い方は、「基盤地図対応GDAL/OGR」の「バイナリ」をクリックし、Zipファイルをダウンロードします。解凍した「mext_gdal」フォルダをCドライブの直下に置くだけです。このままだと使い方がよくわかりませんね。こちらのサイトが参考になります。このサイトはOpenStreetMapの作成に基盤地図情報を利用するための解説ページですが、こちらも「基盤地図対応GDAL/OGR」を使用しています。ページ中段の「GDAL/OGRのダウンロード」で「QGIS_MEXT.zip」をダウンロードし、解凍した「QGIS_MEXT」フォルダをCドライブの直下に置きます。Cドライブの直下に置くかというと、処理を実行するためのバッチファイル(.bat)の中に記載されているシステムの場所が、「C:\QGIS_MEXT\mext_gdal」となっているためです。「QGIS_MEXT」フォルダの中には、「road2shp2tiff_wgs84.bat」と「xml2tif.bat」の2つのバッチファイルが入っています。下画像のリンクをクリックし、「road2shp2tiff_wgs84.bat」をダウンロードして、「QGIS_MEXT」のフォルダに移動します。「road2shp2tiff_wgs84.bat」は、基盤地図情報のXMLファイルをShapefileに変換します。GeoTiffファイルも作成します。「xml2tif.bat」は基盤地図情報の標高DEMをGeoTiffファイルに変換します。このまま使うと、ものすごい時間がかかる場合があるので、すこし編集します。バッチファイルを右クリックし、「編集」でテキストエディタで開きます。100」に変更します。この部分は、tiffファイルの1ピクセルを何mにするかの設定です。ちなみに北海道の行政界を1ピクセル=100mで作成すると、14MB程度でした。この部分は、地図の大きさによって、適宜調整してください。基盤地図情報を「JPGIS形式」でダウンロードし、1つのShapefileにしたい場合は、1つのフォルダにまとめます。(「JPGIS(GML)形式は変換できませんので、必ず「JPGIS形式」で行なってください。)ベクタデータの場合は、「road2shp2tiff_wgs84.bat」のバッチファイルを、まとめたXMLファイルと同じフォルダにコピーし、ダブルクリックで変換が始まります。変換されたShapefileとGeotiffファイルをQGISなどで開いて確認してみてください。<バッチファイルの中身の説明>バッチファイルの中身は、次のような命令が記載されています。ROOTDIR=C:\QGIS_MEXT\mext_gdal↑ここはGDALシステムのあるフォルダを指定しています。PATH=%ROOTDIR%\lib;%ROOTDIR%\bin;%PATH%↑ここは特にイジる必要はありません。%%i↑ここは、フォルダ内のXMLファイルをすべてShapeFileに変換する命令です。「ENCODING=CP932」とすると、属性データの文字コードがShift-JISになります。「EPSG:」で変換後の測地系を設定できます。平面直角座標系JGD2000の12系なら「EPSG:2454」とします。これはQGISの座標参照系と同じですので、番号はQGISで確認するといいと思います。%%i↑一つ上の命令で作成した個別のShapefileを結合して、「Shp1.shp」として保存します。ファイル名は自由に設定できます。shp1.tif↑作成した「Shp1.shp」を「Shp1.tif」というGeotiffファイルで保存します。「-tr」のあとの数字は、1ピクセルを何mにするかという設定です。小さくすると、ものすごい大きなファイルができるので注意してください。Geotiffファイルを作成する必要がない場合には、命令の先頭に「#」を入力するとコメントになりますので、命令が実行されません。tiffファイルが作成されない場合があります。(原因不明)例えば、文字コードShift-JISで、平面直角座標系JGD200の12系でShapefileを作成し、Geotiffファイルが不要な場合には、次のように命令を記載します。shp1_JGD2000_12.tif=====================================また、このパッケージには、基盤地図情報を直接読み込めるQGISも付属しています。「C:\QGIS_MEXT\mext_gdal」の中の「qgis1.9.bat」をダブルクリックすると、QGISが開きます。ROOTDIR=C:\QGIS_MEXT\mext_gdal」となっているか確認してください。QGISで基盤地図情報を直接読み込むと、処理が遅いのであまり実用的ではないなという感想です。パワーのあるパソコンならいいのかもしれません。|この記事へのコメントこんにちは。わかりやすい解説を作って頂き,ありがとうございます!kitaさんはおわかりだと思うのですが,ちょっとコメントを。0.00001」はご指摘の通り作成したGeoTiff画像の解像度です。このバッチプログラム,初期の設定では変換後の座標系が「EPSG:4326」つまりWGS84の緯度経度になっているので,解像度も度で計算されています。つまる,「0.00001」度で,およそ南北方向で1mぐらいになります。100」とすると解像度100mのGeoTiff画像ができる事になります。100」にすると解像度が100mではなく100度の画像ができる事になります。2.0業務で使う林業QGIS徹底使いこなしガイドブログ主が書きましたQGISのマニュアル本です。2017年6月発売!「林業」とありますが、林業だけでなくQGISの基本から応用まで幅広く網羅しています。御覧ください!Amazonでも販売中!下のリンクから、本で利用できるサンプルファイルをダウンロードできます。業務で使う林業QGIS徹底使いこなしガイド本の紹介ページ(全林協)公開しているテキスト当ブログ管理人が作成したQuantumGISのテキストなどをダウンロード出来ます。北海道Shapeファイル北海道全域の基盤地図情報をShapeに変換しました。)、「道路縁」があります。(04/19)ファンメッセージを送るこのブログの読者になる更新情報をチェックするブックマークする友達に教えるカテゴリ日記(14)その他(2)CAD(1)GIMP(1)林道(1)QGIS(50)GPS(4)リンク(1)DBF編集(2)GISその他(22)魚道(0)治山ダム(1)Android(10)カシミール3D(6)vb.net(2)Web地図(5)過去ログ2019年10月(2)2019年02月(1)2018年09月(3)2018年08月(1)2018年03月(1)2017年12月(1)2017年08月(1)2017年07月(3)2017年06月(2)2017年05月(2)2017年02月(1)2017年01月(1)2016年12月(1)2016年11月(4)2016年09月(1)2016年08月(3)2016年07月(3)2016年03月(1)2015年12月(1)2015年09月(2)Seesaaブログ。