Как сопоставить координаты (UTM и геометрию) этих объектов df / sp?

Я был бы очень рад, если бы вы помогли мне с этой проблемой. Я хочу geom_point df daa_84 в файл shp_5. После просмотра нескольких связанных вопросов о stackoverflow и тестирования их ответов (например, создания объекта sp из daa_84 и преобразования координат UTM в соответствие с координатами shp_5) я получаю только что-то вроде графика. Кроме того, я знаю, что зона UTM (19S) и код EPSG, связанный с моей страной (32719) системы координат (WGS84), нужны для чего-то, ха-ха. Любые идеи?

> head(daa_84)
            # A tibble: 6 x 2
              utm_este utm_norte
                 <dbl>     <dbl>
            1   201787   6364077
            2   244958   6247258
            3   245947   6246281
            4   246100   6247804
            5   246358   6242918
            6   246470   6332356
            
    
    > head(shp_5)
            Simple feature collection with 6 features and 1 field
            geometry type:  MULTIPOLYGON
            dimension:      XY
            bbox:           xmin: -7973587 ymin: -3976507 xmax: -7838155 ymax: -3766040
            projected CRS:  WGS 84 / Pseudo-Mercator
                 Comuna                       geometry
            1 Rinconada MULTIPOLYGON (((-7871440 -3...
            2   Cabildo MULTIPOLYGON (((-7842610 -3...
            3   Petorca MULTIPOLYGON (((-7873622 -3...
            4 Panquehue MULTIPOLYGON (((-7874932 -3...
            5     Olmué MULTIPOLYGON (((-7916865 -3...
            6 Cartagena MULTIPOLYGON (((-7973501 -3...
            
        
        ggplot() + geom_sf(data = shp_5, aes()) + 
              geom_point(data = daa_84, aes(x= "utm_este", "utm_norte"),
                         alpha = 0.05, size = 0.5) + 
              labs(x = "Latitude", y = "Longitude")+
              theme_bw()

мой прогресс на данный момент

ИЗМЕНИТЬ

в дополнение к ответу william3031 этот код также работает

library(sf)
daa_84 = tribble(~utm_este, ~utm_norte,
201787,   6364077,
244958,   6247258,
245947,   6246281,
246100,   6247804,
246358,   6242918,
246470,   6332356)


daa_84 = st_as_sf(daa_84, 
                  coords=c('utm_este', 'utm_norte'), 
                  crs=st_crs(32719)) %>% 
      st_transform(st_crs(shp_5))

sf coordinates utm
person MPZ    schedule 08.02.2021    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Это сработает для вас. Я использовал другой набор данных для Южной Америки, поскольку вы не предоставили воспроизводимый пример.

library(tidyverse)
library(sf)
library(spData) # just for the 'world' dataset

# original
daa_84 <- data.frame(
  utm_este = c(201787L, 244958L, 245947L, 246100L, 246358L, 246470L),
  utm_norte = c(6364077L, 6247258L, 6246281L, 6247804L, 6242918L, 6332356L)
)

# converted
daa_84_sf <- st_as_sf(daa_84, coords = c("utm_este", "utm_norte"), crs = 32719)


# load world to get South America
data("world")
sam <- world %>% 
  filter(continent == "South America")

# plot
ggplot() +
  geom_sf(data = sam) + 
  geom_sf(data = daa_84_sf)

введите описание изображения здесь

person william3031    schedule 08.02.2021