¿Extracción de elevaciones de DEM ráster a tablas de atributos vectoriales en ArcGIS Desktop?

Tengo un proyecto ArcMap 10 con múltiples grupos de capas. Una capa es un DEM ráster, otro grupo de capas es un montón de archivos de formas vectoriales (incluidos puntos, líneas y polígonos). El cliente quiere que todos los datos vectoriales contengan valores Z.

Rellenar las tablas de atributos para los datos vectoriales con valores X e Y es bastante fácil, pero ¿cómo relleno los valores Z?

¿Hay alguna manera de usar la Calculadora de campo para "tomar" los valores de altura de la capa DEM para cada entidad vectorial?

Debo agregar que no tenemos una licencia para las extensiones 3D Analyst o Spatial Analyst.

En ArcGIS, solo he hecho esto con vectores de puntos; no estoy seguro de cómo lo haría con polígonos o entidades de línea. Sin embargo, lamentablemente sí necesita Spatial Analyst. Hay una herramienta en Spatial Analyst Tools \ Extracción \ Extraer valores múltiples a puntos y Extraer valores a puntos . Extraer valores múltiples a puntos es mejor porque le permite cambiar el nombre del campo z al que está extrayendo sus valores. Extraer valores a puntos no le permite hacer esto. Tengo problemas para hacer esto en ArcMap (la herramienta se bloquea), pero funciona bien en ArcCatalog.

La extensión (gratuita) Hawths Analysis Tools para ArcMap tenía una herramienta de punto de intersección, pero el desarrollo de la extensión se suspendió hace un par de años. Ha sido reemplazado por el entorno de modelado geoespacial , que nunca he usado.

SAGA-GIS (software independiente gratuito) tiene dos funciones en Shapes: Cuadrícula \ Agregar valores de cuadrícula a puntos y formas - Cuadrícula \ Agregar valores de cuadrícula a formas que probablemente puedan hacer lo que desee. Esta es la descripción de la herramienta Agregar valores de cuadrícula a formas:

Unión espacial: recupera información de las cuadrículas seleccionadas en las posiciones de las formas de la capa de formas seleccionada y la agrega a la capa de formas resultante. Para los puntos, esto es similar al módulo 'Agregar valores de cuadrícula a puntos'. Para líneas y polígonos, los valores promedio se calcularán a partir de celdas de cuadrícula interferentes. Para los polígonos, el módulo 'Estadísticas de cuadrícula para polígonos' ofrece opciones más avanzadas.

SAGA es tu mejor apuesta, en mi opinión.

¡Espero que esto ayude!

No he probado esto para polígonos, pero funciona muy bien para puntos. En Arcmap 10, si tiene acceso a una herramienta de análisis 3D, puede buscar la Forma de interpolación de herramientas. Luego, selecciona su archivo ráster y su capa vectorial ajusta los valores z si desea elevaciones en pies y su DEM está en metros y deja que se ejecute . Una vez que esté completo, busque la herramienta AGREGAR XY y seleccione su capa recién creada y luego calculará sus valores X, Y, Z en la tabla de atributos de esa capa.

tal vez podría usar gdallocationinfo para recuperar todos los valores Z pasándole los valores X, Y. una vez que tenga los valores Z, use la calculadora de campo para agregar el valor Z a los puntos que conforman cada característica. Habría algunas secuencias de comandos involucradas, pero en teoría esto debería funcionar.

Dependiendo de sus unidades o resolución requerida, siempre puede exportar el ráster DEM a un archivo de forma de polígono, directamente desde la caja de herramientas, y luego hacer una unión espacial; esto es más fácil para los puntos. Si necesita valores Z para líneas o polígonos, obviamente necesitaría segmentos definidos o subdivisiones de polígonos en cada cambio de elevación según lo informado por la cuadrícula DEM vectorizada. Eso se puede hacer a través de herramientas de geoprocesamiento estándar y uniones espaciales, pero debe prestar mucha atención a sus requisitos de resolución y la cantidad de procesamiento que consumirá un gran conjunto de datos.

Recomiendo obtener una licencia de 3D Analyst y usar Agregar información de superficie . O use Interpolate Shape seguido de Add Z Information, que permite que algunas funciones devuelvan valores NA sin bloquear todo el programa. De lo contrario, hay dos métodos gratuitos probados que requieren más trabajo preliminar. Es posible implementarlos desde QGIS , de modo que no tenga que abandonar sus cómodos entornos SIG.

La herramienta ESRI que desea utilizar es la "Agregar información Z". Quizás pueda obtener una versión de prueba y ver si la herramienta le sería de mayor utilidad. Si obtener la extensión 3D Analyst no es una opción, lo único que se le viene a la mente (aparte de una solución de programación) es extraer sus datos de ESRI y usar una de las opciones gratuitas como SAGA-GIS / qGIS su "Agregar cuadrícula El valor de las formas "o" Estadísticas de cuadrícula para polígonos "no es tan flexible como el de los ESRI, pero hace lo mismo si busca la altura promedio a lo largo del polígono / línea.

Sé cómo hacer esto por puntos. Las líneas y los polígonos son diferentes porque cruzan varias celdas de la cuadrícula del DEM. Puede hacer puntos a lo largo de las líneas o algo así y extraer valores para ellos. En Arc, necesitará un analista espacial para obtener los valores, pero puede obtenerlos a través de Qgis (gratis), aunque no es tan útil. Deberá ser hábil para trabajar con archivos de Excel y convertirlos de un lado a otro de los archivos de forma. Qgis generará un archivo de forma que solo tenga las elevaciones de puntos extraídas, por lo que deberá copiarlo y pegarlo en su tabla original como un archivo de Excel (con los valores X, Y, para que pueda volver a un archivo de forma con "mostrar datos x, y"). Deberá realizar un seguimiento de sus proyecciones para asegurarse de elegir la correcta al convertir de nuevo a shapefile.