Soy muy nuevo en el uso de datos SIG y solo tengo una modesta experiencia con R. He estado leyendo sobre cómo analizar datos espaciales utilizando el libro PDF de espacial-analyst.net, así que no estoy completamente perdido, pero pensé que podría describirlo. mi problema y la gente pueden sugerir ideas.
Tengo un conjunto de datos con aproximadamente 2000 mediciones en diferentes coordenadas lat / long, aunque probablemente subdividiré este conjunto de datos ya que los datos se recopilaron durante 3 años y las condiciones cambiaron con el tiempo. Llamemos a la variable que se mide "IP".
Quiero crear un mapa de IP en el área completa en cuestión utilizando Kriging o algún otro método de interpolación en los datos de la muestra. Luego quiero crear un histograma que mida la cantidad de tierra en varios segmentos de IP. También necesitaré crear un histograma que muestre el número de muestras en cada cubeta (tenga en cuenta que una muestra podría tener una IP real más alta o más baja de lo que predice kriging para su tierra).
Sigo cómo cargar los datos en un SpatialPointsDataFrame y ejecuto un análisis de kriging, donde tengo problemas es cómo convertir esos datos en un marco de datos en cuadrícula para poder hacer el análisis de histograma.
¿Alguna sugerencia para convertir puntos en cuadrículas?
pts
está en mi ejemplo anterior. ¡Simplemente ejecute el código en su objeto SpatialPointsDataFrame!r
) para usar en unhist
procedimiento similar, que es simplemente una cuestión de expresión similarhist(getValues(r))
.El paquete plotKML tiene una función llamada
vect2rast
. Esta función básicamente extiende larasterize
función disponible en el paquete ráster. La ventaja devect2rast
; sin embargo, es que no requiere entrada del lado del usuario, es decir, determina automáticamente el tamaño de celda de la cuadrícula y el cuadro delimitador en función de las propiedades del conjunto de datos de entrada. El tamaño de la celda de la cuadrícula se estima en función de la densidad / tamaño de las características en el mapa (nndist
función en el paquete spatstat).fuente