tag:blogger.com,1999:blog-77517099620873670402024-03-14T00:04:02.283+09:00spatiohackTakayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.comBlogger50125tag:blogger.com,1999:blog-7751709962087367040.post-72105732299507066662018-04-16T13:21:00.000+09:002018-04-16T13:21:10.694+09:00ESRI ShapefileからGeoPackageへの移行QGIS 3.0より新規ベクターデータのデフォルト保存形式がESRI Shapefile形式からGeoPackage形式に変わったようです。<br />
<br />
ESRI Shapefile形式は<a href="https://blog.chizuwota.net/2016/12/geopackage-osm.html">複数のファイルから構成されるとか属性名の日本語の文字数制約とか</a>、いろいろと問題のある形式でしたが、先日の<a href="https://qgis.connpass.com/event/82120/">QGIS 3.0 Hackfest</a>でいろいろと使ってみたところGeoPackage形式ではそれらの問題は大丈夫そうな感じがしました。<br />
<br />
とはいえ自分の持っている全てのESRI ShapefileをGeoPackageにスクリプトで一括変換するのはまだちょっと怖いので、重要性の低いデータ(Webで入手可能とか授業で例として適当に作ったやつとか)からちょっとづつ変換して行こうと思っています。<br />
<br />
QGIS上で名前を付けて保存でも良いのですが、もうちょっとだけ楽にするためにUbuntu/Debianのnautilus scriptという機能を使って対象ファイルの右クリックからESRI Shapefile -> GeoPackage変換ができるようにしました。<br />
<br />
<br />
まず以下のようなbashスクリプト(ファイル名: shp2gpkg.sh)を作成します<br />
<br />
<br />
#!/bin/bash<br />
<br />
shp_filename=$1<br />
root_filename=${shp_filename%.shp}<br />
ogr2ogr -f GPKG ${root_filename}.gpkg ${shp_filename}<br />
<br />
<br />
それをホームディレクトリの、.local/share/nautilus/script に実行可能にして配置します(シンボリックリンクでもOK)。<br />
<br />
使用する際は、ESRI Shapefile上で右クリック -> scripts -> shp2gpkg.sh と選択します。<br />
変換元のファイルは一応残すようにしています。<br />
<br />
<br />
<br />
先日のQGIS 3.0 Hackfestでの共有ドキュメントはどなたでも参考になれます。GeoPackageの検証以外にも様々なQGIS3.0の機能を多くの参加者で試して得られた知見があります。<br />
<br />
<a href="https://hackmd.io/Y5z0jxdxS7Wxw3D57-116g?view">QGIS 3.0 Hackmd</a><br />
<br />
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-24682023671278360222018-01-09T14:49:00.001+09:002018-01-09T14:49:28.047+09:00Windows 10のWSL(Windows Subsystem for Linux)にJupyter notebookをインストール<h3>
使用した環境:</h3>
Windows 10 Pro (version 1709) 64 bit<br />
<br />
<br />
<h3>
WSLのセットアップ:</h3>
<a href="http://www.atmarkit.co.jp/ait/articles/1608/08/news039.html">こちらのサイト</a>などを参考に。今回はDistributionはUbuntuを選択。<br />
<br />
<br />
<br />
<h3>
Jupyter notebookのインストール(今回はPython 3での環境構築)</h3>
<br />
# まずpipというPythonoパッケージマネージャをインストール<br />
sudo apt install python3-pip<br />
<br />
# Python3 なので pip3 と数字の3がつく。まずはパッケージマネージャ自身の更新<br />
pip3 install -U pip<br />
<br />
# Jupyter notebookのインストール(管理者権限が必要)<br />
sudo pip3 install jupyter<br />
<br />
# Jupyter notebookのカーネル起動<br />
jupyter notebook<br />
<br />
<br />
<br />
しかしWSLから直接Windowsのブラウザ起動できないっぽいのでターミナルに表示される<br />
localhost:8888なんとか<br />
というアドレスをコピーしてブラウザURLに貼り付けると開ける。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjcfrZRJQ42SZAee5aL6h8lMim0F0220V0pPZanCv18tgsCxefkNJgRrRlqiIghBB-fYHH67PZI0p5GBXspqcRcsQCKBAZS-l_z6szD1KDgKpBtLLNLfCtMaDRWZkOI3D8LduXZzhWX8jP9/s1600/jupyter_on_wsl.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="763" data-original-width="1213" height="250" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjcfrZRJQ42SZAee5aL6h8lMim0F0220V0pPZanCv18tgsCxefkNJgRrRlqiIghBB-fYHH67PZI0p5GBXspqcRcsQCKBAZS-l_z6szD1KDgKpBtLLNLfCtMaDRWZkOI3D8LduXZzhWX8jP9/s400/jupyter_on_wsl.png" width="400" /></a></div>
<br />
<br />
<br />
参考サイト<br />
<a href="https://qiita.com/koara-local/items/8d48d0306c85c61aed97">[Windows][Python3] Windowsにpython3とJupyter Notebook(旧: ipython notebook)をインストールする</a><br />
<br />
<br />
<br />
<br />
<br />
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-2528291283580691812017-09-20T09:58:00.005+09:002017-09-20T09:59:33.886+09:00ネパールヒマラヤ、リルン氷河の変動の論文がでましたー<a href="http://www.sciencedirect.com/science/journal/10406182?sdc=1">Quaternary International</a>にてネパールヒマラヤ、リルン氷河の変動の論文がでました。<br />
2008年の現地観測データと各種のリモートセンシングデータをまとめたものです。時間がかかってしまいました…。<br />
<br />
Elsevier系での論文は初めてだったのですが、科研費ゼロの身としては出版料がタダなのは助かりました。(本当はオープンアクセスにしたかったけれども…)<br />
<br />
ただ正式出版後50日間(2017/11/08まで)は論文が無料でダウンロードできるとのことですので、もしよろしければ以下のリンクからダウンロードしてください。<br />
<br />
<a data-saferedirecturl="https://www.google.com/url?hl=en&q=https://authors.elsevier.com/a/1VlBj3ic-F8L7F&source=gmail&ust=1505953381308000&usg=AFQjCNF-0-IPmqeKEgcQj87vqDZZhQn8nw" href="https://authors.elsevier.com/a/1VlBj3ic-F8L7F" rel="noreferrer" style="background-color: white; color: #1155cc; font-family: arial, sans-serif; font-size: 12.8px;" target="_blank">https://authors.elsevier.com/<wbr></wbr>a/1VlBj3ic-F8L7F</a><br />
<br />
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-54492777050031299812017-01-19T13:32:00.000+09:002017-03-15T15:06:01.902+09:00GMTによる散布図作成の基礎GMTというと地図作成ソフトで知られますが、普通のグラフも作ることができます。<br />
今日は簡単な散布図作成につい紹介します。<br />
<br />
<a href="http://spatiohack.blogspot.jp/2016/12/gmt-531-cygwin.html">GMTのインストールは既に完了しているものとします。</a><br />
<br />
<br />
<h3>
1. ファイルの説明</h3>
以下のファイルをダウンロードして同じフォルダにおきます。<br />
<br />
<ul>
<li><a href="https://raw.githubusercontent.com/nuimura/gmt_xyplot_basic/master/sample_xy.sh">sample_xy.sh</a></li>
<li><a href="https://raw.githubusercontent.com/nuimura/gmt_xyplot_basic/master/data1.csv">data1.csv</a></li>
<li><a href="https://raw.githubusercontent.com/nuimura/gmt_xyplot_basic/master/data2.csv">data2.csv</a></li>
</ul>
<br />
sample_xy.sh はGMTのコマンドが並べられたシェルスクリプトで、メモ帳などのエディタで開くことができます。ただし改行コードがUnix(LF)となっているためWindowsのメモ帳では1行に表示されてしまいます。メモ帳ではなく<a href="https://notepad-plus-plus.org/">Notepad++</a>などの正しく様々な改行コードを認識できるエディタを使うことを進めます。<br />
<br />
data1.csv、data2.csv は適当に作成した2変数のスペース区切りデータで、以下のような内容です。<br />
0 1000<br />
-5 900<br />
-18 800<br />
-14 700<br />
-15 600<br />
<br />
<br />
<h3>
2. GMTコマンドの説明</h3>
sample_xy.shの中身は以下のとおりです。<br />
<br />
<script src="https://gist-it.appspot.com/github/nuimura/gmt_xyplot_basic/raw/master/sample_xy.sh"></script>
順番に説明をすると、
<br />
1行目はシェルスクリプト(bash)のおまじない※無くても良い<br />
<br />
4行目はグラフの枠とdata1.csvのプロット、それぞれのオプションの意味は<br />
-Jはグラフの横と縦の長さ、縦をマイナスにすると数値の大小が逆向きになる。<br />
-Rは値の範囲、X軸は-30から30、Y軸は500から1000<br />
-BxはX軸の目盛間隔とラベル<br />
-ByはY軸の目盛間隔とラベル<br />
3つめの-Bは目盛をふる側(WとSだけ大文字なので左と下の目盛をふる)とタイトル<br />
-Wは線の太さと色(ここでは黒色)と種類(実線)<br />
-Vはターミナルに情報表示<br />
-Pはグラフの向きを縦に<br />
-KはGMT描画コマンド(psxy,pscoast,psbasemapなど、psconvertは関係ない)が2つ以上の場合、最後のコマンド以外につける必要がある。「まだ続きがあるよ」という意味<br />
<br />
6行目はdata2.csvのプロット、それぞれのオプションの意味は<br />
-Jは2回目なので省略可能<br />
-Rも2回目なので省略可能<br />
-Wで線の太さと色(ここでは赤色)と種類(破線)<br />
-Vはターミナルに情報表示<br />
-Pはグラフの向きを縦に<br />
-OはGMT描画コマンド(psxy,pscoast,psbasemapなど、psconvertは関係ない)が2つ以上の場合、最初のコマンド以外につける必要がある。「前のコマンドの続きだよ」という意味<br />
-Kは先程と同じく「まだ続きがあるよ」という意味<br />
<br />
8-14行目で凡例の指定<br />
<br />
17行目でpsファイルをPDFファイルに変換<br />
18行目でpsファイルをPNGファイルに変換<br />
<br />
<br />
3. 実行<br />
ファイルのある場所に移動してから、<br />
bash sample_xy.sh<br />
<br />
とターミナルにコマンドを入力します。<br />
sample_xy.shをテキストエディタで開いて一行一行ターミナルにコピペしてもできますが、上のようにコマンドを入力したほうが簡単です。<br />
<br />
<br />
出来上がったグラフはこんな感じ<br />
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjlnjm4kixVnUsDXo9f0lHfZilPMPjydn2d28ZqgydqqdNDzZYKkNZ5X7dNNpxg3zZ7XzeEf4hbnSVKvIFEwOMy4rlDTvTqq9Fudw2gJq1b5f4v5qt4tZxt_DpPzt355wBTCkmHX_nx_RuS/s1600/figure01.png" imageanchor="1"><img border="0" height="320" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjlnjm4kixVnUsDXo9f0lHfZilPMPjydn2d28ZqgydqqdNDzZYKkNZ5X7dNNpxg3zZ7XzeEf4hbnSVKvIFEwOMy4rlDTvTqq9Fudw2gJq1b5f4v5qt4tZxt_DpPzt355wBTCkmHX_nx_RuS/s320/figure01.png" width="305" /></a></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<span style="color: red;"><br /></span>
<span style="color: red;">※2017/1/21 線の種類指定の追加と凡例のプロットを追加</span><br />
<div>
<span style="color: red;">※2017/2/6 ℃記号を修正しました(^_^;)。</span></div>
Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-64073372904613712712016-12-22T17:50:00.002+09:002016-12-29T21:42:11.875+09:00GMTによる地図作成の基礎毎年この時期になるとGMTやLaTeX関連の記事のアクセス数が増加してきて、卒論・修論・D論が架橋に差し掛かっているのだなぁと季節を感じます。<span style="font-size: x-small;">(いや、D論はもう方がついていないとやばいですね)</span><br />
<br />
<br />
<h3>
0. 前準備</h3>
まず、<a href="http://gmt.soest.hawaii.edu/">GMT</a>をインストールします。<a href="http://spatiohack.blogspot.jp/2016/09/gmt-521bash-on-ubuntu-on-windows.html">Bash on Windows</a>でも<a href="http://spatiohack.blogspot.jp/2016/12/gmt-531-cygwin.html">Cygwin</a>でもLinuxでも構いません。この記事ではとりあえずWindows10でCygwinという環境ということで説明していきます。<br />
<br />
データの前処理として<a href="http://qgis.org/ja/site/">QGIS</a>も使いますのでそちらもインストールしてください。<br />
<br />
<br />
<h3>
1. データダウンロード</h3>
今回はUSGS(米国地質調査所)のSRTM3データ(3秒メッシュ≒100 mメッシュ)を用います。<a href="https://dds.cr.usgs.gov/srtm/version2_1/">USGSのダウンロードページ</a>から、SRTM3/Eurasia とたどり、以下のファイルをダウンロード&解凍して同じフォルダ(例. srtm といったわかりやすい名前のフォルダ)に入れておいてください。<br />
<br />
N35E139.hgt.zip<br />
N36E139.hgt.zip<br />
N35E140.hgt.zip<br />
N36E140.hgt.zip<br />
<br />
<br />
2. QGISでデータの変換<br />
先ほど解凍して同じフォルダに入れたhgtファイルをQGIS上で結合して、GMTで読み込める形式(NetCDF形式)に変換します。<br />
<br />
メニューの「ラスター」→「その他」→「結合」と選択して、以下の結合Windowを立ち上げます。<br />
<br />
「ファイルの代わりに入力ディレクトリを選択する」にチェックを入れて、入力ディレクトリ(hgtファイルが入っているフォルダ)を指定、出力ファイル名(GeoTIFF形式で)を指定します。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhyFk0IbLC3LeiiqTe2-m0oELDQOS2yza45RRTeJPx_RaLfmBvjF5e5JX1FsF3U7uGjC86naOCNDIMqyhgHkn3PgqUfuCrVcDmMcJIKF55yyAlRjSC5TQDaRY-LS2D90R99AqJOCpxgotyw/s1600/merge_window.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="320" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhyFk0IbLC3LeiiqTe2-m0oELDQOS2yza45RRTeJPx_RaLfmBvjF5e5JX1FsF3U7uGjC86naOCNDIMqyhgHkn3PgqUfuCrVcDmMcJIKF55yyAlRjSC5TQDaRY-LS2D90R99AqJOCpxgotyw/s320/merge_window.png" width="279" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
図1. 結合Window</div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjrQ0diSJvO_RSPoKfq3TmTaM87qSNJLEmyySwxybMTmsfdWYnZaChRXnikrZ54sdelqQRTF4gYWidgAZgbdNzY1n9kiy4tMVaAHD7cdc7Mq1W3agWxPC4NZGmrGkoyBPLc4bh2TplZaOfl/s1600/dem_output.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="200" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjrQ0diSJvO_RSPoKfq3TmTaM87qSNJLEmyySwxybMTmsfdWYnZaChRXnikrZ54sdelqQRTF4gYWidgAZgbdNzY1n9kiy4tMVaAHD7cdc7Mq1W3agWxPC4NZGmrGkoyBPLc4bh2TplZaOfl/s320/dem_output.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
図2. 結合結果(色は見やすく変えてあります。0の値を透過設定)</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
メニューの「ラスター」→「変換」→「変換(形式変換)」を選択し、出力形式に「GMT NetCDF Grid Format」を選択して保存します。ここではsrtm2_dem.nc というファイル名で保存をしました。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<br />
<h3>
3. GMTで地形表示</h3>
<br />
Cygwinを立ち上げて、データの置いてあるディレクトリ(=フォルダ)に移動します。<br />
<br />
今回使用する地形データは海域に0の値が入っていますので、まずは以下のコマンドでNo data値が正しく設定されたNetCDFファイルを作成します。<br />
<br />
gmt grdmath srtm3_dem.nc 0 NAN = srtm3_dem_withnodata.nc<br />
<br />
上のコマンドはファイル変換に使うだけですので、最初1回だけ実行すればOKです。<br />
srtm3_dem_withnodata.nc というファイルが作成されます。<br />
<br />
<br />
GMT上で下記の3行のコマンドを実行します。<br />
<br />
<b><span style="font-family: "courier new" , "courier" , monospace;">gmt makecpt -Cdem1 -T0/3000/100 -Z > sample.cpt</span></b><br />
<b><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R139/141/35/37 -Jm1:2000000 -Ba1f1g1 -K > output.ps</span></b><br />
<b><span style="font-family: "courier new" , "courier" , monospace;">gmt grdimage srtm3_dem_withnodata.nc -Csample.cpt -Q -R -J -O >> output.ps</span></b><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjyT8Kw1r2cC4WSHnpq-R1Kbi1q3tQ0N8R5SHj9s_HOznM-rlAkGJMVdC1pz8fgp9kZCNV7DIitGSNQIHduxXjinzJKFKuo5fFUlN0Z0jxcjDuGvN60B4wV2-hGhNoaoZnZC_dswMDr6kq0/s1600/gmt1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="320" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjyT8Kw1r2cC4WSHnpq-R1Kbi1q3tQ0N8R5SHj9s_HOznM-rlAkGJMVdC1pz8fgp9kZCNV7DIitGSNQIHduxXjinzJKFKuo5fFUlN0Z0jxcjDuGvN60B4wV2-hGhNoaoZnZC_dswMDr6kq0/s320/gmt1.png" width="225" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
図3. 出力結果</div>
<div>
<br /></div>
<br />
<h3>
4. GMTで地形に他の情報の重ね合わせ</h3>
緯度経度の列を持つ下記のようなデータを用意します。(今回は適当に作りました)<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhzBs2NOtgNepFYAd2Pgd12BRkTj6WdmoYNOPUMm_6fmhVk8vSE3N72YTYabL7BIHWTSjEUoYlvw4WpvVLpTJd7zJ0_qN34kwpeaR9mb15hfyAN0CfmnyI1_Q1USkbJWm7bg-ro6HoMYSi4/s1600/excel_sample.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="304" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhzBs2NOtgNepFYAd2Pgd12BRkTj6WdmoYNOPUMm_6fmhVk8vSE3N72YTYabL7BIHWTSjEUoYlvw4WpvVLpTJd7zJ0_qN34kwpeaR9mb15hfyAN0CfmnyI1_Q1USkbJWm7bg-ro6HoMYSi4/s320/excel_sample.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
図4. ダミーデータ</div>
<br />
このようなデータをCSV形式で保存します。<br />
<br />
<br />
GMT上で下記の4行のコマンドを実行します。<br />
<b><span style="font-family: "courier new" , "courier" , monospace;">gmt makecpt -Cdem1 -T0/3000/100 -Z > sample.cpt</span></b><br />
<b><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R139/141/35/37 -Jm1:2000000 -Ba1f1g1 -K > output.ps</span></b><br />
<b><span style="font-family: "courier new" , "courier" , monospace;">gmt grdimage srtm3_dem_withnodata.nc -Csample.cpt -Q -R -J -O -K >> output.ps</span></b><br />
<b><span style="font-family: "courier new" , "courier" , monospace;">gmt psxy data.csv -Sc0.6 -W1,black -Gred -R -J -O >> output.ps</span></b><br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh2N53L7g72V-tEM6hytwyUrmuWrkgz9okSxJpy6TTG9TOtk08ByNh1OnZRT324LmOZSD3o1UbhnwkYqQ6Djrz_e6pyIZwO3ws1_wSAISvJEw5FKwrgzsHuTvBnjIbi1oOKLSLHCiQjy8QW/s1600/gmt2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="320" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh2N53L7g72V-tEM6hytwyUrmuWrkgz9okSxJpy6TTG9TOtk08ByNh1OnZRT324LmOZSD3o1UbhnwkYqQ6Djrz_e6pyIZwO3ws1_wSAISvJEw5FKwrgzsHuTvBnjIbi1oOKLSLHCiQjy8QW/s320/gmt2.png" width="225" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
図3. 出力結果</div>
<br />
するとこのような図が作成できます。<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-88930692591399527942016-12-22T16:34:00.001+09:002016-12-22T16:39:53.217+09:00GMT 5.3.1のインストール (Cygwin編)<span style="color: #cc0000; font-size: x-small;">2012年の</span><span style="color: #cc0000; font-size: x-small;"><a href="http://spatiohack.blogspot.jp/2012/03/gmt-457.html">GMT 4.5.7のインストール (Cygwin編) </a>に多くのアクセスがありますが、現在のGMTの最新バージョンは5系列です。こちらは過去の記事の内容を最新バージョン用にアップデートした記事となります。</span><br />
<br />
<a href="http://gmt.soest.hawaii.edu/">GMT (Generic Mapping Tools)</a>とはコマンドラインベースの地図作成ツールで、地球科学分野でよく用いられるツールで、綺麗な空間分布図を描画することができます。クロスプラットフォームなソフトなのでWindows、Mac、Linuxで使うことができます。<br />
<br />
<br />
今回はWindows 10でのGMTのインストール(&初期設定)について説明します。<br />
<br />
<hr />
<h4>
1. はじめに</h4>
<br />
WindowsでGMTを使う場合は、いくつかの選択肢があります。<br />
<br />
<ol>
<li><a href="http://spatiohack.blogspot.jp/2016/09/gmt-521bash-on-ubuntu-on-windows.html">Bash on Windows でGMTをコンパイル・インストール</a></li>
<li>昔ながらのCygwinでGMTを使用</li>
</ol>
<br />
今回は2番目のCygwinを用いた方法を紹介します。<br />
以下のツールをインストールする必要があります。<br />
<br />
<br />
<b><span style="color: #073763;">必須</span></b><br />
<div style="background-color: #e8e8ff; border-radius: 10px;">
<dl>
<dt><a href="http://pages.cs.wisc.edu/~ghost/gsview/get50.htm">GSView</a></dt>
<dd>GMTで出力したps/epsファイルの表示や変換処理に必要 (※実際に変換しているのはGhostscriptやpstoedit、GSViewはインターフェース)。</dd><dd></dd>
<dt><a href="https://www.sumatrapdfreader.org/free-pdf-reader.html">SumatraPDF</a></dt>
<dd>GMTで出力したps/eps/pdfファイルの表示に便利。Adobe readerで表示すると排他アクセス処理が働くので、こちらのほうが良い。</dd></dl>
</div>
<span style="color: #073763;"><b>推奨</b></span><br />
<div style="background-color: #e8e8ff; border-radius: 10px;">
<dl>
<dt><a href="http://www.cygwin.com/">Cygwin</a></dt>
<dd>シェルスクリプトでGMTコマンドを使う場合必要。WindowsバッチでもGMTは扱えるがCygwinの方が多機能で便利。</dd>
<dt><a href="http://www.ghostscript.com/download/gsdnld.html">Ghostscript</a></dt>
<dd>GMTで出力したps/epsファイルをビットマップ画像に変換するのに必要。</dd>
<dt><a href="http://www.pstoedit.net/">PStoedit</a></dt>
<dd>GMTで出力したps/epsファイルをベクター形式のままSVG形式に変換できる。Inkscapeなどのイラストソフトを使っている人には便利。</dd>
<dt>改行コードやエンコーディングが変更できるエディタ</dt>
<dd><span style="font-size: x-small;">(<a href="http://www.geocities.jp/npp_euc_jp/">Notopad++EUC-JP対応版</a>がおすすめ、自分は<a href="http://vgoulet.act.ulaval.ca/en/emacs/">emacs modified</a>を使ってるけど…)</span><br />
異なるOSや他人とのファイルの移動を行う場合は、環境の違いによって日本語コメントなどが文字化けする場合がある。<br />
受け渡し先に合わせたエンコーディングや改行コードに変換したり、人からもらったスクリプトファイルのエンコーディング&改行コードが自分の環境に合わない場合などにあると便利<span style="font-size: x-small;">(でも本当はコメントは英語で書いておくのが無難)</span>。</dd></dl>
</div>
<hr />
<h4>
2. インストール</h4>
<a href="http://gmt.soest.hawaii.edu/projects/gmt/wiki/Download">GMTサイトのwindows用ダウンロードのページ</a>から、最新バージョンのインストーラをダウンロードして、インストールを行います。<br />
<br />
ちなみに2016/12/22現在の最新バージョンは5.3.1なので<br />
<ul>
<li>Windows (64 bit)の人は、 <i>gmt-5.3.1_install64.exe</i></li>
<li>Windows (32 bit)の人は、 <i>gmt-<i>5.3.1</i>_install32.exe</i></li>
</ul>
<br />
<br />
GMTの実行ファイルは、 <i>C:\programs\gmt5 </i>フォルダにあります。 昔のGMTではインストール後(というかインストーラが無かったので、圧縮された実行ファイルをダウンロード後、解凍して任意の場所に置くだけ)、パスの設定を行う必要がありましたが、最近のGMTでは自動でパスの設定をしてくれます。<br />
<br />
なのでコマンドプロンプトやCygwinターミナルで何かGMTのコマンド、例えば、<br />
<br />
gmt psxy<br />
<br />
<span style="color: red; font-size: x-small;">※GMT4まではpsxyなどのコマンドだけで良かったのですが、GMT5よりgmtを最初に付ける必要があります。</span><br />
<br />
を打ち込むと、大量のメッセージが出てきたらインストール成功です。<br />
もし「'gmt psxy'は、内部コマンド~(中略)~として認識されていません」とか「gmt psxy:コマンドが見つかりません」などと出てきたらインストール失敗です。そのようなことは通常無いと思うのですが(少なくとも自分は見たことないです)、万が一そのようなことが起こった場合は自分で環境変数の設定をしてパスを通す必要があります。 必要なユーザー環境変数は、以下のとおりです。<br />
<br />
<div align="center">
<table border="0" cellpadding="3" cellspacing="0" style="text-align: center;"><tbody>
<tr><th style="background-color: #cacaff;">変数</th><th style="background-color: #cacaff;">値</th></tr>
<tr></tr>
<tr><td>GMT_SHAREDIR</td><td><i>C:\programs\gmt5\share</i></td></tr>
<tr><td>path</td><td><i>C:\programs\</i><i>gmt5</i><i>\bin</i></td></tr>
</tbody></table>
</div>
<br />
※設定後は再起動が必要です(たぶん)。<br />
<br />
<br />
上記のGMT&周辺ツールのインストールが終わったら、Cygwinターミナルで試しに<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><b>gmt pscoast -Jm1:30000000 -R120/150/20/50 -Ba10f5g5 -Gtan > sample.ps</b></span><br />
<br />
と打ってみて下さい。<br />
<span style="color: #cc0000; font-size: x-small;">※表示の都合で2行にまたがっているかもしれませんが、1行で入力して下さい</span><br />
<span style="color: #cc0000; font-size: x-small;"><br /></span>カレントディレクトリにsample.psというPSファイルができているはずです。このPSファイルをWindowsで見るために、SumatraPDFやGsviewが必要となります。<br />
<br />
ちなみにPSファイルをPDFファイルに変換したい場合は<br />
<b><span style="font-family: Courier New, Courier, monospace;">gmt psconvert -A -Tf world.ps</span></b><br />
<br />
ちなみにPSファイルをPNGファイルに変換したい場合は<br />
<b><span style="font-family: Courier New, Courier, monospace;">gmt psconvert -A -Tg world.ps</span></b><br />
<br />
を実行するとそれぞれ変換できます。<br />
<br />
<br />
<table style="width: auto;"><tbody>
<tr><td><div class="separator" style="clear: both; text-align: center;">
<a href="https://picasaweb.google.com/lh/photo/RgEelmIAJlKABQddJQS3ntMTjNZETYmyPJy0liipFm0?feat=embedwebsite" style="margin-left: 1em; margin-right: 1em;"><img height="400" src="https://lh6.googleusercontent.com/-MHYKd7eT2v0/T3PKy4QnpGI/AAAAAAABk1g/NeCvDKXH6eQ/s400/sample.png" width="345" /></a></div>
</td></tr>
<tr><td style="font-family: arial,sans-serif; font-size: 11px; text-align: right;">From <a href="https://picasaweb.google.com/104578933986350101585/Research_public?authuser=0&feat=embedwebsite">Research_public</a></td></tr>
</tbody></table>
pscoastコマンドで作成した図Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-25689765825152829022016-12-19T11:59:00.002+09:002016-12-19T11:59:29.597+09:00QGIS processing でのRアルゴリズムの自作例 (気象庁の気象観測データのGISへのインポート)この記事は<a href="http://qiita.com/advent-calendar/2016/foss4g">FOSS4G Advent Calendar 2016</a> の18日目の記事として書きました。<br />
<br />
色んな場所のQGISハンズオンにて、QGIS processing機能による他のソフトウェアとの連携処理 (R言語とかGRASS GISとか)のお話を最近ちょこちょこしています。<br />
<br />
Slideshareに資料をあげているものだと<br />
<br />
<ul>
<li>FOSS4G 2014 Tokyo「QGIS&Rハンズオン」 (2014/10/31) @東京大学 (<a href="http://www.slideshare.net/nuimura/141031-qgisr-handsonpresentation">Slideshare</a>)</li>
<li><a href="http://www.chikyu.ac.jp/publicity/events/etc/2015/0827-29.html">地球環境GIS講習会2015</a>「QGISで環境リモートセンシング」(2015/08/28) @地球研 (<a href="http://www.slideshare.net/nuimura/150828-rihn-gisworkshophandsonpresentation">Slideshare</a>)</li>
<li><a href="http://www.osgeo.jp/events/foss4g-2016/foss4g-2016-tokyo/foss4g-2016-tokyo-handson">FOSS4G 2016 Tokyo</a>「QGIS中級+プロセッシングハンズオン」(2016/11/14) @東京大学 (<a href="http://www.slideshare.net/nuimura/161104-foss4-gtokyoqgishandsonpresentation">Slideshare</a>)</li>
</ul>
<br />
※他にもちょこちょこ<br />
<br />
<br />
Rの初期設定やRSXファイルとは何かについては<a href="http://www.slideshare.net/nuimura/161104-foss4-gtokyoqgishandsonpresentation">Sliedshare</a>を参照していただくとして、ここでは気象庁で公開されている気象観測データをGISに読み込むRアルゴリズムを作ってみたので、RSXの自作の仕方の例として順を追って紹介していきます。<br />
<br />
<br />
<h3>
使用データ(気象庁の気象観測データ)について</h3>
<a href="http://www.data.jma.go.jp/obd/stats/etrn/index.php">過去の気象データ検索</a>から日本全国の気象台・アメダス、日時を指定することで該当の気象観測データのページが表示されます。(<a href="http://www.data.jma.go.jp/obd/stats/etrn/view/hourly_s1.php?prec_no=45&block_no=47648&year=2016&month=12&day=18&view=">例. 銚子の2016/12/18の1時間ごとのデータ</a>)<br />
<br />
似たような<a href="http://www.data.jma.go.jp/gmd/risk/obsdl/index.php#">過去の気象データ・ダウンロード</a>とは異なり、URLに規則性があったり、10分毎の観測値があったりとプログラミング上、使い勝手が良いので今回は過去の気象データ検索を使っています。<br />
<br />
<br />
<h3>
まずは普通にRスクリプトの作成
</h3>
<script src="https://gist.github.com/nuimura/57f188445747f6656bed2aa9e937564a.js"></script>
<br />
ざっくりとした処理の流れは、<br />
<ol>
<li>Webページのテーブル情報の読み取り</li>
<li>予め調べておいた観測所ごとの緯度経度の紐付け</li>
<li>SpatialPointsDataFrameクラスに変換</li>
<li>GISデータとして出力</li>
</ol>
となります。Step2の観測所ごとの緯度経度を調べたり、Webページのテーブルの列数が気象台とアメダス(アメダス内でも観測項目の数はいろいろ)で異なり、場合分け処理の記述がめんどいので、今回は関東地方の気象台のみの観測データに対応させてあります。<br />
<br />
<br />
<br />
<br />
<br />
<h3>
RSXファイル独自の記述</h3>
(詳細はQGIS Training Manual <a href="http://17.0.0.31/">17.31</a>, <a href="http://17.0.0.32/">17.32</a> を参照)<br />
<br />
<br />
<script src="https://gist.github.com/nuimura/31f24c8f6d8a78830438b52acd042e08.js"></script>
Rスクリプトの冒頭に以下の記述をしています。<br />
スクリプトでは指定した年月日時の観測値を読み取りGISdataに変換しているので、年月日時を入力するように設定しています。<br />
<br />
year の2016は初期入力値で変更可能です。<br />
month,day,hour は以下のような記述でドロップボックスから選択するようにすることができます。<br />
<br />
##Data downloading=group<br />
##year=number 2016<br />
##month=selection;1;2;3;4;5;6;7;8;9;10;11;12<br />
##day=selection;1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21;22;23;24;25;26;27;28;29;30;31<br />
##hour=selection;1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21;22;23;24<br />
##Output=output vector<br />
<br />
<br />
また、スクリプトの最後ではWriteOGR関数によるGISデータの出力はせずに、以下のように冒頭で指定したOutputという名前の変数(SpatialPointsDataFrameクラス)を作成して終わっています。<br />
<br />
Output <- SpatialPointsDataFrame(sp,data=sp.dfm)<br />
<br />
<br />
<h3>
出力結果</h3>
出力されるポイントデータは属性情報として、気象庁の観測データの各項目が入っています(気温、降水量、気圧や風速、風向など)。<br />
※風向は北を0とした時計回りの角度に変換してます。<br />
<br />
以下の図はOpenStreetMapの背景図に気温分布のラスターデータと観測地点のシンボルを矢印に変えて風向の値に応じて回転させたものが載せてあります。<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgKn6aTwsMVvLzqBCrW4hUWKXlIOtayfrP1Wrccw3aSwM_sw7Cmbvp7laD-LVBgg5gsXEGuxbAK1z35ByLQMndzX5FP8cbSLFLjJzXfsYUTVwS4M5eCS1LBYqhRjcskmPsUwGLKB_xTI0xv/s1600/qgis_output.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="434" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgKn6aTwsMVvLzqBCrW4hUWKXlIOtayfrP1Wrccw3aSwM_sw7Cmbvp7laD-LVBgg5gsXEGuxbAK1z35ByLQMndzX5FP8cbSLFLjJzXfsYUTVwS4M5eCS1LBYqhRjcskmPsUwGLKB_xTI0xv/s640/qgis_output.png" width="640" /></a></div>
<br />
<br />
とまあRでやっている処理をQGISでGUIで操作できるようにできるというのがprocessingのメリットです。<br />
<br />
ぶっちゃけRのプログラミングができる人はそちらで完結したほうが楽なのですが…。同じ処理をプログラミングのできない人とも共有できるのがpurocessingのメリットかなと思います。<br />
<br />
<br />
<br />
<br />
<br />
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-42827246513893210342016-09-02T15:48:00.001+09:002017-03-15T15:14:34.469+09:00GMT 5.2.1のインストール(Bash on Ubuntu on Windows 編)Bash on Ubuntu on Windows (長い…、以下BoUoW)の登場でWindows上でのUnix系の処理のやり方が大きく変わりつつあるのでGMTのインストール方法の情報を更新。<br />
(以前掲載した<a href="http://spatiohack.blogspot.jp/2012/03/gmt-457.html">GMT 4.5.7のインストールのしかた</a>をベースに作成しています。)<br />
<br />
<br />
<a href="http://gmt.soest.hawaii.edu/">GMT (Generic Mapping Tools)</a>とはコマンドラインベースの地図作成ツールで、地球科学分野でよく用いられるツールで、綺麗な空間分布図を描画することができます。クロスプラットフォームなソフトなのでWindows、Mac、Linuxで使うことができます。<br />
<span style="color: #cc0000; font-size: x-small;">追記:2016年9月2日現在GMT本体の最新バージョンは5.2.1です。海岸線データ(GSHHG)の最新バージョンは2.3.6です。後に出てくるインストールスクリプトのバージョン番号を自分が使用するバージョン番号に置き換えてください。</span><br />
<br />
今回は主にWindows 10でのGMTのインストール(&初期設定)について説明します。<br />
<br />
<hr />
<h4>
1. はじめに</h4>
Bash on Ubuntu on Windows のインストールは下記サイトを参照<br />
<a href="http://qiita.com/Aruneko/items/c79810b0b015bebf30bb">Bash on Ubuntu on Windowsをインストールしてみよう!</a><br />
<br />
作成したps/pdfファイルの表示は<a href="http://www.sumatrapdfreader.org/free-pdf-reader-ja.html">SumatraPDF</a>がおすすめ。<br />
Adobe Readerとかだとファイルを開くとロックが掛かり更新ができない。<br />
<br />
また、改行コードやエンコーディングが変更できるエディタ(<a href="https://notepad-plus-plus.org/">Notepad++</a>など)も必要。<br />
<br />
また、GMTインストールにあたって以下のパッケージが必要となりますので、apt-get でインストールしてください。<br />
<br />
build-essential<br />
make<br />
cmake<br />
netcdf-bin<br />
libnetcdf-dev<br />
gdal-bin<br />
libgdal-dev<br />
<br />
<div>
<hr />
<h4>
2. GMTのインストール</h4>
BoUoWにて、下記スクリプトを実行してください。</div>
<br />
<script src="https://gist-it.appspot.com/github/nuimura/gmt4BoUoW/raw/master/compile_gmt.sh"></script>
<br />
<br />
処理内容の説明<br />
<br />
Line 4-5:<br />
バージョン番号は最新のものに置き換えてください。<br />
<br />
Line 10-14:<br />
ファイルをダウンロードして展開。<br />
<br />
Line 19-22:<br />
海岸線データをインストールするための処理<br />
<br />
Line 25-30:<br />
インストール処理(時間かかるかも)<br />
<br />
<br />
<div>
<hr />
<h4>
3. GMTの動作確認</h4>
</div>
上記のGMTのインストールが終わったら、BoUoWで試しに<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><b>gmt pscoast -R0/360/-80/80 -Jm1:300000000 -Ba30df15dWeSn -A10000 -Dl -W0.5 -P > world.ps</b></span><br />
<div>
<br /></div>
と打ってみて下さい。<br />
<span style="color: #cc0000; font-size: x-small;">※表示の都合で2行にまたがっているかもしれませんが、1行で入力して下さい</span><br />
<br />
world.psというPSファイルがカレントディレクトリ※できているはずです。このPSファイルをWindowsで見るために、SumatraPDFやGsviewが必要となります。<br />
<br />
ちなみにPSファイルをPDFファイルに変換したい場合は<br />
<span style="font-family: "courier new" , "courier" , monospace;"><b>gmt </b><b>psconvert -A -Tf </b></span><b style="font-family: "courier new", courier, monospace;">world.ps</b><br />
<b style="font-family: "courier new", courier, monospace;"><br /></b>
ちなみにPSファイルをPNGファイルに変換したい場合は<br />
<span style="font-family: "courier new" , "courier" , monospace;"><b>gmt </b><b>psconvert -A -Tg </b></span><b style="font-family: "courier new", courier, monospace;">world.ps</b><br />
<br />
を実行するとそれぞれ変換できます。<br />
<br />
<br />
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh2csW341DSrXPEV83pZb9wYxxz9hTU6LkpXL9fdZv9qGpHd_qxW9IaHL5SCRWQQ-IlqvNp1fVaYFo2b3CmB6akaDbUeCwZe4Xx1W_cukpV-NeM1illq9d_hnWVWh5OvVh6TIhI4WQ0zHvf/s1600/world.png" imageanchor="1"><img border="0" height="243" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEh2csW341DSrXPEV83pZb9wYxxz9hTU6LkpXL9fdZv9qGpHd_qxW9IaHL5SCRWQQ-IlqvNp1fVaYFo2b3CmB6akaDbUeCwZe4Xx1W_cukpV-NeM1illq9d_hnWVWh5OvVh6TIhI4WQ0zHvf/s320/world.png" width="320" /></a><br />
pscoastコマンドで作成した図<br />
<br />
<br />
※BoUoWのホームディレクトリは以下の場所となります。<br />
(nuimuraをあなたのユーザー名に置き換えてください。)<br />
C:\Users\nuimura\AppData\Local\lxss\home\nuimura<br />
<br />
<br />
また、別のやり方としてはsources.list を強制的にxenialに書き換えて、16.04にアップグレードしてからソースからではなくapt-getからgmt5をインストールするやり方もあるそうですが、そのやり方は私もまだ試したことがありません。試される場合は自己責任でお願いします。<br />
<br />
<br />
※GMTコマンドや説明の誤り、16.04にアップグレードの情報などを @Say_no さんにいただき修正・追記いたしました。どうもありがとうございましたー。<br />
<br />
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-30729393761659965712016-02-03T14:06:00.001+09:002016-02-03T14:07:02.125+09:00SRTM GL1 (1秒メッシュ) の欠損値処理他大の学生さんから解析手法について問い合わせがあったのでついでにブログ投稿。<br />
<br />
<br />
地形データを利用する様々な分野の多くの研究者によって使われているSRTM DEMの前処理についてです。<br />
<br />
2000年2月のスペースシャトルエンデバーにより計測されたSRTM DEM は3秒メッシュ (約90 m) 解像度 (<a href="https://lpdaac.usgs.gov/dataset_discovery/measures/measures_products_table/srtmgl3_v003">SRTM3</a>) については当初から全世界のDEMが公開されていましたが、1秒メッシュ(約30 m) 解像度は米国内の領域しか公開されていませんでした。<br />
<br />
2014年には米国以外の領域の1秒メッシュ解像度DEMも公開され (<a href="https://lpdaac.usgs.gov/dataset_discovery/measures/measures_products_table/srtmgl1">SRTM GL1</a>) ダウンロードができるようになっていますが、利用者の利便性を考えてなのか、本来計測できていない部分の標高値も他のDEM (<a href="http://www.jspacesystems.or.jp/ersdac/GDEM/J/3.html">ASTER GDEM</a>、GMTED、NED) を使って補間してしまっています。<br />
<br />
手っ取り早く高解像度のDEMを使いたい人にとっては問題ないのですが、地形や氷河体積の時間変化などを明らかにしたい人にとっては、2000年時点 (スナップショット) での標高値ではない標高 (例えばASTER GDEM はグリッドによって2000年~2011年と異なる) が混ざるため困ります。<br />
<br />
そこで標高の時間変化計算に用いるためにはSRTM GL1から外来DEMのグリッドを除去する処理が不可欠となります。SRTM GL1には対応するタイルごとに<a href="https://lpdaac.usgs.gov/dataset_discovery/measures/measures_products_table/srtmgl1n">SRTMGL1N</a>というフラグデータがあり、SRTM GL1のグリッドごとに、外来DEMを使ったかどうか使った場合はどの種類かがわかる情報があります。この情報を元にオリジナル のSRTMを使っているグリッドのみを抽出すればOKです。<br />
<br />
<h3>
QGISでの読み込み</h3>
以上の処理は<a href="https://cran.r-project.org/">R</a>やMATLABなどのプログラミングでももちろんできますが、ここでは<a href="http://qgis.org/">QGIS</a>を使ってやり方を紹介します。<br />
<br />
<a href="https://lpdaac.usgs.gov/dataset_discovery/measures/measures_products_table/srtmgl1n_v003">SRTM GL1Nの仕様のページ</a>を見ると、値は8 bit (0-255) で格納されており、201以上の値を使えば大丈夫そうです。SRTM GL1のフォーマット (*.hgt) は単純な raw binary形式ですが、QGISはSRTMのファイル名(ex. NaaEbbb.hgtの場合は北緯aa度東経bbb度を北西端)から自動で認識して正しい位置に表示してくれます。<br />
<br />
SRTM GL1Nも単純な raw binary 形式となっていますがこちらはQGISが自動認識はしてくれないためヘッダーファイル (*.hdr) を用意して位置情報などのデータに関する上をQGISに教えてやる必要があります。<br />
<br />
<h3>
<br />hdrファイルとは</h3>
テキスト形式でデータについての情報を記述したものです。 raw binaryファイルと拡張子より前の部分の名前を揃える必要があります (ex. データがN27E086.numの場合はN27E086.hdr)。<br />
<br />
詳しくはこちらのページのEHdr -- ESRI .hdr labelled を参照<br />
<a href="http://www.gdal.org/frmt_various.html">http://www.gdal.org/frmt_various.html</a><br />
<br />
<br />
N27E086.num の場合は、N27E086.hdrという名前で以下の内容が記述してあるテキストファイルを用意すればOKです。<br />
<br />
<br />
BYTEORDER I<br />
LAYOUT BIL<br />
NROWS 3601<br />
NCOLS 3601<br />
NBANDS 1<br />
NBITS 8<br />
BANDROWBYTES 7202<br />
TOTALROWBYTES 7202<br />
PIXELTYPE UNSIGNEDINT<br />
ULXMAP 86<br />
ULYMAP 28<br />
XDIM 0.000277777777777778<br />
YDIM 0.000277777777777778<br />
NODATA 0<br />
<br />
<br />
一応こちらからダウンロードできます。適当なエディタで開いてみてください。<br />
<a href="https://dl.dropboxusercontent.com/u/870568/blogger/N27E086.hdr">https://dl.dropboxusercontent.com/u/870568/blogger/N27E086.hdr</a><br />
<br />
<br />
準備ができたら N27E086.num をQGISにドラッグ&ドロップすれば表示されるはずです。あとはQGISのラスターメニューの変換コマンド (または右クリックで別名保存) でGeoTIFF形式に変換すると他のソフトウェアでも簡単に読み込むことができます。<br />
SRTM GL1のファイル名と同じだとややこしいので N27E086_num.tifとして保存します。<br />
<br />
<br />
<h3>
SRTM GL1からオリジナルのSRTMのみの抽出</h3>
QGISにSRTM GL1 (N27E086.hgt) とSRTM GL1N (N27E086_num.tif) の2つを読み込んで置いて下さい。<br />
<br />
まずラスタ計算機 (Raster calculater) で以下の内容を図のようにウィンドウ下部の演算欄に入力します。<br />
<br />
"N27E086@1" * ("N27E086_num@1" >= 201)<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiSob4MZehfiqzdzAyii_m2_qYE3HQa1orYYzv6QwUSm0-Mf5BqSl8xQZLSW-wJ32gRTXdE_YuDyXWnzF6aonIfdBdx5-w6kG_rcGyPdl5jlFIonN3naKsmXpUMrJrOF_yIJICf3M-vymlY/s1600/raster_calc.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="346" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiSob4MZehfiqzdzAyii_m2_qYE3HQa1orYYzv6QwUSm0-Mf5BqSl8xQZLSW-wJ32gRTXdE_YuDyXWnzF6aonIfdBdx5-w6kG_rcGyPdl5jlFIonN3naKsmXpUMrJrOF_yIJICf3M-vymlY/s400/raster_calc.png" width="400" /></a></div>
<br />
<br />
フラグデータで201以上の値が入っている場所のSRTM GL1のみを抽出するという意味です。ウィンドウ右上の出力ファイル名も指定しておいて下さい。<br />
出力ファイル名を N27E086_temp.tif とします。<br />
<br />
変換して出力されたファイルは以下の図のようにSRTMじゃない部分 (つまり欠損値) に数値の0が入っています。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhzEwV26BNO9RmC0GnIYfc5NWvADqY4WOUEFGGBfREq5IippXVg5LXSKOj9KNCthTaILVeFTTdirBK0oywi6pfAG0RumNfD9ABzgKUzA3VZ0Vqt-HZfDS1F0e5nfnaQM0kyg_umaTDLPDD1/s1600/result.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="300" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhzEwV26BNO9RmC0GnIYfc5NWvADqY4WOUEFGGBfREq5IippXVg5LXSKOj9KNCthTaILVeFTTdirBK0oywi6pfAG0RumNfD9ABzgKUzA3VZ0Vqt-HZfDS1F0e5nfnaQM0kyg_umaTDLPDD1/s400/result.png" width="400" /></a></div>
<br />
<br />
<br />
標高変化の計算のためには0じゃなくてNo Dataだと認識されて欲しいのでもうひと手間かけます。<br />
<br />
N27E086_temp.tif を右クリックで別名保存で、ウィンドウ下部のNo data valueで0を指定します。出力ファイル名を指定してOKをクリックします。<br />
ここでは出力ファイル名をN27E086_void.tif とします。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjGGsEdSxcereiL9DnID_p_pRqiX6_irjz_n-ksJQIGKvEjaw5albCTTfBvNw6Ff9kE4n_9PhxNr694Q3OIk_pYwgOQ0LWRZrCUumvNDvz-cYOQdrujqPp63j0SPF5YNl9tJFGUfhJu2pAT/s1600/saveas.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="400" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjGGsEdSxcereiL9DnID_p_pRqiX6_irjz_n-ksJQIGKvEjaw5albCTTfBvNw6Ff9kE4n_9PhxNr694Q3OIk_pYwgOQ0LWRZrCUumvNDvz-cYOQdrujqPp63j0SPF5YNl9tJFGUfhJu2pAT/s400/saveas.png" width="397" /></a></div>
<br />
<br />
<br />
以下のような感じでできあがりです。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgQcghYL1FMQ2nrDkr7lovdkik6xA_H9haDlfjXl2Uh0tr0mtjKUxtz2s0xmBJaRr5nvGqTv853v1Rgc2xOv0VLYQN2r-YKdCnDAQ25sVBBlbhrXqwvsfnVkD1IkyPZUfR9H3A0mo9KOSUl/s1600/result2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="300" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgQcghYL1FMQ2nrDkr7lovdkik6xA_H9haDlfjXl2Uh0tr0mtjKUxtz2s0xmBJaRr5nvGqTv853v1Rgc2xOv0VLYQN2r-YKdCnDAQ25sVBBlbhrXqwvsfnVkD1IkyPZUfR9H3A0mo9KOSUl/s400/result2.png" width="400" /></a></div>
<div>
<br /></div>
Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-37931364701324999592016-02-02T16:49:00.001+09:002017-03-15T15:16:51.129+09:00R言語(ggplot2)による気温変動グラフの作成毎年(といっても去年から)後期授業も一段落した2月ぐらいに、ゼミの学生たち+α対象に勉強会を開いています。<br />
<div>
<br /></div>
<div>
去年は<a href="http://env-crisis-cis.blogspot.jp/2015/02/gis.html">GISの勉強会</a>と<a href="http://env-crisis-cis.blogspot.jp/2015/02/blog-post_25.html">統計(Excel)の勉強会</a>を行いましたが、今年はそれに加えてプログラミング(R言語)の勉強会も追加しようと考えています。</div>
<div>
<br /></div>
<div>
でいろいろネタ探しをしていたら、気象庁のWebページで以下のような報道発表があり、グラフに使われたデータもありました(CSVなどではなくWebページ上の表ですが)。</div>
<div>
<br /></div>
<div>
<a href="http://www.jma.go.jp/jma/press/1602/01b/worldtemp2015_0201.html">世界の年平均気温がこれまでの最高値を更新 ~2015年(平成27年)の世界と日本の年平均気温(確定値)~</a></div>
<div>
<br />
<br /></div>
<div>
</div>
報道発表のリンク先にある<a href="http://www.data.jma.go.jp/cpdinfo/temp/an_wld.html">世界の年平均気温の偏差のページ</a>にあるような図をちょっとRのggplot2ライブラリで作成してみました。<br />
<div style="-webkit-text-stroke-width: 0px; color: black; font-family: メイリオ; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: 1; word-spacing: 0px;">
<div style="margin: 0px;">
<br />
<br />
<br />
<br />
<br /></div>
</div>
<script src="https://gist-it.appspot.com/github/nuimura/globalTemperature/raw/master/globalTemperature.R"></script>
以下のような図が出来上がります。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEg4-YCm-CPDPH9ASX5iFNU5hRPd7v8epunsQvaVbTv4-jAl9yVD3cfCG6bOzAENn00N_PY0QcDFezdiJIFsKYWLfOpjfNy61aOI4cXbkc9ZsKwPuCA8omXhwSKH796TjQhwVgQZhxTDQjJS/s1600/all_detail.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="250" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEg4-YCm-CPDPH9ASX5iFNU5hRPd7v8epunsQvaVbTv4-jAl9yVD3cfCG6bOzAENn00N_PY0QcDFezdiJIFsKYWLfOpjfNy61aOI4cXbkc9ZsKwPuCA8omXhwSKH796TjQhwVgQZhxTDQjJS/s400/all_detail.png" width="400" /></a></div>
<div style="text-align: center;">
気象庁のページに似せて作った図。</div>
<div style="text-align: center;">
世界全体での年平均値(黒丸)、回帰直線(青)および5年移動平均(赤)</div>
<div style="text-align: center;">
<br /></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgMa3zQNJMcjM5dU7kijLlLFqt_XxnCovQFomzJY_SxQPLVpmjp7x4bZn1euRYTYBtMt5coCWWX5wQDoVvGI3vwGM895k1enN2GkQVtq6w0-JoT-PbXcBmFoxshymr9CtkMlvHFhgGI2hMR/s1600/compare_mv.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="250" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgMa3zQNJMcjM5dU7kijLlLFqt_XxnCovQFomzJY_SxQPLVpmjp7x4bZn1euRYTYBtMt5coCWWX5wQDoVvGI3vwGM895k1enN2GkQVtq6w0-JoT-PbXcBmFoxshymr9CtkMlvHFhgGI2hMR/s400/compare_mv.png" width="400" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
世界全体(赤)、北半球(緑)および南半球(青)の5年移動平均</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
Githubリポジトリはこちら</div>
<div class="separator" style="clear: both; text-align: left;">
<a href="https://github.com/nuimura/globalTemperature">https://github.com/nuimura/globalTemperature</a></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-67691145115466852342015-12-02T16:44:00.000+09:002015-12-03T12:40:54.065+09:00Earthの風データをFOSS4Gで可視化<a href="http://qiita.com/advent-calendar/2015/foss4g">FOSS4G Advent Calendar 2015</a>の12/2の記事です。<br />
2013年、2014年は参加してなかったので<a href="http://spatiohack.blogspot.jp/2012/12/agu2012-fall-meetingosgeo-workshop.html">2012年</a>以来となります。<br />
<br />
<br />
ご存じの方も多いと思うのですが、気象データの美しい可視化サイトで有名なEarthというサイトがあります。<br />
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjVdudmFfeQIxuI7_nxUoFdE9FfsIlHawWNweAgSlx4TL8reZE9co4TDaoqYP_uyQawXijPCE1obLvR4r1HFcQaICisgLho3MOflkfsCzAsTdk1_o54fQz34jI_A2_Vok0sR7AwLZsfrWhG/s1600/earth_web.png" imageanchor="1"><img border="0" height="252" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjVdudmFfeQIxuI7_nxUoFdE9FfsIlHawWNweAgSlx4TL8reZE9co4TDaoqYP_uyQawXijPCE1obLvR4r1HFcQaICisgLho3MOflkfsCzAsTdk1_o54fQz34jI_A2_Vok0sR7AwLZsfrWhG/s320/earth_web.png" width="320" /></a></div>
<div style="text-align: center;">
<a href="http://earth.nullschool.net/jp/">http://earth.nullschool.net/jp/</a></div>
<div style="text-align: center;">
<br /></div>
<br />
こちらはCameron Beccarioさんという方がNOAAやNASAなどの研究機関が公開している気象データからD3などのツールを使って可視化を行っており、大気や海洋の様子をわかりやすいアニメーションで見ることができます。<br />
<br />
もちろんこのまま表示をいじっていろいろ見るだけでも十分楽しいのですが、研究者的には元データってどんなデータなんだろう?とか、どこからダウンロードして手元でいじれるのかなどと考えてしまうわけです。<br />
<br />
そこでこの記事では、EARTHで使われている気象データのうち、風データのみに着目してデータ配布サイトから入手、フォーマット変換をしてからQGISで可視化するまでを紹介したいと思います。<br />
<br />
<br />
<h3>
データの入手</h3>
一応<a href="http://earth.nullschool.net/jp/about.html">EARTHの説明ページ</a>を見るとデータの種類ごとに、元データの配布機関へのリンクは張っています。風データの場合はNational Centers for Environmental Prediction (NCEP) により提供されている気象予測モデルGlobal Forecast System (GFS)の気象データを使っているようです。<br />
<br />
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgzHVtfu2iqONC8e02PuRsbttQme9D72hPqQQz69GaJ4A-mllVskghoCfQY66R0Ck6BgTjxu2yJYqRNL9JI8cfwp-vh55xo3doZStBbpc6eVsyT1opqA1MDl5rznxt1pMm6m23irsB1jmH_/s1600/gfs_web.png" imageanchor="1"><img border="0" height="196" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgzHVtfu2iqONC8e02PuRsbttQme9D72hPqQQz69GaJ4A-mllVskghoCfQY66R0Ck6BgTjxu2yJYqRNL9JI8cfwp-vh55xo3doZStBbpc6eVsyT1opqA1MDl5rznxt1pMm6m23irsB1jmH_/s320/gfs_web.png" width="320" /></a></div>
<div style="text-align: center;">
<a href="https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs">GFSのサイト</a></div>
<br />
ぱっと見ただけではそれらしいデータへのリンクがたくさんあってよくわかりませんorz.<br />
そこで、<a href="https://github.com/cambecc/earth">githubに公開されているEARTHのソースコード</a>を確認してみます。<br />
<br />
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiTIjGjUSiQa-3-p4Etkn25phEBfQCgo0Yl_84b3lR46gwhWd9zlqLjmBZ8By5VB481fKjxd7njwAmxL8eVapjnX_tDJOhReIF9f-hCj1MUlEEtEnku-YCh34Jo9wvA9OGOfzr2cUvRIE9O/s1600/earth_github1.png" imageanchor="1"><img border="0" height="240" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiTIjGjUSiQa-3-p4Etkn25phEBfQCgo0Yl_84b3lR46gwhWd9zlqLjmBZ8By5VB481fKjxd7njwAmxL8eVapjnX_tDJOhReIF9f-hCj1MUlEEtEnku-YCh34Jo9wvA9OGOfzr2cUvRIE9O/s320/earth_github1.png" width="320" /></a></div>
<div style="text-align: center;">
<a href="https://github.com/cambecc/earth">https://github.com/cambecc/earth</a></div>
<div style="text-align: center;">
<br /></div>
<div style="text-align: left;">
ページの下の方に行くと <a href="https://github.com/cambecc/earth#getting-weather-data">"getting weather data"</a> というそれらしい説明部分がありました。<a href="http://nomads.ncep.noaa.gov/">NOAA Operational Model Archive and Distribution System(NOMADS)</a>というデータ配布システムからダウンロードできるみたいです。</div>
<div style="text-align: left;">
※そういえばRにも <a href="https://cran.r-project.org/web/packages/rNOMADS/index.html">rNOMADS</a> というこのサイトのデータへアクセスできるライブラリがありますね。</div>
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjV_s4og_0p66p626qxBQQpoq6DZvFDeEhVXUBVqJGhYXG908ZOZnXpm6nciExw64ElhAGlbTHbuCtG4Y-WH7oldzHAy-os1YJvrsU8BlRiigptMXfvh7eBKJVGJ4JNDIFfpd_kLuwdmWfI/s1600/earth_github2.png" imageanchor="1"><img border="0" height="167" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjV_s4og_0p66p626qxBQQpoq6DZvFDeEhVXUBVqJGhYXG908ZOZnXpm6nciExw64ElhAGlbTHbuCtG4Y-WH7oldzHAy-os1YJvrsU8BlRiigptMXfvh7eBKJVGJ4JNDIFfpd_kLuwdmWfI/s320/earth_github2.png" width="320" /></a></div>
<br />
でNOMADSのページから "GFS 0.25 Degree"の <a href="http://nomads.ncep.noaa.gov/cgi-bin/filter_gfs_0p25.pl">"grib filter"というリンク</a>をたどり。<br />
ダウンロードしたい日時(GMT)を選びます。<a href="http://nomads.ncep.noaa.gov/cgi-bin/filter_gfs_0p25.pl">ここでは2015/12/01 00時のデータ</a>を使いました。<br />
(↑だんだんWebページのスクリーンショットを貼るのが面倒になってきましたので文字ベースで)<br />
<br />
<br />
GRIB Filterのページでは以下の指定を行います。<br />
<br />
<ul>
<li>ドロップダウンリストで "gfs.t00z.pgrb2.0p25.f000" (fのあとはその日時から何時間後の予測データなのかを示す。f003だと3時間後)</li>
<li>levelとして "10 m above ground" にチェック</li>
<li>variable(変数)として "UGRD" と "VGRD" を選択してダウンロードを開始</li>
</ul>
<br />
<div>
※とりあえず地上の風(1000 hPa)を見たいので10 mの高さで、風の東西成分(UGRD)、南北成分(VGRD)を指定しています。</div>
<div>
<br /></div>
<div>
ダウンロードされたデータは、 "gfs.t00z.pgrb2.0p25.f000" という名前のGRIB2形式のデータです。GRIB2形式はGDALでも対応しているのでQGISで開くことができます。</div>
<div>
<br /></div>
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEg125efmEmi43XkZ44tQioRLx5o4kQ7cIlfCBFxRx3KHDgLyXaS7NX6NvfMabJooDpDa00Wwti_HJFUt3f38dXEm-9ZF7avqU0fSHcqH6G3g4wC53pw8lg_jO-t86pjbzxvraFDfScuJwQR/s1600/qgis1.png" imageanchor="1"><img border="0" height="191" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEg125efmEmi43XkZ44tQioRLx5o4kQ7cIlfCBFxRx3KHDgLyXaS7NX6NvfMabJooDpDa00Wwti_HJFUt3f38dXEm-9ZF7avqU0fSHcqH6G3g4wC53pw8lg_jO-t86pjbzxvraFDfScuJwQR/s320/qgis1.png" width="320" /></a></div>
<div style="text-align: center;">
Band1のU(東西成分)が赤色、Band2のV(南北成分)が緑色で表示されているみたいですが、このままではよくわかんないですね…。</div>
<div style="text-align: center;">
<br /></div>
<div>
<br /></div>
<h3>
データの変換</h3>
この2バンドのラスターGRIB2を以下の2種類のデータに変換します。<br />
<br />
<ol>
<li>風向・風速を示すベクターデータ</li>
<li>風速を示すラスターデータ</li>
</ol>
<br />
<br />
<h4>
風向・風速を示すベクターデータへの変換</h4>
<div>
WindowsにQGIS 2.12.1(スタンドアロン版)がインストールされているとして話を進めていきます。<br />
<br />
先ほどのGRIB2データがDドライブの直下に置いてあるとして、MSYS (QGISと一緒にインストールされるUNIXライクなコンソール)を起動して以下のコマンド入力していきます。#で始まる行はコマンドじゃなく説明です。<br />
<span style="font-family: "courier new" , "courier" , monospace;"><b><br /></b></span>
<span style="font-family: "courier new" , "courier" , monospace;"><b> cd d</b></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><b><br /></b></span>
<span style="font-family: inherit;">#gdalinfo コマンドでGRIB2データの情報を確認</span><br />
<span style="font-family: inherit;"><br /></span>
<span style="font-family: "courier new" , "courier" , monospace;"><b> gdalinfo gfs.t00z.pgrb2.0p25.f000</b></span><br />
<span style="font-family: inherit;"><br /></span>
#GRIB2形式のband1(UGRD)とband2(VGRD)をそれぞれCSVに変換<br />
<span style="font-family: inherit;"><br /></span>
<span style="font-family: "courier new" , "courier" , monospace;"><b> gdal2xyz.py -band 1 -csv gfs.t00z.pgrb2.0p25.f000 wind_x.csv</b></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><b> gdal2xyz.py -band 2 -csv gfs.t00z.pgrb2.0p25.f000 wind_y.csv</b></span><br />
<div>
<span style="font-family: "courier new" , "courier" , monospace;"><b><br /></b></span></div>
#改行コードをUNIX形式に変換<br />
<br />
<b style="font-family: 'Courier New', Courier, monospace;"> </b><span style="font-family: "courier new" , "courier" , monospace;"><b>dos2unix wind_x.csv</b></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><b> dos2unix wind_y.csv</b></span><br />
<div>
<br /></div>
#1つのCSVファイルに結合<br />
<br />
<b style="font-family: 'Courier New', Courier, monospace;"> </b><span style="font-family: "courier new" , "courier" , monospace;"><b>echo "lon,lat,u,v" > wind.csv</b></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><b> paste -d"," wind_x.csv wind_y.csv | awk -F, '{print $1","$2","$3","$6}' >> wind.csv</b></span><br />
<div>
<br /></div>
<div>
これで以下のように、longitude、latitude、U成分、V成分のCSVファイルが出来上がります。</div>
<div>
<div>
lon,lat,u,v</div>
<div>
0.000,90.000,-2.79,-1.51</div>
<div>
0.250,90.000,-2.8,-1.5</div>
<div>
0.500,90.000,-2.8,-1.48</div>
</div>
<br />
あとはQGISでCSVファイルを読み込み、Shapefile形式で別名で保存を行って下さい。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhUnyiPTsSVwj5qlMradwzqUMLHz_9EaqvdTnanPIBSkLebH_2ns5RM6m6AEAswMiT9JoO4Q3MRGkPwbdAT0gaNqNMBeLH-5VV_WJF3FS8hYaCB3xQKVv8jqMwbAgFCZJ2-IRbivn-i5Sgi/s1600/qgis2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="242" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhUnyiPTsSVwj5qlMradwzqUMLHz_9EaqvdTnanPIBSkLebH_2ns5RM6m6AEAswMiT9JoO4Q3MRGkPwbdAT0gaNqNMBeLH-5VV_WJF3FS8hYaCB3xQKVv8jqMwbAgFCZJ2-IRbivn-i5Sgi/s320/qgis2.png" width="320" /></a></div>
<div style="text-align: center;">
CSVファイル読み込みウィンドウ</div>
<div style="text-align: center;">
X fieldにlon、Y fieldにlatを指定</div>
<br />
<h4>
風速を示すラスターデータへの変換</h4>
#再びMSYS上で以下のようにコマンドを入力して、GRIB2形式のファイルを風速のラスターデータ(GeoTIFF)に変換します。<br />
<b style="font-family: 'Courier New', Courier, monospace;"> </b><span style="font-family: "courier new" , "courier" , monospace;"><b>gdal_calc.py -A gfs.t00z.pgrb2.0p25.f000 --A_band=1 -B gfs.t00z.pgrb2.0p25.f000 --B_band=2 --calc="sqrt(A*A+B*B)" --outfile=scalar.tif</b></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><b><br /></b></span>
<span style="font-family: "courier new" , "courier" , monospace;"><b><br /></b></span>
<span style="font-family: "courier new" , "courier" , monospace;">変換したデータを読み込んで色を設定すると以下の様になります。</span><br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhH0GuxU4Y-EkxBTfdOwIRRym8K-asKyzNwuQOnp7Xda-ET9VbI2aTv_YuyvghDrWdnBw4uHVsTKAFgFQZ9uI0Gt9v_RNkEFErMBBIckMM6BpJeb5pLHgk58IQ2TUcEv7qhQlsVz_6PairG/s1600/qgis3.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em; text-align: center;"><img border="0" height="183" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhH0GuxU4Y-EkxBTfdOwIRRym8K-asKyzNwuQOnp7Xda-ET9VbI2aTv_YuyvghDrWdnBw4uHVsTKAFgFQZ9uI0Gt9v_RNkEFErMBBIckMM6BpJeb5pLHgk58IQ2TUcEv7qhQlsVz_6PairG/s320/qgis3.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
風速が強いところが赤、弱いところが青</div>
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><b><br /></b></span>
<br />
<h3>
QGISでの操作</h3>
<div class="separator" style="clear: both; text-align: left;">
風のベクトルプロットをするにあたって追加でプラグインが必要になります。</div>
<div class="separator" style="clear: both; text-align: left;">
QGISのプラグインメニューより "Vector field renderer" プラグインをインストールして下さい。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
その後、風速のラスターデータの上に、風向・風速のベクターデータを読み込んで重ねて下さい。ポイントの数が多いので、読み込んだ直後は以下の図のように真っ黒に見えます。日本周辺などに拡大してみましょう。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhcXk49eiRsrl1lp-MTTsPiE-iaEXdF6QTpCVeWVU7TmcfuEvxlpEsHmV_3z7ag-LzISSyJohVLRPsASbDyYJ5C2LB9sQstQmvyYUM3iLVx3IcbZHCZ2ewJZI_hR_xps0gpTE69Ucd2HIQL/s1600/qgis4.png" imageanchor="1"><img border="0" height="184" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhcXk49eiRsrl1lp-MTTsPiE-iaEXdF6QTpCVeWVU7TmcfuEvxlpEsHmV_3z7ag-LzISSyJohVLRPsASbDyYJ5C2LB9sQstQmvyYUM3iLVx3IcbZHCZ2ewJZI_hR_xps0gpTE69Ucd2HIQL/s320/qgis4.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
データ全体表示</div>
<div class="separator" style="clear: both; text-align: center;">
よくわからない…。</div>
<br />
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjs0K4bWwKJer9BxN2GuCSUGIgfC0zOgPVfBg-wjcFwbKEOgXnH6HdNS4YqNK_oDTAYB4EHFCxBJazLUAFbH7R_ljNCakdRtfoXtXkj_OUg9UwAWacBvyqwM6zH3XdFMs0wQoU0p9JxyIjA/s1600/qgis5.png" imageanchor="1"><img border="0" height="184" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjs0K4bWwKJer9BxN2GuCSUGIgfC0zOgPVfBg-wjcFwbKEOgXnH6HdNS4YqNK_oDTAYB4EHFCxBJazLUAFbH7R_ljNCakdRtfoXtXkj_OUg9UwAWacBvyqwM6zH3XdFMs0wQoU0p9JxyIjA/s320/qgis5.png" width="320" /></a></div>
<div style="text-align: center;">
日本の北海道あたりに拡大するとポイントデータが格子状に並んでいるのがわかります。</div>
<br />
<br />
このそれぞれのポイントデータが風の東西成分と南北成分を持っていいるのでその値を元にそれぞれのポイントにおけるベクトルを表示していきます。<br />
<br />
ベクトルデータのプロパティをみると、スタイルタブで "ベクターフィールド" という項目が選べるようになっています。<br />
そこで以下のように "X attribute" に "u"、 "Y attribute" に"v" を指定して、"Base size" を "0" に指定してOKをクリックします。<br />
<br />
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgDZJ27tUQrpX-csS5_d8hKtc14K3b58g_3HPUKiRAgvBZRyABZKXG6yZvJlLgyAuTj_V-x9JqZfALHdaj_pGdU5XoHRD2c_I98RQ90ff0zBrwdm8o6QzKCWrZ4UmIH8ETsE9VRxKRhGGJf/s1600/qgis6.png" imageanchor="1"><img border="0" height="268" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgDZJ27tUQrpX-csS5_d8hKtc14K3b58g_3HPUKiRAgvBZRyABZKXG6yZvJlLgyAuTj_V-x9JqZfALHdaj_pGdU5XoHRD2c_I98RQ90ff0zBrwdm8o6QzKCWrZ4UmIH8ETsE9VRxKRhGGJf/s320/qgis6.png" width="320" /></a></div>
<div style="text-align: center;">
”Vector field” の設定</div>
<br />
<br />
すると最終的に以下のように風向・風速の分布を表示することができます。<br />
<div style="text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgJ7qdOl_DALBbtMNq9SWfbc8vMSrqyiGpkllpWU0iZrvUS7yb1H_CWwDQgvP2sgS2uuyo086i3Dyt5_xb7KXyPyWGLERgnWKXjlzr3C0QuXfCQb7-ZOLZfL-13ede-af4LcNSc229nORPk/s1600/qgis7.png" imageanchor="1"><img border="0" height="230" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgJ7qdOl_DALBbtMNq9SWfbc8vMSrqyiGpkllpWU0iZrvUS7yb1H_CWwDQgvP2sgS2uuyo086i3Dyt5_xb7KXyPyWGLERgnWKXjlzr3C0QuXfCQb7-ZOLZfL-13ede-af4LcNSc229nORPk/s400/qgis7.png" width="400" /></a></div>
<br />
<br />
まあ実際解析する際はこのようなスナップショットではなく、一定期間ごとに統計処理を行って解析するのですが、そこまでやる場合はプログラミングでの処理が必要になってきますね。そこら辺は追々Rを使った処理についてまたブログに書きたいと思います。<br />
<br />
ちょっと学生への手順書としてもかねての記事なので説明が冗長な部分もありますが…。<br />
こんなところで終わります。<br />
<br /></div>
Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-54331855016083159952015-06-03T10:42:00.000+09:002015-06-03T10:42:04.714+09:00GRASS GIS 7 を使ったShapefileのエラーチェック&修正GRASS GISではデータを処理する際、一般的に多く使われているShapefileやGeoTIFFなどのファイル形式からGRASS GIS独自のファイル形式に変換をします。<br />
<br />
ベクターデータの変換の際にはトポロジーなどのデータのエラーを自動的に検知して修正してくれるため、GRASS GISにベクターデータをインポートして、そのままエクスポートするだけで、エラーのないきれいなShapefileを作成することができます。<br />
<br />
デジタイズなどでShapefileを編集して作成する場合、途中切り取ったり穴あきポリゴンにしたり、マルチポリゴンにしたりなどの作業をしていくうちにエラーの含まれたShapefileとなってしまい、様々なソフトウェアでジオメトリを計算する際に正しい値が算出されない場合があります(ドーナツポリゴンの面積が正しくないなど)。<br />
そのためにもShapefileのエラーチェックは必要となってきます。<br />
<br />
以下が作業手順です。<br />
GRASS GIS 7で"gcs"という名前の地理座標系のLocationを作成し、"nuimura"という名前のマップセットを作成したとして説明していきます。<br />
マップセットの作成などの手順は<a href="http://spatiohack.blogspot.jp/2015/06/grass-gis-7-bat.html">以前の投稿</a>を参考にして下さい。<br />
<br />
<h3>
ベクターデータ(Shapefile) のインポート</h3>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhmSdn6lOieD1v9h-AV41R4HHRfgnAKLIYYo2wAIuZA3v2D-UKP6EOYh57ryc97zN6Zb42Fj38qkKvmLZbl_HugKLMF9tF7qiXYfhXxZjyB3zRU4LFQ_jIduS_OK8VwVW6uUOaeM8N0Qjcg/s1600/aa.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="256" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhmSdn6lOieD1v9h-AV41R4HHRfgnAKLIYYo2wAIuZA3v2D-UKP6EOYh57ryc97zN6Zb42Fj38qkKvmLZbl_HugKLMF9tF7qiXYfhXxZjyB3zRU4LFQ_jIduS_OK8VwVW6uUOaeM8N0Qjcg/s320/aa.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: left;">
こちらの図のように、「File」-> 「Import vector data」 -> 「Common import formats [v.in.ogr]」を選択します。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgbosESYvdWVNoN0PSU-RRGO2JU95LKB8XR2UWq5AglcGedP73zewEz5mkQpVG_GIAPDmugz0oZbCAetoxNGZNCVimkCT0M9CVA-rwk_JEAreuB3OFXN1AHnbe9zPdBWJgpwCJEYKe09xD2/s1600/bb.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="320" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgbosESYvdWVNoN0PSU-RRGO2JU95LKB8XR2UWq5AglcGedP73zewEz5mkQpVG_GIAPDmugz0oZbCAetoxNGZNCVimkCT0M9CVA-rwk_JEAreuB3OFXN1AHnbe9zPdBWJgpwCJEYKe09xD2/s320/bb.png" width="307" /></a></div>
<div class="separator" style="clear: both; text-align: left;">
「Source settings」で「Browse」から読み込みたいShapefileを指定して、「Import」をクリックします。(こちらの例では、Dドライブ直下にあるaaa.shp)</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiRyBpPi1mgAI9XkW5MJMSjCPDcy_QUeajtHkmeeoSpGh0ex6EGgMRPrBGF5juCWm6o8guadAyh_FwapavIVyymnHj-3o1PjRTFd3BsW_mo_2jszlVLCCuc6nLIRLQPlfyV2VivIYCwFnTm/s1600/ccc.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="268" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiRyBpPi1mgAI9XkW5MJMSjCPDcy_QUeajtHkmeeoSpGh0ex6EGgMRPrBGF5juCWm6o8guadAyh_FwapavIVyymnHj-3o1PjRTFd3BsW_mo_2jszlVLCCuc6nLIRLQPlfyV2VivIYCwFnTm/s320/ccc.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: left;">
その後、「File」-> 「Export vector map」 -> 「Common export formats [v.out.ogr]」を選択します。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhRPyTKeRMKvVNgv2mn9-lGrOIR7R8gK5Wed6t9cHwHLiCii8fpqtdNshnNB_1RXDxAEdqs5vqZUo01SIqI1x7su9RRK0r_7VvITx4oMzW9G8RkwWEn4aeFvYK9y4meujSMXfcbS0O38YyK/s1600/ddd.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="240" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhRPyTKeRMKvVNgv2mn9-lGrOIR7R8gK5Wed6t9cHwHLiCii8fpqtdNshnNB_1RXDxAEdqs5vqZUo01SIqI1x7su9RRK0r_7VvITx4oMzW9G8RkwWEn4aeFvYK9y4meujSMXfcbS0O38YyK/s320/ddd.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: left;">
「Name of input vector map to export」にファイル名@マップセット名を選び、「Name of output OGR datasource」に出力先のパスを指定して「Run」をクリックします。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
これでエラーチェックしたきれいなShapefileは出力出来ました。</div>
<div class="separator" style="clear: both; text-align: left;">
最後にGRASS GIS内に残ったファイルを削除して終了です。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEihjbPEPUaub8OlHMMRm9lbtB1bghTIUFpwIvckH64-4ROM3qM6vVm6LeNmg7xszx1fIllgo7iYPF_ZDi85qVlLM-hbrVvexCXpKmabFKUQGrB3jgRpofbkNiEQHJaw9aiAgS6qmk8M6PJt/s1600/eeee.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="320" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEihjbPEPUaub8OlHMMRm9lbtB1bghTIUFpwIvckH64-4ROM3qM6vVm6LeNmg7xszx1fIllgo7iYPF_ZDi85qVlLM-hbrVvexCXpKmabFKUQGrB3jgRpofbkNiEQHJaw9aiAgS6qmk8M6PJt/s320/eeee.png" width="276" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
<span style="text-align: left;">「File」-> 「Manage maps」 -> 「Delete [g.remove]」を選択します。</span></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
「Data type」にvector mapを選択して「Force removal」にチェックを入れて</div>
<div class="separator" style="clear: both; text-align: left;">
消したいファイル名を選択してRunをクリックすると削除されます。</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-6021702093490834692015-06-02T10:28:00.000+09:002015-06-02T10:28:05.848+09:00GRASS GIS 7 でのbatファイル実行までのメモ共同研究者の方への説明メモを作成したので、ついでにブログにも投稿します。<br />
<br />
<b id="docs-internal-guid-5d933739-b1da-d468-7c8c-8fcdecfa592d" style="font-weight: normal;"><br /></b><br />
<h3 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 16px; font-style: normal; font-variant: normal; font-weight: bold; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS 7のインストール(Windows)</span></h3>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GISのサイトよりインストーラをダウンロード&インストール</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><a href="http://grass.osgeo.org/grass70/binary/mswindows/native/">http://grass.osgeo.org/grass70/binary/mswindows/native/</a></span></div>
<b style="font-weight: normal;"><br /></b>
<h3 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 16px; font-style: normal; font-variant: normal; font-weight: bold; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">初期設定</span></h3>
<h4 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: underline; vertical-align: baseline; white-space: pre-wrap;">GRASS GISのデータ置き場</span></h4>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GISでは一時的なデータ置き場を指定する必要があるため、予めWindows上でフォルダを作成しておく必要がある。今回の場合はDドライブ直下にgrassdataという名前のフォルダを作成して進める。</span></div>
<b style="font-weight: normal;"><br /></b>
<h3 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 16px; font-style: normal; font-variant: normal; font-weight: bold; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS上での初期設定</span></h3>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS 7を起動すると以下のようなスタート画面とコマンドプロンプトの2つが立ち上がる。</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="283px;" src="https://lh3.googleusercontent.com/f4N7dTferRkbWHKA5fhbGU7k9PqSeI-hu-27P-yiO3yHDHbGoXhaMBCUTOlrvcTiW96amXS-UqmaZgn7zyHkh1f0vjJyVvwsnwuk9UvV0SLOHi19cwwiKFyiPeECw6e2rM7omOQ" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="289px;" /></span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<br /></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">ここではまず、GRASS GISのデータ置き場を上記のように入力する(今回の例ではD:\grassdata)。</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">その後、GRASS Locationを設定する。GRASS Locationでは座標系を設定するので自分が使用するデータに応じて設定をする必要がある。</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"> </span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">以下に地理座標系の場合の設定方法を記す。</span></div>
<b style="font-weight: normal;"><br /></b>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">“Selectg GRASS Location”で「New」をクリック</span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">“Project Location”で適当な名前を指定する(スペースを空けずに英数字のみが無難)今回は地理座標系を扱うので”gcs”(←Geographic Coordinate System)としておく</span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">「Next」をクリック</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="193px;" src="https://lh3.googleusercontent.com/ToWxiWAR7P1nU0q1ClEjuLK5krvwBs_4aXD2x2piTtnhwL4mRRnUXNEE0RCwQJvOgRRLSqyZJK6dWIOq-ULokRwUVxpINCInxvzAziEtHHrfsoJMt7Xu2oT3AX-j-O_5e3dNsIs" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="312px;" /></span></div>
<b style="font-weight: normal;"><br /></b><br />
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">“Select coordinate system ~”を選択して「Next」をクリック</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="187px;" src="https://lh4.googleusercontent.com/pplSjEpsqs2pCN57cCF9sBLEi5SSjC1Ow1bkqXEC64wCIa4EI7ullaYo40HqdpobsY3bHrEhhPs0ABAp1DF2Y9YVxer1R3QvuWTiZSiM9uoj6PCWBt3cki71YAI9OmSuiACcChg" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="305px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">“Projection code”に ll と入力(小文字のLが2つ:Longitude Latitudeの略で地理座標系の場合はこれを指定)して、「Next」をクリック</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="198px;" src="https://lh4.googleusercontent.com/lmssXSe8V6hq6K7csh_7ouU6KwJ2wl9juaQzPd1aj_6fZmdHqu0Jg-IZzA0BP-mswE_e8PG4uktzFtYXyyjJfpcBSRRMFDu2nROhqD2tuhNAxOo2xYp88ewcmjrPxPR3Q1AoPDI" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="319px;" /></span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"> </span></div>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">ここでは何も変更せず「Next」をクリック</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="201px;" src="https://lh3.googleusercontent.com/_ug5oCZG0Ao0Qy1elxEW-1g_YX-rpu3MhWBIMpa7_y0X_EhDcVCW9QdXkkqUlCOitnLk0Ce0O5uQztKORBJjp-jmVE7eNWFUi5wb00s1TMzStV7owZB_ixZOpOYTROYqSQVwL9I" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="326px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">ここでは楕円体の指定をする必要があるので、”wgs84” と入力して「Next」をクリック</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="205px;" src="https://lh5.googleusercontent.com/1IB_noVm7YjWbgH9gLj9aT2JOo01LzKXMnAsB8kvAdbBylWcBkaRySQy3no07DrMM5c1ZtX9PLkV2hOx0dnxM__yLJx7Llf2SxXKGLhPKyNwMMcoAxM0cPIiksEwP6mkYVICDPc" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="331px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">ここでは何も変更せず「OK」をクリック</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="170px;" src="https://lh4.googleusercontent.com/spUJEmNqKavgSmQCmPUZvFRcvjosWvprKuSZ4CVeDgtI9uJAWj1Ko3jFuIC-PiRxjv5HT3NhFL-YdSGe92Bm7fzIDrqItDfA9cIrIrBnhLIRVF5MdWAUYfzLwE7PMyf2zPOT8rA" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="332px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"> </span></div>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">「Finish」をクリック</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="227px;" src="https://lh3.googleusercontent.com/aEH_aQ9LF3rxEYKs6vZ9ZSnVAZrSy1W1Dt8LiHp2o7OOppLaH-H3e_OedLSK9qE7oIM_xyMPd4lULgWwfkI-9RET8D_oL7Fu2DpV8M4UD-0g-ErPv-N5QdwqnW5DgT-18xVO-gM" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="368px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">その後、デフォルトの範囲と解像度を指定するかと聞いてきますが、こちらは後で変更できるので指定しないでかまいません(NoをクリックしてOK)。こちらの解像度指定はラスターデータの処理の場合に重要となりますが、ベクターデータの処理のみの場合は関係がありません。Yesをクリックしてしまった場合でも何もせずに「Set region」をクリックして終了させてかまいません。</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="185px;" src="https://lh4.googleusercontent.com/zYVdgWaeqbjqARvRJcUWe6OqpXMufnHQHm1JqGPkpXK-gesPSZSX7NRUfzXzllhQBY6JB2GAirYudO0cWbS0fRObO8N1ymQbcZxP4QCJwU8QA_mKrI7Kn0qycMaNYk0iqZDK0ss" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="466px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">最後にマップセットの名前を決めるように行ってきますが、こちらは好きな名前でかまいません。一つのPCを複数ユーザーで使う場合にユーザーごとに分けるためにある機能だと思います。おそらく。今回の例では “nuimura” としておきます。</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="151px;" src="https://lh3.googleusercontent.com/D0lh8cTFeQSnE--LHCIto1mrRPpCZuJnLf8ZTZ6QPRlEAu9SGSW4meD5JCDjQ3Ir5XXF5crmJsScYAPUcG0PD0x4SVyI72-PGCkTb5f7y81BPmJkVUr5mKN6EnZsZxN4qy9xA6Q" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="336px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"> </span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">ここまでの設定が済むと、以下の画面のようにLocationには”gcs”、mapsetには”nuimura”と”PERMANENT”の2つが表示されていると思います。mapsetの”nuimura”の部分をクリックして選択された状態にしてからウィンドウ下部の”Start GRASS session”をクリックすると、GRASSが起動します。</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="399px;" src="https://lh4.googleusercontent.com/Qq6BHO7cuySIEetVzQfFh4d2UJWc-EgrG7ABNYLqnf-Q8xaKkmyCUsEDF8jIR44uDLtvpd89yR-N3SNDA7D227rw58b_kFnwIGfO15zNwtQhvAIZkmdYVKBPfs_RMV9V2Ga1RTk" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="406px;" /></span></div>
<b style="font-weight: normal;"><br /></b><br />
<h3 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 16px; font-style: normal; font-variant: normal; font-weight: bold; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"> </span></h3>
<h3 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 16px; font-style: normal; font-variant: normal; font-weight: bold; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS 7 起動後</span></h3>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">正しく起動すると以下の3つのウィンドウがたちあがります。左から「コマンドプロンプト」、「レイヤマネージャー」、「マップディスプレイ」です。使用するPCのディスプレイサイズによっていは「コマンドプロンプト」ウィンドウは他のウィンドウの裏側に隠れているかもしれません。</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><img height="352px;" src="https://lh4.googleusercontent.com/MAWMd4pntgFlW7OFvzmMENmF5NwutpSYiyhwzrV3uSB7AN9hj8K9pLe7CzrwXK9KCUgXEjAGHfPQ7OOBPoFa0wV3PewzGKw84Q73R0CuvwfVN-uO8QfwPYZMpQW14i2opUwR__M" style="-webkit-transform: rotate(0.00rad); border: none; transform: rotate(0.00rad);" width="624px;" /></span></div>
<b style="font-weight: normal;"><br /></b>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS 7 では、「レイヤマネージャー」からはメニューをクリックすることで様々なGIS処理を行うことができます。「コマンドプロンプト」からはGRASS のコマンドを入力することで様々なGIS処理を行うことができます。「マップディスプレイ」にはそれらの処理結果を表示して確認することができます。</span></div>
<b style="font-weight: normal;"><br /></b>
<h3 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 16px; font-style: normal; font-variant: normal; font-weight: bold; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS 7でのスクリプト処理</span></h3>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS 7でも以前のバージョン同様にスクリプトを使った自動処理を行うことができます。</span></div>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">スクリプトの種類としては、</span></div>
<ul style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">bash</span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">bat</span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: disc; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">python</span></div>
</li>
</ul>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">などがあり、”bash”は「コマンドプロンプト」で”bash”と入力してbashモードに入ると使用でき、”bat”は「コマンドプロンプト」からそのまま使用可能で、”python”は「レイヤマネージャー」の”Python shell”タブから入力ができます。今回はすでにGRASS の処理を記述した”bat”ファイルがあるとして説明をしていきます。</span></div>
<b style="font-weight: normal;"><br /></b>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"> </span></div>
<h3 dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 8pt;">
<span style="background-color: transparent; color: #666666; font-family: 'Trebuchet MS'; font-size: 16px; font-style: normal; font-variant: normal; font-weight: bold; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">GRASS GIS の処理を記述したbatファイルの読み込み</span></h3>
<ol style="margin-bottom: 0pt; margin-top: 0pt;">
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: decimal; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">batファイル(ファイル名:cleaningShapefile.bat)が、以下の場所にあるとします。</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">D:\gamdam</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: decimal; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">コマンドプロンプトでは最初は以下の場所にいますので、batファイルのある場所に移動する必要があります。</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">C:\Users\[ユーザー名]</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">(今回の例では、C:\Users\nuimura)</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: decimal; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">Dドライブに移動するために以下のコマンドを入力</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">D:</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: decimal; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">Dドライブ直下のgamdamという名前のフォルダに移動するので以下のコマンドを入力</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">cd gamdam</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: decimal; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">移動後に dir というコマンドを入力するとそのフォルダ内のファイルがリストアップされます。使用したいbatファイルもあるか確認をして下さい。</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">dir</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span></div>
</li>
<li dir="ltr" style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; list-style-type: decimal; text-decoration: none; vertical-align: baseline;"><div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">その後、batファイルの実行のためにbatファイル名を入力します。</span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;"><br class="kix-line-break" /></span><span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">cleaningShapefile.bat</span></div>
</li>
</ol>
<b style="font-weight: normal;"><br /></b><br />
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: transparent; color: black; font-family: Arial; font-size: 15px; font-style: normal; font-variant: normal; font-weight: normal; text-decoration: none; vertical-align: baseline; white-space: pre-wrap;">以上で一連の手順は終わりです。</span></div>
<b style="font-weight: normal;"><br /></b>
<div dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;">
<br /></div>
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-73441215159253047462014-06-16T13:44:00.003+09:002014-06-19T09:07:18.992+09:00動画投稿テスト(Google Drive埋め込み)うちで飼っているサバトラ模様の猫さんです。<br />
<br />
<br />
<iframe src="https://docs.google.com/file/d/0B_6nHp3JYTaRbDZpSzVEUDJwaUk/preview" width="436" height="256"></iframe>
<br />
この動画の転載・配布は禁止いたします。Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-22422360183628118152012-12-27T15:20:00.000+09:002012-12-27T15:20:58.257+09:00GDAL/OGRコマンドをWindowsでもLinuxライクに使う手順<span style="text-align: -webkit-auto;">研究室で購入した衛星データを他のメンバーに使いやすくするために(HDF->GeoTIFF変換や各種派生データ(コンポジット、地形指標など)の作成)、その処理の多くをGDAL/OGRコマンドを使ったシェルスクリプトに記述しています。</span><br />
<div>
<div style="text-align: -webkit-auto;">
それらの処理を研究室の後輩への引き継ぐために、一連の手順の中でも汎用的な部分 (Cygwinのインストール&初期設定) をブログに載せます。</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<ol style="text-align: -webkit-auto;">
<li>Cygwinのインストール</li>
<li>OSGeo4Wのインストール</li>
<li>Cygwin.batの編集</li>
</ol>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
<span style="color: #0b5394;"><b>1. Cygwinのインストール</b></span></div>
<div style="text-align: -webkit-auto;">
CygwinはWindows上でLinuxの各種コマンドを使えるようにするツールです。後述のOSGeo4WでインストールされるGDAL/OGRコマンドはCygwinをインストールしなくてもWindowsコマンドプロンプトからでも使用出来るのですが、非常に使いづらいので(※個人の感想です)、Cygwinをインストールします。</div>
<div style="text-align: -webkit-auto;">
<br clear="none" /></div>
<div style="text-align: -webkit-auto;">
以下はインストールの手順です。</div>
<div style="text-align: -webkit-auto;">
Cygwinのサイト ( <a href="http://www.cygwin.com/" shape="rect">http://www.cygwin.com/</a> ) からインストーラをダウンロードして実行します。</div>
<div style="text-align: -webkit-auto;">
<div class="separator" style="clear: both; text-align: center;">
<a href="http://3.bp.blogspot.com/-G4ASXlUetzI/UNvlOY5kVAI/AAAAAAACYX0/6a3FaSIbO_4/s1600/Image.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="253" src="http://3.bp.blogspot.com/-G4ASXlUetzI/UNvlOY5kVAI/AAAAAAACYX0/6a3FaSIbO_4/s400/Image.png" width="400" /></a></div>
<br /></div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image.png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
インターネット経由でダウンロードが選択されているのを確認して次へ(つまりインストール時にはインタネット環境が必要)。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="http://2.bp.blogspot.com/-LlhR289mmT8/UNvlMevS1fI/AAAAAAACYW8/eGqcQYPaWho/s1600/Image+%255B1%255D.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="253" src="http://2.bp.blogspot.com/-LlhR289mmT8/UNvlMevS1fI/AAAAAAACYW8/eGqcQYPaWho/s400/Image+%255B1%255D.png" width="400" /></a></div>
<br />
<br /></div>
<div style="text-align: -webkit-auto;">
インストールディレクトリを変更することもできます。が特別な理由がなければデフォルトのままでいいでしょう。日本語フォルダをインストール先にするなどもってのほか。<br />
<span style="color: #cc0000; font-size: x-small;">※Windowsのユーザー名に日本語やスペースが入っていると以前はいろいろエラーが出たらしいです。最近はどうだか知りません。</span><br />
<div class="separator" style="clear: both; text-align: center;">
<a href="http://2.bp.blogspot.com/-rDo_Q7e_Psk/UNvlMRHxM3I/AAAAAAACYXE/mJRqbFpmiw4/s1600/Image+%255B2%255D.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="253" src="http://2.bp.blogspot.com/-rDo_Q7e_Psk/UNvlMRHxM3I/AAAAAAACYXE/mJRqbFpmiw4/s400/Image+%255B2%255D.png" width="400" /></a></div>
<br />
<br /></div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image(2).png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
インターネットからダウンロードしたインストーラの一時的な置き場所。インストール完了後はこのファイルは削除してもOK。このファイルを残しておくと次からはインターネット環境がなくてもこのファイルを呼んでCygwinインストールができる。まあよくわかんない人は消してOK。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="http://2.bp.blogspot.com/-CF3tsxpM2ig/UNvlMYGKEJI/AAAAAAACYXA/tUxYegLdiT4/s1600/Image+%255B3%255D.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="253" src="http://2.bp.blogspot.com/-CF3tsxpM2ig/UNvlMYGKEJI/AAAAAAACYXA/tUxYegLdiT4/s400/Image+%255B3%255D.png" width="400" /></a></div>
<br />
<br /></div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image(3).png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
で次へ次へと進んでいくとダウンロード元のサーバーを選択する画面がでます。適当に好きなとこを選んで下さい。日本国内のサーバーだとなんとなく早いような気がします(※個人の感想です)。選んだサーバーの状態が良くない場合は別の所を選んでと言われるので、その場合は別の所を選択して次へ進んで下さい。</div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image(4).png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<div class="separator" style="clear: both; text-align: center;">
<a href="http://1.bp.blogspot.com/--rh120PIni0/UNvlNHlvrSI/AAAAAAACYXI/hbHTKg1A1mA/s1600/Image+%255B4%255D.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="253" src="http://1.bp.blogspot.com/--rh120PIni0/UNvlNHlvrSI/AAAAAAACYXI/hbHTKg1A1mA/s400/Image+%255B4%255D.png" width="400" /></a></div>
<br />
<br />
<br /></div>
<div style="text-align: -webkit-auto;">
次にインストールしたいソフトウェアを選択する場面になります。カテゴリ別に多くのソフトウェアが含まれているのでインストールしたいソフトウェアがあれば"Skip"となっているところをクリックして"Install"に変更してください。デフォルトの状態で基本的なソフトウェアはすでに選択されているので特に追加インストールしたいソフトウェアが無ければそのまま次へをクリックして下さい。</div>
<div style="text-align: -webkit-auto;">
<span style="color: #cc0000; font-size: x-small;">※大量のソフトウェアがあり探すのが大変なので検索窓にソフトウェア名を打ち込んで探すと楽だと思います。例えば Imagemagick と打ち込むとGraphicカテゴリだけ表示されます。今回はためしに Imagemagick を追加でインストールしてみましょう。</span><br />
<div class="separator" style="clear: both; text-align: center;">
<a href="http://2.bp.blogspot.com/-6yZnlafOy6o/UNvlNDAzz_I/AAAAAAACYXQ/N4a8oU8e_mE/s1600/Image+%255B5%255D.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="271" src="http://2.bp.blogspot.com/-6yZnlafOy6o/UNvlNDAzz_I/AAAAAAACYXQ/N4a8oU8e_mE/s400/Image+%255B5%255D.png" width="400" /></a></div>
<br /></div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image(5).png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
インストールが完了すると下のような画面が出てきます。アイコンをデスクトップに置きたい場合はチェックを付けてください。<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="http://3.bp.blogspot.com/-95q-yraHTrQ/UNvlNS1XSRI/AAAAAAACYXU/wXeQ3uOblc8/s1600/Image+%255B6%255D.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="253" src="http://3.bp.blogspot.com/-95q-yraHTrQ/UNvlNS1XSRI/AAAAAAACYXU/wXeQ3uOblc8/s400/Image+%255B6%255D.png" width="400" /></a></div>
<br /></div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image(6).png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
<b><span style="color: #0b5394;">2. OSGeo4Wのインストール</span></b></div>
<div style="text-align: -webkit-auto;">
OSGeo4Wのサイト ( <a href="http://trac.osgeo.org/osgeo4w/">http://trac.osgeo.org/osgeo4w/</a> ) からインストーラをダウンロードして実行します。</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
インストールモードの選択画面が出てくるので、よくわかんない人は「エクスプレスデスクトップインストール」が選択されているのを確認して次へ進みましょう。</div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image(7).png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<div class="separator" style="clear: both; text-align: center;">
<a href="http://1.bp.blogspot.com/-J9bNesfJuUg/UNvlN8MrdyI/AAAAAAACYXg/7t2VFE1ddjM/s1600/Image+%255B7%255D.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="266" src="http://1.bp.blogspot.com/-J9bNesfJuUg/UNvlN8MrdyI/AAAAAAACYXg/7t2VFE1ddjM/s400/Image+%255B7%255D.png" width="400" /></a></div>
<br />
<br />
<br /></div>
<div style="text-align: -webkit-auto;">
インストールするソフトウェアを選択する画面が出てきます。よくわかんない人は全て選択して次へ進んで下さい。</div>
<div style="text-align: -webkit-auto;">
<span style="color: #cc0000; font-size: x-small;">※GDAL/OGRコマンドだけを使いたい場合はGDALだけにチェックを入れればOKなのですが、QuantumGISはあると何かと便利ですし、QuantumGISを入れると依存関係にあるGRASS GISもインストールされます。uDigとOpenEVはあまり使う機会はないかも(※個人の感想です)。</span><br />
<div class="separator" style="clear: both; text-align: center;">
<img border="0" height="266" src="http://3.bp.blogspot.com/-Cm_vYPChFck/UNvlOdlaf4I/AAAAAAACYXo/ohyT6A58Fjk/s400/Image+%255B8%255D.png" width="400" /></div>
<br /></div>
<div style="text-align: -webkit-auto;">
<img src="file:///R:/TEMP/enhtmlclip/Image(8).png" style="cursor: default;" /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
<b><span style="color: #0b5394;">3. Cygwin.batの編集</span></b></div>
<div style="text-align: -webkit-auto;">
実は先ほどまでの説明はググればいくらでも説明しているサイトはあります。ここからが重要な部分です。</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
Cygwinのアイコンのプロパティを開いてリンク先を見るとわかるのですが、Cygwinの起動はC:\cygwinにあるCygwin.batというバッチファイルで行われています。</div>
<div style="text-align: -webkit-auto;">
CygwinのなかでGDAL/OGRコマンドを使う場合はそのバッチファイル中にパスを通すコマンドを書く必要があります。</div>
<div style="text-align: -webkit-auto;">
<span style="color: #cc0000; font-size: x-small;">※環境変数をいじってもいいのですが別のソフトウェアに悪影響を与える可能性もあるので(特にpythonまわり)、バッチファイルに記述しておくのが安全です。</span></div>
<div style="text-align: -webkit-auto;">
<br />
<br /></div>
<div style="text-align: -webkit-auto;">
まずCドライブのcygwinフォルダに移動して、Cygwin.batというバッチファイルで右クリックをして"編集"をクリックします。</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
バッチファイルには以下のように記述されていると思います。</div>
<div style="background-color: #e8e8ff; border-radius: 10px;">
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
@echo off</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
C:</div>
<div style="text-align: -webkit-auto;">
chdir C:\cygwin\bin</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
bash --login -i</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
<br /></div>
</div>
<div style="text-align: -webkit-auto;">
<br />
ここに以下の3行を@echoのあとに加えます。</div>
<div style="background-color: #e8e8ff; border-radius: 10px;">
<div style="text-align: -webkit-auto;">
<span style="text-align: -webkit-auto;">set OSGEO4W_ROOT=C:\OSGeo4W</span></div>
<div style="text-align: -webkit-auto;">
PATH=%OSGEO4W_ROOT%\bin;%PATH%<br />
for %%f in (%OSGEO4W_ROOT%\etc\ini\*.bat) do call %%f</div>
</div>
<div style="text-align: -webkit-auto;">
<br />
追加後は以下のような感じになると思います。</div>
<div style="text-align: -webkit-auto;">
※追加行は赤で表示。</div>
<div style="background-color: #e8e8ff; border-radius: 10px;">
<div style="text-align: -webkit-auto;">
<span style="text-align: -webkit-auto;">@echo off</span></div>
<div style="text-align: -webkit-auto;">
<br />
<span style="color: red;">set OSGEO4W_ROOT=C:\OSGeo4W<br />
PATH=%OSGEO4W_ROOT%\bin;%PATH%<br />
for %%f in (%OSGEO4W_ROOT%\etc\ini\*.bat) do call %%f</span><br />
<br />
C:<br />
chdir C:\cygwin\bin<br />
<br />
bash --login -i</div>
</div>
<div style="text-align: -webkit-auto;">
<br /></div>
<div style="text-align: -webkit-auto;">
で保存すれば、次にCygwinを起動したときからGDAL/OGRのコマンドが使えるようになります。<br />
<br /></div>
</div>
Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-54034904990933144762012-12-13T00:01:00.000+09:002012-12-13T00:01:18.358+09:00AGU2012 Fall Meetingで行われたOSGeo Workshop参加報告FOSS4G Advent Calendar 2012 の13日目の記事です。<br />
<a href="http://atnd.org/events/34052">http://atnd.org/events/34052</a><br />
<br />
今年のFOSS4G Advent Calendarは今のところなかなかの猛者揃いで、個人的には満足です。自分も当初は「クロスプラットフォームなスクリプト作成法(GDAL、GRASS、GMT、R)」というネタを密かに考えていたのですが、AGUのOSGeoの<a href="http://fallmeeting.agu.org/2012/events/introduction-to-the-web-enabled-landsat-data-weld-products-using-open-source-software/">Workshop</a>が面白かったのでその紹介をします。当初予定していたネタはまた別の機会に。
<br />
<br />
<ol>
<li>はじめに </li>
<li>ハンズオンの概要 </li>
<li>WELD (Web-Enabled Landsat Data)について </li>
<li>ハンズオンの内容 </li>
<li>おわりに </li>
</ol>
<br />
<h3>
<span style="color: #0b5394;">1. はじめに</span></h3>
12月3日(月)から7日(金)の間サンフランシスコで行われた<a href="http://fallmeeting.agu.org/2012/">AGU (American Geophysical Union: 米国地球惑星科学連合) Fall Meeting</a>という、地球惑星科学系では最大の研究集会(参加者数20000名以上)に参加してきました。自分はそちらでポスター発表を行なってきたのですが、今回は個人的に注目すべきイベントとしてアメリカのOSGeoメンバーによる「オープンソースソフトウェアを用いたLandsatデータ解析」というハンズオンがありました。通常のセッション修了後の6:00 pm~10:00 pmと遅い時間だったのですが、早めの夕食を取り後輩と2人で参加してきました。<br />
<br />
ハンズオンの会場は、AGUのメイン会場であるMoscone Centerではなく、すぐ近くのMarriott Marquisという高級そうなホテルの会議室で行われました。講師を務めたのはサウスダコタ大学の人(<a href="http://globalmonitoring.sdstate.edu/faculty/roy/roy.html">David Roy</a>と<a href="http://globalmonitoring.sdstate.edu/people.php?view=6&a=show&id=116">Giuseppe Amatulli</a>)で、アシスタントは4名、参加者は意外と多く20数名はいました。WELDの紹介はDavid Roy教授が行なっていましたが、その他は全てDr. Giuseppe Amatulli(←じーだる発音派)が行なっていました。<br />
<br />
<hr />
<h3>
<span style="color: #0b5394;">2. ハンズオンの概要</span></h3>
このハンズオンは以下のような構成になっていました。<br />
<br />
<プレゼンテーション><br />
6:00 pm - 6:10 pm ワークショップの概要説明とUSBドライブでのデータ(VirtualBoxイメージとインストーラ(Windows、Mac))配布。<br />
6:10 pm - 6:30 pm WELDプロジェクトの説明 (Roy)<br />
6.30 pm - 7.00 pm VirtualBoxのインストールと機能の説明、Linuxのコマンドの説明 (Amatulli)<br />
7.00 pm - 7.30 pm 使用するオープンソースソフトウェアの紹介:GDAL/OGR、R、OFT、PKTOOLS (Amatulli)<br />
7.30 pm - 7.45 pm 休憩<br />
<br />
<ハンズオン><br />
7.45 pm - 8.15 pm OpenEVでWELDデータの可視化 (Amatulli)<br />
8.15 pm - 8.30 pm WELDデータのグリッド単位での時系列変化について (Amatulli)<br />
8.30 pm - 9.30 pm WELDデータの森林火災エリアのマッピング:自動分類や変化の検出 (Amatulli)<br />
9.30 pm - 9.45 pm WELDデータの大気補正前と後の比較 (Amatulli)<br />
9.45 pm - 10.00 pm オープンディスカッション (Roy) <br />
<br />
会場に入るとまず<a href="https://www.virtualbox.org/wiki/Downloads">VirtualBox</a>インストーラー&イメージの入ったUSBディスクのコピーをするように言われます。VirtualBoxイメージはzip圧縮された状態で約3GB(解凍後は約8GB)で、このイメージを起動すると今回使用するオープンソースソフトウェア、使用データ、マニュアルの含まれたXubuntuが起動します。自分はこのzipファイルをWindows 7 上で解凍する際に何故か「600PB(ペタバイト)以上の空き領域が必要です」とのエラーメッセージが出て解凍できず(無理w)。デュアルブートのUbuntuで解凍して事なきを得ました。その他にも小さい危機(xubuntuのキーボードがUS設定なので記号入力に最初戸惑うなど)がありましたが、無事最後まで脱落せずについていくことができました。 <br />
<br />
前半(プレゼンテーション)では、<a href="http://landsat.usgs.gov/WELD.php">WELD</a>プロダクトの説明や、使用するオープンソースソフトウェアの紹介をしていました。使用したソフトウェアについては、<a href="http://www.gdal.org/">GDAL/OGR</a>と<a href="http://www.r-project.org/">R</a>についてはよく知っていたのですが、<a href="http://km.fao.org/OFwiki/index.php/Open_Foris_Geospatial_Toolkit">OFT</a>と<a href="http://www.spatial-ecology.net/dokuwiki/doku.php?id=wiki:pk_tools">PKTOOLS</a>は今回のハンズオンで初めて知りました。また、後半のハンズオンではLinuxのコマンドを多く使用するため、初歩的なLinuxコマンドの紹介(ls、less、head、awkなど)もしていました。さらりとawkを紹介している辺りに後半は荒れるなという予感がしました。 <br />
<br />
後半(ハンズオン)の大まかな流れは、WELDデータを使って森林火災エリアを含むシーンの教師付き分類をSVM(サポートベクターマシン)法とRandomforest法で行い、各プロセスが終わるごとに<a href="http://openev.sourceforge.net/">OpenEV</a>あるいはQGISで結果の確認を行い、最後におまけとしてWELD version 1.5(大気補正前)とversion 2.0(大気補正後)の違いを紹介していました。
<br />
<br />
<hr />
<h3>
<span style="color: #0b5394;">3. WELD (Web-Enabled Landsat Data)について</span></h3>
<a href="http://landsat.usgs.gov/WELD.php">WELD</a>はLandsat 7 ETM+ オルソ補正画像をモザイクしてシームレスに作成した画像プロダクトで、土地被覆分類や、様々な地表面解析を目的として、USGSとサウスダコタ大学の共同プロジェクトで作成されたプロダクトです。対象領域はアメリカ本土とアラスカのみで、30m解像度で提供されています。<br />
<br />
複数時期のモザイクしたプロダクトのため、各グリッドに用いられた元データの時期がわかる補助データも付随しています。特定の時間での地表面状態を知りたいという研究にはあまり向いていないのですが、地表面被覆の大まかな分布がわかれば良い研究分野にとってはハンドリングが楽なデータであるといえます。US限定ですが。
<br />
<br />
<hr />
<h3>
<span style="color: #0b5394;">4. ハンズオンの内容</span></h3>
<span style="color: #0b5394;"><a href="http://globalmonitoring.sdstate.edu/projects/weld/training/weld_visual_fire.html">4-1. OpenEVでWELDデータの可視化</a></span><br />
gdal_mergeによりコンポジット画像を作成しOpenEVでヒストグラムの調整などをして可視化。その後はgdal_calcでログスケール変換した各バンドを同様にコンポジット→可視化。<br />
<br />
gdal_calcはほんとうに便利で簡単なバンド演算はほとんどこれで可能。自分が従来GRASSで行なっていた単純なバンド演算のいくつかも既にgdal_calcに移行しています。
<br />
<br />
<span style="color: #0b5394;"><a href="http://globalmonitoring.sdstate.edu/projects/weld/training/time_series.html">4-2. WELDデータのグリッド単位での時系列変化について</a></span><br />
WELDではサイト上で緯度経度で指定したグリッドにおける放射量の時系列情報を得ることもできます。ハンズオンでは、CSVデータとして用意されたあるグリッドにおける時系列データをheadコマンドで確認をしつつ、awkによるデータ整形、Rによる時系列プロットを行なっていました。自分はRは普段良く使っているので問題なかったのですが、Rの未経験者にとってはなかなかしんどそうなところでした。一回目の山場です。<br />
<br />
<span style="color: #0b5394;"><a href="http://globalmonitoring.sdstate.edu/projects/weld/training/openforis_burned.html">4-3. WELDデータの森林火災エリアのマッピング:自動分類や変化の検出</a></span><br />
このステップは今回の二回目の山場で、脱落者も続出していました。隣席のお姉さんも途中で「オー・・・。」とかいって大きなため息をついて脱落してました。自分もついていくのが精一杯で助ける余裕はありませんでした。<br />
<br />
まずはPKTOOLSを使って多バンドのコンポジット作成(pkcropコマンド)とマスク画像の作成(pkgetmask)。gdal_rasterizeで教師データとしての森林火災域(とおもわれる領域)のポリゴンをラスタライズを前準備として行いました。ちなみにPTOOLSはGDALライブラリを使ってC++で実装されたコマンド群です。PKTOOLSコマンドの多くはすでにあるような処理ばかりでしたが…、SVM分類などが簡単にできるのは便利だと思いました。 <br />
<br />
その後1つ目の教師付き分類(SVM)の処理に入ります。ラスタライズした教師領域の0.01%を処理時間節約のためにランダム抽出して教師データとして利用してpkclassify_svmコマンドでSVM分類。できたものはOpenEVで可視化して確認を行いました。 <br />
<br />
続いて2つ目の教師付き分類(Randomforest)の処理に入ります。ここが非常に長くて複雑。流れとしてはまずOFTを使ってK-means NN教師なし分類でセグメンテーション分類したうえで、教師データを用いてR上でRandomforestの教師付き分類を行います。ただひとつ残念なのがRでのRandomforest分類の結果は直接ラスターイメージとして出力されるのではなく、初めのセグメンテーション分類IDとRandomforest分類の対応を示すテキストファイルとして出力されるため、そこからラスターイメージに持っていくために処理が複雑化しているという部分です。この部分の個別の処理を知りたい方はリンク先の説明を参照下さい。<br />
<br />
<span style="color: #0b5394;"><a href="http://globalmonitoring.sdstate.edu/projects/weld/training/weld_visual.html">4-4. WELDデータの大気補正前と後の比較</a></span><br />
ここではPKTOOLSのpkinfoコマンドで待機補正前と後のRGBバンドの画素値を、それぞれ個別のテキストファイルに出力して、それをgnuplotにて可視化して比較を行いました。<br />
<br />
<span style="color: #0b5394;">4-5. オープンディスカッション</span><br />
最後にこのハンズオンセッションについて、難易度や作業量はどうだったか参加者の意見を聞いていました。自分が聞き取れた範囲では、少し作業量が多杉という意見がちらほら出ていました。<br />
<br />
<hr />
<h3>
<span style="color: #0b5394;">5. おわりに</span></h3>
このハンズオンでは多くのオープンソースソフトウェアを使用するため、VirtualBoxを使ってあらかじめ主催者側が環境設定したイメージを利用しました。おかげで複数のソフトウェアを組み合わせて一連の解析を行うという実践的な処理が可能となり、個人的には非常に勉強になるハンズオンだと感じました。 <br />
<br />
ちなみに今週土曜日に中部大学で第3回FOSS4Gツール勉強会@名古屋を行います。こちらは全然スパルタではないので安心して気軽にご参加ください。<br />
<a href="https://sites.google.com/site/foss4gnagoya/foss4gnagoya20121215"><span style="font-size: x-small;">https://sites.google.com/site/foss4gnagoya/foss4gnagoya20121215</span></a><br />
<br />
あとAGUの通常のセッションの参加報告を研究室のブログにあげていますのでこちらも興味がありましたらどうぞ。<br />
<a href="http://hello.ap.teacup.com/snowman/1602.html">http://hello.ap.teacup.com/snowman/1602.html</a>Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-59408179116621340312012-10-31T12:53:00.001+09:002012-10-31T12:53:40.876+09:00FOSS4G2012(Tokyo/Osaka)のお知らせオープンソースGISソフトウェアの祭典、FOSS4G2012が11月に東京と大阪で行われます。
昨年はいち参加者として、初めてFOSS4G2011 Tokyoに行きましたが、今年は運営側として初めての参加となります。<br />
<br />
<h3>
<a href="http://www.osgeo.jp/%E3%82%A4%E3%83%99%E3%83%B3%E3%83%88/foss4g2012tokyo/">
FOSS4G Tokyo (11/3, 4, 5)
</a>
</h3>
<ul>
<li>11/3 <a href="http://www.osgeo.jp/%E3%82%A4%E3%83%99%E3%83%B3%E3%83%88/foss4g2012tokyo/foss4g2012tokyoprogram/">コミュニティデイ</a></li>
<li>11/4 <a href="http://www.osgeo.jp/%E3%82%A4%E3%83%99%E3%83%B3%E3%83%88/foss4g2012tokyo/foss4g2012tokyoprogram/">コアデイ</a></li>
<li><a href="http://www.osgeo.jp/%E3%82%A4%E3%83%99%E3%83%B3%E3%83%88/foss4g2012tokyo/foss4g2012tokyohandson/">11/5 ハンズオン</a></li>
</ul>
<h3>
<a href="http://www.osgeo.jp/%E3%82%A4%E3%83%99%E3%83%B3%E3%83%88/foss4g2012osaka/">
FOSS4G Osaka (11/7, 8)
</a>
</h3>
<ul>
<li>11/7 <a href="http://www.osgeo.jp/%E3%82%A4%E3%83%99%E3%83%B3%E3%83%88/foss4g2012osaka/foss4g2012osakahandson/">ハンズオン</a></li>
<li>11/8 <a href="http://www.osgeo.jp/%E3%82%A4%E3%83%99%E3%83%B3%E3%83%88/foss4g2012osaka/foss4g2012osakaprogram/">コアデイ</a></li>
</ul>
<br />
2010年度までは自分の博士論文の執筆に追われて余裕がなく、このようなコミュニティの活動にも興味を持ってはいたのですが、まずはD論書かなきゃヤバイ(汗)ということで、密かに?FOSS4G活動(といってもGRASS GISやGDAL、QGIS、Rで解析処理をするだけでアウトプットは皆無)をしていました。<br />
<br />
<br />
そんなこんなで蓄積した自分の知識・経験をコミュニティに還元しつつ、それを武器として人脈を広げてみようと思い、<a href="http://corpus-study.info/nagoyar/">Nagoya.R</a>というR言語の勉強会が名古屋大学で不定期に行われていたので、そちらでRによるGISやリモートセンシング解析の話をしたりWebに発表資料の公開などを始めました。するとそのWebの資料を見つけて、以前衛星ワークショップや学内勉強会などで少しだけ知り合いだった平松さんにRをつかったGISとリモセン解析に興味を持って連絡をいただき、各種FOSS4Gツールについての勉強会面白そうだねということで<a href="https://sites.google.com/site/foss4gnagoya/">FOSS4Gツール勉強会@名古屋</a>の話が立ち上がりました。<br />
<br />
<br />
そして立ち上げにあたって、FOSS4Gの本丸にていろんな方に宣伝したりご助言などを頂けたらと思い、FOSS4G2011 Tokyoへと参加して、気がついたら中の人になってました。FOSS4G2011 Tokyoに参加して以来、多くの方々と知り合うことができたため、この1年間は自分にとっては実りのあるとしでした(プライベートでも節目の年でした)。<br />
<br />
<br />
最後に宣伝を、11/4と11/8には様々なFOSS4Gソフトウェアを使った地形解析処理の紹介の話をします。また、11/5にはQGISのハンズオンの講師を、11/8にはRでリモートセンシングハンズオンの講師をつとめます。これを機会に研究&業務利用のGIS及びリモートセンシング環境をオープンソース環境への切り替えを!Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-25443733114667194092012-06-22T14:38:00.004+09:002012-06-22T20:25:35.932+09:00Octaveを使った地上開度および地下開度の計算<h3>
<span style="color: #0b5394;">
地上開度/地下開度という地形量</span></h3>
DEM(Digital Elevation Model:デジタル標高モデル)から計算できる地形量には様々なものがあります。有名なものとしては「傾斜」、「斜面方位」などが挙げられ、他にも斜面が凸か凹かの具合を示す「曲率」などもあります。<br />
今回とりあげる地上開度というのは、ある地点で空を見上げたとき、どれくらい開けているかを示すものです。谷底では地上開度は小さくなり、尾根上では大き
くなります。また、同じ谷底地形の中でも、V字谷のようにすぐそばに斜面が迫っている谷では小さく、U字谷のように谷底に広い平坦地が広がっている谷では
V字谷と比べて大きくなります。<br />
<br />
<h3>
<span style="color: #0b5394;">
地上開度/地下開度のアルゴリズム</span></h3>
Prima et al. (2006)で地上開度の計算アルゴリズムが説明されています。細かい説明はこの論文をご覧いただくとして、ざっくり説明すると以下のとおりです。<br />
<div>
まず計算するにあたって探索半径というものを決めます。そのあとDEMの各グリッドにおいて、東西南北とそれらの中間の計8方位の探索半径内のグリッドの
標高値と中心グリッドの標高値から仰角を計算して8方位それぞれにおける最大仰角を計算します。そしてその8つの値から平均最大仰角を求めて、最終的に天
頂から平均最大仰角を引いたものが地上開度となります。</div>
<div>
ちなみに地下開度の場合は仰角を俯角に、天頂を天底(←この日本語初めて知りました、英語ではNadir)に置き換えます。<br />
<br /></div>
<h3>
<span style="color: #0b5394;">
Octaveで計算してみよう</span></h3>
<div>
前置きが長くなりましたが、Octaveで地上開度および地下開度を計算するOctaveプログラムをこちらに公開します。</div>
<ul>
<li><a href="https://dl.dropbox.com/u/870568/blogger/openness_calc.m">openness_calc.m</a> : <span style="background-color: white;">地上開度/地下開度のメインの計算部分</span></li>
<li><a href="https://dl.dropbox.com/u/870568/blogger/gdalread.m">gdalread.m</a> : <span style="background-color: white;">GeoTIFFなどのGDALライブラリで対応しているラスターデータをOctaveに読み込む関数プログラム</span></li>
<li><a href="https://dl.dropbox.com/u/870568/blogger/gtiffwrite.m" style="background-color: white;">gtiffwrite.m</a><span style="background-color: white;"> : </span><span style="background-color: white;">計算結果をGeoTIFFに書きだす関数プログラム</span></li>
<div>
</div>
</ul>
<div>
<span style="color: #cc0000; font-size: x-small;">※2012-6-22 17:10 追記</span><br />
<span style="color: #cc0000; font-size: x-small;">gdalreadとgtiffwriteはGDAL/OGRのgdal_translateというプログラムを利用しているため、OSGeo4Wというインストーラを利用してGDAL/OGRをインストールしておく必要があります。</span><br />
<span style="color: #cc0000; font-size: x-small;">OSGeo4W (</span><span style="font-size: x-small;"><a href="http://trac.osgeo.org/osgeo4w/" style="background-color: white;">http://trac.osgeo.org/osgeo4w/</a><span style="background-color: white; color: #cc0000;">)</span></span><br />
<br />
<br /></div>
<div>
使用方法:</div>
<ol>
<li>上記3つのプログラムと計算したいDEM(例. aster2003_dem30.tif)ファイルを同じディレクトリにおきます。※DEMは投影座標じゃなきゃダメ。</li>
<li>Octaveを起動して、そのディレクトリに移動します。</li>
<li>Octaveコンソールで、</li>
</ol>
<pre><code>openness_calc('aster2003_dem30.tif', 300, 'up')
</code>
</pre>
<div>
というようにコマンドを打ちます。上のオプションは1つめがファイル名、2つめが任意の探索半径、3つめが地上開度(up)か地下開度(down)かの選択フラグです。</div>
<div>
<br clear="none" /></div>
<div>
結構計算時間はかかります。フリーソフトのSAGA GISでも開度計算はできるらしいのでそちらの方が楽かも(汗)。<span style="background-color: white;">上記ソースコードをもとに計算アルゴリズムをチューニングして遊びたいという方はこのOctaveプログラムは参考になるかもしれません。</span></div>
<div>
<br clear="none" /></div>
<div>
<br /></div>
<div>
引用文献</div>
<div>
Prima,
O. D. A., Echigo, A., Yokoyama, R., & Yoshida, T. (2006).
Supervised landform classification of Northeast Honshu from DEM-derived
thematic maps. Geomorphology, 78(3-4), 373-386. Elsevier Science Bv.
doi:10.1016/j.geomorph.2006.02.005</div>
<div style="font-weight: bold;">
<br /></div>Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com2tag:blogger.com,1999:blog-7751709962087367040.post-38971201222755629712012-03-29T12:16:00.000+09:002016-12-22T16:40:21.682+09:00GMT 4.5.7のインストール (Cygwin編)<span style="color: #cc0000; font-size: xx-small;">2016/12/22現在の最新バージョンはGMT 5.3.1です。<a href="http://spatiohack.blogspot.jp/2016/12/gmt-531-cygwin.html">そちらに対応させた記事</a>を書いていますので、この古い記事ではなくそちらを参考にすることをおすすめします。</span><br />
<br />
<a href="http://gmt.soest.hawaii.edu/">GMT (Generic Mapping Tools)</a>とはコマンドラインベースの地図作成ツールで、地球科学分野でよく用いられるツールで、綺麗な空間分布図を描画することができます。クロスプラットフォームなソフトなのでWindows、Mac、Linuxで使うことができます。
<br />
<br />
<br />
<br />
今回は主にWindows 7でのGMTのインストール(&初期設定)について説明します。<br />
<br />
<hr />
<h4>
1. はじめに</h4>
<br />
WindowsでGMTを使う場合は、GMTそのもの以外にも以下のツールが必要となってきます。<br />
<br />
<br />
<b><span style="color: #073763;">必須
</span></b>
<br />
<div style="background-color: #e8e8ff; border-radius: 10px;">
<dl>
<dt><a href="http://pages.cs.wisc.edu/~ghost/gsview/get50.htm">GSView</a></dt>
<dd>GMTで出力したps/epsファイルの表示や変換処理に必要
(※実際に変換しているのはGhostscriptやpstoedit、GSViewはインターフェース)。</dd>
</dl>
</div>
<span style="color: #073763;"><b>推奨
</b></span>
<br />
<div style="background-color: #e8e8ff; border-radius: 10px;">
<dl>
<dt><a href="http://www.cygwin.com/">Cygwin</a></dt>
<dd>シェルスクリプトでGMTコマンドを使う場合必要。WindowsバッチでもGMTは扱えるがCygwinの方が多機能で便利。</dd>
<dt><a href="http://www.ghostscript.com/download/gsdnld.html">Ghostscript</a></dt>
<dd>GMTで出力したps/epsファイルをビットマップ画像に変換するのに必要。</dd>
<dt><a href="http://www.pstoedit.net/">PStoedit</a></dt>
<dd>GMTで出力したps/epsファイルをベクター形式のままSVG形式に変換できる。Inkscapeなどのイラストソフトを使っている人には便利。</dd>
<dt>改行コードやエンコーディングが変更できるエディタ</dt>
<dd><span style="font-size: x-small;">(<a href="http://www.geocities.jp/npp_euc_jp/">Notopad++EUC-JP対応版</a>がおすすめ、自分は<a href="http://vgoulet.act.ulaval.ca/en/emacs/">emacs modified</a>を使ってるけど…)</span><br />
異なるOSや他人とのファイルの移動を行う場合は、環境の違いによって日本語コメントなどが文字化けする場合がある。<br />
受け渡し先に合わせたエンコーディングや改行コードに変換したり、人からもらったスクリプトファイルのエンコーディング&改行コードが自分の環境に合わない場合などにあると便利<span style="font-size: x-small;">(でも本当はコメントは英語で書いておくのが無難)</span>。</dd>
</dl>
</div>
<hr />
<h4>
2. インストール</h4>
GMTサイトの<a href="ftp://ftp.scc.u-tokai.ac.jp/pub/gmt/windows/">windows用ダウンロードのページ</a>から、<br />
<br />
<ul>
<li>Windows (64 bit)の人は、
<i>gmt-4.5.7_install64.exe</i></li>
<li>Windows (32 bit)の人は、
<i>gmt-4.5.7_install32.exe</i></li>
</ul>
<br />
をダウンロードして、インストールを行います。<br />
<br />
ちなみにGMTの実行ファイルは、
<i>C:\programs\GMT4 </i>フォルダにあります。
昔のGMTではインストール後(というかインストーラが無かったので、圧縮された実行ファイルをダウンロード後、解凍して任意の場所に置くだけ)、パスの設定を行う必要がありましたが、最近のGMTでは自動でパスの設定をしてくれます。<br />
<br />
なのでコマンドプロンプトやCygwinターミナルで何かGMTのコマンド(例えば、psxy)を打ち込むと、大量のメッセージが出てきたらインストール成功です。<br />
もし「'psxy'は、内部コマンド~(中略)~として認識されていません」とか「psxy:コマンドが見つかりません」などと出てきたらインストール失敗です。そのようなことは通常無いと思うのですが(少なくとも自分は見たことないです)、万が一そのようなことが起こった場合は自分で環境変数の設定をしてパスを通す必要があります。
必要なユーザー環境変数は、以下のとおりです。
<br />
<br />
<div align="center">
<table border="0" cellpadding="3" cellspacing="0" style="text-align: center;">
<tbody>
<tr><th style="background-color: #cacaff;">変数</th><th style="background-color: #cacaff;">値</th></tr>
<tr>
</tr>
<tr><td>GMT_SHAREDIR</td><td><i>C:\programs\GMT4\share</i></td></tr>
<tr><td>path</td><td><i>C:\programs\GMT4\bin</i></td></tr>
</tbody></table>
</div>
<br />
※設定後は再起動が必要です(たぶん)。<br />
<br />
<br />
上記のGMT&周辺ツールのインストールが終わったら、Cygwinターミナルで試しに<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><b>pscoast -Jm1:30000000 -R120/150/20/50 -Ba10f5g5 -Gtan > sample.eps</b></span><br />
<br />
と打ってみて下さい。<br />
<span style="color: #cc0000; font-size: x-small;">※表示の都合で2行にまたがっているかもしれませんが、1行で入力して下さい</span>
<br />
<span style="color: #cc0000; font-size: x-small;"><br /></span>
カレントディレクトリにsample.epsというEPSファイルができているはずです。このEPSファイルをWindowsで見るために、Gsviewが必要となります。<br />
<br />
すでにGsviewをインストールしていれば、EPSファイルに関連付けられているはずなのでダブルクリックで作成した図を見ることができると思います。
Adobe Illustratorをインストールしている場合はそちらに関連付けられているかもしれません。<br />
<br />
その場合は図を見るだけならGSview、編集する場合はIllustratorと使いわけると良いと思います(Illustratorの起動はGSviewより重いため)。<br />
<br />
<table style="width: auto;"><tbody>
<tr><td><div class="separator" style="clear: both; text-align: center;">
<a href="https://picasaweb.google.com/lh/photo/RgEelmIAJlKABQddJQS3ntMTjNZETYmyPJy0liipFm0?feat=embedwebsite" style="margin-left: 1em; margin-right: 1em;"><img height="400" src="https://lh6.googleusercontent.com/-MHYKd7eT2v0/T3PKy4QnpGI/AAAAAAABk1g/NeCvDKXH6eQ/s400/sample.png" width="345" /></a></div>
</td></tr>
<tr><td style="font-family: arial,sans-serif; font-size: 11px; text-align: right;">From <a href="https://picasaweb.google.com/104578933986350101585/Research_public?authuser=0&feat=embedwebsite">Research_public</a></td></tr>
</tbody></table>
pscoastコマンドで作成した図Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-67056394553327650412011-12-14T09:10:00.000+09:002011-12-25T10:06:51.236+09:00FOSS4G勉強会@名古屋の開催報告FOSS4G勉強会@名古屋の開催報告
かなり久々のブログ更新w。<br />
<br />
FOSS4G Advent Calendarの企画です。<br />
<a href="http://atnd.org/events/23085">http://atnd.org/events/23085</a><br />
<br />
昨日までのFOSS4G Advent Calendarではハイレベルな技術的な話が続いたので、ここでいったんハードルを下げときます。先日第1回FOSS4G勉強会@名古屋を開催致しましたので、その報告を。<br />
<br />
日時:2011/12/11(日)13:00-17:00<br />
場所:名古屋大学 環境総合館 3F 講義室2<br />
<br />
<br />
初回にも関わらず、40名弱の参加者数&10件の発表件数とみなさんに多くの関心を持っていただきありがとうございました。参加者も東海地方のみならず、西は奈良、東は東京からもお越しいただきました。<br />
<br />
<br />
今回の以下の10件の発表をいただきましたので、順に紹介していきます(敬称略)。<br />
<br />
<hr />
<span class="Apple-style-span" style="color: #0b5394;"><b>1. 縫村崇行(名古屋大学)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">初心者の方向けにQGISの機能の紹介
</span><br />
<div id="__ss_10554263" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/nuimura/111211-foss4g-nagoyapresentation1-10554263" target="_blank" title="111211 foss4g nagoya_presentation1">111211 foss4g nagoya_presentation1</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10554263" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/nuimura" target="_blank">Takayuki Nuimura</a> </div>
</div>
私の発表です。
発表しておいてなんですが、実はQGISはちょっとしたデータのビューアーとしてしか使っておらず。解析はほとんどしたことがありません(笑)。なので、今回は超基本的なところをさらっとしかお話できませんでしたが、今回の話でインストールしただけで一度も使っていないという初心者の方の後押しができたらなと(笑)。今回は他の方々の超高度なQGISの使いこなしに圧倒されまくりでした。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>2. 平松孝晋(アジア航測)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">QGIS+GRASS GIS関連や地形解析(SAGA GIS)</span><br />
<span class="Apple-style-span" style="color: #0b5394;"><br /></span><br />
特にすごかったのがSAGAで多様な地形解析をサクッとデモンストレーションされたところ。デモでの操作も非常に手際がよくプロフェッショナルを感じました。SAGAはコマンドラインでも使えるので、Linuxでbashで自動処理もできるそうです。時代はSAGAかも。Ubuntuでもパッケージにあるので簡単にインストールできるし。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>3. 山本一清(名古屋大学)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">学部教育と研究におけるGRASS GISの利用例の紹介</span><br />
<span class="Apple-style-span" style="color: #0b5394;"><br /></span><br />
<div id="__ss_10683433" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/nuimura/foss4-g2011" target="_blank" title="Foss4 g勉強会2011">Foss4 g勉強会2011</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10683433" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/nuimura" target="_blank">Takayuki Nuimura</a> </div>
</div>
学部での実習授業に関してとrubyを使ったGRASSの解析のお話、GRASSのバッチ処理には多くのやり方があり(自分はbash派)、山本さんはrubyで行なっているそうです。実習授業ではUbuntuサーバーに学生がVNCでアクセスして授業を行うとのことで、名大生命農学研究科の学生さんが羨ましくなりました。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>4. 福山 薫(三重大学)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">小学生がQGISで地域の防災マップ作り</span><br />
<span class="Apple-style-span" style="color: #0b5394;"><br /></span><br />
<div id="__ss_10619618" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/FukuyamaKaoru/qgis-10619618" target="_blank" title="小学生がQGISで地域の防災マップ作り_名古屋FOSS4G勉強会">小学生がQGISで地域の防災マップ作り_名古屋FOSS4G勉強会</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10619618" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/FukuyamaKaoru" target="_blank">FukuyamaKaoru</a> </div>
</div>
オープンソースGIS歴20年の実力を見せつけられました。本当に勉強になりました。小学生たちの教育ではGMTで出力したマップを用いて、街の防災調査を行い、情報をQGIS上に位置情報&写真付きで集約。最終的にQGISのeVisプラグインで地図上に位置づけられた写真が簡単に見られたり、geoPDFに出力したりと、途中かなり高度な処理が多くありましたがさらっと話していらっしゃいました。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>5. 嘉山陽一(OSGeo.JP/朝日航洋)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">QGIS日本語化の現状と課題
</span><br />
<div id="__ss_10554544" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/Yoichi_gmf/current-status-and-issuesjapaneselocalizationofqgis" target="_blank" title="QGIS日本語化の現状と課題_名古屋FOSS4G勉強会">QGIS日本語化の現状と課題_名古屋FOSS4G勉強会</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10554544" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/Yoichi_gmf" target="_blank">Yoichi Kayama</a> </div>
</div>
QGIS日本語化の第一人者の嘉山さんによる日本語化のお話。日本語化の仕組み。日本語化すべき項目(インターフェース、ヘルプ、Web)と参加方法などを丁寧にお話して頂きました。私もフリーライダーだけでなくコントリビュートしなくては…と反省しました。日本語マニュアルのLateXの処理でエラーが出て困ってらっしゃるとのことです。エンコーディングをEUC-JPに指定して(英語だけならutf-8で大丈夫だけど)、platexからの、dvipdfmxで解決するとかいう単純な問題ではないのかな?<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>6. 岩崎亘典(OSGeo.JP/農業環境技術研)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">FOSS4Gで衛星データがを使いやすくーFOSS4G を活用した衛星データ利用のためのオープン・リソースの構築ー
</span><br />
<div id="__ss_10553897" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/wata909/foss4g-10553897" target="_blank" title="FOSS4Gツール勉強会@名古屋発表資料(岩崎)">FOSS4Gツール勉強会@名古屋発表資料(岩崎)</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10553897" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/wata909" target="_blank">wata909</a> </div>
</div>
文科省プロジェクト(FOSS4G を活用した衛星データ利用のためのオープン・リソースの構築)、チュートリアル・e-learningのご紹介など。このe-learningとチュートリアルの内容は非常に豊富で(もちろん日本語)、もうオープンソースGISは商用GISに比べてドキュメントが少ないとは言えなくなりつつあると感じました。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>7. 尾崎 剛(カナエジオマチックス)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">QGISを利用したシステム開発~業務支援システム構築事例~</span><br />
<span class="Apple-style-span" style="color: #0b5394;"><br /></span><br />
QGISは多機能で、既存のGISユーザーには好評だけど、全くの初心者の方にとってはその多機能性が足かせとなるとのこと。QGISのAPIを用いた独自の実装により、ニーズに合わせた必要最小限の機能に絞ることで初心者の方にとっても使いやすくなる。開発にはPython、C++を利用されているとのことでした。やはりGeo系プログラマーにはPythonは必須なのかも。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>8. 東 修作(OSMFJ)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">OpenStreetMapの概要と利用例
</span><br />
<div id="__ss_10547453" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/higa4/ss-10547453" target="_blank" title="自由な地図を作るオープンストリートマップの活動">自由な地図を作るオープンストリートマップの活動</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10547453" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/higa4" target="_blank">Shu Higashi</a> </div>
</div>
OpenStreetMapの概要や活用例の紹介。バリアフリーマップを始め様々な活用例を紹介して頂きました。海外の大学ではOpenStreetMapとpgRoutingを使った、バリアフリーな経路検索サイトを実現しているところもあるそうです。私自身もOpenStreetMapマッパーなのですが(←最近サボリ気味です、ごめんなさい)、OpenStreetMapの多くの可能性を感じられました。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>9. 東 修作(OSMFJ)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">10分で始められるCrowdMap/Ushahidiのデプロイ
</span><br />
<div id="__ss_10547456" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/higa4/crowd-map" target="_blank" title="10分でできるCrowd mapのデプロイ">10分でできるCrowd mapのデプロイ</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10547456" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/higa4" target="_blank">Shu Higashi</a> </div>
</div>
Ushahidiによるマッピングの実装のお話。クラウドのマップサイトの作成は専門家だけのものでなく、CroudMap.comを使えば素人の方でも簡単にできるそうです。今回WiMaXの調子が悪くデモンストレーションができなかったのは残念でしたが(ネット環境の事前準備ができていなかった私の責任です。申し訳ありません(泣))。このようなマッピングサービスに慣れておくのはGeo系研究者としては必須なのかもと感じました。<br />
<br />
<br />
<br />
<span class="Apple-style-span" style="color: #0b5394;"><b>10. 縫村崇行(名古屋大学)</b></span><br />
<span class="Apple-style-span" style="color: #0b5394;">FOSS4Gツール(GRASS、R、GMT)を用いた氷河の解析事例
</span><br />
<div id="__ss_10546545" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/nuimura/111211-foss4g-nagoyapresentation2" target="_blank" title="111211 foss4g nagoya_presentation2">111211 foss4g nagoya_presentation2</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/10546545" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/nuimura" target="_blank">Takayuki Nuimura</a> </div>
</div>
私の本業のお話です。専門的な解釈や実は複雑なアルゴリズムについては完全にはしょり、FOSS4Gツールをどんな使い方をしたのかに焦点をあててお話ししました。実はいまリバイス中の論文の内容(使用している図は投稿版とレイアウトなど変えてあります)です。早くリバイス終わらせなければ(汗)。空間分布図はGMTで全て作成しています。<br />
<br />
<br />
<br />
次回はGMTのお話とかしようかなぁ。
次回の勉強会は民間企業の繁忙期がすぎる来年4月ごろを考えています。
運営に協力してくれる方、絶賛募集中です。
<br />
<br />Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com1tag:blogger.com,1999:blog-7751709962087367040.post-832053861880057912011-08-26T16:24:00.000+09:002011-08-26T17:41:17.104+09:00Nagoya.R#6で使用するファイルのリンク8/27 (土)のNagoya.R#6で使用するファイルのリンク<br />
<br />
<br />
<ul>
<li><a href="http://dl.dropbox.com/u/870568/nagoyar6/csv2kml.r">Rスクリプトファイル</a></li>
<li><a href="http://dl.dropbox.com/u/870568/nagoyar6/eqs7day-M2.5.txt">地震の CSV データ: USGS (c)</a></li>
<li><a href="http://dl.dropbox.com/u/870568/nagoyar6/eqs7day-M2.5.kml">変換後の KML データ : USGS (c)</a></li>
<li><a href="http://dl.dropbox.com/u/870568/nagoyar6/110827_nagoyar_presentation.pdf">スライドファイル(高解像度版)</a></li>
</ul>
<div>
<br /></div>
<div id="__ss_9019114" style="width: 425px;">
<strong style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/nuimura/110827-nagoyar-presentationslideshare-9019114" target="_blank" title="110827 nagoyar presentation_slideshare">110827 nagoyar presentation_slideshare</a></strong> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/9019114" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">
View more <a href="http://www.slideshare.net/" target="_blank">presentations</a> from <a href="http://www.slideshare.net/nuimura" target="_blank">Takayuki Nuimura</a> </div>
</div>
Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-81901456897245252392011-06-15T16:11:00.000+09:002011-06-15T16:11:06.632+09:00研究室でArcGIS講習会をやりました研究室のGISをこれまで使ったことのないメンバーを対象に、ArcGIS講習会を行いました。<br />
<br />
<br />
こちらが講習会で使用したスライドです。<br />
<div id="__ss_8300181" style="width: 425px;"><b style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/nuimura/arcgis" title="研究室内でのArcGIS講習会の資料">研究室内でのArcGIS講習会の資料</a></b> <iframe frameborder="0" height="355" marginheight="0" marginwidth="0" scrolling="no" src="http://www.slideshare.net/slideshow/embed_code/8300181" width="425"></iframe> <br />
<div style="padding: 5px 0 12px;">View more <a href="http://www.slideshare.net/">presentations</a> from <a href="http://www.slideshare.net/nuimura">Takayuki Nuimura</a> </div></div><br />
約5年前にこの研究室に来たときは、GISユーザーは僕だけでしたが、現在、研究室のメンバーのArcGIS使用率は14名中9名になってしまいました。<br />
個人的にはQGISやGRASS GISなどのFOSS4Gツールを広めてベンダーロックインから抜け出したいのですが、まだまだ難しそうです。<br />
<br />
ニーズがあればQGISやGRASS GIS講習会のほうをやってみたい。Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-34794072319377651672011-04-04T12:51:00.000+09:002011-04-04T12:51:02.040+09:00Nagoya.R#5に参加3/26にプログラミング言語「R」の勉強会のNagoya.R#5に参加してきました。<br />
Nagoya.R#4に続いて2回目の参加で、発表は今回が初めて。<br />
<br />
発表テーマはR言語のrgdalというパッケージを使ったリモートセンシング解析の話。<br />
参加者は言語研究をされている方がほとんどで、自然科学系はわずか(自分ともう一人)でした。専門外の方にプレゼンする貴重な機会だったので、研究のこむずかしい話は減らして、R言語での処理に重点をおいて話してみました。<br />
<br />
GISやGPSに関して興味を持ってくださった方もいて、質問を結構してもらえたのは発表してみてよかったと思いました。<br />
<br />
<br />
発表に使用したスライド<br />
<div id="__ss_7383482" style="width: 425px;"><b style="display: block; margin: 12px 0 4px;"><a href="http://www.slideshare.net/nuimura/110326-nagoyar-presentationslidshare" title="110326 nagoyar presentation_slidshare">110326 nagoyar presentation_slidshare</a></b> <object height="355" id="__sse7383482" width="425"> <param name="movie" value="http://static.slidesharecdn.com/swf/ssplayer2.swf?doc=110326nagoyarpresentationslidshare-110325041032-phpapp01&stripped_title=110326-nagoyar-presentationslidshare&userName=nuimura" /><param name="allowFullScreen" value="true"/><param name="allowScriptAccess" value="always"/><embed name="__sse7383482" src="http://static.slidesharecdn.com/swf/ssplayer2.swf?doc=110326nagoyarpresentationslidshare-110325041032-phpapp01&stripped_title=110326-nagoyar-presentationslidshare&userName=nuimura" type="application/x-shockwave-flash" allowscriptaccess="always" allowfullscreen="true" width="425" height="355"></embed> </object> <br />
<div style="padding: 5px 0 12px;">View more <a href="http://www.slideshare.net/">presentations</a> from <a href="http://www.slideshare.net/nuimura">Takayuki Nuimura</a> </div></div><br />
SlideshareにPDFを直接アップロードすると、日本語の部分が文字化けor非表示になってしまったため、PDFを画像化してアップロードしています。Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-46465560926375364132011-02-11T20:42:00.000+09:002011-02-11T20:42:14.274+09:00spatiohack平和公園1万歩コースをジョギング、コースの外は雪が積もっているところも見られた。<br /><br /><iframe src="http://www.jogtracker.com/HistoryDetailsEmbed.aspx?SyncID=309285d3-b79f-4337-9a3c-c5540f5593cc&UserName=nekogahora" style="width: 590px; height: 600px;" frameborder="0" scrolling="no">Embedding not supported. View jog details at &amp;amp;lt;a href="http://www.jogtracker.com/HistoryDetails.aspx?SyncID=309285d3-b79f-4337-9a3c-c5540f5593cc&amp;amp;amp;UserName=nekogahora"&amp;amp;gt;JogTracker&amp;amp;lt;/a&amp;amp;gt;</iframe><br /><br /><a href="http://goo.gl/photos/S5aGv0Zkbd" imageanchor="1" style="clear: right; margin-bottom: 1em; margin-left: 1em;"><img src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjR0rE3d-lc6yQdpMQgf-nvZHmvqKP9nKeF_j23QK1tVYJxm_AyY6FSjD81kZCcApCFW8obU2XEJmy0DN3Q4HWAGwHwojzI_GvPctHuaoQiMtEmTmSFeEFXzsYgRTOHEkP8VAI4mOM8ukC_/s512/-405828789.jpg" border="0" /></a>Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0tag:blogger.com,1999:blog-7751709962087367040.post-6932044138933021702011-02-02T18:05:00.000+09:002011-02-02T18:05:42.141+09:00latex2rtfで変換するときの注意点texソースコードをword形式に変換する際のやり方の一つとして、以前紹介したlatex2rtfでrtf(リッチテキスト)形式を経由する方法がある。<br />
<br />
<a href="http://spatiohack.blogspot.com/2010/12/latex_14.html">LaTeX原稿のスペルチェック (Dec. 14, 2010)</a><br />
<br />
<br />
latex2rtfで変換する際には、\graphicpathコマンドで画像のパスを指定すると、出力されるrtfファイルに画像が含まれない。<br />
<br />
\graphicpathコマンドを使うと、画像を同じ名前のディレクトリ(例.image)に入れている場合は、\includegraphics{image/filename1.png}<br />
\includegraphics{image/filename2.png}<br />
\includegraphics{image/filename3.png}<br />
と書かずに、<br />
<br />
\graphicpath{image/}<br />
\includegraphics{filename1.png}<br />
\includegraphics{filename2.png}<br />
\includegraphics{filename3.png}<br />
とすることができるようになって便利だが、このやり方ではlatex2rtfコマンドで図を読み込んでくれないらしい。<br />
<br />
<br />
texlipseでlatexを書いている場合のコマンド例(texファイル名がdocument.texの場合)は以下の通り。<br />
<br />
latex2rtf -a ./tmp/document.aux -b ./tmp/document.bbl -o output.rtf document.tex<br />
<br />
linenoやmultirowなど一部の環境はサポートされていないが、文献情報や図の取り込みはうまく出来ている。<br />
また、見出し情報もちゃんと含まれている。<br />
<br />
ただ、Ubuntu上でlatex2rtfを使った場合は大丈夫だったが、Windows上のw32texに含まれているlatex2rtfでは画像がページをはみ出すぐらい大きくなっていた。<br />
<br />
画像形式がpngなのに、bb情報を含めなかったためかもしれない。Takayuki Nuimurahttp://www.blogger.com/profile/05195923173300610556noreply@blogger.com0