¿Generando contornos desde DEM usando GDAL?

16

Acabamos de recibir un gran conjunto de DEM en el trabajo y me gustaría generar contornos a partir de ellos. Los DEM tienen una resolución de 1m y un tamaño de 1kmx1km.

Salida de gdalinfo:

Driver: AAIGrid/Arc/Info ASCII Grid
Files: 380000_6888000_1k_1m_DEM_ESRI.asc
Size is 1000, 1000
Coordinate System is ''
Origin = (380000.000000000000000,6888000.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Corner Coordinates:
Upper Left  (  380000.000, 6888000.000)
Lower Left  (  380000.000, 6887000.000)
Upper Right (  381000.000, 6888000.000)
Lower Right (  381000.000, 6887000.000)
Center      (  380500.000, 6887500.000)
Band 1 Block=1000x1 Type=Float32, ColorInterp=Undefined
  NoData Value=-9999

Sé que puedo usar gdal_contour para generar los contornos ( mi publicación de blog sobre el tema ) pero me pregunto cuáles son algunas de las mejores prácticas para generar contornos.

¿Hay alguna regla que deba seguir para aprovechar al máximo los contornos pero no inventar o perder demasiada información?

Diga que quiero generar tres conjuntos de contornos:

  • 250 mm
  • 1m
  • 5m

¿Debo hacer algo con el DEM antes de cada serie?

¿El suavizado posterior de las líneas es una buena manera de avanzar o el suavizado de la trama es una mejor opción?

    
pregunta Nathan W 28.09.2011 - 00:34

3 respuestas

14

Las reglas cartográficas para representar el relieve como contornos se presentan en El famoso Imhof Libro sobre representación en relieve, capítulo C . Algunas de estas reglas se encuentran en esta página de wikipedia . La recomendación principal al simplificar los contornos para preservar las características principales del terreno.

Suavizar los contornos de forma independiente no evita que se superpongan: se recomienda suavizar el DTM primero. Un alisado gaussiano tradicional con un parámetro adecuado en función de la resolución objetivo permite borrar los pequeños detalles. Un inconveniente es que el suavizado DTM llena los valles y depresiones, y aplana las crestas y picos. El uso del algoritmo de filtro douglass-peucker como en este documento puede ser una solución. También hay número de métodos basados en el uso de algoritmos de suavizado de contorno restringidos por la red de drenaje o un esqueleto . Finalmente, para evitar que los contornos se superpongan en partes descuidadas, es posible borrarlas localmente o también deformarlos .

    
respondido por el julien 28.09.2011 - 09:34
1

Quiero secundar el comentario de @ whuber. El análisis cuantitativo siempre es mejor desde un DEM directamente y el análisis visual a menudo es mejor cuando se realiza desde un sombreado en lugar de contornos.

Para responder a la pregunta directamente: En ArcGIS usaría Focal Statistics o Aggregate [Spatial Analyst Toolbox] para suavizar las líneas de contorno resultantes. Debido a que los contornos son una función de análisis visual, la cantidad de suavizado variará según su necesidad. Así que tendrás que experimentar y ver qué funciona mejor para tu proyecto.

El suavizado de líneas después de generarlas funciona, pero es un poco torpe en comparación con la modificación del ráster primero. Una limpieza de generación de contorno posterior que a menudo hago es seleccionar las líneas de cierta longitud (por ejemplo, < 10 'de largo) y eliminarlas. Esto elimina los datos de pequeños bits "ruidosos" de contornos cerrados (es decir, pequeños bucles) que probablemente no reflejen el carácter real de la superficie que están siendo modelados por los contornos.

Otras cosas que podría ver [usuarios de ArcGIS] son:

Hacer un mosaico ráster y / o un mosaico ráster de estadísticas focales y usar un modelo para generar contornos para todo el conjunto de datos.

Realización de contornos 3D para el uso de AutoCAD.

Tenga cuidado con las áreas de datos inusuales de DEM (por ejemplo, grandes extensiones de terreno bajo que necesitan diferentes intervalos de contorno para representar con precisión la superficie, áreas de vegetación densa que dan valores falsos de Bare Earth, áreas de relieve vertical - acantilados, etc.).

    
respondido por el user23715 08.11.2013 - 17:32
0

Hay una manera fácil de usar gdal_contour. Después de configurar todas las opciones en la ventana de diálogo, puede editar la línea de comandos y en su lugar el "-i intervalo" puede usar niveles fijos "-fl niveles". Como muestra la imagen a continuación. Puede consultar otras opciones aquí enlace

    
respondido por el christian.gobel 14.02.2017 - 18:02

Lea otras preguntas en las etiquetas