Хранение метаданных при использовании gIntersection из пакета rgeos в R - PullRequest
Новая
       7

Хранение метаданных при использовании gIntersection из пакета rgeos в R

6 голосов
/ 10 сентября 2011

У меня есть набор координат в виде объекта SpatialPointsDataFrame в R, и я обрезаю эти точки, используя многоугольник, чтобы получить только те точки, которые находятся внутри этого многоугольника.Я делаю это с помощью функции gIntersection из пакета rgeos.Моя проблема в том, что функция возвращает только координаты этих точек, а не метаданные, связанные с ними.Есть ли способ, чтобы gIntersection передавал все данные в результат, а не только координаты?

Вот пример:

SpatialPointsDataFrame, представляющий точки появления видов: 

> spexample
            coordinates SpAbbr InstitutionCode CatalogNumberText
1   (-76.8727, 3.66282) BanRot       EBIRD_COL       OBS81997559
2   (-76.9749, 3.71683) BanRot          AUDCLO       OBS89767945
3    (-76.884, 3.61609) BanRot          AUDCLO       OBS89769896
4   (-77.5167, 5.51667) BanRot            AMNH       Skin-123476
5   (-76.0334, 4.86669) BanRot            LACM             34848
6   (-78.4333, 1.43333) BanRot           LSUMZ             38939
7        (-78.55, 0.95) BanRot            ANSP            182799
8  (-79.2139, 0.471944) BanRot          AUDCLO       OBS58485973
9  (-78.5104, 0.895349) BanRot          AUDCLO       OBS84822747
10  (-78.3781, 1.51028) BanRot          AUDCLO       OBS67916517
11       (-75.15, 7.07) BanRot   8110002317-09         4743-5160

Обрезать эти точки до многоугольника: 

> gIntersection(spexample,bufferclip)
SpatialPoints:
          x         y
1 -78.55000 0.9500000
1 -78.51036 0.8953493
1 -78.43333 1.4333333
1 -78.37810 1.5102800
1 -76.97495 3.7168289
1 -76.88397 3.6160872
1 -76.87271 3.6628163
1 -76.03337 4.8666900
Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84
+ellps=WGS84 +towgs84=0,0,0

1 Ответ

6 голосов
/ 10 сентября 2011

Используйте gIntersects (не gIntersection) с byid = TRUE, чтобы получить вектор TRUE / FALSE, точки которого находятся в вашем буфере.Затем установите подмножество ваших данных.Работа выполнена.

...