27

注意：这个问题是针对映射的，但我希望在标准笛卡尔坐标系中绘图时能够使用它。

我喜欢基本图形，但也喜欢 ggplot2 在很多方面。我最常用的微调图形的基本函数之一是 locator(n) ，但这会在 ggplot2 中产生错误。

library(ggplot2) 
county_df <- map_data('county')  #mappings of counties by state
ny <- subset(county_df, region=="new york")   #subset just for NYS
ny$county <- ny$subregion

ggplot(ny, aes(long, lat, group=group)) +  geom_polygon(colour='black', fill=NA)
locator(1)

现在grid.locator()正如 Dason 在talkstats.com 上向我指出的那样（这里）可以返回一些东西。我只是不知道如何使用那个东西来获取地图坐标。

> grid.locator()
$x
[1] 286native

$y
[1] 133native

单位似乎没有帮助，因为它们不是地图坐标。也许我需要某种转换。

先感谢您。

编辑：（基于 DWin 的回应）

Dwin 的想法是正确的，但转换系数有点偏离。帮助将不胜感激。在下面的示例中，我在坐标 (x = -73 & y = 40.855) 处有一个带有红点的地图。我将 Dwin 的响应放入一个返回坐标的函数中。我希望结果是我输入的坐标，但事实并非如此。

想法？

require(maps); library(ggplot2); require(grid)

county_df <- map_data('county')  #mappings of counties by state
ny <- subset(county_df, region=="new york")   #subset just for NYS
ny$county <- ny$subregion


NY <- ggplot(ny, aes(long, lat)) +  
          geom_polygon(aes(group=group), colour='black', fill=NA) +
          coord_map() + geom_point(aes(-73, 40.855, colour="red"))
NY  

gglocator <- function(object){
    require(maps); require(grid)
    z <- grid.locator("npc")
    y <- sapply(z, function(x) as.numeric(substring(x, 1, nchar(x))))
    locatedX <- min(object$data$long) + y[1]*diff(range(object$data$long))
    locatedy <- min(object$data$lat)  + y[2]*diff(range(object$data$lat))
    return(c(locatedX, locatedy))
}

#click on the red dot
gglocator(NY)  #I expect the results to be x = -73 & y = 40.855

编辑 2：（脱离 Baptise 的回答）

在那里

NY <- ggplot(ny, aes(long, lat)) +  
          geom_polygon(aes(group=group), colour='black', fill=NA) +
          coord_map() + geom_point(aes(-73, 40.855, colour="red")) +
          scale_x_continuous(expand=c(0,0)) + scale_y_continuous(expand=c(0,0))


NY 
x <- grid.ls()[[1]][grep("panel-", grid.ls()[[1]])] #locate the panel
seekViewport(x)
y <-  grid.locator("npc")
y <- as.numeric(substring(y, 1, nchar(y)-3))

locatedX <- min(NY$data$long) + y[1]*diff(range(NY$data$long))
locatedy <- min(NY$data$lat) + y[2]*diff(range(NY$data$lat))
locatedX; locatedy

更新： ggmap 包的gglocator功能现在包含此功能。

4

5 回答 5

10

需要使用有意义的单位系统并将信息保存在 ggplot 对象中，以便您可以从“npc”单位转换为地图单位：

require(maps)
require(grid)
NY <- ggplot(ny, aes(long, lat, group=group)) +  geom_polygon(colour='black', fill=NA)
 grid.locator("npc")
# clicked in middle of NY State:

#$x
#[1] 0.493649231346082npc
#
#$y
#[1] 0.556430446194226npc
 range(NY$data$long)
#[1] -79.76718 -71.87756
 range(NY$data$lat)
#[1] 40.48520 45.01157
 locatedX <- min(NY$data$long) + 0.493649231346082*diff(range(NY$data$long))
 locatedX
#[1] -75.87247
locatedY <- min(NY$data$lat) +  0.556430446194226*diff(range(NY$data$lat))
locatedY
#[1] 43.00381
于 2012-02-26T17:54:40.027 回答
7

如果我添加scale_x_continuous(expand=c(0,0)) + scale_y_continuous(expand=c(0,0))到情节中，我会得到正确的结果，seekViewport("panel-3-4")之前grid.locator()

于 2012-02-26T20:08:40.927 回答
5

我写信给 ggplot 帮助列表，并收到了来自 David Kahle 的非常有帮助的回复，他恰好对同样的问题感兴趣。他的功能很棒：

1）您不必在绘图中添加比例 y 和比例 x

2）它可以一次找到多个点并将它们作为数据框返回

