R'de bir SpatialPointsDataFrame nesnesi olarak bir koordinat kümesi var ve yalnızca bu çokgen içinde bulunan noktaları almak için bu noktaları bir çokgen kullanarak kırpıyorum. Bunu, rgeos paketindeki gsntersection fonksiyonu ile yapıyorum. Benim problemim, fonksiyonun sadece bu noktaların koordinatlarını değil, onlarla ilişkili meta verileri döndürmesidir. Gingersection'ın tüm verileri yalnızca koordinatlardan ziyade sonuca iletmesinin bir yolu var mı?rGeos paketinden gIntersection kullanıldığında meta verileri tutuluyor R
SpatialPointsDataFrame temsil türleri oluşumu puanı:
> 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
bir çokgen bu noktaları kırpma: byid = ile
> 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
Bahşiş için teşekkürler! Bu harika bir çalışma. – Pascal