GUIA PARA OBTENER, CLASIFICAR, GENERALIZAR LAS PENDIENTES EN ARCGIS – PRO
Datos de partida:
- DEM – en este caso con un GRID de 3x3 de resolución espacial (SIG-Tierras)
- Un marco (polígono) que será la máscara de recorte del RASTER
OBTENER PENDIENTES
1. Verificar el sistema de referencia de los datos de partida, deben ser el mismo; en este caso WGS84-UTM
Zona 17 SUR.
2. Cargamos el fichero Raster del DEM y el polígono de recorte
El área de interés es el polígono rojo, pero se debe hacer un recorte con una cierta tolerancia mayor para
que las pendientes de los bordes de nuestro polígono no se vean afectados por cortes abruptos; para ello
hacemos un marco (polígono verde) para establecer nuestra área de trabajo
Crear el nuevo polígono
Para recortar aplicamos la herramienta Extract by Mask, como se muestra a continuación, el Imput raster
es el raster que vamos a cortar; el Input raster of feature mask data es el polígono Marco.
Obtenemos el área de interés, se ve modificado los tonos del DEM, eso es por que cambiaron los valores
máximos y mínimos.
Trabajamos sobre el Raster cortado en el paso anterior, al que le aplicaremos un SLOPE (pendientes), le
definimos un nombre de salida y podemos elegir entre medir la pendiente en Grados o Porcentaje, en este
caso lo haremos con porcentaje, para ajustarnos al formato de FAO que pide en %.
RECLASIFICACIÓN
Obtenemos un las pendientes en porcentaje, pero debemos hacer una re-clasificación para ajustarnos a lo
que nos dice FAO
Cuadro de FAO clasificación pendientes
Le damos click derecho sobre la capa creada y vamos a Simbology
Se nos habilita un panel Symbology, aquí primero seleccionamos en el método Natural Breaks (Jenks)
Para luego modificar, el número de Clases (en este caso 6) y los valores mayores de cada clase en relación a
nuestra referencia de la FAO, como se muestra a continuación.
Pero ahí no esta realizada la clasificación aún, ahora vamos ArcToolBox y buscamos la herramienta
Reclassify, seleccionamos o arrastramos el Raster en el que estamos trabajando y se nos carga la
herramienta con los datos que ya hemos establecido en el paso anterior.
Obtenemos el Raster con las categorías o clases establecidas
GENERALIZACIÓN
Para eso hay que considerar el área mínima cartografiable
El área mínima cartografiable (AMC) es una de las premisas a tener en cuenta a la hora de digitalizar bajo
una escala fija identificando elementos que definan una unidad mínima de superficie. Un concepto que no
sólo te ayudará a gestionar los elementos potencialmente cartografiables sino que también servirá de guía
para digitalización técnicamente coherente.
Las relaciones para identificar el área mínima cartografiable a partir de una escala se toman partiendo de
una delimitación óptima legible que permita obtener un píxel de 4×4 mm. Aquí tienes una referencia para
la identificación del área mínima cartografiable a partir de una escala fija de trabajo definida en base a
autores como Rossiter o Salitchev.
Primero se ejecuta el filtro “Eliminate” para generalizar áreas no significativas para el área de trabajo de
acuerdo al área mínima cartografiable que para la escala de trabajo 1:4000, AMC = (4×4) × (4×4) = 256 m2
como indica Rossiter.
En este caso por la escala 1:4000 que estamos trabajando y aplicando la formula, nos resulta un área mínima
cartografiable de 256 m2.
El programa lo que hará es fundir (unir) a los polígonos más grandes a los polígonos pequeños que cumplan
con la condición del área mínima cartografiable.
Se lo puede hacer con FILTROS RASTER, en ese caso lo vamos a hacer convirtiendo el Raster a Vector
Poligonal. El Geoproceso es Raster to Polygon
Obtenemos el siguiente resultado
Para generalizar el mapa tenemos primero hacer un corte para delimitar el área de trabajo, que sería el
polígono Rojo.
Para eso realizaremos un CLIP,
Ahora para generalizar lo que tenemos, que calcular el área de los elementos del SHP, abrimos la tabla de
atributos y le damos en Add Field, le colocamos un nombre y el tipo de dato float. Se daclick derecho sobre
la tabla que se está editando y se guarda los cambios.
Luego calculamos el área dando click derecho sobre la cabecera de la columna y seleccionando Calculate
Geometry y se coloca los siguientes parámetros.
Resultaran algunos polígonos el área mínima cartografiable de 256 m2, se debe hacer una consulta por
atributos:
Le podemos cambiar para que se pueda hacer la consulta SQL, para poder variar más fácilmente el valor del
área de ser necesario
Ahora si se aplica la GENERALIZACIÓN, para lo cual abrimos la herramienta Eliminate, en el cual cargamos en
el Input la capa con los elementos seleccionados.
Se fusionan los polígonos más pequeños en los más grandes aledaños, actualizamos nuevamente el campo
de área hectárea del nuevo shape creado, hacemos la consulta otra vez y ejecutamos las veces que sean
necesarias este proceso para que se eliminen todos los polígonos con el menores o iguales al área mínima
cartografiable.
Ahora tenemos que hacer un DISOLVE, para fusionar los elementos en función del campo gridcode – que son
las clases que definimos en pasos anteriores.
En la tabla de atributos del shp podemos ver que solamente existen tres elementos que son la fusión de
todos los pequeños polígonos en cada una de las clases que les corresponde. Le añadimos nuevamente el
campo área_ha para conocer la superficie y tambien la clase en texto , para que sea más detallada.
Se debe aplicar una simbología dando clic derecho al shp Symbology Unique Values y se selecciona
rampa de colores que más se ajuste a la temática.
PERFIL TOPOGRÁFICO
Un gráfico de perfil es uno de los gráficos disponibles en ArcGIS Pro. Los gráficos de perfil muestran el
cambio en la elevación de una superficie a lo largo de una o varias líneas. Los gráficos de perfil pueden servir
de ayuda para evaluar la dificultad de un camino o la viabilidad de ubicar una línea ferroviaria a lo largo de
una ruta determinada. Puede crear gráficos de perfil en ArcGIS Pro a partir de clases de entidad de línea 3D o
a partir de entidades de línea 3D digitalizadas de una superficie en un mapa.
Creamos primero una polilinea, y habilitamos Has Z, para que tenga un componente en 3D y tome los datos
de elevación del DEM de elevaciones.
Una vez creado el shp tipo polilinea, lo editamos. Se selecciona el elemento (línea) y en Edit, seleccionamos
Create, elegimos la línea
graficamos sobre el DEM que contenga las elevaciones que queremos representar en el perfil topográfico.
Finalizamos la edición
Ahora debemos asignarle los valores de elevación a la línea, para lo cual aplicamos Interpolate Shape
y damos click derecho el fichero de salida create chart Profile Chart
Y así obtenemos el perfil topográfico o de elevaciones de la zona que no interese, la podemos editar a
nuestro gusto.