0

何千もの地理的に局在化された観測データを含む R のデータセットがあり、それらをマップにプロットする必要があります。私はこのコードを使用して単一の点をプロットすることができました:

originale<-read.table("file.txt", header=TRUE,sep=";")
require(ggplot2)
require(ggmap)
map <- get_map(location = c(lon=13.781693, lat=45.623124), zoom = 14, maptype = "terrain",source = "google")
p <- ggmap(map)
p_punti <- p + geom_point(data=originale, aes(x=lon, y=lat),size=5)
plot(p_punti)

ここで、そのポイントにある観測の数に基づいて、5 ～ 6 色のセットを使用してそれらをプロットしたいと思います。CartoDB で作成したこのようなもの: http://i.stack.imgur.com/Be1DL.png

実際、CartoDB でプロットされたマップはおおまかすぎます。色ごとに独自の範囲を設定する必要があります。

4

1 に答える 1

0

新しいデータフレームの観測ごとに繰り返しをカウントし、geom_point と scale_colour_discrete および break を組み合わせてマップをプロットする実用的なソリューションを見つけました。

提案してくれた Josh と Ista に感謝します。

コードは次のとおりです。

require(plyr)
df <-read.table("input.txt", header=TRUE,sep=";")

##selection
original <- df[(df$Year==2013),] #

##count recurrence 
conteggi2 <- ddply(original, .(original$lon, original$lat), nrow)
names(conteggi2) <- c("lon", "lat", "Freq")

##sort
conteggi2 <- conteggi2[with(conteggi2, order(Freq)), ]

##Classification
verdi <- conteggi2$Freq <= 5
conteggi2$classe[verdi] <- "1-5"
gialli <- conteggi2$Freq > 5 & conteggi2$Freq <= 20
conteggi2$classe[gialli] <- "6-20"
rossi <- conteggi2$Freq > 20 & conteggi2$Freq <= 50
conteggi2$classe[rossi] <- "21-50"
neri <- conteggi2$Freq > 50 & conteggi2$Freq <= 500
conteggi2$classe[neri] <- "51-500"

##plot map
require(ggplot2)
require(ggmap)
library(reshape2)
Trieste <- qmap(location = c(lon=13.781693, lat=45.623124), zoom = 13,
               legend = "topleft", color="bw")#maptype = "terrain",
Trieste1 <- Trieste +
        geom_point(aes(x = lon, y = lat, colour = classe, 
                       title='Year 2013'),
                        data = conteggi2) +
        guides(size=FALSE) + #remove useless label
        scale_colour_discrete(name="Legend Title",
                              breaks=c("1-5","6-20","21-50", "51-500"))
plot(Trieste1)
于 2015-03-02T20:58:31.620 に答える