Conteo de puntos sobre polígonos en R

En el post anterior se creó una Bounding Box de un ráster para generar 300 puntos aleatorios. En éste vamos a cargar un archivo vectorial, con 4 features, que cubre la misma zona del ráster y vamos a contar cuántos puntos caen sobre cada feature individual.

Las instrucciones para cargar el archivo vectorial y visualizarlo son las siguientes:

m_pol <- shapefile('multi_polygon.shp') #require(raster)
plot(m_pol, add=T) #agregado al gráfico del post anterior

obteniéndose el resultado a continuación:

R7

El objetivo, tal como se mencionó anteriormente, es realizar un contaje de cuantos puntos están comprendidos dentro de cada feature de m_pol e incorporar el resultado como una nueva columna de su Data Frame. Para ello se va a emplear el método over de la librería sp.

res <- over(pts, m_pol)
class(res)
[1] "data.frame"

Este objeto res (que es un Data Frame) se convierte en una columna de datos con la cuenta por id con:

count <- as.integer(table(res$id))

La incorporamos al Data Frame de m_pol simplemente con:

m_pol$count <- count
m_pol@data
  id     area    length count
0  1 20732.06  556.4997    19
1  2 26616.03  711.5804    19
2  3 25169.60  688.7720    17
3  4 60059.51 1329.1579    47
Esta entrada fue publicada en Lenguaje R. Guarda el enlace permanente.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión /  Cambiar )

Google photo

Estás comentando usando tu cuenta de Google. Cerrar sesión /  Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión /  Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión /  Cambiar )

Conectando a %s