16 votos

¿Qué técnica de interpolación es adecuada para una batimetría de un pequeño lago?

Se trata de un lago de aproximadamente 13 ha con 81 puntos de profundidad muestreados en unas 10 líneas de transecto:

Sampled points

Anteriormente en ARCgis con un lago de 50 ha y unos 100 puntos de muestreo de profundidad, obtuve una salida decente con el TopotoRaster herramienta. Sin embargo, no parece existir un equivalente exacto con QGIS o el software de código abierto.

En QGIS, obtuve el raster interpolado de abajo, utilizando el método TIN en el plugin de Interpolación Raster. Sin embargo, esto parecía ser un resultado afortunado, cuando traté de repetir, obtuve resultados diferentes (insatisfactorios) a pesar de no variar realmente ningún parámetro.

Interpolated raster

Se trata más de depósitos urbanos (construidos) que de lagos, por lo que su perfil de lecho es bastante regular, de ahí la baja densidad de muestreo puntual. Algunos artículos parecen sugerir la ponderación inversa de la distancia (IDW) como la técnica más adecuada, pero parece que obtuve los peores resultados con ella (posiblemente debido a la baja densidad de puntos).

¿Existe alguna heurística que pueda emplearse aquí teniendo en cuenta el tamaño del lago, su perfil regular y la densidad de muestreo de puntos para llegar a qué método de interpolación es apto? (entre Kriging, IDW, Bilinear, Convolución cúbica, TIN o spline)

¿O es siempre un poco de ensayo y error?

0 votos

Intuitivamente, creo que la interpolación con el vecino más cercano es lo más lógico, ya que las diferencias de profundidad son relativamente regulares. Sería estupendo si alguien pudiera sugerir un recurso en el que se describan algunos procedimientos/directrices para asegurarse de que se consigue un resultado coherente con respecto a la interpolación de este tipo y también en otros casos.

0 votos

Reefmaster crea contornos y relieves sombreados de buen aspecto con TIN reefmaster.com.au/index.php/forum/support/ .

2 votos

Mirando la imagen QGIS no utiliza el límite automáticamente. Añadiría un montón de puntos con profundidad cero a lo largo del límite.

10voto

Andre Silva Puntos 2910

Una alternativa es interpolación spline como se sugiere en el post relacionado: Interpolación de la batimetría multihaz .

Desde QGIS, utilice la herramienta GRASS v.surf.rst :

Realiza la interpolación de superficies a partir de un mapa de puntos vectoriales mediante splines.

A continuación, puede probar diferentes tipos de parametrización disponibles en la herramienta. Existe la opción de aplicar una validación cruzada de un solo uso para la optimización de los parámetros, que variará gradualmente la tensión, el suavizado, etc., con el fin de encontrar el mínimo error de predicción.

Como sugiere user30184, asegúrese de que el conjunto de datos de entrada utiliza puntos a nivel del suelo en los límites del lago, para evitar la extrapolación.

0 votos

Spline es lo que utilicé en un proyecto de batimetría fluvial que hice como contratista para el USGS. Nunca pregunté exactamente por qué, pero quedamos satisfechos con los resultados. Estuve a punto de hacer una respuesta similar pero no estaba seguro de la justificación. Buena respuesta.

0 votos

Gracias André, te agradezco mucho tu respuesta rápida y nítida. Tuve algunos problemas tratando de ejecutar v.surf.rst, se negó dentro de QGIS para darme cualquier salida raster adecuada con los valores. Sospecho que hay problemas con la forma en que he proyectado los datos. Cuando lo resuelva pondré los resultados en las ediciones de la pregunta.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X