Использование R GIS для пересечения шейп-файлов: буфер и данные точек
У меня есть два шейп-файла, которые я бы хотел пересечь. Один является буфером для границы вокруг Германии, а другой - шейп-файлом с точками широты и долготы, каждый из которых содержит столбец для "засухи".
Вот сводка для каждого из данных:
summary(Buffer)
Object of class SpatialPolygonsDataFrame
Coordinates:
min max
x 4.948681 12.81217
y 48.765813 52.51086
Is projected: FALSE
proj4string :
[+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0]
Data attributes:
id
Min. :1
1st Qu.:1
Median :1
Mean :1
3rd Qu.:1
Max. :1
И для шейп-файла данных о засухе:
summary(Shapefile)
Object of class SpatialPointsDataFrame
Coordinates:
min max
coords.x1 -11.75 44.75
coords.x2 27.25 70.75
Is projected: FALSE
proj4string :
[+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0]
Number of points: 3062
Data attributes:
lon lat pdsi_440
Min. :-11.75 Min. :27.25 Min. :-5.6550
1st Qu.: 6.75 1st Qu.:47.25 1st Qu.:-0.3078
Median : 17.25 Median :54.25 Median : 0.9800
Mean : 17.13 Mean :54.31 Mean : 0.7636
3rd Qu.: 26.75 3rd Qu.:62.25 3rd Qu.: 1.9207
Max. : 44.75 Max. :70.75 Max. : 5.2610
Столбец pdsi_440 предназначен для значений засухи.
Я хотел бы пересечь оба этих шейп-файла, а затем взять среднее значение для значений засухи.
Однако у меня проблемы с пересечением. Вот мой код:
Results <- intersect(Buffer, Shapefile)
Когда я делаю это, это просто дает мой буфер без каких-либо данных о засухе. Что я делаю неправильно? Спасибо за помощь.