Elevación de interpolación precisa de las curvas de nivel en QGIS?

8

Estoy tratando de interpolar elevaciones desde las curvas de nivel (~ 70000 a intervalos de 1 m, que abarcan ~ 20 km ^ 2) en QGIS. 

Cuando interpolo usando Inverse Distance Weighting (IDW)en el Interpolationcomplemento, ¿usa todas las líneas de contorno en lugar de solo los contornos adyacentes?

El uso del GDAL Rastercomplemento (distancia inversa a una potencia) también parece usar todos los contornos y tratar las líneas de contorno no como líneas sino como puntos.

También lo intenté Triangular interpolation (TIN), tratando las líneas de contorno como líneas de estructura, pero bloquea QGIS, probablemente debido a la cantidad y densidad de las líneas de contorno.

Para potencias inferiores a 5, utilizando el complemento de trama GDAL, el DEM resultante es muy ruidoso y los puntos de contorno tienden a convertirse en máximos o mínimos locales.

es decir, para las celdas de cuadrícula que tienen una línea de contorno que las atraviesa, pero no contienen los puntos que definen la línea de contorno, las elevaciones interpoladas se desplazan hacia la elevación promedio de todos los datos de contorno, en lugar de que se les asigne la elevación de la línea de contorno en ejecución a través de ellos; Esto se debe a que la interpolación explica innecesariamente las elevaciones de cada punto de contorno (aunque ponderado por la distancia inversa).

Si bien este problema se reduce cuando uso IDW con potencias más altas, las regiones relativamente planas comienzan a parecerse a las interpolaciones vecinas más cercanas, como pasos, en lugar de como una región grande con una pequeña inclinación.

De otras publicaciones, parece que el r.surf.contouralgoritmo GRASS tiene problemas similares.

Usar un radio de búsqueda es mejor, pero aún suaviza las características más pequeñas que el radio de búsqueda, y aún permite que los contornos no relevantes para una celda de la cuadrícula afecten su elevación interpolada. 

¿Hay alguna manera en QGIS de interpolar usando solo las líneas de contorno inmediatamente adyacentes en todas las direcciones?

¿Cuál es la mejor práctica para interpolar datos de contorno en QGIS de tal manera que los datos de contorno se traten como contornos, en lugar de como un conjunto grande de puntos o líneas? 

jlnzbk
fuente
Por lo que vale, encontré r.surf.contourque tenía un aspecto mucho más natural de lo v.surf.rstque puede ver mis esfuerzos en gwprojects.org/forum/viewtopic.php?f=3&t=2084 y la cartilla que solía comenzar
Mr Purple

Respuestas:

3

Vengo del mundo de AutoCAD / ESRI y no estoy familiarizado con las herramientas de QGIS. Sin embargo, he creado muchas superficies a partir de contornos.

Para crear una superficie a partir de contornos de superficie, se requiere un procesamiento adicional. Primero debe desmalezar y complementar los puntos de contorno Consulte https://knowledge.autodesk.com/support/autocad-civil-3d/learn-explore/caas/CloudHelp/cloudhelp/2017/ENU/Civil3D-UserGuide/files/GUID- 523368FB-BF6D-4E7C-8E38-7BA9A2B7E38A-htm.html

En segundo lugar, debe crear líneas de ruptura para el fondo de los valles, líneas de cresta y otras rupturas en la superficie para evitar áreas planas artificiales.

Luego, el TIN proporciona los mejores resultados para la topografía utilizando los puntos de contorno y las líneas de ruptura. Alisar los contornos resultantes los hace más bonitos pero no más precisos.

Explicación de ESRI de contornear datos estadísticos y gráficos: https://www.esri.com/news/arcuser/0704/files/interpolating.pdf

alemán
fuente