R：在地图上用颜色强度绘制点答案

【问题标题】：R: draw points on map with colour intensityR：在地图上用颜色强度绘制点
【发布时间】：2015-04-27 19:02:32
【问题描述】：

我在 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 绘制的地图确实过于近似：我需要为每种颜色设置自己的范围。

【问题讨论】：

Here's what I came up with 作为类似问题的答案。
您的问题并没有真正提供足够的信息来回答。如果您已经有一个包含计数的列，您可以使用geom_point(data=originale, aes(x=lon, y=lat, color = counts))。你可以控制scale_color_manual使用的颜色。
对不起@Ista，我正在尝试更具体。数据集的有趣部分是这样制作的：'number longitude latitude; 1 45.654 13.645； 2 44.876 12.987； 3 45.654 13.645。如您所见，几行具有相同的坐标，因此我没有计数。我认为没有计数以便每次都对我想要的数据帧部分进行子集化。但是根据您的观察，我会尝试这种方式。
很好，请使用此附加信息更新问题，包括数据样本。

标签： r google-maps plot cartodb

【解决方案1】：

我找到了一个可行的解决方案，计算新数据框中每个观察的重复次数，然后使用 geom_point 结合 scale_colour_discrete 和中断来绘制地图。

感谢乔希和伊斯塔的建议！

代码如下：

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)

【讨论】：