スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

岩国~青海島

install.packages("XML")
library(XML)
doc <- xmlParse("06Apr2014.gpx")
r <- xmlRoot(doc)

library(maptools)
library(shapefiles)
pref <- 35
shapefilename <- paste("F:/GMLdata/行政区域/N03-140401_", pref, "_GML/N03-14_", pref, "_140401", sep="")
datashp <- readShapePoly(shapefilename,proj4string=CRS("+proj=longlat"))
pdf("map16Aug2014.pdf",height=5.5,width=6)
plot(datashp, col="#F1E6D0", bg="#A8D7FF", lwd=0.025)
dev.off()
for(i in c(2:15)){
data <- unlist(xmlToList(r[[i]][[2]]))
lat <- data[seq(8,length(data),9)]
lon <- data[seq(9,length(data),9)]
lines(lon,lat, col="red", lwd=2)
}

MapTools
・市町村名を入れるにはどうしたらよいか
・縦横比が変わる件
スポンサーサイト
Profile

zoa

Author:zoa
Metropolis (1927)

Calendar
07 | 2014/08 | 09
- - - - - 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 - - - - - -
Labels
Admin
Previous Posts
Recent Comments
Recent Trackback
Archive
このページのトップへ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。