以前、Google Earthで図を描いたが、今回はアニメーションを作ってみる。図が3枚だけのアニメーションだ。
手順1
Gradsで図を3枚つくる。GradsでGoogle Earth向きの図をつくるやり方は以前の記事を参照のこと。下の例はNCEP/NCARの最解析データを使って、2008年から2010年の8月の表面大気温度のアノマリをプロットしている。colorbarはxcbarを使って図の上に一緒に乗せた。
"reinit" "sdfopen http://www.esrl.noaa.gov/psd/thredds/dodsC//Datasets/ncep.reanalysis.de rived/surface/air.mon.mean.nc" "sdfopen http://www.esrl.noaa.gov/psd/thredds/dodsC//Datasets/ncep.reanalysis.de rived/surface/air.mon.ltm.nc" year=2008 while (year<=2010) "set time aug"year "set x 1 144" "set y 1 73" "set parea 0 11 0 8.5" "set grid off" "set grads off" "set mproj scaled" "set mpdraw off" "set gxout shaded" "color -5 5 0.5 -kind blue->gray->red" "d air.1-air.2(t=8)" "xcbar 4.5 6.5 4 4.2 -edge triangle -fstep 10 " "printim airtemp_"year"aug.png x2048 y1024 -t 0" "c" year=year+1 endwhile
こんな感じの図ができる。
手順2
基になるkmlファイルをGradsで作成する。やり方は以前の記事と同じだ。
*************KML "reinit" "sdfopen http://www.esrl.noaa.gov/psd/thredds/dodsC//Datasets/ncep.reanalysis.de rived/surface/air.mon.mean.nc" "set time aug2010" "set x 1 144" "set y 1 73" "set kml airtemp" "set gxout kml" "d air"
以下のようなkmlファイルが出来る。
<?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://www.opengis.net/kml/2.2"> <GroundOverlay> <name>Monthly Mean Air Temperature</name> <Icon> <href>airtemp.tif</href> </Icon> <LatLonBox> <west> -1.25</west> <east> 358.75</east> <south> -91.25</south> <north> 91.25</north> <rotation>0.0</rotation> </LatLonBox> </GroundOverlay> </kml>
手順3
kmlファイルを以下の様に書き換える。
ポイントは
・GroundOverlayから/GraoundOverlayの部分を図の数だけ繰り返す。
・それぞれ図の名前を指定する。
・図がどの時間を占めるべきかTimeSpan /TimeSpanの部分で指定する。なお、今の場合、図は8月だけの図であるが、アニメーションで図が続く様に1年中同じ図が表示されるよう、1月1日から12月31日をTimeSpanで設定している。
・繰り返しの全体をFolder と/Folderではさむ。
<?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://www.opengis.net/kml/2.2"> <Folder> <GroundOverlay> <name>Monthly Longterm Mean of Air Temperature</name> <TimeSpan> <begin>2008-01-01</begin> <end> 2008-12-31</end> </TimeSpan> <Icon> <href>airtemp_2008aug.png</href> </Icon> <LatLonBox> <west> -1.25</west> <east> 358.75</east> <south> -91.25</south> <north> 91.25</north> <rotation>0.0</rotation> </LatLonBox> </GroundOverlay> <GroundOverlay> <name>Monthly Longterm Mean of Air Temperature</name> <TimeSpan> <begin>2009-01-01</begin> <end> 2009-12-31</end> </TimeSpan> <Icon> <href>airtemp_2009aug.png</href> </Icon> <LatLonBox> <west> -1.25</west> <east> 358.75</east> <south> -91.25</south> <north> 91.25</north> <rotation>0.0</rotation> </LatLonBox> </GroundOverlay> <GroundOverlay> <name>Monthly Longterm Mean of Air Temperature</name> <TimeSpan> <begin>2010-01-01</begin> <end> 2010-12-31</end> </TimeSpan> <Icon> <href>airtemp_2010aug.png</href> </Icon> <LatLonBox> <west> -1.25</west> <east> 358.75</east> <south> -91.25</south> <north> 91.25</north> <rotation>0.0</rotation> </LatLonBox> </GroundOverlay> </Folder> </kml>
出来たファイルをkmzにまとめたものをこちらに置いた。kmzファイルはkmlと図をまとめたもので、Google Earthで一度kmlファイルを開いて、「名前をつけて保存」で作成できる。
Google Earthで開くとこんな感じ。
0 件のコメント:
コメントを投稿