Estoy planeando utilizar el método de interpolación utilizando la herramienta de análisis espacial y he encontrado varios métodos como spline NN y kriging. Sé que puedo hacer la interpolación utilizando el analista geoestadístico, pero algunos métodos de interpolación no están disponibles en el analista geoestadístico, como spline y NN, aunque la medición de errores se puede obtener del analista geoestadístico. Mi pregunta es cómo obtener los resultados de los errores, como los errores cuadráticos medios y los datos observados y previstos del modelo de interpolación utilizando Spatial Analyst en ArcMap.
Respuestas
¿Demasiados anuncios?
Antti Huima
Puntos
15465
Dependiendo del tipo de spline, lo más probable es que se trate de un caso especial de interpolación de función de base radial, pero RDF es teóricamente equivalente a una de las formas de kriging. Con RBF se obtiene una función interpoladora, con kriging la función interpoladora está implícita pero no se da explícitamente.
Con kriging se puede utilizar un vecindario de búsqueda móvil pero no con RBF (aunque M. Powell ha considerado la aproximación secuencial que es análoga a vecindarios de búsqueda cada vez mayores).
RBF no genera los "errores de interpolación", ya que también es un interpolador exacto, por lo que habría que utilizar una variación de la validación cruzada.
Véase 1992, Myers, D.E., Selection of a radial basis function for data interpolation en Advances in Computer Methods for Partial Differential Eq. VII, R. Vichnevetsky, D. Knight y G.Richter (eds), IMACS, 553-558.
Puede descargarse en www.u.arizona.edu/~donaldm
lauren
Puntos
9
Puede calcular fácilmente el error cuadrático medio (RMSE) de cualquier método de interpolación de Spatial Analyst. Así es como yo lo haría:
- Añade dos campos extra a tu capa de puntos, y llámalos algo así como: interpolación y SqDeviation
- Ejecute el interpolador de su elección para crear la nueva superficie (IDW, Kriging, Vecino más próximo, etc.).
- Ejecute la herramienta "Extraer valores a puntos" utilizando sus valores de puntos originales y la nueva superficie.
- Una la salida a sus puntos originales para rellenar la columna "interpolación".
- Rellene la columna "SqDeviation" utilizando la calculadora de campo para restar los valores interpolados de los valores de los puntos originales y elevarlos al cuadrado. (P1 - P2)2
- Ahora que has calculado todas las desviaciones al cuadrado, súmalas, divídelas por n y saca la raíz cuadrada.
Los pasos 5 y 6 deben satisfacer esta ecuación para el RMSE: 
Además, recomiendo encarecidamente la lectura de esta pequeña hoja informativa de 4 páginas sobre los métodos de interpolación de Spatial Analyst, para entender lo que hace cada herramienta de ArcGIS: http://www.esri.com/news/arcuser/0704/files/interpolating.pdf .
Por ejemplo, con Spline, la superficie se fuerza a través de los puntos originales, por lo que el RMSE básicamente carece de sentido porque cada una de las desviaciones al cuadrado es igual a 0. Es bueno tener en cuenta que el RMSE es sólo una métrica para evaluar los resultados, pero no debe ser la prueba de todo para la toma de decisiones. Hay muchos artículos en Internet sobre cómo el RMSE puede ser engañoso. En general, utilizo spline para hacer superficies bastante suaves para las presentaciones, pero en realidad utilizo EBK en Geostatistical Analyst como mi interpolador caballo de batalla. Miro el RMSE para comparar los ajustes de mis diversas ejecuciones de EBK para ayudarme a decidir sobre sus ajustes y opciones en mis datos subyacentes, pero siempre tomo el RMSE con un grano de sal.
