¡Descarga SIG (Arcmap) y más Ejercicios en PDF de Topografía solo en Docsity!
CLASE 10 – VECTORIAL (EJERCICIO)
- Reclasificar el mapa de vegetación que se proporciona simplificándolo de acuerdo a las unidades que se muestran en la siguiente tabla: Vegetación antigua Vegetación nueva Pinus sylvestris Pinares Melojares Quejigares Frondosas P. sylvestris + Q. pyrenaica Mezcla de pinos y frondosas Vegetación riparia Vegetación riparia Piornales Matorral acidófilo (escobas) Matorrales Lastonares Cervunales Pastizales mesófilos Pastos xerofíticos Pastos reticulares Pastos Afloramientos rocosos Afloramientos rocosos Urbanizados Urbanizados La reclasificación se hará de tres maneras diferentes: a) A mano, editando Primero añado un campo con add field (vege_n), escribo en cada campo uno a uno la información Selecciono por atributos Pinus sylvestris, después en el campo Vege_n -> Field calculator -> “Pinares” Selecciono por atributos melojares y quejigares: "VEGETACIóN" = 'Melojares' OR "VEGETACIóN" = 'Quejigares' después en el campo Vege_n->Field calculator-> “Frondosas” Así con cada campo b) Generando la tabla en Excel y vinculándola a la capa Herramienta ‘Join’, la tabla de Excel tiene que cumplir tener un campo en común con la tabla de atributos a la que la queremos unir. ver taba subida a Moodle vege vege_n Pinus sylvestris Pinares Melojares frondosas Quejigares frondosas P. sylvestris + Q. pyrenaica mezcla de pinos y frondosas Vegetación riparia vegetacion riparia Piornales matorrales Matorral acidófilo (escobas) matorrales Lastonares pastos Cervunales pastos Pastizales mesófilos pastos Pastos xerofíticos pastos Pastos reticulares pastos Afloramientos rocosos afloramientos rocosos Urbanizados Urbanizados
No se puede combinar las celdas de la tabla como hacemos en Word porque no funcionaría en arcgis Los elementos del campo en común tienen que estar escritos exactamente igual que en Argis, si no no funciona. Abrimos la hoja de Excel en arcgis Botón derecho sobre la capa vegetación, joins and relates-> join… ; Elegimos el campo de cada capa mediante la cual queremos hacer el Join (Vegetación en la capa vegetación y vege en la capa Hoja1). Join no es permanente, si quiero hacerlo permamente tengo que exportar la campa como una nueva. Botón derecho sobre nombre de la capa vegetación->data->export data. Guardar en formato shapefile. [Caerá en el examen] c) Generando una tabla simplificada en ArcGIS Hago la tabla reclasificadora que hemos hecho en Excel desde arcgis utilizando la herramienta summarize -> Cuando utilizo esta herramienta obtengo una tabla. Doy clic con el botón derecho sobre el campo mediante el cual voy a agrupar. Guardo la tabla en formato dBASE Table (siempre guardo las tablas en este formato). Añado la tabla a argis y hago join con la tabla de atributos de la capa de vegetación (¿)
2. Generar unidades paisajísticas como combinación entre las cuencas hidrográficas (zonas con intervisibilidad) y las unidades de vegetación simplificadas anteriormente. Paisaje = cuenca + vegetación, tengo que utilizar intersect o unión, me valen ambos porque el contorno de ambas capas es el mismo. Sin tener nada seleccionado->Geoprocessing->Union-> seleccionar las dos capas a unir (cuencas y Vegetasimple) y las arrastro a Features. Tengo toda la información de vegetación y cuencas en una sola capa, puede representarlos por ambas formas. Si abrimos la tabla de atributos hay polígonos que no tienen vegetación asociada, eso es imposible, estos polígonos son errores. Selecciono los polígonos errores, puedo cargármelos o unirlos al de al lado (mejor esta opción)- >search for tools-> eliminate Voy a añadir un campo que combine cuenca y vegetación para poder representar por tipos de paisaje. Boton derecho nombre del campo creado ->field calculator-> [VEGETACIóN] &"-" & [CUENCA] "-" representa campos variables -> pinares- a) ¿cuantas unidades paisajísticas diferentes hay? Miro la tabla de atributos, abajo del todo pone que tengo 66 elementos en total. b) ¿Qué superficie tiene cada una? Para agrupar clases por un campo (una clase para queda tipo de vegetación) utilizo summarize o dissolve, para que se repita una vez cada tipo de paisaje. Primero tengo que calcular la superficie de cada polígono. (botón derecho sobre el campo nuevo sum, calculate greometry, sum) Despues hago un summarize dando botón derecho sobre el campo paisaje y en las opciones seleccionando campo sum. Me crea una tabla nueva
Compruebo que tengo los campos que necesito, la calidad ambiental está en el campo valor y para saber a que línea de ave corresponden tengo el campo trazado, me falta la longitud, la calculo añadiendo un campo y calculate geometry Calculo el impacto multimplicando la columna valor por la columna longitud, para eso añado un campo nuevo de tipo doble (I_CA)-> botón derecho sobre nombre del campo-> fiel calculator-> [VALOR]* [long] Ahora quiero un valor de I_CA para cada trazado, lo hago mediante summarize->botón derecho sobre trazado->summarize->despliego I_CA y le doy a Sum->lo guardo en formato dbaseTable El que me da menor valor de impacto es la línea de ave 3 (trazado 3) Para decidir la mejor alternativa teniendo en cuenta los dos casos, puedo sumar los dos impactos. Antes de eso tengo que unificar escalas, no es lo mismo sumar longitudes que longitudes multiplicadas por calidad ambiental. Hay muchas técnicas para esto, una de ellas es dividir por el valor más alto en cada uno de los casos, tendríamos un rango entre cero y uno. Teniendo en cuenta los puntos de interior geológico Lo hago calculando el área de influencia. Puede hacernos un área de influencia para todas las líneas junas o una para cada línea. Escogemos la de una para cada línea porque es lo que nos interesa en este caso para poder decidir. Geoprocessing->buffer->linear unit:250->dissolve type:none me ha salido demasiado pequeña el área de influencia porque no he tenido en cuenta la distancia a los puntos de interés. Vuelvo a hacer buffer con un área de influencia de 6000m Si quiero saber cuántos puntos coinciden con cada punto de influencia tengo que hacer intersección. Geoprocessing->intersect->entre ave_6000 y pigs->se me crea una capa nueva Para saber cuantos puntos tengo por trazado (en la capa nueva que se me ha creado) hacemos un summarize por trazado (botón derecho sobre trazado) y no elijo nada mas porque no quiero sumar nada solo quiero un conteo. Se me crea una tabla nueva donde sale cuantos puntos corresponden a cada trazado.
CLASE 12 – VECTORIAL (EJERCICIO)
En este ejercicio se pretende analizar la vulnerabilidad de los acuíferos de agua subterránea en la comunidad autónoma de la Rioja. Para ello se va a utilizar una simplificación del modelo DRASTIC, que llamaremos modelo DAT. Este modelo considera los siguientes parámetros: Parámetros Ponderación D Profundidad del nivel del agua
A Tipo de suelo 3 T Topografía (pendientes) 1 Cada uno de estos parámetros tiene un rango de valores en función de su vulnerabilidad a la contaminación. Vulnerabilidad = Dv ∙ Dp + Av ∙ A (^) p + T (^) v ∙ T (^) p Siendo v el valor del parámetro y p la ponderación. La vulnerabilidad será mayor cuanto mayor sea el número resultante en la ecuación anterior. [Dp=5; Ap=3; Tp=1; El resto van a variar] Información de partida: Capas: La_Rioja.shp : límite administrativo de La Rioja Niveles_Agua.shp : niveles de agua subterráneos de La Rioja Geologia_CHebro.shp : geología de la cuenca del Ebro Pendientes.shp: capa con polígonos de pendiente similar Dom_HG_Ebro: capa con los dominios hidrogeológicos Tablas: Nivel_agua.xls : tabla con los valores que adopta el modelo para cada clase de nivel de agua. Tipo_suelo: tabla con los valores que toma el modelo en función de los materiales. Topografia.xls : tabla con los valores que toma el parámetro topografía en función de las clases de pendientes. Se pide:
1. Realizar el análisis DAT de La Rioja. Nivel de agua (Dv) Dv=Niveles de agua; Queremos quedarnos con la información de niveles agua con el contorno de la rioja clip->input: niveles agua; clip: la rioja. Quiero relacionar una capa y una tabla-> join; pero no tengo campo común así que o creo un campo común para poder hacer join y añado los datos a mano porque son pocos. (lo hago por la segunda opción) Creo nuevo campo Dv y le añado el campo Valor de la tabla nivelagua.xls hoja1$. Este valor pondera la variale en función de la profundiad (nivel de agua=profundidad) -> "NIVEL_AGUA"
1,5 AND "NIVEL_AGUA" <4,5 así para todos los valores de profundidad asigno un valor que luego ponderará el parámetro Nivel del agua. El primer sumando el vulnerabilidad es 5xDv, para tener este primer sumando tendríamos que multiplicar en la capa NivelesAgua_Clip 5 por el campo Dv Tipo de Suelo (Av)
RASTER. Herramientas fundamentales
Lo primero es conocer el tamaño del pixel. El tamaño depende de la precisión (precisión (mm)=0,2*escala). Siendo 0,2 el límite de percepción visual.
MDE: modelo digital de elevaciones. Representación continúa del terreno en altitud. Nos va
a servir para hacer otros mapas. Pendientes (slope) Orientaciones (aspect) Sombras (hillshade) Cuencas visuales Volver de un MDE a curvas de nivel (contour) Para pasar de curvas de nivel a MDE tenemos dos comandos: Create Tin y Topo to raster Recordatorio: no se pueden utilizar carpetas con espacios en sus nombres, espacios, tildes, etc
Pasar de curvas de nivel a MDE -> argis interpola para hacer esto de dos diferentes maneras
Create Tin : Triangula entre las diferentes coordenadas creando una red de muchos triángulos que conforman la superficie. (geoprocessing->search for tolos->Create Tin) [ A veces la herramientas para trabajar con raster no están actives, hay que ir a Customize-
extensions->activamos 3D Analyst ] Input feature Class: añadimos las capas de las cuales va a sacar la información (pint_Merge y curvas_Merge). En Height Field seleccionamos el campo del cual tiene que sacar la información (el campo que tenga la altitud) en este caso GGroup. Topo to raster : parte de una capa de pixel, cada pixel tiene un nivel de elevación. Va a interpolar con otro algoritmo Input feature Class: añadimos las capas de las cuales va a sacar la información (pint_Merge y curvas_Merge). En Field seleccionamos el campo del cual tiene que sacar la información (el campo que tenga la altitud) en este caso GGroup. En type en la capa de puntos seleccionamos PointElevation. Suponemos que estacas capas están a 1:100.000 luego En Output Cell size ponemos (0,2*100.000) 20. Leer los datos en una capa raster. Ver datos de la distribución de los datos: botón derecho sobre nombre de la capa->propiedades->Source; Podemos ver: Tamaño de pixel, extensión, min-max, media (mean). La opción de abrir tabla de atributos está desactivada. Si queremos saber el valor completo de un pixel le damos al botón de información. y clic sobre el pixel.
Mapa de pendientes : (geoprocessing->search for tolos->slope)
Input raster: la capa que hemos creado antes con topo to raster o create tin (modelo digital de elevaciones). La salida podemos ponerla en porcentaje o en grados. Aparece la capa pueda representada en unos intervalos determinados pero podemos cambiarlo. Boton derecho sobre la capa->propiedades->simbology->Classify: puedo cambiar cada cuanto hago los intervalos.
Mapa de orientaciones : (geoprocessing->search for tolos->aspect) Vamos a obtener una capa
donde figura la orientación, la orientación se mide en grados. Input raster: modelo digital de elevaciones Flat (-1) = Todos los vientos Tenemos 8 intervalos, para ello el norte está dividido en dos secciones, porque nos queda una parte de la sección “norte” a la derecha del norte geográfico y otra parte a la izquierda, por lo que arcgis tiene que representarlo en dos trocitos. Si queremos dividirlo en 5 intervalos los ángulos donde “corta” serían: 45, 135, 225, 315
Mapa de sombras : (geoprocessing->search for tolos->hillshade) Nos representa en una
gradación entre 0 y 255 la luz que recibe cada pixel teniendo en cuenta la altitud (diferenciar zonas de humbría y solana) Input raster: modelo digital de elevaciones
Volver de MDE a capas de nivel : (geoprocessing->search for tolos->contour)
En intervalo meter el valor del intervalo entre curcas de nivel (en este caso 20, se mira en el mapa)
Mapa de cuencas visuales: (geoprocessing->search for tolos->viewshed). Una Cuenca visual
es donde alcanza la vista desde un punto. Se puede hacer desde puntos o desde polilineas. Nosotros en este caso vamos a calcularlo desde puntos. Necesito tener un punto de referencia, añado la capa punto_origen Input raster: modelo digital de elevación. Observadores: capa con el punto, puntos, polilinea o polilineas de referencia Nos da las zonas visibles o no visibles desde esos puntos Abrimos la tabla de atributos de la nueva capa (CuencasVisu) el valor 0 es la zona no visible. Los valores 1, 2 y 3 es el número de puntos desde los cuales es visible un pixel. La herramienta no me permite distinguir desde cuales de los puntos es visible el pixel, para saberlo tendría que hacer el mapa de cuencas visuales punto a punto. El conteo (count) que aparece en la tabla de atributos es el número de pixeles que contiene cada valor, nos está dando la superfice = Nºpixels*tamaño pixel (20x20 en este caso) Aplicaciones: prevención de incendios, evaluación de impacto visual. Puedo gestionar si estoy mirando desde una torre o si quiero saber si veo un elemento que está más alto que el nivel del suelo. No necesito saber de memoria como puedo mirar la ayuda del programa (search->viewshed->Using Viewshed and observer points for visibility analysis) Puedo limitar el radio de la visual y modificar otros parámetros específicos del estudio que estemos haciendo. Para poder cambiar cualquiera de estos campos los tengo que añadir abriendo la tabla de atributos en la capa de puntos de origen. Añado un campo a la tabla (por ejemplo Radio (Radius2), le pongo un valor (ej 5000 -m-) y vuelvo a calcular la cuenca visual. Añadimos la capa pend_baja a argis. La capa no tiene tabla de atributos. Podemos clasificar mediente la simbología de la capa los intervalos que queremos. Botón derecho->propiedades-
simbología->Classified->classify-> pongo que haya 4 intervalos.
Reclassify (geoprocessing->search for tolos->reclassify) Seleccino la capa de entrada y el campo
que quiero reclasificar. Como no hay tabla de atributos así que reclasificamos los valores. En classify->break values->pongo los puntos por los que quiero que corte cada intervalo de valores En este caso lo que he hecho es convertir una capa de pendientes, a una capa de riesgo, donde las pendientes de 0-5 son riesgo 1, de 5-10 riesgo 2… etc. Como en la capa Riesgo no tengo rangos, sino un mismo valor para cada apartado, si que tengo tabla de atributos.
Raster to polygon : pasar de raster a vectorial. ((geoprocessing->search for tolos->Raster to
polygon) Input raster: la capa que quiero transformer. No podemos pasar nunca de raster a vectorial una capa que tenga un valor por cada pixel, necesitamos reclasificar primero para que tenga una capa de atrbutos. Field: El valor que queremos conservar de la tabla de atributos, normalmente es el valor del pixel. (VALUE, en este caso corresponde al valor del riesgo) Tras hacer esta operación el campo value pasa a llamarse GRIDCODE Cargo la capa monte, tiene un campo que se llama rodal en su tabla de atributos. Quiero sabe r cual es la pendiente media de cada rodal. Necesito la capa que delimita los rodales en polígonos y la capa de pendientes de la zona. Hago una media de las pendientes que corresponden a cada rodal.
