2011-09-09 18 views
5

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 

cevap

6

kullanımı gIntersects (değil gIntersection) Burada

bir örnektir DOĞRU/FALSE vektörünü, tamponunuzda hangi noktaların bulunduğunu TRUE. Ardından, puan veri çerçevenizi bölümlere ayırın. İş bitmiş.

+0

Bahşiş için teşekkürler! Bu harika bir çalışma. – Pascal

İlgili konular