3

KML データ ソースから ggplot2 を使用してコロプレスまたはテーマ マップをプロットするにはどうすればよいですか?

KML の例: https://dl.dropbox.com/u/1156404/nhs_pct.kml

データ例: https://dl.dropbox.com/u/1156404/nhs_dent_stat_pct.csv

ここに私がこれまでに持っているものがあります：

install.packages("rgdal")
library(rgdal)
library(ggplot2)

fn='nhs_pct.kml'

#Look up the list of layers
ogrListLayers(fn)

#The KML file was originally grabbed from Google Fusion Tables
#There's only one layer...but we still need to identify it
kml=readOGR(fn,layer='Fusiontables folder')

#This seems to work for plotting boundaries:
plot(kml)

#And this:
kk=fortify(kml)
ggplot(kk, aes(x=long, y=lat,group=group))+ geom_polygon()

#Add some data into the mix
nhs <- read.csv("nhs_dent_stat_pct.csv")

kml@data=merge(kml@data,nhs,by.x='Name',by.y='PCT.ONS.CODE')

#I think I can plot against this data using plot()?
plot(kml,col=gray(kml@data$A.30.Sep.2012/100))
#But is that actually doing what I think it's doing?!
#And if so, how can experiment using other colour palettes?

#But the real question is: HOW DO I DO COLOUR PLOTS USING gggplot?
ggplot(kk, aes(x=long, y=lat,group=group)) #+ ????

だから私の質問は: たとえば kml@data$A.30.Sep.2012 の値を使用して地域に色を付けるにはどうすればよいですか?

補足的な質問として、もう一度 ggplot のコンテキストで、さまざまなカラー パレットをどのように試すことができますか?

4

1 に答える 1

5

R でマップをプロットするのは非常に面倒です。https://github.com/hadley/ggplot2/wiki/plotting-polygon-shapefilesの Hadley のチュートリアルに大きく従う回答を次に示します。

library(maptools)
library(rgdal)
library(ggplot2)
library(plyr)

fn='nhs_pct.kml'
nhs <- read.csv("nhs_dent_stat_pct.csv")

kml <- readOGR(fn, layer="Fusiontables folder")

注: オーファン ホールに関するメッセージを受け取りました。https://stat.ethz.ch/pipermail/r-help/2011-July/283281.htmlを読んだ後、次の行を含めました

slot(kml, "polygons") <- lapply(slot(kml, "polygons"), checkPolygonsHoles)

## The rest more or less follows Hadley's tutorial 
kml.points = fortify(kml, region="Name")
kml.df = merge(kml.points, kml@data, by.x="id",by.y="Name",sort=FALSE)
kml.df <- merge(kml.df,nhs,by.x="id",by.y="PCT.ONS.CODE",sort=FALSE,all.x=T,all.y=F)

## Order matters for geom_path!
kml.df <- kml.df[order(kml.df$order),]

nhs.plot <- ggplot(kml.df, aes(long,lat,group=group,fill=A.30.Sep.2012)) + 
  geom_polygon() +
  geom_path(color="gray") +
  coord_equal() +
  scale_fill_gradient("The outcome") + 
    scale_x_continuous("") + scale_y_continuous("") +   theme_bw()
于 2012-12-07T20:51:20.563 に答える