Преобразование на север / восток в десятичные градусы в R - PullRequest
Новая
       116

Преобразование на север / восток в десятичные градусы в R

0 голосов
/ 05 февраля 2020

Я пытаюсь преобразовать векторы, содержащие данные о северном и восточном направлениях, в десятичные градусы, используя R. До сих пор я был в состоянии выполнить пакет rgdal для выполнения проекции, но у меня возникла проблема. Аргументом для зоны UTM должна быть строка из одного символа, но в моих данных есть несколько зон UTM. Есть ли способ использовать функцию rgdal :: SpatialPoints () для установки нескольких зон UTM? Я попытался ввести набор символов и использовать вектор в функции CRS (), но получил предупреждение о том, что использовалась только первая зона. 


df<- data.frame(X = c(774869, 771437, 1051883, 524468), 
                Y = c(414498, 403790, 184967, 779682), 
                WGS.1984.UTM.Zone = c("57N", "57N", "54N", "59N"))

 library(rgdal)
  sputm <- SpatialPoints(df[c("X", "Y")], proj4string=CRS("+proj=utm +zone=57N +datum=WGS84"))  
  spgeo <- spTransform(sputm, CRS("+proj=longlat +datum=WGS84"))


#resultant coordinates are all transformed relative to UTM zone 57N

df$projargs<- paste("+proj=utm +zone=", df$WGS.1984.UTM.Zone, " +datum=WGS84", sep = "")

#trying again with set projargs character strings

spUTM<- SpatialPoints(df[c("X", "Y")], proj4string=CRS(df$projargs))

Предупреждающие сообщения:

1 : In if (! Is.na (projargs)) {: условие имеет длину> 1, и будет использоваться только первый элемент

2: In if (! Is.na (projargs)) {: the условие имеет длину> 1, и будет использоваться только первый элемент

3: В случае если (is.na (projargs)) uprojargs <- projargs else uprojargs <- paste (unique (unlist (strsplit (projargs,: условие имеет длину> 1, и будет использоваться только первый элемент

Ответы [ 2 ]

1 голос
/ 05 февраля 2020

Вот еще одно решение с пакетом sf.

Создайте объект sf с st_as_sf(). Обратите внимание на новый столбец geometry. 

library(sf)
library(dplyr)

df_sf <- st_as_sf(df,
                  coords = c("X", "Y"),
                  remove = FALSE,
                  crs = projargs)

df_sf

Simple feature collection with 4 features and 4 fields
geometry type:  POINT
dimension:      XY
bbox:           xmin: 524468 ymin: 184967 xmax: 1051883 ymax: 779682
epsg (SRID):    32657
proj4string:    +proj=utm +zone=57 +datum=WGS84 +units=m +no_defs
        X      Y WGS.1984.UTM.Zone                         projargs               geometry
1  774869 414498               57N +proj=utm +zone=57N +datum=WGS84  POINT (774869 414498)
2  771437 403790               57N +proj=utm +zone=57N +datum=WGS84  POINT (771437 403790)
3 1051883 184967               54N +proj=utm +zone=54N +datum=WGS84 POINT (1051883 184967)
4  524468 779682               59N +proj=utm +zone=59N +datum=WGS84  POINT (524468 779682)

Перепроектируйте с помощью st_transform(), который обновляет столбец геометрии. Мы можем вытянуть их в новые столбцы с помощью st_coordinates(), который возвращает матрицу из двух столбцов для каждой строки. 

df_sf_longLat <- st_transform(df_sf, crs = "+proj=longlat +datum=WGS84") %>% 
  mutate(long_degrees = st_coordinates(.)[,1],
         lat_degrees = st_coordinates(.)[,2])

df_sf_longLat

Simple feature collection with 4 features and 6 fields
geometry type:  POINT
dimension:      XY
bbox:           xmin: 159.2216 ymin: 1.667156 xmax: 163.9555 ymax: 7.053618
epsg (SRID):    4326
proj4string:    +proj=longlat +datum=WGS84 +no_defs
        X      Y WGS.1984.UTM.Zone                         projargs                  geometry long_degrees lat_degrees
1  774869 414498               57N +proj=utm +zone=57N +datum=WGS84 POINT (161.4747 3.746529)     161.4747    3.746529
2  771437 403790               57N +proj=utm +zone=57N +datum=WGS84  POINT (161.4435 3.64983)     161.4435    3.649830
3 1051883 184967               54N +proj=utm +zone=54N +datum=WGS84 POINT (163.9555 1.667156)     163.9555    1.667156
4  524468 779682               59N +proj=utm +zone=59N +datum=WGS84 POINT (159.2216 7.053618)     159.2216    7.053618
0 голосов
/ 05 февраля 2020

Функция не векторизована. Поэтому нам нужно l oop

library(rgdal)
spUTM <- lapply(df$projargs, function(x)  
          SpatialPoints(df[c("X", "Y")], proj4string=CRS(x)))

str(spUTM)
#List of 4
# $ :Formal class 'SpatialPoints' [package "sp"] with 3 slots
#  .. ..@ coords     : num [1:4, 1:2] 774869 771437 1051883 524468 414498 ...
#  .. .. ..- attr(*, "dimnames")=List of 2
#  .. .. .. ..$ : NULL
#  .. .. .. ..$ : chr [1:2] "X" "Y"
#  .. ..@ bbox       : num [1:2, 1:2] 524468 184967 1051883 779682
#  .. .. ..- attr(*, "dimnames")=List of 2
#  .. .. .. ..$ : chr [1:2] "X" "Y"
# ...

или использовать map из purrr

library(dplyr)
library(purrr)
df %>% 
     mutate(spUTM = map(projargs, ~ 
          SpatialPoints(df[c('X', 'Y')], proj4string = CRS(.x))))
...