Estoy tratando de analizar cómo los conjuntos de rásteres se relacionan entre sí utilizando algunas técnicas estadísticas. Como no tengo mucha experiencia usando las herramientas de estadísticas espaciales en ArcGIS , estaba exportando mis rásteres como archivos Ascii y analizándolos usando R (específicamente el maptools
paquete y readAsciiGrid()
). Esto ha estado funcionando bien (pero como los conjuntos de datos tienen 90,000 puntos, es lento para ejecutar el análisis), pero no sé si estoy recreando en R , la funcionalidad existente en ArcGIS .
Por ejemplo, quiero realizar regresiones entre cada uno de estos rásteres utilizando algunas transformaciones diferentes (logarítmica, exponencial, etc.). ¿Se puede hacer esto dentro de ArcGIS ? Una segunda pregunta más amplia es si hay métodos estadísticos estándar para examinar este tipo de datos.
Cada par de ráster tiene valores coincidentes de datos / sin datos y todos los parámetros son idénticos, aparte del valor de la celda de cuadrícula.
Mire el paquete ráster R , fue diseñado con este tipo específico de problema en mente. Intenta mantener la menor cantidad de ráster en la memoria y realiza una serie de operaciones espaciales básicas: a través de GDAL puede manejar archivos de cuadrícula ESRI nativos, entre muchos otros formatos. La viñeta proporciona algunos buenos ejemplos de su uso.
fuente
Acabo de leer sobre un paquete de ArcGIS 10 que es un complemento R. No he tenido la oportunidad de explorarlo por completo, pero podría ser posible modificar esto para hacer lo que describo anteriormente.
fuente