convert to kml - グリッドデータや測点データをkmlファイルにする：GoogleEarthに持って行く　

---

  　使うコマンド　 makecpt, grd2kml, (gmt)2kml



bash スクリプト例
#  convert data files to KML

# convert grid data to kml
grdfile=JTOPO1_30.grd                   prefix=JTOPO
cptfile=haxby_grad.cpt                  #
gmt makecpt -Chaxby -T-8000/4000/1000 -Z > $cptfile
gmt grd2kml $grdfile -N$prefix -C$cptfile
#
# convert table data to kml
ifile=DSDP.list
kmlfile=DSDP_SB.kml
#
gmt 2kml $ifile -Ag0 -Gblue+f -Fs -N > $kmlfile


Tips

• スクリプト前半は，グリッドファイルをカラーテーブルを指定してkmlに変換している．入力のグリッドファイルのヘッダーがcartesian coordinateとなっているとエラーになる．nearneighborやsurfaceなどで-fgオプションをつけずにグリッドを作成すると，境界やグリッド間隔は一見geographicで，grdimage等の図化する場合はなんの問題もないのだが，ヘッダーに[Cartesian Coordinate]と記載されてしまい，それがエラーになる．ヘッダはgmt grdinfoで確認できる．-fgなしで作成したグリッドの場合は，grdeditでヘッダだけ書き換えればOK.　
• 　スクリプト後半は，点データをkmlに変換している．-Nオプションをつけているので各点にラベルが付いていない．細かい指定はマニュアル参照のこと．