3）它适用于任何类型的ggplot，而不仅仅是地图

gglocator <- function(n = 1, object = last_plot(), 
    message = FALSE, xexpand = c(.05,0), yexpand = c(.05, 0)){ 

  #compliments of David Kahle
  if(n > 1){
    df <- NULL
    for(k in 1:n){
      df <- rbind(df, gglocator(object = object, message = message, 
        xexpand = xexpand, yexpand = yexpand))
    }
    return(df)
  }

  x <- grid.ls(print = message)[[1]]
  x <- x[ grep("panel-", grid.ls(print=message)[[1]]) ] #locate the panel
  seekViewport(x)
  loc <-  as.numeric(grid.locator("npc"))

  xrng <- with(object, range(data[,deparse(mapping$x)]))
  yrng <- with(object, range(data[,deparse(mapping$y)]))    

  xrng <- expand_range(range = xrng, mul = xexpand[1], add = xexpand[2])
  yrng <- expand_range(range = yrng, mul = yexpand[1], add = yexpand[2])    

  point <- data.frame(xrng[1] + loc[1]*diff(xrng), yrng[1] + loc[2]*diff(yrng))
  names(point) <- with(object, c(deparse(mapping$x), deparse(mapping$y)))
  point
}

#Example 1
require(maps); library(ggplot2); require(grid)
county_df <- map_data('county')  #mappings of counties by state
ny <- subset(county_df, region=="new york")   #subset just for NYS
ny$county <- ny$subregion


NY <- ggplot(ny, aes(long, lat)) +  
          geom_polygon(aes(group=group), colour='black', fill=NA) +
          coord_map() + geom_point(aes(c(-78, -73), c(41, 40.855), 
          colour=c("blue", "red"))) + opts(legend.position = "none") 


NY 
gglocator(2)

#Example 2
df <- data.frame(xvar = 2:10, yvar = 2:10)
ggplot(df, aes(xvar, yvar)) + geom_point() + geom_point(aes(x = 3, y = 6))
gglocator()

更新： ggmap 包的gglocator功能现在包含此功能。

于 2012-02-27T18:44:51.147 回答
5

这是使用 DWin 和 Baptise 给我的所有东西打包成一个函数的最终结果。我还查询了 ggplot 帮助列表，并将在此处报告和其他信息。

require(maps); require(ggplot2); require(grid)

ny <- map_data('county', 'new york')

NY1 <- ggplot(ny, aes(long, lat)) +  
          geom_polygon(aes(group=group), colour='black', fill=NA) +
          coord_map() + geom_point(aes(c(-78, -73), c(41, 40.855), 
          colour=c("blue", "red"))) + opts(legend.position = "none") 

NY <- NY1 + scale_x_continuous(expand=c(0,0)) + 
          scale_y_continuous(expand=c(0,0))
          #the scale x and y have to be added to the plot

NY 

ggmap.loc <- function(object){
    x <- grid.ls()[[1]][grep("panel-", grid.ls()[[1]])] #locate the panel
    seekViewport(x)
    y <-  as.numeric(grid.locator("npc"))
    locatedX <- min(object$data$long) + y[1]*diff(range(object$data$long))
    locatedy <- min(object$data$lat) + y[2]*diff(range(object$data$lat))
    return(c(locatedX, locatedy))
}

ggmap.loc(NY)
于 2012-02-27T06:36:09.447 回答
3

这些帖子非常有帮助，但是已经有几年了，所以有些东西坏了。这是一些对我有用的新代码。找到正确视口的代码不起作用，所以我过去current.vpTree()手动搜索正确的视口，然后将其复制到seekViewport(). 请注意我的视口是'panel.6-4-6-4'而不是旧样式panel-*。最后，在 rstudio 中渲染时我没有得到正确的答案，而是我不得不使用x11().

这是一个完整的例子。希望这会有所帮助。

library(ggplot2)
library(grid)

object<-ggplot(dat=data.frame(x=1:5,y=1:5),aes(x=x,y=y)) + 
  geom_point()  +
  scale_x_continuous(expand=c(0,0)) +
  scale_y_continuous(expand=c(0,0))
x11()
print(object)
formatVPTree(current.vpTree()) #https://www.stat.auckland.ac.nz/~paul/useR2015-grid/formatVPTree.R
seekViewport('panel.6-4-6-4')
y <-  as.numeric(grid.locator("npc"))
locatedX <- min(object$data$x) + y[1]*diff(range(object$data$x))
locatedY <- min(object$data$y) + y[2]*diff(range(object$data$y))
locatedX; locatedY
于 2017-09-22T19:51:38.553 回答