Я новичок в SF и собираюсь, надеюсь, мой вопрос достаточно ясен.Мне удалось создать набор линий, соединяющих 1 точку с набором точек по всей территории США.Я могу читать округа США в мультиполигоны.Моя цель - найти и геолокации всех точек, где линии, которые я создал, пересекают границы округа.
До сих пор мне удавалось создавать линии из точек:
points_to_lines <- dt %>%
st_as_sf(coords = c("lon", "lat"), crs = 4326) %>%
group_by(lineid) %>%
summarize(do_union = FALSE) %>% lineid
st_cast("LINESTRING")
Этозаголовок линий
Simple feature collection with 1628 features and 1 field
geometry type: LINESTRING
dimension: XY
bbox: xmin: 30.1127 ymin: -91.32484 xmax: 37.23671 ymax: -82.31262
epsg (SRID): 4326
proj4string: +proj=longlat +datum=WGS84 +no_defs
# A tibble: 1,628 x 2
lineid geometry
<int> <LINESTRING [°]>
1 1 (33.51859 -86.81036, 36.16266 -86.7816)
2 2 (33.51859 -86.81036, 34.61845 -82.47791)
Это заголовок набора данных округа.
Reading layer `US_county_1930_conflated' from data source `~/county_gis/1930' using driver `ESRI Shapefile'
Simple feature collection with 3110 features and 18 fields
geometry type: MULTIPOLYGON
dimension: XY
bbox: xmin: -7115608 ymin: -1337505 xmax: 2258244 ymax: 4591848
epsg (SRID): NA
proj4string: +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=37.5 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs
Очень наивно я попытался дать им оба одинаковых набора координат, а затем st_intersectsих.Не разреженные матричные швы говорят о том, что все линии пересекаются только с одним округом.
gis1930_p <- st_set_crs(gis1930, 4326) %>% st_transform(4326)
st_intersects(points, gis1930_p, sparse=FALSE)
Я также рисую линии поверх округов, но отображается только карта округов США.
plot(gis1930_p[0], reset = FALSE)
plot(points[0], add = TRUE)
Любая помощь будет принята с благодарностью, и, пожалуйста, дайте мне знать, если я могу предоставить какие-либо дополнительные сведения.