R拉/圆形拉伸圆（半径（r米））（epsg 4326）

Question

我想在我的数据集中创建一个大约100个点的缓冲区，其长度为r并且在点之间变化（基本上它是客户愿意前往点/商店的第80个百分点距离）

我的数据包含以下列： lat_centre，lon_centre，radius_km

这是一个1公里的例子：

使用Google地图的Javascript。

   var draw_circle = new google.maps.Circle({
        center: point_i,
        radius: 1000,
        strokeColor: col,
        strokeOpacity: 0.15,
        strokeWeight: 2,
        fillColor: col,
        fillOpacity: 0.15,
        map: map
    });

在左侧绘制图像：

但是，以下R-attempt在右侧绘制圆圈（您可以看到更小）：

library(dismo)
emory <- gmap("Bishopsgate, London", zoom = 14, scale = 2)
d <- data.frame(lat = c(51.51594), lon = c(-0.08248))
coordinates(d) <- ~ lon + lat
projection(d) <- "+init=epsg:4326"

mm <- "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"
d_mrc <- spTransform(d, CRS = CRS(mm))

# Buffer creation
d_mrc_bff <- gBuffer(d_mrc, width = 1000)
library(scales) # for `alpha()` function

plot(emory)
plot(d_mrc_bff, col = alpha("red", .35), add = TRUE)
points(d_mrc, cex = 2, pch = 20)

Answer 1

我已经编写了googleway软件包来访问Google Maps Javascript API，因此您可以复制图片。

为此，您需要有效的Google API密钥

library(googleway)

## your api key 
map_key <- "your_api_key_goes_here"

d <- data.frame(lat = c(51.51594), lon = c(-0.08248))

google_map(data = d, key = map_key, height = 800) %>%
    add_markers() %>%
    add_circles(radius = 1000)