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 Interpolation
complemento, ¿usa todas las líneas de contorno en lugar de solo los contornos adyacentes?
El uso del GDAL Raster
complemento (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.contour
algoritmo 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?
r.surf.contour
que tenía un aspecto mucho más natural de lov.surf.rst
que puede ver mis esfuerzos en gwprojects.org/forum/viewtopic.php?f=3&t=2084 y la cartilla que solía comenzar