後方流跡線をRで描いてみる・続

後方流跡線のデータを読み込んだら

あとはこっちのもの（なにが？）

地図を描いて、後方流跡線を乗っけるだけ。

今回、地図はNature Earthのものを使用しました。

荒めの国境データですが、今回にはこれで十分かと。

shapeファイルを読み込んで、ggplot2で描けるように

変換しました。

library(ggplot2)

library(maptools)

worldCountry10_shp <-  readShapeLines("ne_10m_admin_0_countries.shp")
worldCountry10_map <- fortify(worldCountry10_shp)

ggmap1 <-
  ggplot() +
  geom_path(data = worldCountry10_map, aes(x = long, y = lat, group = group),                                    colour="#999999") +
  xlim(c(80, 150)) + ylim(c(20, 60)) +
  coord_fixed(ratio=1.35/1)

coord_fixedは目視で設定した、個人的な趣味の数字です。

正解は、何なんだろう？

あと、中国は土地勘が無いので適当な都市を地図で表示。

pointdata <- read.csv("pointdata.csv", header=T)

> pointdata
name lat lon
1 Beijing 39.90449 116.3915
2 Shanghai 31.24787 121.4727
3 Hong_Kong 22.27838 114.1743
4 Ulan_Bator 47.92138 106.9055
5 Shenyang 41.80030 123.4338
6 Chongqing 29.57000 106.5800

ggmap2 <-
  ggmap1 +
  geom_point(data = pointdata,aes(x = lon, y = lat), size=4, shape=16, colour="blue") +
  geom_text(data = pointdata, aes(x= lon+1, y = lat-1, label = name), size=5)

この段階で地図を描くと、こんな感じ。

あとは、これに後方流跡線のデータを載せるだけ。

 それはまた、次に。