r - 将纬度/经度转换为正确的格式

Question

我有一些看起来像这样的数据:

     long_bnk                 lat_bnk                  
[1,] "3<U+00B0> 52' 30.1\" W" "40<U+00B0> 44' 3.7\" N" 
[2,] "2<U+00B0> 44' 54.4\" E" "42<U+00B0> 7' 18.1\" N" 
[3,] NA                       NA                       
[4,] "2<U+00B0> 7' 54.2\" E"  "41<U+00B0> 31' 21.9\" N"
[5,] "0<U+00B0> 1' 54.5\" W"  "39<U+00B0> 58' 59.3\" N"
[6,] "3<U+00B0> 41' 15.5\" W" "40<U+00B0> 27' 47.2\" N"

我正在尝试将数据放入正确的纬度/经度格式。我正在运行以下内容:

pts_bnk[pts_bnk==""] <- NA 
pts_bnk <- pts_bnk[complete.cases(pts_bnk),]
pts_bnk <- matrix(as.numeric(sp::char2dms(as.vector(pts_bnk), "°")), ncol=2)

但是，我不断得到:

Error in if (any(abs(object@deg) > 90)) return("abs(degree) > 90") :
missing value where TRUE/FALSE needed In addition: Warning message: In asMethod(object) : NAs introduced by coercion

在转换为正确的纬度/经度格式时我哪里出错了？ 数据:

pts_bnk <- structure(c("3<U+00B0> 52' 30.1\" W", "2<U+00B0> 44' 54.4\" E", 
NA, "2<U+00B0> 7' 54.2\" E", "0<U+00B0> 1' 54.5\" W", "3<U+00B0> 41' 15.5\" W", 
"40<U+00B0> 44' 3.7\" N", "42<U+00B0> 7' 18.1\" N", NA, "41<U+00B0> 31' 21.9\" N", 
"39<U+00B0> 58' 59.3\" N", "40<U+00B0> 27' 47.2\" N"), .Dim = c(6L, 
2L), .Dimnames = list(NULL, c("long_bnk", "lat_bnk")))

编辑:

基本上我想使用以下方法绘制数据:

library(ggrepel)
library(ggmap)
register_google(key = "MyKey")
spain <- get_map("Spain", zoom = 6)
ggmap(spain, extent = "normal") +
geom_point()

编辑 2:

我拥有的(有效的)原始数据如下:

dms_lat <- readLines(n=5)
1  40° 25' 35.8" N
2  40° 26' 28.4" N
3  40° 28' 39.8" N
4                 
5  38° 59' 15.0" N
dms_long <-readLines(n=5)
1   3° 41' 19.9" W
2   3° 47' 42.2" W
3   3° 41' 11.7" W
4                 
5   3° 55' 29.6" W
pts <- cbind(dms_long, dms_lat)
pts <- sub("^\\d+\\s+", "", pts)
pts[pts==""] <- NA
pts <- pts[complete.cases(pts),]
pts <- matrix(as.numeric(sp::char2dms(as.vector(pts), "°")), ncol=2)
library(rworldmap)
plot(subset(getMap(resolution = "low"), NAME=="Spain"))
points(pts[,1], pts[,2], col = "red", pch=3, cex = 0.6)

我目前的数据(不起作用)是:

x <- structure(c("3<U+00B0> 52' 30.1\" W", "2<U+00B0> 44' 54.4\" E", 
NA, "2<U+00B0> 7' 54.2\" E", "0<U+00B0> 1' 54.5\" W", "3<U+00B0> 41' 15.5\" W", 
"40<U+00B0> 44' 3.7\" N", "42<U+00B0> 7' 18.1\" N", NA, "41<U+00B0> 31' 21.9\" N", 
"39<U+00B0> 58' 59.3\" N", "40<U+00B0> 27' 47.2\" N"), .Dim = c(6L, 
2L), .Dimnames = list(NULL, c("long_bnk", "lat_bnk")))

x %>%
  data.frame() %>% 
  mutate(
    lat = sub("<U\\+00B0>", "\u00B0", lat_bnk),
    long = sub("<U\\+00B0>", "\u00B0", long_bnk)
  ) %>% 
  select(lat, long) %>% 
  drop_na()

我正在尝试使第二个数据与第一个数据相等，以便可以使用 ggmap() 绘制它。

Answer 1

根据我们的评论，让我知道这是否适合您:

library(rworldmap)
library(sp)
library(dplyr)

pts <- x %>%
  data.frame() %>% 
  mutate(
    lat = sub("<U\\+00B0>", "d", lat_bnk),
    long = sub("<U\\+00B0>", "d", long_bnk)
  ) %>% 
  select(lat, long) %>% 
  drop_na()

pts_long <- as.numeric(char2dms(pts[["long"]]))
pts_lat <- as.numeric(char2dms(pts[["lat"]]))

plot(subset(getMap(resolution = "low"), NAME=="Spain"))
points(pts_long, pts_lat, col = "red", pch=3, cex = 0.6)

请注意，char2dms 的用法如下:

char2dms(from, chd = "d", chm = "'", chs = "\"")

默认的度数字符终止符是字母 d(如果未指定 chd，则为度数符号的替代)。

使用 ggmap 可以在 geom_point 中传递经度和纬度:

library(ggrepel)
library(ggmap)
library(ggthemes)

pts_data <- data.frame(pts_long, pts_lat)

# Note requires Google key
spain <- ggmap::get_map("Madrid, Spain", zoom = 6)
ggmap(spain, extent = "normal") +
  geom_point(data = pts_data, aes(x = pts_long, y = pts_lat)) +
  theme_map()