SIG (Arcmap), Ejercicios de Topografía. Universidad Politécnica de Madrid (UPM)
lieblan
lieblan

SIG (Arcmap), Ejercicios de Topografía. Universidad Politécnica de Madrid (UPM)

PDF (212 KB)
9 páginas
1Número de descargas
2Número de visitas
Descripción
Asignatura: Topografía, Sistemas de Información Geográfica y Teledetección, Profesor: Lozar Lozar, Carrera: Ingeniería Forestal, Universidad: UPM
20 Puntos
Puntos necesarios para descargar
este documento
Descarga el documento
Vista previa3 páginas / 9
Esta solo es una vista previa
3 páginas mostradas de 9 páginas totales
Descarga el documento
Esta solo es una vista previa
3 páginas mostradas de 9 páginas totales
Descarga el documento
Esta solo es una vista previa
3 páginas mostradas de 9 páginas totales
Descarga el documento
Esta solo es una vista previa
3 páginas mostradas de 9 páginas totales
Descarga el documento

CLASE 10 – VECTORIAL (EJERCICIO)

1. 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-1

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

CLASE 11 – VECTORIAL (EJERCICIO) Objetivo: Determinar el trazado de una línea de Alta Velocidad que provoque un impacto ambiental menor. Para calcular el impacto se han de tener en cuenta los criterios siguientes: Afección a lugares de interés comunitario (LICs) Afección a puntos de interés geológico (PIGs) Impacto sobre la vegetación en función de su calidad

Información disponible: Lics del Mediterráneo Puntos de interés geológico del sureste de la Comunidad de Madrid Hojas del mapa forestal español (mfe) de la zona Tabla reclasificadora que relaciona la especie vegetal del mfe con la calidad ambiental Alternativas para el trazado del tren de alta velocidad. Teniendo en cuenta lics Que alternativa de la capa AVE (tres líneas) afecta menos a los lics?-> Tenemos que decidirlo fijándonos en cuanta superficie de cada línea de ave se superpone más o menos con los lics (Cuanta longitud). Para ello rompo cada una de las alternativas por los puntos de intersección con la capa lics para saber cuánto superpone cada una, operación->intersección Geoprocessin->intersecction->input features: ave, lics Abro la table de atributos de la nueva capa, añadir campo->tipo double->botón derecho sobre el campo->calculate geometry Agrupo por el campo trazado, botón derecho sobre campo trazado->summarize->long->sum (Suma las superfices que corresponden a cada línea de ave) Viendo esta tabla vemos que la peor es el trazado 3 y la mejor el trazado 1. Ahora tengo en cuenta el impacto de la vegetación en función de la calidad Voy a ponderar la longitud multiplicándola por la calidad ambiental EN la tabla reclasificadores se le asigna un valor de calidad ambiental a cada unidad ambiental. Cuando mayor sea la calidad mas perjudicial será el impacto de las líneas del ave. Tenemos que unir una tabla a una capa -> join. Para que join funcione necesitamos un campo en común. El campo Rotulo1 en los mapas coincide con el campo Rotulo en la tabla. El mapa está dividido en seis partes asi que primero unimos las capas de capas para tener solo una. Herramienta->merge; geoprocessing->search for tools->merge-> Input datasets: selecciono todas las capas de mapa. En las zonas donde se juntan las capas queda una línea, frontera artificial que no quiero tener, para quitar esa línea y unir las zonas que son las misma, utilizo DISSOLVE (es lo mismo que summarize en las tablas pero con capas) Geoprocessing->dissolve->input: mfe->dissolve fields: Rotulo 1-> quitar tick de créate multipart features. Para calcular el impacto sobre la vegetación en función de la calidad hago la suma para alternativa (cada línea de ave) de la longitud por la calidad de la zona. Si atravieso zonas con calidad más alta el impacto será mayor. [Impacto= sumatorio (longitudxcalidad)] Botón derecho capa mfe1->join->uno mfe1 campo rotulo1 con tabla reclasificadora campo rotulo Al abrir la tabla de atributos tras hacer la unión salen campo <Null>, esto es porque a las zonas urbanas no se les han atribuido valores. En estas zonas el campo valor va a ser cero, luego el impacto en estas zonas (urbanas) el impacto será cero. Guardo la capa nueva con el join hecho (para que el join sea permante), botón derecho sobre nombre de capa mfe1->data->export data->guardar como shapefile Quiero romper las líneas de ave por las intersecciones con los distintos tipos de vegetación, esto lo hago con la herramienta INTERSECCION Geoprocessing->intersect->Input features: ave y calidad->se crea una capa nueva (Ave_Calidad)

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 5

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)

