Преобразование lat / long в правильный формат

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"

Я пытаюсь разместить данные в правильном формате lat / long. Я выполняю следующее:

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)

Тем не менее, я продолжаю получать:

Ошибка в if (any (abs (object@deg)> 90)) return ("abs (градус)> 90 "):
пропущенное значение, где требуется ИСТИНА / ЛОЖЬ. Дополнительно: предупреждающее сообщение: в asMethod (объект): NA, введенные с помощью принуждения

Где я иду неправильно в преобразование в правильный формат lat / long? Данные:

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 you может передавать долготу и широту в 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()

Answer 2

Мы можем добавить mutate_at в конце

library(dplyr)
x %>%
  data.frame() %>% 
  mutate(
    lat = sub("<U\\+00B0>", "\u00B0", lat_bnk),
    long = sub("<U\\+00B0>", "\u00B0", long_bnk)
  ) %>% 
  select(lat, long) %>% 
  drop_na()%>% 
  mutate_at(vars(matches('^(lat|long)')), ~ as.numeric(sp::char2dms(., "°")))