¿Kriging con barreras en ArcGIS Desktop?

10

Estoy tratando de hacer una interpolación en una zona costera con puntos de datos en varios ríos que desembocan en una bahía y en lados opuestos de las penínsulas.

He creado una capa de costa para usarla como barrera, pero no puedo entender cómo hacerlo en ArcGIS 10.

No veo ningún objeto de barrera en ninguna parte de ninguna de las herramientas de kriging.

Intenté simplemente hacer el kriging y luego recortarlo en la tierra, pero obtengo resultados poco realistas donde los puntos de datos cercanos pero separados por tierra crean problemas.

Veo documentación en línea para hacer esto en scripts de geoprocesamiento en 9.3 pero nada similar para 10.

elayne927
fuente
Solo para agregar algunos enlaces para ver los cambios, compare la documentación de ArcGIS 9.3 versus ArcGIS 10
Mike T

Respuestas:

3

¿Los ríos corren generalmente en una dirección específica? Si es así, puede usar la agrupación direccional a lo largo de las líneas centrales de los ríos para restringir los semivariogramas empíricos solo a los pares de puntos que corren a lo largo de los ríos. (Consulte aquí para obtener más detalles: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Exploring_directional_autocorrelation/0031000000m2000000/ )

Otra opción sería dividir las muestras de puntos en conjuntos de datos separados y ejecutar su kriging en cada conjunto por separado. De todos modos, esto podría ser más apropiado, ya que cada río puede tener diferentes modelos estadísticos para ellos; pero no funcionaría si el tamaño de la muestra es pequeño para cada río.

¿Necesita usar un método geoestadístico o funcionaría un método determinista? Creo que puede lograr lo que desea con la Interpolación de difusión con barreras. Podrías usar tu tierra como una barrera de costos extremadamente altos. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00310000002s000000.htm

Además, aquí está la Interpolación de difusión con barreras absolutas como herramienta de script de geoprocesamiento. Su costa sería una barrera absoluta perfecta. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Diffusion_Interpolation_With_Barriers/003000000005000000/

castillo-blord
fuente
0

¿Has probado la herramienta Topo to Raster? Puede especificar el tipo de entidades de entrada (Elevación de punto, Arroyo, sumidero, contorno, límite, lago). Puede intentar configurar las entidades de capa de línea de costa en "TopoLake"

Además, el resultado de las herramientas Topo to Raster es un ráster "hidológicamente correcto"

Citando de la Ayuda: TopoLake Una clase de entidad poligonal que especifica la ubicación de los lagos. Todas las celdas ráster de salida dentro de un lago se asignarán al valor mínimo de elevación de todas las celdas a lo largo de la costa.

Jakub Sisak GeoGraphics
fuente
Creo que, en lugar de construir una representación de la topografía, el OP quiere crear una superficie de mediciones, como la calidad del agua, realizadas dentro de los cuerpos de agua. Topo2Raster no hará eso. (Tampoco Krige.)
whuber
Sí, estoy tratando de hacer exactamente eso.
elayne927
Tienes toda la razón Bill, Topo to Raster no permite kriging. Acabo de tener la impresión de que Elayne se refería a la interpolación en general. (A menudo escucho a la gente agrupar toda interpolación bajo kriging.)
Jakub Sisak GeoGraphics
Elayne, creo que sé lo que quieres decir con resultados poco realistas. Puede intentar convertir sus lagos en una cuadrícula de puntos espaciados regularmente con el valor de elevación del lago asignado a cada punto y agregarlo como entrada a la herramienta Topo to Raster como "Elevación de puntos". Esto tampoco dará un resultado perfecto, especialmente a lo largo de la costa pero parecerá más realista. ¿Has probado TIN y líneas de ruptura y luego has convertido el TIN resultante en ráster? ¿Puede publicar el enlace a los scripts de geoprocesamiento a los que se refiere?
Jakub Sisak GeoGraphics
0

Soy un adoptante tardío, por lo que solo puedo decirte qué funciona en 9.3 con la extensión Spatial Analyst: publicar solo en caso de que el método siga funcionando en 10. Primero, crea una cadena lineal cerrada que represente tu área de interés (suena como el inverso de tu costa aquí). Especifique su "contenedor" en el cuadro de diálogo de opciones de Spatial Analyst en la pestaña de extensión. Realice su análisis seleccionando Kriging de Interpolar a Ráster en el menú de Spatial Analyst.

¡Buena suerte!

rec.thegeom
fuente