Tenemos una capa y una tabla (con el valor para ponderar), nos gustaría hacer un join, y hay uncampo común así que hacemos un join (campo común MATERIALES). Boton derecho sobre nombre de la capa-> joins and relates-> join. Join no es permanente así que para hacerlo permamente data->export data o añadir un campo nuevo en la cama, copiar la columna que he añadido con el join y después desacer el join (elegimos esta opción). Para copiar los valores de la columna valor en la nueva columna Av. Botón derecho sobre columna que quiero rellenar (Av) -> field calculator-> = [tipo_suelo.VALOR] (para que la columna Av sea igual a la columna Valor). Deshago el join (joins anda relates-> remove joins. Se me ha quedado un nombre muy largo que era lo que idenficaba durante el join que ese campo pertenecia a la capa de geología. Una vez deshecho el join para evitar tener ese nombre tal largo descargo la capa (remove sobre la capa) y la vuelvo a cargar (con deshacer) entonces me aparecerá el nombre corto que es el que quería poner al inicio. Topografía (Tv) Tenemos que relacionar la capa topografía con la tabla rectas_pendientes . Lo hago de la misma manera que en el caso DV. Creo un campo Tv y para cada intervalo de pendientes le añado un valor al campo Tv, según me muestra la tabla de Rectas_pendientes Tengo tres capas cada una con uno de los sumandos que necesito. Si quiero poder hacer la suma necesito tenerlo todo en una sola capa. Tengo tres opciones: unión, intersect y merge. En este caso no hay diferencia entre unión e intersección porque el contorno es el mismo. Cuando creo una unión se quedan todas las combinaciones posibles en la geografía y en la tabla (hago esto). No hago merge porque merge no hace combinaciones, junta las capas pero no las combina (no divide las cosas). Una vez he unido las tres capas en una añado un nuevo campo V (para indicar la vulnerabilidad). Botón derecho encima de V->field calculator-> 5* [Dv]+3* [Av]+ 1*[Tv]

2. Obtener la vulnerabilidad media de los acuíferos para los dominios hidrogeológicos

Calcular una media de la vulnerabilidad para los dominios hidrogeológicos. Cargamos la capa de los dominios Primero tenemos que hacer un intersect entre la capa de dominios y la capa de DAT. Para simplificar podemos usar summarize o dissolve (no me da una tabla y otro una capa con su tabla) En mi tabla tiene que aparecer vulnerabilidad y dominio, necesito que estén las dos juntas. Utilizo un intersect porque las capas no tienen la misma extensión por lo que no da lo mismo que utilizar unión. (se crea capa dom_DAT) Como quiero agrupar por dominio (por nombre) le doy con el botón derecho para hacer summarize sobre el campo nombre->V->Average (media); Se crea una tabla.

Si nos pidiera cuanta superficie tiene DAT>50

Primero tengo que calcular la supercifie creando un nuevo campo->Calculate geometry- >Area Selecciono los polígonos que tienen V>50 Ahora puedo utilizar summarize y statistics. Con statistics: botón derecho encima de campo superficie->statistics->me da las estadisiticas de lo que tengo seleccionado. Uno de los valores que me da es la suma, cojo ese valor. [Casi seguro que cae en examen] Con summarize: necesito un campo que nos agrupe por vulnerabilidad normal y extrema, aun no lo tenemos pero lo podemos añadir. Añado campo Riesgo. Botón derecho sobre nombre del campo->field calculator-> “Alto”; invierto selección y con field calculator añado campo “normal”. Desselecciono, botón derecho sobre el campo riesgo-> summarize->Area- >sum

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.

Zonal statistics as table (geoprocessing->search for tolos->zonal statistics as table) Input raster: La zona para la cual queremos calcular las pendientes (monte) Zone field: rodal Input value raster: capa de pendientes Nos sale una tabla con las estadísticas de cada rodal.

Cargo capa punto_origen Vamos a hacer una capa raster de distancias a los puntos pertenecientes a la capa que acabo de añadir. Alrededor de cada punto el ordenador hace una malla y a cada pixel le da un valor de la distancia a ese punto. Euclidean distance (geoprocessing->search for tolos->Euclidean distance)

Input rater: capa a la cual pertenecen los puntos Input cell size: 20 en este caso Para controlar hasta que distancia quiero llegar. geoprocessing->search for tolos->Euclidean

distance->Environment (esto lo tienen todas las herramientas raster)->Processing Extent->Extent(le digo hasta donde quiero que calcule, en este caso si selecciono la capa de pendientes le estoy diciendo que calcule hasta donde llega la extensión de la capa pendientes)

geoprocessing->search for tolos->Euclidean distance->Environment (esto lo tienen todas las herramientas raster)->Raster Analysis (cambiar el tamaño de celdilla)->[En este caso: cell size-> as specified below:20, que es lo mismo que hemos hecho antes] Imaginemos que la capa de puntos son nucleos urbanos y nos piden calcular el riesgo de incendio mediante la fórmula Rincendio=R*dist. Puedo hacer un mapa nuevo donde sus pixeles tengan como valor R*dist Raster calculator: me permite hacer operaciones con las capas (geoprocessing->search for tolos->Raster calculator)

Multiplico la capa riesgo por la capa distancia: "Riesgo" * "distance1” Si multiplico dos capas de diferente tamaño siempre la resultante va a tener la extensión de

la capa más pequeña Hacer un “clip” en raster: Raster calculator->Enviroment->Raster analysis->mask: monte (la capa que quiero utilizar de molde)-> En extent también pongo la capa que quiero utilizar de molde (monte). Al calcular decimos que la capa resultante que queremos es igual a la capa mediante la cual queremos recortar.

comentarios (0)
No hay comentarios
¡Escribe tú el primero!
Esta solo es una vista previa
3 páginas mostradas de 9 páginas totales
Descarga el documento