ミネソタ州のすべての郡に関するデータ セットがあり、列の 1 つはその形状です。各郡については、次のようになります。
エイトキン郡の場合:
<Polygon><outerBoundaryIs><LinearRing><coordinates>-93.051956,46.15767700000001,0 -93.434006,46.15313,0 -93.43261,46.240253,0 -93.80480900000001,46.23817100000001,0 -93.80933400000001,46.580681,0 -93.77426199999999,46.59050400000001,0 -93.77412400000001,46.802605,0 -93.77500100000002,47.030445,0 -93.058258,47.022362,0 -93.05964600000001,46.766071,0 -93.05208600000002,46.417576,0 -93.051956,46.15767700000001,0</coordinates></LinearRing></outerBoundaryIs></Polygon>
私はRを初めて使用し、Google API、HTMLなどについて何も知りません.ggplot2およびmapsパッケージを使用して、ミネソタ州のすべての郡のさまざまな側面の強度マップを作成しようとしています. これらの座標をそのまま使用して郡のレイヤーを作成する方法はありますか、それとも何か他のことをする必要がありますか?
これまでのコードは次のとおりです。
ミネソタ州の地図:
library(maps)
library(ggplot2)
all_states <- map_data("state")
mn<-subset(all_states, region %in% c("minnesota"))
p<-ggplot()
p<-p+geom_polygon(data=mn, aes(x=long, y=lat, group=group), colour="black", fill="white")
p
そして、私の計画は、これらのポリゴンを取得したら、次を変更して各郡に適用することです。
dataset <- data.frame(region=states,val=runif(49, 0,1))
us_state_map <- map_data('state')
map_data <- merge(us_state_map, dataset, by='region', all=T)
map_data <- map_data[order(map_data$order), ]
(qplot(long, lat, data=map_data, geom="polygon", group=group, fill=val)
+ theme_bw() + labs(x="", y="", fill="")
+ scale_fill_gradient(low='#EEEEEE', high='darkgreen')
+ opts(title="Title",
legend.position="bottom", legend.direction="horizontal"))
どんな提案でも大歓迎です！