- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我需要计算栅格质心(长纬度坐标)和引用点之间的旅行时间(例如乘车)。然而,不幸的是,许多栅格质心不在 googlemaps 的覆盖范围内,我想知道是否有一种解决方法可以让我以最快的方式到达引用点(例如结合步行到街道然后开车到引用点),如果是这样,我如何有效地做到这一点?
问题是,当我只是寻找下一个定位点时,我会产生偏差,将到该点的距离作为步行距离,然后从那里到引用点的 googlemaps 行驶距离,因为最近的点可能不是最快路线..
我试图使用 gmapsdistance-function 来做到这一点。下面是一个例子:
library(gmapsdistance)
#While this can be located and works well
> gmapsdistance(origin = "5.600451+-0.202553",
+ destination = "5.591622+-0.187677",
+ mode = "walking")
$Time
[1] 2101
$Distance
[1] 2667
$Status
[1] "OK"
#Changing the origin to other points often does not work as the points cannot be located and results in an NA-output
> gmapsdistance(origin = "7.9254948+-0.6283887",
+ destination = "5.591622+-0.187677",
+ mode = "walking")
$Time
[1] NA
$Distance
[1] NA
$Status
[1] "ROUTE_NOT_FOUND"
最佳答案
我认为解决这个问题的一种方法是获取加纳道路的形状文件并执行地理空间操作以找到最近的道路。从那里你应该能够使用谷歌的 API 来获得驾驶距离
此解决方案中的步骤是
library(sf) ## geospatial operations
library(googleway) ## plotting and google API
googleway
是我的包,您需要一个 Google API key 才能使用它
## Setting API keys I'll be using, one for the maps and one for directions
set_key("my_map_key"), api = "map")
set_key("my_other_api_key")
## Ghana roads
ghanaRoads <- sf::st_read("~/Downloads/gha_roads_dcw/GHA_rds_1m_dcw.shp")
## origin piont
df <- data.frame(lat = 7.9254948, lon = -0.6283887)
## destination point
dest <- data.frame(lat = 5.591622, lon = -0.187677)
google_map() %>%
add_markers(data = df) %>%
add_markers(data = dest) %>%
add_polylines(ghanaRoads)
## convert origin into an 'sf' object
sf_origin <- sf::st_sf(geometry = sf::st_sfc(sf::st_point(x = c(-0.6283887, 7.9254948))))
## create a line between the origin and destination
m <- matrix(c(-0.6283887, 7.9254948, -0.187677, 5.591622), ncol = 2, byrow = T)
sf_line <- sf::st_sf(geometry = sf::st_sfc(sf::st_linestring(x = m)))
## The coordinate reference system needs to match between the two for
## spatial operations
sf::st_crs(sf_line) <- sf::st_crs(ghanaRoads)
sf::st_crs(sf_origin) <- sf::st_crs(ghanaRoads)
## find all the intersecting points
sf_intersections <- sf::st_intersection(ghanaRoads, sf_line)
google_map() %>%
add_markers(data = df) %>%
add_markers(data = dest) %>%
add_polylines(data = ghanaRoads) %>%
add_markers(data = sf_intersections)
sf
找到最近的方法是使用
sf::st_distance
,但这取决于
lwgeom
已安装在您的系统上,但我在安装时遇到问题,因此我必须使用不同的方法
library(data.table)
coords <- matrix(unlist(sf_intersections$geometry), ncol = 2, byrow = T)
## Taking a fucntion I wrote for this answer
## https://stackoverflow.com/a/42014364/5977215
dt.haversine <- function(lat_from, lon_from, lat_to, lon_to, r = 6378137){
radians <- pi/180
lat_to <- lat_to * radians
lat_from <- lat_from * radians
lon_to <- lon_to * radians
lon_from <- lon_from * radians
dLat <- (lat_to - lat_from)
dLon <- (lon_to - lon_from)
a <- (sin(dLat/2)^2) + (cos(lat_from) * cos(lat_to)) * (sin(dLon/2)^2)
return(2 * atan2(sqrt(a), sqrt(1 - a)) * r)
}
dt <- as.data.table(coords)
dt[, `:=`(origin_lon = df$lon, origin_lat = df$lat)]
dt[, distance := dt.haversine(origin_lat, origin_lon, V1, V2)]
## min distance
sf_nearest <- dt[order(-distance)][1, .(lon = V1, lat = V2)]
sf_nearest <- sf::st_point(c(sf_nearest$lon, sf_nearest$lat))
sf_nearest <- sf::st_sf(geometry = sf::st_sfc(sf_nearest))
sf_nearest$colour <- "green"
google_map() %>%
add_markers(data = df) %>%
add_markers(data = dest) %>%
add_markers(data = sf_nearest, colour = "colour")
orig <- sf_nearest$geometry[[1]]
orig <- as.numeric(orig)
df_orig <- data.frame(lat = orig[2], lon = orig[1])
google_map() %>%
add_markers(df_orig)
res <- google_directions(origin = df_orig,
destination = dest)
## all the api results are now stored in the 'res' object.
direction_legs(res)$distance
# text value
# 1 397 km 396829
## you can look at the route through the polygoin
df_route <- data.frame(polyline = direction_polyline(res))
google_map() %>%
add_markers(data = df_orig) %>%
add_markers(data = dest) %>%
add_polylines(data = df_route, polyline = "polyline")
dt[order(-distance)]
给我们从原点到我们新原点的距离,
dt[order(-distance)][1, distance]
# [1] 1329904
sf_buffer <- sf::st_buffer(sf_origin, dist = 0.5)
sf::st_crs(sf_buffer) <- sf::st_crs(ghanaRoads)
google_map() %>%
add_polylines(ghanaRoads) %>%
add_polygons(sf_buffer)
sf_intersection <- sf::st_intersection(sf_buffer, ghanaRoads)
google_map() %>%
add_markers(data = df) %>%
add_polylines(sf_intersection)
sf_intersection
距离计算中的对象。
关于r - 是否有解决方法来获得 googlemaps 未涵盖的点(长纬度)之间的旅行时间(使用 gmapsdistance)?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/49854447/
我需要修复 getLineNumberFor 方法,以便如果 lastName 的第一个字符位于 A 和 M 之间,则返回 1;如果它位于 N 和 Z 之间,则返回 2。 在我看来听起来很简单,但我不
您好,感谢您的帮助!我有这个: 0 我必须在每次点击后增加“pinli
Javascript 中是否有一种方法可以在不使用 if 语句的情况下通过 switch case 结构将一个整数与另一个整数进行比较? 例如。 switch(integer) { case
我有一列是“日期”类型的。如何在自定义选项中使用“之间”选项? 最佳答案 请注意,您有2个盒子。 between(在SQL中)包含所有内容,因此将框1设置为:DATE >= startdate,将框2
我有一个表,其中包含年、月和一些数字列 Year Month Total 2011 10 100 2011 11 150 2011 12 100 20
这个问题已经有答案了: Extract a substring between double quotes with regular expression in Java (2 个回答) how to
我有一个带有类别的边栏。正如你在这里看到的:http://kees.een-site-bouwen.nl/ url 中类别的 ID。带有 uri 段(3)当您单击其中一个类别时,例如网页设计。显示了一
这个问题在这里已经有了答案: My regex is matching too much. How do I make it stop? [duplicate] (5 个答案) 关闭 4 年前。 我
我很不会写正则表达式。 我正在尝试获取括号“()”之间的值。像下面这样的东西...... $a = "POLYGON((1 1,2 2,3 3,1 1))"; preg_match_all("/\((
我必须添加一个叠加层 (ImageView),以便它稍微移动到包含布局的左边界的左侧。 执行此操作的最佳方法是什么? 尝试了一些简单的方法,比如将 ImageView 放在布局中并使用负边距 andr
Rx 中是否有一些扩展方法来完成下面的场景? 我有一个开始泵送的值(绿色圆圈)和其他停止泵送的值(簧片圆圈),蓝色圆圈应该是预期值,我不希望这个命令被取消并重新创建(即“TakeUntil”和“Ski
我有一个看起来像这样的数据框(Dataframe X): id number found 1 5225 NA 2 2222 NA 3 3121 NA 我有另一个看起来
所以,我正在尝试制作正则表达式,它将解析存储在对象中的所有全局函数声明,例如,像这样 const a = () => {} 我做了这样的事情: /(?:const|let|var)\s*([A-z0-
我正在尝试从 Intellivision 重新创建 Astro-Smash,我想让桶保持在两个 Angular 之间。我只是想不出在哪里以及如何让这个东西停留在两者之间。 我已经以各种方式交换了函数,
到处检查但找不到答案。 我有这个页面,我使用 INNER JOIN 将两个表连接在一起,获取它们的值并显示它们。我有这个表格,用来获取变量(例如开始日期、结束日期和卡号),这些变量将作为从表中调用值的
我陷入了两个不同的问题/错误之间,无法想出一个合适的解决方案。任何帮助将不胜感激 上下文、FFI 和调用大量 C 函数,并将 C 类型包装在 rust 结构中。 第一个问题是ICE: this pat
我在 MySQL 中有一个用户列表,在订阅时,时间戳是使用 CURRENT_TIMESTAMP 在数据库中设置的。 现在我想从此表中选择订阅日期介于第 X 天和第 Y 天之间的表我尝试了几个查询,但不
我的输入是开始日期和结束日期。我想检查它是在 12 月 1 日到 3 月 31 日之间。(年份可以更改,并且只有在此期间内或之外的日期)。 到目前为止,我还没有找到任何关于 Joda-time 的解决
我正在努力了解线程与 CPU 使用率的关系。有很多关于线程与多处理的讨论(一个很好的概述是 this answer )所以我决定通过在运行 Windows 10、Python 3.4 的 8 CPU
我正在尝试编写 PHP 代码来循环遍历数组以创建 HTML 表格。我一直在尝试做类似的事情: fetchAll(PDO::FETCH_ASSOC); ?>
我是一名优秀的程序员,十分优秀!