Как я могу отобразить широту / долготу в ближайшую точку моего полигона / шейп-файла в R - PullRequest
Новая
       46

Как я могу отобразить широту / долготу в ближайшую точку моего полигона / шейп-файла в R

2 голосов
/ 08 апреля 2020

Я пытаюсь сопоставить набор латов / длин с моим файлом формы, чтобы получить информацию о CSDUID / разделах переписи.

Однако, есть некоторые выбросы или широта / долгота, которые находятся за пределами шейп-файла. Мне было интересно, можно ли сопоставить эти точки с ближайшей точкой шейп-файла, чтобы я мог получить информацию о делении для последующей части моего анализа?

Воспроизводимый пример

Допустим, У меня есть шейп-файл, подобный gr

nc <- st_transform(st_read(system.file("shape/nc.shp", package="sf")), 2264)                
gr = st_sf(
    label = apply(expand.grid(1:10, LETTERS[10:1])[,2:1], 1, paste0, collapse = " "),
    geom = st_make_grid(nc))

# gr
# Simple feature collection with 100 features and 1 field
# geometry type:  POLYGON
# dimension:      XY
# bbox:           xmin: 406262.2 ymin: 48374.87 xmax: 3052887 ymax: 1044158
# epsg (SRID):    2264
# proj4string:    +proj=lcc +lat_1=36.16666666666666 +lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.2192024384 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs
# First 10 features:
#    label                           geom
# 1   J  1 POLYGON ((406262.2 48374.87...
# 2   J  2 POLYGON ((670924.7 48374.87...
# 3   J  3 POLYGON ((935587.2 48374.87...
# 4   J  4 POLYGON ((1200250 48374.87,...
# 5   J  5 POLYGON ((1464912 48374.87,...
# 6   J  6 POLYGON ((1729575 48374.87,...
# 7   J  7 POLYGON ((1994237 48374.87,...
# 8   J  8 POLYGON ((2258900 48374.87,...
# 9   J  9 POLYGON ((2523562 48374.87,...
# 10  J 10 POLYGON ((2788225 48374.87,...

И у меня есть точка, которая находится за пределами ограничительной рамки: 

point = st_point(c(106160,28370))

library(ggplot2)
ggplot() + geom_sf(data = gr) + geom_sf(data = point)

enter image description here

Как мне сопоставить точки с ближайшим полигоном?

1 Ответ

1 голос
/ 08 апреля 2020

У вас есть отличные функции для пространственного анализа, особенно пространственные объединения, в пакете sf. Предполагая, что ваши точечные данные имеют имя point, а полигоны имеют имя .... polygons (оба объекта sf), вы можете использовать sf::st_join для пространственных объединений.

st_join принимает разные типы пространственного картирования. Я думаю, что вы заинтересованы в st_nearest_feature слиянии (сначала удалите точки, которые попадают в ваши полигоны, например, первым st_join, используя st_within) 

library(sf)

st_join(
  points,
  polygons,
  join = st_nearest_feature
)

В отредактированном примере это даст: 

point = st_sf(st_sfc(point), crs = st_crs(gr))
st_join(point, gr, join = st_nearest_feature)

см. ?st_join до c.

...