Memoria De Trabajo

PROYECTO SIG DE LA MONTAÑA ORIENTAL

CONSTRUCCIÓN DE UN MODELO DIGITAL DE ELEVACIONES MEMORIA DE TRABAJO

Fundación Naturaleza y Hombre Autor Fernando Gaviria Stewart, biólogo Año 2003

1

AGRADECIMIENTO Quiero agradecer muy sinceramente a Enrique Pascual la amabilidad que ha tenido al facilitarme las Demos de ArcView 8.2, gracias a lo cual he podido hacer mi trabajo sin ningún problema y sin interrupción. Le estoy doblemente agradecido porque he aprendido y he disfrutado mucho haciéndolo. En nombre de la Fundación Naturaleza y Hombre, para la que he trabajado, quiero agradecerle también que, como Director de la zona norte de la empresa ESRI GEOSISTEMAS, haya tenido la amabilidad de hacer las gestiones necesarias para que le fuera donada a la Fundación una licencia de ArcView. Por último el agradecimiento a ESRI GEOSISTEMAS y a ESRI (USA) por habernos tramitado y concedido la licencia de ArcView. Esperamos hacer un buen uso de tan valiosa ayuda.

2

Índice Pág. Agradecimiento

2

Presentación

4

Introducción

5

Objetivos

11

Materiales y métodos Cartografía Hardware Software Procedimiento 1. Escaneado de la cartografía original en papel

11 11 12 12 12 12

2. Georreferenciación de las imágenes raster 3. Digitalización manual 4. Creación del TIN y el Raster

13 14 15

Resultados 1. Obtención y tratamiento de las imágenes raster La resolución espacial La georreferenciación El solapamiento de imágenes 2. Digitalización manual 3. Creación del MDE El TIN El Raster

16 16 17 18 20 20 22 22 23

Conclusiones

27

Bibliografía

29

Algunos recursos en Internet

29

Figuras

31

3

Presentación La finalidad de esta memoria es exponer la metodología que he seguido y los resultados que he obtenido, pensando sobre todo en la persona o personas que deban seguir adelante con el proyecto, puesto que no ha quedado terminado aún y presumiblemente yo no podré llevarlo a término. He procurado explicar con detalle todo lo que he considerado importante en cada paso que he dado, intentando dejar lo más claro posible cuál era su propósito y los problemas que han ido apareciendo. He querido hacer especial énfasis en los errores que sé que he cometido (seguro que hay más que ignoro), porque creo que la mayoría pueden y deben corregirse. Creo que el trabajo en general ha estado correctamente planteado y que puede tener una buena continuidad, pero quizá puedan introducirse cambios si se consideran necesarios. Para realizar el trabajo me he apoyado en dos fuentes bibliográficas, los capítulos 2 y 3 del curso sobre Modelos Digitales del Terreno de A. M. Felicísimo y el manual de ayuda que viene incluido con ArcView. Me han sido de mucha ayuda y recomiendo su consulta para proseguir con el proyecto. He querido llegar hasta el final y generar el modelo raster, aún dejando cosas a medias por el camino (fases del proceso sin terminar y errores sin corregir), para poder ver los resultados de todo el proceso. Los análisis de error y su posterior corrección son una parte importante del trabajo que queda por hacer y que no he abordado. En esta memoria he incluído algunos ejemplos de estos errores y algunas sugerencias sobre cómo creo que podrían corregirse. Otras dos fases que tampoco he tocado son las de validación (estimación de la precisión del modelo) y la de creación de los metadatos. Ambas son muy importantes y creo que deben hacerse. He incluido en el CD los archivos creados a lo largo del trabajo y el mapa que se ha construido con ellos (está en “Arc Map document” y debe ser abierto con ArcView), también lo he incluido como “Publisher Map”, para que pueda ser leído con ArcReader (por algún error en la gestión de los archivos algunos de ellos son inaccesibles en este documento, entre los que lamentablemente está el TIN). La memoria está ilustrada con algunas imágenes tomadas de este mapa. Se notará que entre algunas de ellas, pese a figurar escrita la misma escala, parecen estar a una escala distinta; la razón de esto es que algunas imágenes las he copiado usando la pantalla completa de ArcMap y otras no (al pasar de un tipo de pantalla al oto el tamaño de la imagen cambia). Al final de la memoria cito la bibliografía que he utilizado, parte de la cual procede de Internet y puede ser consultada accediendo a ella mediante el vínculo correspondiente. También he añadido algunas páginas relacionadas con la aplicación de los SIG en la Conservación, que igualmente pueden ser consultadas pinchando en el vínculo. Espero que el esfuerzo haya merecido la pena y se lleve adelante con éxito este ambicioso e importante proyecto. Suerte.

4

Introducción1 La Naturaleza es un tema que parece atraer y suscitar interés a un número cada vez mayor de personas y sin embargo, paralelamente, su conservación se está convirtiendo en un problema cada vez más acuciante y más complicado de abordar, aún contando con la consideración, al menos sobre el papel, de distintas administraciones. Al margen de las causas específicas para cada caso particular, el fondo del problema es que la conservación de la Naturaleza se contrapone en muchos casos a otro tipo de asuntos que suscitan un mayor interés y consideración, o porque se consideran más importantes, o porque rinden más beneficios económicos. Si queremos preservar el patrimonio natural, aunque sea una parte de él, es necesario tomarse el problema de su conservación en serio y otorgarle la importancia y el respeto que merece antes de emprender cualquier acción que pueda dañarlo. La solución pasa por adquirir el compromiso de buscar en todo momento el equilibrio entre conservación y desarrollo. Para garantizar el éxito de cualquier proyecto de conservación o de recuperación que se afronte, ya sea referido a una especie o población animal o vegetal, a un hábitat, a un paisaje, o a un territorio entero, es necesario tener un conocimiento previo profundo y detallado de la realidad del medio natural y humano (entendido medio como el conjunto de todos los elementos y factores bióticos y abióticos que participan y de las interacciones que se dan entre ellos), así como de su fragilidad y de los factores de amenaza presentes y futuros. Este conocimiento, que debería hacerse extensible a todos, implica el manejo de un buen número de datos y de información procedente de muy diversos campos y referente a muy diferentes temas, unas veces obtenido a partir de fuentes externas y otras como fruto de estudios propios. Información y datos que deben ser tratados de un modo integrado para que sea posible tener una visión adecuada de la realidad del problema que se aborda y se puedan plantear políticas eficaces para resolverlo. Frecuentemente esta visión global no puede tenerse debido a lo difícil que resulta el tratamiento y gestión de todos los datos que se recopilan y que muy a menudo quedan finalmente almacenados en archivos tan grandes como inútiles. Lo más habitual es emprender proyectos partiendo de una percepción sólo parcial del problema, dejando de lado asuntos importantes, a veces clave, que están implicados en él (por ejemplo el respaldo de la gente). En auxilio de la causa, los medios informáticos con los que se cuenta hoy en día proporcionan poderosas herramientas de gestión y tratamiento de datos que hacen posible enfocar los problemas de conservación de un modo más holístico y realista y a partir de ahí, diseñar políticas de actuación más eficaces para resolverlos. Dentro de estas herramientas y dado el componente geográfico y espacial intrínseco de la información y de los datos que se manejan, los Sistemas de Información Geográfica (SIG) aparecen como una herramienta idónea. El papel que los SIG pueden desempeñar en la conservación es, para decirlo de un modo sencillo, impresionante. Para tener una visión general de ello recomiendo leer el capítulo “The role of GIS in Conservation”, dentro del documento “ESRI Conservation Program" Program”, y si se quiere tener un conocimiento más "ESRI Conservation profundo y detallado, con abundante información y gran cantidad de ejemplos, ir a “GIS Behavioural Environmental Biology. index html \ web-based Conservation Unit” (http://149.170.202.120/teaching.htm). Ambos están incluidos teching en el CD. Beh.l Env. Biol.\ Web-based teaching Dado que los SIG se han popularizado recientemente y todavía pueden dar lugar a una mala interpretación de su finalidad y aplicación, quiero incluir algunos fragmentos del capítulo 2, “Los modelos como herramientas de gestión de sistemas complejos”, del libro de F. Javier 1

Con esta introducción he querido situar el objeto del proyecto, la construcción de un MDE, dentro del contexto del proyecto de creación de un SIG de la Montaña Oriental y dado que el trabajo que presento en esta memoria es el primer paso que se da para su creación, he creído conveniente hacer una introducción general de todo el proyecto. 5

Moldes “Tecnología de los SIG”, porque expone bastante bien la función principal de los sistemas de información geográfica en el contexto en el que nos movemos: Un modelo es una representación simplificada de la realidad. Los modelos nos permiten una aproximación al conocimiento de la realidad, la cual está formada por sistemas complejos, no siempre comprensibles en toda su complejidad. Los modelos simplifican la realidad y nos presentan aquellos aspectos en los que estamos interesados, mediante objetos e imágenes (por ejemplo, un mapa de Europa nos permite concebir distintos aspectos de su geografía, régimen político, etc.). Estos modelos se denominan modelos instrumentales porque son la base para construir nuestro modelo intelectual de la realidad. Los modelos instrumentales serán todas aquellas representaciones que, teniendo un soporte real, tal como un mapa, una ecuación, un histograma, etc., nos permiten generar modelos conceptuales, es decir, ideas concebidas o adquiridas por la mente que permiten organizar percepciones y conocimientos. Cuando se pretende alcanzar un objetivo, tal como reducir el número de siniestros en la carretera o lograr la protección del águila imperial, planteamos lo que se llama una política, es decir, un conjunto de acciones encaminadas a lograr el objetivo propuesto. Este conjunto de acciones se tiene que basar en un profundo conocimiento de la realidad si queremos alcanzar nuestro objetivo. Este conocimiento de la realidad sólo es posible a través del uso de un conjunto adecuado de modelos instrumentales que nos permita generar un modelo intelectual acertado de la realidad.

En la figura de arriba se recoge un esquema en el que se refleja el proceso cíclico de la gestión de un sistema complejo, En la parte superior está el modelo conceptual, que se apoya para su elaboración en los modelos instrumentales. Abajo está la realidad con toda su complejidad. A la derecha

6

está la política como instrumento de modificación de la realidad a través de un análisis del modelo, la toma de decisiones y por último un conjunto de normas e inversiones. Y a la izquierda están los datos que alimentan los modelos. Nótese que se almacenan datos y no información, ya que ésta se elabora precisamente a través de los modelos instrumentales. Los sistemas de información según este esquema lo constituyen las técnicas e instrumentos que posibilitan la parte izquierda del ciclo. Es difícil discernir entre información y modelo, sobre todo cuando la información está elaborada y organizada en un sistema de información. Un sistema de información es la base para la creación de modelos instrumentales e incluso el propio sistema de información es un modelo instrumental. En un sistema de información caben los modelos matemáticos, ecológicos, estadísticos, etc., constituyendo éstas las partes del modelo instrumental global que es el sistema de información. Una información es el contenido de un mensaje que modifica el estado de conocimiento del receptor; es decir, modifica el modelo conceptual del receptor sobre el objeto real en cuestión. Un conjunto de datos no constituyen un modelo, a no ser que se organicen, se elaboren o se interpreten por parte del receptor. Un conjunto de datos que deja de estar actualizado no permite crear un modelo de la realidad actual, lo será de la realidad del día en que estaba actualizado; así pues, un sistema de información no actualizado deja de constituir un modelo, y por lo tanto deja de ser un instrumento adecuado para la toma de decisiones Los SIG deben ser diseñados como modelos de la realidad que se pretende gestionar. Por lo tanto, deben evitarse aquellos planteamientos en los que el SIG se conciba como sistema de gestión de cartografía y sistema de emisión de gráficos y esquemas, ya que para tal propósito hay herramientas más específicas y baratas, tales como los sistemas de cartografía automática o sistemas multimedia, que nos permitirán obtener estos productos con más calidad y menor coste. La implantación de un SIG supone una nueva metodología de trabajo que invalida una parte de los métodos tradicionales de la toma de decisiones, creando una nueva doctrina metodológica que en pocos casos se ha aplicado y de cuya eficacia no tengo la menor duda. Dentro de la idea de modelo instrumental podemos hacer una distinción, atendiendo a su finalidad, en dos tipos principales (según se explica en el manual de ayuda de ArcView): aquellos que representan o describen los objetos en el paisaje (modelos de representación) y aquellos que tratan de simular procesos en el paisaje (modelos de procesos). ƒ

Los modelos de representación, o descriptivos, tratan de describir los objetos en el paisaje, tales como edificios, arroyos o bosques. El modo en que estos modelos son creados en un SIG es a través de un conjunto de capas temáticas. Estas capas de datos pueden definir localizaciones (raster), o entidades u objetos (vectoriales). Las capas raster se representan por una malla rectangular (grid) y cada localización en cada capa es

7

representada por una celda que tiene un valor. Las celdas de varias capas se apilan unas obre otras, describiendo muchos atributos de cada localización. Las capas vectoriales se representan por puntos, líneas o polígonos que describen la forma y definen la posición de los objetos. Los modelos de representación tratan de capturar las relaciones espaciales dentro de un objeto (la forma del edificio, por ejemplo) y entre los otros objetos en el paisaje (la distribución de edificios). Además el modelo de representación también modela los atributos de los objetos (a quién pertenece el edificio). ƒ

Los modelos de procesos tratan de describir la interacción de los objetos que son modelados en el modelo descriptivo. Las relaciones son modeladas usando herramientas de análisis espacial. Los modelos de procesos pueden ser usados para describir procesos y también parara predecir acontecimientos en función de determinados factores.

La creación de estos modelos referidos a la Montaña Oriental tendría gran cantidad de aplicaciones; entre otras cosas nos ayudarían a: -

Identificar, localizar y valorar2 los distintos hábitats (o biotopos) y comunidades vegetales. Dónde se localizan. Qué hábitat o vegetación hay en determinado sitio. Cómo es el relieve o el tipo de suelo donde se encuentran. Qué extensión tienen y cómo de fragmentados están. Qué valor natural o cultural tienen. En qué estado de conservación se encuentran actualmente. Qué tipo de amenaza pesa sobre ellos.

-

Estudiar la biodiversidad y su distribución dentro del territorio. Qué especies hay en un determinado sitio (qué especies de pícidos podemos encontrar en determinado bosque). Dónde están localizadas (dónde podemos encontrar nidos de alimoche). Qué características definen el sitio donde están.

-

Analizar las interacciones que se producen entre las poblaciones animales y vegetales, el medio físico y el medio humano. Qué selección de hábitat hace dentro de su territorio una determinada población animal (qué vegetación está asociada con el área de campeo de la población x de rebecos; dentro de qué rango de altitud vive). Qué tipo de relaciones intra o extraespecíficas se dan entre determinadas especies dentro de las comunidades donde viven (¿hay ganado dentro de las áreas de campeo de la población de rebecos? Dónde se han producido ataques de lobos).

-

Predecir la distribución potencial de determinadas especies o poblaciones animales y vegetales. De la información obtenida en los puntos anteriores se puede estimar el patrón de distribución que caracteriza a una determinada especie o población y cuáles son sus requerimientos mesológicos. Localizando en qué zonas del territorio se dan éstos, se pueden por ejemplo diseñar áreas de reintroducción o corredores biológicos para conectar poblaciones aisladas.

-

Estudiar la evolución espacio-temporal de las poblaciones animales y vegetales y de su medio. Cuánto ha retrocedido un determinado tipo de comunidad vegetal en los últimos años.

2

En la figura figura 11, tomada de “Metodología y práctica de la Biogeografía”, capítulo 4: “Valoración de la vegetación”. En este capítulo se expone una metodología para la diagnosis del estado actual de la vegetación y para su evaluación como recurso y patrimonio natural con fines conservacionistas. El cuadro representa los distintos parámetros que se usan para la diagnosis y evaluación. Éstos descansan en dos conceptos valorativos diferenciados: El Interés de conservación (interés natural e interés cultural) y la Prioridad de Conservación (basada en la amenaza). 8

Cómo ha variado la distribución de los nidos de alimoche y qué ha cambiado en el paisaje que pueda explicar esta variación. -

Simular alteraciones en el medio y analizar sus posibles consecuencias. Qué efecto macroecológico tendría la erradicación de los incendios provocados en los prados montanos y la consecuente recuperación de la vegetación natural. Qué efecto barrera o efecto borde produciría sobre un ecosistema una determinada infraestructura humana.

-

Presentar la información de un modo claro, detallado y atractivo. Un proyecto de conservación que se vaya a llevar a cabo puede ser presentado al público y al organismo del que se busque apoyo, haciendo uso de mapas temáticos interactivos a distintas escalas, gráficos, fotografías, etc.; mostrando cuál es la finalidad del proyecto, el área donde se pretende realizar, la situación actual del problema, los factores que intervienen en él, los objetivos que se pretenden conseguir y los resultados, cuando se obtengan. Una presentación así ayudaría a que nuestros argumentos fueran mejor entendidos y posiblemente compartidos, obteniendo con ello el respaldo por parte de todos y garantizando el éxito del proyecto.

La primera etapa en la creación del modelo, una vez que se ha diseñado a grandes rasgos, sería la obtención de las distintas capas temáticas que formarán parte de él: relieve, vegetación, suelos, hidrología, red viaria, propiedad y usos del suelo, son algunas de ellas. El relieve (o la topografía) del terreno es una capa básica dentro del modelo, ya que tiene una gran influencia en el resto. El relieve, a través de sus componentes de altitud, pendiente y exposición (orientación), determina de manera importante las condiciones del medio físico, influyendo directa o indirectamente en la distribución y en la actividad no sólo de la vegetación y de la fauna, sino también humana. La altitud y la exposición determinan en gran medida el ambiente climático temperatura, humedad, insolación, presencia de nieve, régimen de viento-, que para la mayoría de las especies marca el límite de distribución: la variación altitudinal de la temperatura y de la pluviometría, a la que se suma el efecto producido por la orientación, determina la aparición de los pisos de vegetación y la distribución cliserial de gran parte de la fauna (un ejemplo de esto puede verse en la Figura 2). La pendiente del terreno condiciona directamente la aparición de determinadas especies vegetales y animales, como es el caso del buitre común y del alimoche, que localizan sus nidos en paredes verticales de farallones, o del rebeco, que establece sus áreas de campeo donde encuentra fuertes pendientes que usa como zonas de escape. Indirectamente, al estar ligada a procesos de erosión, deslizamientos, retención de agua, propagación de incendios, etc., la pendiente también ejerce influencia sobre la vegetación y la fauna. El relieve puede originar la aparición de barreras geográficas que afectan a la distribución de algunas especies animales, como en el caso de las sierras, que limitan los desplazamientos de la nutria; o la aparición de corredores utilizados por muchos animales en sus desplazamientos, como los collados, usados por los rebecos en sus migraciones. Por último el relieve determina también las actividades humanas de aprovechamiento agrario y ganadero, el uso lúdico y deportivo, la implantación de infraestructuras, etc. La representación del relieve en un SIG se obtiene mediante el Modelo Digital de Elevaciones (MDE), que puede definirse como “una estructura numérica de datos que representan la distribución espacial de la altitud de la superficie del terreno” (en Felicísimo: “Modelos Terreno,tema tema2"22” se ofrece una descripción clara y detallada sobre lo que "ModelosDigitales Digitales del del Terreno, tema son, los tipos que hay y cómo se construyen los MDE). Un MDE permite obtener información sobre la altitud, la pendiente y la orientación del terreno, además de ofrecer una visión bastante

9

realista de su relieve. Los dos tipos de modelos más importantes son el TIN (vectorial) y el de Matrices Regulares (raster). ƒ

El modelo TIN está basado en entidades geométricas, consistentes en triángulos irregulares adosados, definidas por las coordenadas de sus nodos o vértices. El terreno queda representado por el conjunto de superficies planas (triángulos) que se ajustan a un conjunto previo de puntos con atributos de altitud.

ƒ

El modelo de Matrices Regulares está basado en localizaciones espaciales (obtenidas al superponer una retícula sobre el área representada) definidas por su posición en la matriz, una vez se han definido el origen y el valor del intervalo entre sus filas y columnas. A cada localización se le asigna el valor medio de altitud para la unidad elemental de superficie, la celda.

En función del modo en que estructura sus datos, cada modelo tiene unas cualidades específicas (viene muy bien explicado en el manual de ayuda de ArcView) El TIN tiene la ventaja de que: ƒ

Representa el relieve con más precisión y detalle que el raster, ya que los vértices que forman su estructura conservan intactos los valores de altitud reales y se puede adaptar a la complejidad del terreno incorporando más vértices, así como elementos estructurales que definen su forma (crestas, fondos de valle, etc.)

ƒ

Permite representar entidades lineales y puntuales (ríos, cimas, etc.)

ƒ

Permite medir distancias lineales (euclídeas) y superficies planas de forma más precisa que en el raster

ƒ

Permite identificar y tratar cada elemento o entidad individualmente.

El raster tiene la ventaja de que: ƒ

Proporciona unas condiciones mucho mejores para la creación y uso de modelos, especialmente cuando están involucradas variables continuas

ƒ

Permite trabajar con distintos mapas temáticos (superposición de capas) de un modo más flexible y rápido

ƒ

Permite realizar potentes operaciones de análisis espacial a partir de las cuales se deriva nueva información Permite la integración de información procedente de sensores remotos (fotografías aéreas e imágenes de satélites)

ƒ

Para nuestro propósito de usar el modelo como instrumento de trabajo, no cabe duda de que el modelo de datos raster es el idóneo. La capacidad del TIN para representar más fielmente la topografía del terreno lo hace útil para capturar la información y generar a partir de él el raster y también para tenerlo de referencia.

10

Objetivos El objetivo propuesto es la obtención de un Modelo Digital de Elevaciones (MDE) TIN y otro Raster del área LIC de la Montaña Oriental de Cantabria y su zona de influencia a escala 1:25.000, que represente lo más fielmente posible la topografía del territorio. Su finalidad es incorporarlos posteriormente en modelos descriptivos (“representation models”) y en modelos de procesos (“process models”) como herramienta de trabajo para proyectos de conservación en la Montaña Oriental. La consecución del objetivo supone obtener un modelo: 9 Georreferenciado correctamente 9 Preciso planimétrica y altimétricamente 9 Validado El grado de precisión planimétrica, altimétrica y de georreferenciación, requeridas para los MDE, teniendo en cuenta sus futuras aplicaciones, no se ha determinado exactamente. No se han establecido por tanto los límites de error que pueden ser aceptados para considerar válidos los modelos. Se tiene propuesto no establecer límites muy estrechos de error a priori, en función de consideraciones externas, aunque sí procurar que éstos sean lo más bajos posible. El planteamiento propuesto es localizar las fuentes de error, cuantificar su valor, estudiar el efecto que tienen en el resultado final y, en la medida en que sea posible, cuando el error sea importante tratar de reducirlo. La validación posterior del modelo y de todos los pasos intermedios permitirá estimar su bondad de ajuste, que será tenido en cuenta en las aplicaciones posteriores.

Materiales y métodos Cartografía: •

Hojas en papel del Mapa Topográfico Nacional del I.G.N. a escala 1:25.000

o o o o o o o o o

Nombre Sarón Arredondo Selaya Veguilla Ramales de la Victoria Concha Vega de Pas Bárcenas Bercedo

Nº de hoja 59-I 59-II 59-III 59-IV 60-I 60-III 84-I 84-II 85-I

Año de edición (1994) (1994) (1994) (1994) (1997) (1997) (1997) (1997) (1997)

La mayoría de las hojas habían sido adquiridas con anterioridad al inicio del trabajo pero en apariencia presentaban buenas condiciones. Tan solo se compraron tres hojas nuevas expresamente para usarlas en este proyecto.

11

•

Archivos cartográficos vectoriales procedentes del CD realizado por CETYMA SL con el título “Red Natura 2000 Cantabria”, que incluye distintas capas de información geográfica de toda la Comunidad de Cantabria3.

•

Hardware: •

Escaner Hewlett Packard Scan Jet 3400 C

Software: •

Sistema operativo Microsoft Windows XP Home Edition versión 2002

•

ArcView 8.2, con la extensión ArcGIS 3D Analyst

•

ArcPublisher, ArcReader, Microsoft Word y Adobe Acrobat (para presentar los resultados y redactar la memoria).

Procedimiento: El procedimiento ha consistido en la digitalización manual en pantalla de las curvas de nivel, los puntos de cota y la red fluvial, utilizando como molde imágenes raster obtenidas por escaneado de mapas topográficos del IGN en papel, posteriormente georreferenciadas. Los resultados de la digitalización se han usado como datos de entrada para la construcción del MDE vectorial con estructura TIN. A partir de éste se ha generado un modelo raster con estructura de Matrices Regulares. A continuación se detallan los pasos seguidos: 1. Escaneado de la cartografía original en papel Las hojas del IGN han sido escaneadas con un escáner de formato A4 (Heulett Packard Scan Jet 3400 C). El tamaño reducido del escaner ha obligado a dividir cada hoja en cuatro partes, dejando zonas de solapamiento entre ellas. La resolución vertical y horizontal para el escaneado no se ha prefijado para la mayoría de las hojas (por descuido), tan solo hay cuatro imágenes raster para las que se ha ajustado a 600 ppp. En cualquier caso se ha considerado que la resolución era aceptable para realizar el trabajo. Las imágenes obtenidas no han sido manipuladas y se han guardado directamente como archivos JPEG para ser usadas con ArcView. 2. Alineación y georreferenciación de las imágenes raster Previamente se habían añadido a ArcMap ficheros vectoriales (shapefile) con información sobre los límites municipales y provinciales, las áreas LIC (lugares de importancia comunitaria), los núcleos de población y la red de carreteras de toda la comunidad de Cantabria 3

Este trabajo ha sido publicado por el Gobierno de Cantabria. Consejería de Ganadería, Agricultura y Pesca. Dirección General de Montes y Conservación de la Naturaleza. Y fue presentado en la última conferencia nacional de ESRI en Madrid. 12

Estos ficheros están georreferenciados con respecto al mismo sistema de proyección que los mapas del IGN de escala 1:25.000 que se iban a utilizar (Datum Europeo 1950 UTM zona 30 T). La finalidad de incorporarlos era la de partir de un mapa digital con buena calidad para crear con él una red de referencia de coordenadas en el SIG. A parte de esto, la información que incluyen es útil para la base de datos del SIG. Estos ficheros se obtuvieron de un CD realizado por la empresa Cetyma, conteniendo información sobre la red Natura 2000 en Cantabria. Para realizar la alineación y georreferenciación de las imágenes raster se han construido dos capas vectoriales, rejilla y malla; la primera contiene los puntos de intersección de los ejes de coordenadas XY que delimitan las cuadrículas UTM de 1 Km de lado en nuestra proyección; la segunda, las líneas que definen estas coordenadas. La edición de esta rejilla se ha hecho incorporando, con ayuda del ratón, uno a uno los puntos en la pantalla, ajustando su posición con la mayor precisión posible, guiándose por los valores X e Y que aparecen en el monitor. Se ha trabajando a una escala de 1:2.500, porque a esta escala la resolución espacial en la pantalla está por debajo de 70 cm (del modelo real), lo que permite obtener errores menores de 1 m en cada punto, más que aceptables para nuestros propósitos, (al aumentar la escala de trabajo lógicamente aumenta la resolución espacial en la pantalla y con ello la precisión espacial de los vértices introducidos). Los puntos de intersección de coordenadas de la rejilla se han usado más tarde como puntos control de referencia (target) para la alineación de las imágenes raster (usando los correspondientes puntos de intersección de coordenadas del mapa original como puntos control homólogos). Las líneas de la malla se han usado para comprobar visualmente en las cuadrículas el error de georreferenciación cometido. La ventaja de usar los vértices de las coordenadas UTM como puntos control es que resulta muy sencillo localizarlos y por otro lado que la alineación se realiza homogéneamente en todo el raster. El único problema se presenta cuando en los márgenes de la imagen hay un espacio grande más allá de las coordenadas, porque no hay otros puntos comunes entre la imagen raster y el mapa de referencia para usar como puntos control. Cada imagen escaneada ha sido añadida a ArcMap y se ha georreferenciado usando su herramienta de georreferenciación. La alineación de la imagen se ha realizado enlazando primeramente los puntos control en las cuatro esquinas, después se han ido creando enlaces (“links”) en puntos interiores críticos separados entre sí. El número de puntos control enlazados ha estado en función de que se estimara que la imagen presentaba ya una alineación suficiente con la rejilla, o que la introducción de más enlaces sólo contribuyera a aumentar el error residual en los otros. El tipo de transformación (orden del polinomio) usado para la alineación del raster no ha sido especificado previamente, de manera que el programa (creo) ha recurrido por defecto a una transformación lineal o de primer orden, que implica una traslación, un giro y un cambio de escala del raster. No se ha recurrido expresamente a una transformación de orden superior, que produce una deformación geométrica de la imagen, pese a que quizá habría sido aconsejable para algunas de ellas. El error cometido en la alineación de cada punto control enlazado -error residual- y el error total -error cuadrático medio- (quedan registrados automáticamente en una tabla) no han sido analizados estadísticamente, por lo que no se ha podido cuantificar la precisión del resultado. Se ha estimado, en cualquier caso que, salvo en algunas imágenes, el error cometido es despreciable cuando se trabaja a escala 1:25.000 ó 1:10.000. En los resultados se comenta este punto.

13

3. Digitalización manual La digitalización de las curvas de nivel y de los otros elementos se ha realizado con las herramientas de edición de ArcMap. Primeramente se había empezado trabajando a escala 1:10.000, pero después se ha aumentado a 1:2.500 (escala 1:10.000 en el mapa y uso de la lupa con un aumento del 400%) para tener una mayor precisión en el trazado. Se ha procurado seguir bien el contorno de las líneas, usando la herramienta de arcos (arc tool) para las líneas curvas e introduciendo el número de vértices necesario para que la polilínea no rebasara los pixels de la línea de referencia del raster. Se ha procurando seguir el criterio de que la distancia entre vértices en una línea sea parecida a la distancia entre líneas contiguas. En cualquier caso, el número final de vértices en cada línea lo establece el propio programa en función del que necesite para ajustar los tramos curvos. Los valores de altitud de cada elemento se han registrado en la correspondiente tabla de atributos. Primeramente se han digitalizado las curvas de nivel maestras de todo el área de estudio (su equidistancia es de 50 m). Las curvas secundarias (equidistancia 10 m) se han digitalizado parcialmente en determinadas zonas y en algunos casos sólo fragmentos de ellas. Se ha escogido un área de ensayo (comprendida en las imágenes “Concha 1” y “Veguilla 2”) donde se ha realizado una digitalización exhaustiva de prácticamente todas las curvas secundarias, con la idea de ver el efecto que tiene incorporar una alta densidad de datos altimétricos a la hora de generar el MDE y comparar el resultado con el que se obtiene en zonas donde la densidad es menor. En otras zonas del área de estudio, hasta que el tiempo lo ha permitido, se han digitalizado fragmentos de curvas secundarias en zonas donde la distancia entre curvas maestras era grande y en zonas con relieve crítico (collados, fondos de valle, líneas de cresta, etc.). En ocasiones ha sido necesario recurrir al mapa de papel, lupa en mano, para distinguir líneas que estaban muy próximas entre sí y que debido a la pequeña resolución de la imagen raster no podían ser discriminadas en la pantalla. En estos casos se ha procurado seguir su contorno cometiendo el mínimo error. A cada curva de nivel se le ha asignado su correspondiente valor de altura en la tabla de atributos. Para diferenciar que una curva representa una depresión en el terreno, como una hondonada o una dolina, se le ha dado un valor de 1 en el correspondiente campo de la tabla; las curvas normales tienen un valor 0 en ese campo. Los resultados están guardados en el archivo “Curvas de nivel”. De los puntos de cota sólo se han digitalizado unos pocos, que al igual que las curvas de nivel, tienen registrado su valor de altitud en la tabla de atributos. Tienen además asignado un código que permite distinguir si el punto corresponde a una cima (1), a un hoyo o dolina (0), o a un collado (2). Están guardados en el archivo “Mass points” (debería llamare “puntos de cota”, lo siento). Los elementos fluviales se han digitalizado como polilíneas sin valores de altura (por error); están almacenados en el archivo “Hidrografía”. Tienen registrado el atributo del nombre del río o arrollo si es conocido. Se ha digitalizado tan solo una polilínea, representando un arrollo, a cuyos vértices se les ha incorporado información altimétrica como coordenada Z. Se encuentra en el archivo “Ríos”.

14

El inconveniente que tiene incorporar datos altimétricos cuando se digitaliza una línea como un río o una carretera, cuya altura varía continuamente, es que es necesario asignar un valor de altitud a todos y cada uno de sus vértices (ya que ésta no se puede registrar como un atributo, como sí ocurre con las curvas de nivel). El problema aparece cuando hay vértices importantes que incorporar a la línea (un recodo en un río, por ejemplo) de los que no se conoce su altura, porque ninguna curva de nivel o punto de cota se cruza con ellos. Cuando esto ha ocurrido en nuestro trabajo, o bien se ha optado por omitir el vértice, o bien se le ha asignado una altura estimada. En el primer caso se comete un error de generalización, que da lugar a un error posicional (localización imprecisa de la trayectoria de la línea) y en el segundo se comete un error atributivo (asignación imprecisa de la altitud). La decisión de elegir una u otra opción se ha tomado ponderando intuitivamente los dos errores posibles y optando por la opción menos mala. En cualquier caso se estima que el error cometido es despreciable. El polígono que delimita el área de estudio se ha fijado de modo que se incluyera dentro de él todo el área LIC de la Montaña Oriental y una porción del territorio circundante hasta llegar más o menos a la divisoria de aguas próxima. Este elemento no tiene ningún tipo de atributo y carece de datos altimétricos. La falta de tiempo no ha permitido (ni mucho menos) concluir el trabajo de digitalización. El proceso de digitalización manual es muy laborioso y ha ocupado la mayor parte del tiempo (unos cuatro meses de trabajo continuado). 4. Creación del MDE El TIN Los elementos usados para crear el TIN han sido las “Curvas de nivel” (maestras y secundarias), los puntos de cota (“Mass points”), que definen cumbres, fondos de dolinas y collados, la “Hidrografía” y los “Ríos”. Los límites del modelo se han fijado con el polígono del “Área de estudio”. Estos elementos se han incorporado como cuatro tipos de datos de entrada: Puntos acotados (Mass Points): los puntos de cota del archivo “Mass points” y los vértices de las “Curvas de nivel”; estos últimos sin haber pasado por un proceso de generalización para reducir su número. Líneas estructurales (Hard Breaklines): el único elemento del archivo “Ríos”. Líneas de inflexión o rotura (Soft Breaklines): la red fluvial del archivo “Hidrografía”. Zonas de recorte (Clip Polygons): el polígono que delimita el “Área de estudio”. El programa procesa automáticamente estos datos para crear la red de triángulos adosados TIN. De los cuatro tipos, sólo los dos primeros tienen información altimétrica y participan activamente en la construcción del TIN. Los puntos acotados (Mass points) y los vértices de las Curvas de nivel proporcionan la información altimétrica puntual y se convierten en los nodos o

15

vértices de la red de triángulos. Las líneas estructurales (Ríos) representan discontinuidades en la pendiente del terreno y contribuyen a dar forma a la superficie, definiendo bordes de triángulos. Las líneas de inflexión o rotura (la Hidrografía) sólo tienen información planimétrica, su efecto en el TIN se limita a superponerse a la superficie creada sin modificar su forma. Las zonas de recorte (el Área de estudio) delimitan la frontera más allá de la cual cualquier dato es ignorado en la interpolación. El Raster El modelo raster (de Matrices Regulares) se ha creado a partir del TIN, seleccionado el atributo elevaciones y ajustando el tamaño de la celda a 25 m (el tamaño que proponía la aplicación por defecto era mayor). El programa automáticamente procesa los datos y genera el modelo. Se ha elegido generar el modelo raster a partir del TIN y no directamente a partir de los datos digitalizados (vértices de curvas de nivel, puntos de cota, etc.), porque el proceso de interpolación que se emplea en el primer caso es específico para esta estructura de datos y permite en teoría aprovechar parcialmente la capacidad del TIN para integrar discontinuidades en la generación del modelo (los métodos de interpolación usados para tratar datos puntuales no discriminan las líneas estructurales ni los puntos singulares, como las cimas)4.

Resultados y discusión Puesto que el trabajo ha quedado interrumpido sin haberse podido terminar, los resultados que se muestran deben considerarse provisionales. 1. Obtención y tratamiento de las imágenes raster El resultado final es un mosaico de veintiuna imágenes raster georreferenciadas con respecto al sistema de coordenadas definido para nuestro trabajo, que cubre todo el área de estudio. Puede verse en la Figura 3. Las imágenes, tal y como han sido obtenidas en el escaneado y se han usado como molde para la posterior digitalización, pueden verse en el fichero imágenes raster incluído en el CD. La posición que ocupa cada imagen en el mosaico aparece en la figura de abajo.

4

No puedo afirmar con certeza que ArcView realice este tipo de interpolación específica para datos TIN; sólo lo supongo. 16

Arredondo 4

Arredondo 3

Ramales 3

Sarón 1

Arredondo 1

Arredondo 2

Ramales 1

Ramales 2

Selaya 1

Veguilla 4

Veguilla 2

Concha 1

Concha 2

Selaya 3

Veguilla 1

Veguilla 3

Concha 3

Vega de Pas 2

Bárcenas 1

Bárcenas 2

Bercedo 1

A continuación se comentan los resultados obtenidos en esta etapa del proceso. La resolución espacial de las imágenes escaneadas Se pretendía que la calidad de las imágenes escaneadas permitiera seguir bien el trazado de las líneas, discriminando entre ellas cuando están muy juntas, a una escala de trabajo de 1:2.500. La resolución espacial de cada imagen, medida en puntos por pulgada, puede verse consultando sus propiedades en la carpeta Imágenes raster. Como ejemplo la figura 4 muestra un fragmento de una imagen raster obtenida ajustando la resolución del escaneado a 600 ppp. La figura siguiente muestra un fragmento de una imagen raster cuya resolución no fue ajustada previamente y quedó determinada por el propio programa en aproximadamente 100 ppp (la línea que se ve a la derecha de la coordenada es un doblez de la hoja original). Si se hubiera querido hacer una digitalización automática, esta última imagen habría dado problemas a la hora de eliminar selectivamente los elementos no deseados, porque es muy difícil discriminarlos por el color del pixel. Para una digitalización manual esta resolución permite trabajar bien salvo en las zonas donde las líneas se aproximan mucho. Según el grado de precisión que se quiera obtener esto puede ser un problema o no. La georreferenciación Cada imagen raster georreferenciada está almacenada en un archivo con la extensión .aux, que contiene la información de sus coordenadas en el mapa. Este archivo puede modificarse permanentemente usando el comando de rectificación (Rectify) de la herramienta de georreferenciación de ArcMap; al hacerlo se le asignan nuevas coordenadas.

17

La georreferenciación de las imágenes raster ha sido el paso más crítico en todo el proceso de construcción del MDE, debido a lo difícil que es conseguir un resultado preciso y a la trascendencia que tiene para la fiabilidad del modelo final cometer un error importante. La calidad de la georreferenciación de una imagen raster en nuestro trabajo depende de dos factores: la precisión espacial con la que se localizan los puntos control en el mapa de referencia (la rejilla) y la precisión con la que éstos se alinean con los puntos control homólogos del raster una vez realizada la transformación en éste. Del grado de fiabilidad de la georreferenciación depende el grado de fiabilidad del modelo digital de elevaciones que se obtenga al final, ya que el error cometido en este punto se transmite a lo largo del proceso de digitalización posterior como un error posicional. La precisión de los puntos que forman la rejilla está en función de la resolución espacial que presenta la pantalla a la escala en la que se ha construído (distancia mínima que se puede registrar entre dos coordenadas X y entre dos coordenadas Y). A la escala 1:2.500 (nuestro caso) la resolución ha permitido localizar los puntos con un error de posición por debajo de los 70 cm., tanto en la coordenada X como en la coordenada Y (algunos puntos tienen un error de mm). Hay que decir que, curiosamente, la resolución no es homogénea a lo largo de la pantalla, no sé porqué. En cualquier caso se ha considerado que la precisión con que está hecha la rejilla es más que aceptable para un modelo a escala territorial como el que pretendemos crear (1:25.000 – 1:10.000). En las imágenes que se muestran como ejemplo los puntos de la rejilla tienen el símbolo +. La alineación de las imágenes raster varía de un caso a otro, dependiendo (es de suponer) del grado de deformación que presentaba el mapa original de papel a partir del cual se ha obtenido cada una. Para algunas imágenes no se ha necesitado más que un mínimo de puntos control (siempre igual o superior a tres) para conseguir un resultado satisfactorio en toda ella (como en el ejemplo que se muestra en la figura 5). En otros casos, pese a haberse utilizado un buen número de puntos, la alineación del raster no puede considerarse aceptable, ya que muchos puntos presentan errores residuales aparentemente importantes, que además afectan al solapamiento de imágenes contiguas. El error cometido en la transformación del raster se mide comparando las coordenadas de los puntos control en la rejilla (target), con las coordenadas de los puntos homólogos en el raster, resultantes de aplicar la transformación. El error cometido en cada punto es el error residual; el error global, aplicable a la totalidad de los puntos utilizados en la transformación, se mide como error cuadrático medio (raíz cuadrada de la media de los cuadrados de los errores residuales). Estos dos errores se expresan separadamente para la coordenada X y para la coordenada Y. La herramienta de georreferenciación de ArcMap registra automáticamente en una tabla los valores de las coordenadas de los puntos control de la rejilla y de los puntos transformados del raster a medida que éstos se van enlazando, compara estos dos valores y proporciona el valor del error residual; también realiza los cálculos para dar el ECM. Teniendo en cuenta que los puntos de la rejilla usados para alinear y transformar el raster tienen un pequeño error de posición (comentado arriba), el error de georreferenciación real (que es el que nos interesa) debe medirse en relación al sistema de coordenadas reales de nuestra proyección UTM y no al de la rejilla (usado sólo como referencia aproximativa):

18

ex = Xtransformado –Xreal ; ey = Ytransformado –Yreal La Figura 6 muestra dos ejemplos de error de alineación (error residual): El primero de ellos corresponde a un fragmento de la imagen “Selaya 3”. El vértice de coordenadas UTM (en metros) que aparece es el (XY) 441.000 4.784.000. En la rejilla creada este vértice (+) tiene los valores 440.999,81 4.784.000,47 (error de posición: X= -0,19m; Y= 0,47m). El valor que toma este punto en la imagen raster después de aplicar la transformación es aproximadamente 441.001,91 4.784.008,88 (error residual X≈2,1m; Y≈8,41m). El error de georreferenciación real es aproximadamente X≈1,91 m; Y≈8,88 m. El segundo ejemplo es un fragmento de la imagen “Concha 1”. El vértice de coordenadas UTM es el 462.000 4.788.000. En la rejilla tiene los valores 461.999,74 4.787.999,96 (error de posición X=-0,26m; Y=-0,04m). Después de aplicar la transformación el valor de ese punto en el raster es aproximadamente 462.005,44 4.788.003,27 (error residual X≈5,7m; Y≈3,31m). El error de georreferenciación real X≈5,44m; Y≈3,27m. Un error de georreferenciación en una imagen raster que va a ser usada como molde para la digitalización de distintos elementos cartográficos (curvas de nivel, puntos de cota, ríos, etc.) se propaga como un error posicional (planimétrico) en la capa vectorial resultante y en todas las capas que se creen a partir de ella. El efecto de la inexacta localización de elementos con datos de cota se traduce en un error de altitud en el modelo, que es proporcional a la tangente de la pendiente del terreno (localizar una cota en el modelo con un error planimétrico de por ejemplo 10 m en un terreno llano no tiene ningún efecto, pero localizarla en el pico de una montaña supone desplazar esta cota 10 m y situarla en un punto donde en el peor de los casos ya no hay montaña). Este tipo de errores puede dar problemas si posteriormente el modelo se va a ver implicado en un proceso de análisis espacial o va a formar parte de un mapa temático, casos en los que se van a superponer otras capas de distinta procedencia (mapa de vegetación, mapa de suelos, etc.). Corregir los errores graves de georreferenciación es fundamental. Para hacerlo, lo primero sería examinar todos los puntos control (+) e identificar aquellos que tienen un error de alineación importante; marcarlos con un color distinguible (rojo, por ejemplo) sería muy práctico para tenerlos localizados. Si toda la imagen raster estuviera mal georreferenciada (cosa que creo, no ocurre en ningún caso), habría que rectificarla de nuevo. Si sólo hay algunos puntos control afectados por errores graves de alineación (como es el caso en alguna imágenes), se podrían hacer dos cosas: ● como la finalidad de las imágenes raster es sólo servir de molde para la digitalización y luego se pueden eliminar, no es necesario que cada una de ellas esté bien georreferenciada en su conjunto (aunque es preferible); basta con que lo esté el fragmento que se está digitalizando en ese momento. De este modo se podrían respetar los elementos que se han digitalizado sobre fragmentos de la imagen comprendidos entre puntos control bien alineados y digitalizar los que falten, luego rectificar los puntos control afectados por un

19

error de alineación (usando el comando de rectificación) y posteriormente digitalizar los elementos del fragmento de la imagen comprendido entre ellos. ● rectificar toda la imagen aplicando una transformación de orden superior. Con esto quizá se ganaría en precisión general de la imagen (reducción del ECM), pero podría perderse precisión en los puntos que originalmente estaban bien alineados. El solapamiento de imágenes Cuando dos imágenes raster contiguas están más o menos bien georreferenciadas, al menos en los bordes que comparten, solapan correctamente entre sí. Esto ha ocurrido en un buen número de casos como los que se muestran en la Figura 7 y en la Figura 8; en la primera el solapamiento es entre dos imágenes obtenidas a partir de hojas distintas del mapa topográfico y en la segunda el solapamiento es entre imágenes procedentes de la misma hoja. En otras ocasiones, los errores de georreferenciación de al menos una de las dos imágenes vecinas en el mosaico se evidencian, a veces de un modo dramático, en la zona de solapamiento. Un ejemplo de solapamiento defectuoso (y no el más grave) se muestra en la figura 9. El problema del mal solapamiento de imágenes en teoría podría solucionarse de una de estas maneras: ● rectificando el error de alineación en el borde conflictivo (con el comando de rectificación) usando los vértices de coordenadas si están cerca del borde, y si no, buscando otros puntos control que sean comunes en las dos imágenes (un punto de cota, un símbolo del mapa) y seleccionando como los de referencia (target) los que correspondan a la imagen que esté mejor georreferenciada. ● rectificando la imagen entera aplicando una transformación de orden superior ● si las dos imágenes comparten una zona de solapamiento más o menos ancha y el problema se presenta sólo en el borde visible, de manera que más allá las dos se superponen correctamente (como ocurre en algunos casos), prescindir de la imagen problemática (la superior) y utilizar la imagen inferior como molde para la digitalización hasta donde el problema de solapamiento desaparezca. 2. Digitalización manual Para construir el modelo TIN lo más fielmente ajustado al modelo real posible, es necesaria la participación de distintas estructuras vectoriales que ayuden a reproducir con el suficiente detalle y precisión los elementos topográficos que definen el relieve. Estas estructuras son: curvas de nivel maestras; curvas de nivel secundarias, en zonas donde la distancia entre las curvas maestras sea demasiado grande5, o donde desempeñan un papel importante para definir el relieve (cimas, fondos de valle, collados, etc.); puntos de cota; red fluvial (ríos, arroyos, lagunas); red viaria (carreteras, caminos, sendas, etc.); líneas no representadas explícitamente en la cartografía, pero que definen implícitamente la forma del relieve, como líneas de cresta, líneas de fondo de barranco y de valle, aristas de ladearas, etc. 5

Este punto es especialmente importante cuando se quiere generar un modelo raster. 20

De todos estos elementos vectoriales mencionados, hasta ahora tan sólo se han digitalizado los siguientes: la totalidad de las curvas de nivel maestras, una parte de las curvas de nivel secundarias (están en el archivo “Curvas de nivel”), un reducido número de puntos de cota (en el archivo “Mass points”) y casi toda la red fluvial (en los archivos “Hidrografía” y “Ríos”). El resultado de la digitalización en toda el área de estudio se muestra en la figura 10. A continuación quiero comentar algunos puntos importantes. Para crear la estructura del TIN (construir los triángulos) sólo son aprovechables los elementos que contienen datos de altura: las curvas de nivel, los puntos de cota y un arroyo (el Barranco del Rosario, en el archivo “Ríos”). El resto de la red hidrográfica (en “Hidrografía”) no cuenta con información altimétrica y tiene que ser digitalizada de nuevo, asignando a cada vértice de las polilíneas su correspondiente altura como valor de z (cuando se creó este shapefile, desde ArcCatalog, no se indicó que las coordenadas contendrían valores z y por tanto no se pueden incluir hora –un fallo lamentable). Las polilíneas creadas de nuevo deben incluirse en el archivo “Ríos”, porque en éste sí se especificó que las coordenadas tendrían valores z. Queda pendiente también digitalizar los demás elementos estructurales. Al igual que en el caso de la red fluvial, se deben incorporar los valores de Z en cada vértice de la polilínea, tomándolos de los datos de altura que aportan las curvas de nivel o los puntos de cota que se cruzan con ellos. Los tramos entre dos datos de altura conocidos tendrán que ser rectos y si esto supone cometer un error de generalización importante (línea recta donde debería ser curva), habrá que optar por introducir un vértice con una altura aproximada. Los elementos digitalizados pueden estar sujetos a tres tipos de error: errores derivados de la inexactitud en la georreferenciación de las imágenes raster, errores cometidos por la imprecisión en el trazado de las líneas al usar una escala de trabajo pequeña (1:10.000) y errores al asignar una altitud incorrecta a un punto de cota, a una curva de nivel, o a un vértice de una polilínea. El primer tipo de error afecta a los puntos y líneas digitalizados a partir de imágenes mal georreferenciadas. Para corregirlo es necesario rectificar primero estas imágenes (del modo que se ha comentado en el apartado de “Georreferenciación” en “Materiales y métodos”) y después, con la imagen ya bien georreferenciada, corregir el trazado de las líneas y puntos, bien moviendo uno a uno los vértices mal localizados, o, en el caso de las líneas, también redibujándolas mediante el comando “Reshape Feature” de la herramienta de edición de ArcMap. La escala de trabajo sería 1:2.500. El error cometido por la imprecisión en el trazado afecta sólo a algunas curvas maestras y a unos pocos puntos de cota que se digitalizaron al principio (un ejemplo puede verse en la figura 14. Se corrigen moviendo los vértices o redibujando las líneas, como he indicado arriba. La asignación errónea de un dato de altura se corrige localizando la curva de nivel, o el punto de cota afectado y cambiando el valor erróneo por el valor real en la tabla de atributos. Si el error afecta a un vértice de una polilínea, se selecciona primero la línea y con el botón derecho del ratón se selecciona el vértice y se rectifica el valor Z.

21

El número de vértices que se han incluido en cada polilínea es quizá excesivo para ser usados como datos de entrada de tipo punto acotado (mass point) en la creación del TIN. En principio cada uno de estos vértices da lugar a un nodo en la red de triángulos, de modo que si hay demasiados vértices habrá un número excesivo de pequeños triángulos, muchos de los cuales no aportan información útil y ocupan memoria innecesariamente. No sé si el programa lleva a cabo automáticamente una generalización de vértices, eliminando los que no son necesarios en el proceso de interpolación; si no es así, quizá sería conveniente hacer esta generalización por nuestra cuenta (previamente habría que hacer una copia del archivo “Curvas de nivel” y trabajar con ella, para no perder los datos del archivo original). En la Figura 11 y en la Figura 13 se muestran ejemplos de digitalización de todas las curvas de nivel, maestras y prácticamente todas las secundarias y los puntos de cota. En la Figura 12 sólo se han digitalizado las maestras y algunos puntos de cota. 3. El Modelo Digital de Elevaciones El TIN El programa ha creado automáticamente el archivo TIN a partir de los datos de entrada que se le han suministrado. El TIN es una estructura vectorial constituida por una red de superficies triangulares adosadas entre sí por los vértices o nodos y los bordes que forman cada triángulo. Los nodos se corresponden con los vértices que se han usado como datos de entrada y conservan intacto en el TIN su valor de altitud y posición. Cada triángulo de la red contiene información de altitud, pendiente, orientación, área y perímetro, de modo que a partir de ellos se puede obtener un mapa de elevaciones, un mapa de pendientes, o un mapa de orientaciones. El mapa que se presenta en esta memoria es un mapa de elevaciones, en el que las alturas se han agrupado en 10 clases de 200 m de amplitud distinguibles con una rampa de color (la clase correspondiente al último intervalo de alturas, 1.800-2.000, no debería figurar, ya que se ha formado debido a un error en la asignación de la altura de un punto de cota; después fue corregido). En la Figura 15 se muestra la superficie TIN de casi todo el área de estudio, primero con los elementos a partir de los cuales se ha creado (curvas de nivel, puntos de cota y red hidrográfica) superpuestos, y después sin ellos. La Figura 16 muestra un fragmento de la superficie a una escala mayor. Una de las características más interesantes que tiene el MDE de estructura TIN es que puede ajustarse cuanto se desee a la complejidad del relieve, incorporando los vértices y las líneas estructurales que sean necesarios. La importancia de las líneas estructurales (hard breaklines) a la hora de modelar un relieve se pone de manifiesto cuando aparecen “zonas críticas”: valles, barrancos, cauces fluviales, dolinas, crestas, aristas de laderas, carreteras y en general cualquier elemento que suponga un cambio abrupto o una discontinuidad en la pendiente de la superficie. En estos casos la misión de una línea estructural es delimitar esta discontinuidad formando una secuencia de bordes de triángulos que se interpone entre secuencias de vértices próximos a ambos lados de la discontinuidad. En ausencia de la línea estructural se formarían triángulos entre estos vértices próximos, que ocultarían la zona crítica; el resultado sería, por ejemplo, que donde debería

22

figurar una línea de cresta aparecería una meseta, o donde debería haber un fondo de valle, un cauce de un río o de un barranco, aparecería una terraza o un “peldaño”. La Figura 17 ilustra un ejemplo de esto. En ella se muestra un tramo del valle de Soba con un tramo del río Gándara discurriendo por su fondo hacia el Norte y cuatro arroyos que confluyen en él por sus dos márgenes. Para construir esta superficie TIN se han usado las curvas de nivel (que aparecen dibujadas en azul) y la red hidrográfica (también en azul). El arrollo que desciende hacia el río Gándara por la parte superior izquierda del mapa es el Barranco del Rosario (el único elemento fluvial que contiene datos de altitud), que se ha incorporado al modelo como línea estructural (hard breakline); por otro lado, el río Gándara y los otros tres arroyos que se le unen se han añadido como líneas de inflexión o rotura (shoft breaklines). En la imagen inferior que acompaña a la figura 15 se muestra el fragmento del mapa topográfico sobre el que se ha realizado la digitalización de estas líneas (corresponde a la imagen “Concha 1”). Puede verse cómo incluyendo el río Gándara, los tres arroyos, las carreteras y las líneas de aristas, como líneas estructurales, el modelo ganaría en realismo. También se mejoraría mucho si se añadieran algunas curvas de nivel secundarias de la ladera situada en la margen derecha. La Figura 18 y la Figura 19 muestran el resultado que se obtiene cuando se incorpora una alta densidad de curvas de nivel y algunos puntos de cota. Puede apreciarse que el relieve queda bastante bien definido sin necesidad, casi, de añadir otro tipo de líneas estructurales en las zonas críticas. Quizá la mejor manera de construir el TIN sea empezar haciendo un modelo sencillo y posteriormente ir incorporando los elementos estructurales que se vayan requiriendo para ajustarlo al modelo real con el detalle que se desee. En este sentido los resultados que se muestran aquí pueden ser útiles para tener una idea de cómo actuar. El Raster El modelo raster de Matrices Regulares se ha generado a modo de prueba, para tener una idea de la precisión y resolución que pueden conseguirse en función de los datos que se aportan y de las características del relieve. La estructura de Matrices Regulares consiste básicamente en una matriz de n filas y m columnas que dan lugar a un conjunto de celdas cuadradas de un tamaño determinado; la matriz se superpone imaginariamente sobre la superficie del terreno como si fuera una retícula. Cada celda (unidad se superficie) representa una porción concreta del terreno, georreferenciada, y describe una cualidad suya determinada mediante un valor que lleva asociado; en nuestro caso esta cualidad es la altitud. El valor de cada celda se obtiene por estimación a partir de una muestra de datos conocidos. Para calcularlo primero se definen las localizaciones de los puntos problema (los puntos para los que se realiza la estimación), que se corresponden con los centros de cada celda, y después se estima por interpolación el valor en cada uno de ellos en función de los datos del entorno, proporcionados en nuestro caso por el TIN. Los valores de altitud que se obtienen no tienen porqué corresponder exactamente con los valores reales en el terreno, porque los valores estimados están influidos por los valores vecinos. Debido a su estructura de datos, el modelo raster no puede representar elementos lineales como ríos, carreteras, puntos singulares (cimas, fondos de dolinas), etc.

23

Nuestro modelo se ha creado a partir del TIN, seleccionando el atributo elevaciones para su conversión a raster y reduciendo el tamaño de celda que por defecto ofrecía la aplicación, a 25 m (no sé si ello afectó al cálculo de los valores de los puntos problema). No estoy seguro de si la conversión de TIN a Raster se realizó satisfactoriamente, porque tras seguir los pasos especificados y pulsar la tecla OK, apareció un mensaje en la pantalla advirtiendo que se había encontrado un determinado número de puntos (unos 2.000, creo recordar) con el mismo valor (¿o con un único valor?); no sé si esto puede suponer un inconveniente a la hora de llevar a cabo el proceso de conversión6. Aceptando el mensaje se mostró un típico raster en escala de grises, que posteriormente se pasó a rampa de color. Como resultado final se obtuvo la imagen coloreada, pero sin el sombreado que se esperaba, de un mapa de elevaciones, donde las alturas aparecen agrupadas (por defecto) en 10 clases, con una amplitud de 163 m. Se muestra en la Figura 20 (el área que cubre el modelo es más pequeña que la que aparece en otras figuras, porque el raster se creó antes de terminar la digitalización. Se ha incluido la hidrográfica como referencia para distinguir fácilmente las zonas bajas de las altas). Se puede tener una idea comparativa del raster con respecto al TIN con más detalle viendo la Figura 16, por un lado y la Figura 18 y Figura 27, por otro. En cuanto a la precisión del modelo, ésta está en función en primer lugar, de la precisión del TIN (porque ha sido generado directamente a partir de sus datos), de modo que cualquier error que aparezca en él se transmite al raster. El TIN puede tener errores de precisión altimétrica en las zonas de la superficie que están comprendidas entre vértices y bordes de triángulos (donde los valores sí son exactos); este tipo de errores se producen en el proceso de interpolación que se sigue para establecer los triángulos y son, supongo, inevitables. Otro tipo de errores, evitables al menos parcialmente, son debidos a la ausencia de elementos estructurales en zonas críticas (como expliqué en el apartado anterior). Así por ejemplo, si en el TIN aparece un aterrazamiento o una meseta, donde debería haber un fondo de valle o una cresta, en el modelo raster las celdas registrarán valores de altitud correspondientes a estas formas erróneas. Ejemplos de esto pueden verse en las dos siguientes figuras: La Figura 22 muestra el mismo tramo del valle de Soba que aparecía en la figura 17. En el TIN se podía apreciar un aterrazamiento en las márgenes del río Gándara (sobre todo en la parte inferior izquierda de la imagen) y en las márgenes de los dos arroyos inferiores, producido por no haberse incorporado los cauces de estos cursos de agua como líneas estructurales. En el modelo raster este aterrazamiento no puede verse por la ausencia de sombreado, pero accediendo (mediante el botón de identificación) a los valores de altitud de las celdas que comprenden estas zonas aterrazadas, se ve que estos valores coinciden con los de las curvas de nivel que las rodean (200 m). La Figura 23 muestra otro ejemplo de aterrazamiento, que aunque sin mucha importancia en cuanto a la fidelidad del modelo, es bastante ilustrativo. En este caso el aterrazamiento se produce porque una curva de nivel dibuja un bucle cóncavo en una ladera, que hace que vértices de la misma curva que están enfrentados formen triángulos entre sí. En el TIN aparece una especie de peldaño en lugar de la hondonada real (no hay imagen de esto, pero en la figura 16 se ven ejemplos similares). En el modelo raster este “peldaño” queda registrado por los valores de altitud de las celdas que están dentro del bucle, que en este caso concreto se puede apreciar muy bien porque la curva de nivel responsable (con una cota de 560 m) marca la división entre dos clases de altura en el mapa: el interior del bucle debería ser de color marrón oscuro y no claro, como aparece. 6

Sería interesante conocer el mecanismo de interpolación a partir de datos TIN (superficies, nodos y bordes de triángulos) 24

Este tipo de errores se pueden corregir localizando las formas erróneas en el TIN y añadiendo los elementos estructurales necesarios. Una vez hecho se puede generar un nuevo raster más ajustado a la realidad. La precisión del modelo raster está también determinada por factores que intervienen en el proceso de interpolación y de creación de la matriz. La topografía del terreno por un lado y el número, altitud y distribución espacial de los nodos, que determinan la estructura TIN por otro, influyen a la hora de calcular el valor de altitud asignado a cada celda. Para calcular este valor cada punto problema debe integrar los datos de altitud de su entorno, lo que será tanto más complejo cuanto mayor sea el número de éstos y más variados sus valores. De este modo si en el entorno de un punto problema (el centro de una celda) hay un número determinado de nodos (o triángulos) cuyos valores de altitud varían mucho entre sí, esto afectará mucho al cálculo de su valor de altitud y el resultado lo reflejará dando un valor más impreciso. El tamaño de celda asignado a la hora de crear la matriz determina la resolución del modelo; cuanto más pequeña sea, mayor será su número y por tanto más puntos problema habrá definidos. Más puntos definidos significa más valores de altitud en el modelo y por tanto, más resolución. La cuestión está en establecer un tamaño de celda apropiado, teniendo en cuenta la precisión que se quiere obtener para el modelo y la precisión de los datos de los que se parte para crearlo. Mediante un correcto análisis del error, realizado sobre una muestra amplia de puntos control representativos del relieve se puede: calcular el error de cada punto (error local), estimar el grado de precisión del modelo (validación) y ver si los errores de precisión detectados son aleatorios, o presentan un patrón de distribución asociado por ejemplo a determinadas formas del relieve, o a la distribución de los nodos o triángulos del TIN. El análisis del error se haría tomando los datos de altitud reales (aportados por el mapa topográfico original) en los puntos control seleccionados y comparándolos con los valores de altitud de esos mismos puntos en el modelo, de modo similar a como se propuso para analizar el error de georreferenciación. El error en cada punto sería la diferencia entre la altitud en el modelo y la altitud real. El cálculo del error medio (EM) permitiría ver si los errores observados son de tipo aleatorio (en cuyo caso el EM sería próximo a cero) o se deben a algún factor que los determina. Para estimar la precisión del modelo no se puede usar el EM, porque los errores positivos y los negativos se anulan entre sí; en su lugar se emplea el error cuadrático medio (ECM). Por último, si los errores obtenidos son reflejados en una tabla de atributos, pueden ser posteriormente representados como una capa superpuesta sobre el modelo raster o el TIN, lo que nos permitirá ver su localización en el mapa y analizar un posible patrón de distribución. La elección del número de puntos control y su localización es muy importante y hay que hacerlo con cuidado. La finalidad del análisis del error es, por un lado detectarlo para corregirlo cuando sea posible y por otro, validar el modelo, proporcionando información sobre el grado de precisión que presenta (que será necesario conocer para sus posteriores aplicaciones). En ambos casos es importante distinguir el tipo de error del que se trata, si es aleatorio, o por el contrario está asociado con algún elemento topográfico, de la estructura de datos del TIN, o con ambos. Por esta razón los puntos control deben localizarse sobre todos los tipos de estructuras topográficas existentes en el terreno: picos de montañas, cimas de montes y colinas, crestas, fondos de dolinas, fondos de valles y barrancos, cursos de agua, laderas empinadas, laderas suaves, collados, zonas llanas, etc. Como lógicamente el punto elegido debe tener una altitud conocida, debe situarse sobre un punto de cota o una curva de nivel. Para ver si la distribución de los datos del TIN tienen un papel en la aparición de errores, estas localizaciones deberán elegirse en zonas con una alta y una baja densidad de vértices o nodos de triángulos. Por otra parte quizá

25

sería conveniente elegir también algunos puntos control fuera de curvas de nivel y puntos de cota, o sea fuera de la influencia de los vértices del TIN. La altura exacta no se conoce en estas zonas, pero se sabe que está dentro de un intervalo de 10 m (la equidistancia de las curvas secundarias. El valor de altura de un punto en el modelo debe estar también dentro de ese intervalo de altura. Por falta de tiempo el análisis del error no se ha llevado a cabo. Para ilustrar un poco cómo se procedería, quiero poner algunos ejemplos donde se calculan los errores en puntos de cota. La Figura 24 muestra un fragmento de la Sierra del Hornijo (“Veguilla 2”) a escala 1:5.000. En la imagen pueden verse las curvas de nivel y seis puntos de cota; también puede apreciarse el tamaño de las celdas (25 m) y cómo éstas “tratan” de ajustarse al contorno de las curvas de nivel (donde están situados los vértices de los triángulos del TIN que han servido como datos para calcular los valores de altura en los puntos problema). La imagen que acompaña al raster en la figura 22 muestra el fragmento correspondiente del mapa topográfico, en el que se pueden ver los datos de altitud de los puntos de cota y de las curvas de nivel (su equidistancia es de 10 m). En la tabla de abajo se representan las altitudes reales y las calculadas por el modelo para los seis puntos de cota de la figura (de izquierda a derecha y de arriba a abajo). Todos los valores están expresados en metros. punto i

Coord. X

Coord. Y

Altura en MDE

Altura real

Error (ei)

1

455.964

4.787.774

1008,8151

1008

0,8151

2

456.228

4.787.862

974,5884

973

1,5884

3

456.055

4.787.298

1.242,3811

1.245

-2,6189

4

456.320

4.787.353

1.253,0746

1.258

-4,9254

5

456.572

4.787.269

1.211,5702

1.212

-0,4298

6

456.873

4.787.432

1.231,0151

1.236

-4,9849

La columna de errores muestra que en los dos primeros puntos el modelo sobreestima ligeramente la altitud real y en los cuatro de abajo la subestima. El error medio para estos seis puntos sería: EM = 1/n ∑ ei = -1,7592 m El error cuadrático medio sería: ECM = √1/n ∑ei = 0.9482 m (la raíz cuadrada afecta a todo el segundo miembro de la igualdad) 2

En Figura 25 los cuatro puntos de los que se van a calcular los errores locales se encuentran localizados en una zona muy llana, las curvas de nivel están muy separadas entre sí y por tanto los vértices de los triángulos del TIN están muy distanciados entre sí también. La imagen corresponde al pueblo de Arredondo (“Arredondo 2”), donde por cierto se encuentra situado el Centro Ecológico de la Montaña Oriental, que gestiona la Fundación Naturaleza y Hombre y donde. Los puntos están ordenados de arriba abajo y de izquierda a derecha.

26

punto i

Coord. X

Coord. Y

Altura en MDE

Altura real

Error (ei)

1

451.072

4.792.010

192,9731

192

0,9731

2

450.816

4.791.886

163,3330

162

1,3330

3

451.562

4.791.706

164,7427

164

0,7427

4

451.612

4.791.444

162,9780

162

0,9780

Todos los errores son pequeños y de sobreestimación de la altura. Quizá el hecho de que los cuatro puntos de cota estén más debajo que las curvas de nivel que los rodean (representan “depresiones” en el terreno) tenga algo que ver. EM = 1/n ∑ ei = 1,0067 m

ECM = √1/n ∑ei = 0,0286 m 2

Por último la Figura 26 muestra un ejemplo de digitalización exhaustiva de prácticamente todas las curvas de nivel y los puntos de cota, donde tan solo faltaría incluir el río que discurre por el fondo del valle y algunas líneas estructurales definiendo fondos de barranco. El aporte de datos para construir el TIN es pues máximo y la densidad de curvas de nivel en algunas zonas es tal que muchas celdas del raster se ven afectadas por más de una y más de dos. Los valores de altitud de los puntos de cota y los errores son los siguientes (de izquierda a derecha y de arriba abajo) punto i

Coord. X

Coord. Y

Altura en MDE

Altura real

Error (ei)

1

460.290

4.787.540

883,1404

885

-1,8596

2

460.520

4.787.587

803,6638

805

-1,3362

3

460.798

4.787.611

904,4124

909

-4,5876

4

461.018

4.787.574

896,6171

899

-2,3829

5

461.126

4.787.561

868,2658

868

0,2658

6

461.221

4.787.598

894,4811

896

-1,5189

7

462.846

4.786.973

717,8712

719

-1,1288

8

461.047

4.786.190

454,8694

455

-0,1306

Todas las alturas calculadas por el modelo, salvo la quinta, subestiman la altura real. EM = 1/n ∑ ei = -1,5848 m

ECM = √1/n ∑ei = 2,1086 m 2

27

Conclusiones El objetivo del proyecto no se ha cumplido, puesto que no se ha obtenido el modelo digital de elevaciones definitivo y todavía quedan pendientes cosas importantes en los pasos previos. Posiblemente el objetivo era demasiado pretencioso, pues construir un MDE a partir de cartografía en papel por el procedimiento seguido es una tarea laboriosa, que requiere mucho tiempo y además ciertos conocimientos técnicos y teóricos con los que en un principio apenas contaba. Siendo realistas podría decirse que se ha conseguido hacer la mitad del trabajo. Visto desde un punto de vista positivo, podría considerarse el trabajo realizado en estos escasos cinco meses como una primera fase del proyecto, en la que se han conseguido unos resultados más o menos satisfactorios, que no son definitivos, pero sí útiles. Un aspecto interesante del trabajo es que se ha pasado por todas las etapas del proceso de construcción del modelo, desde el escaneado de los mapas hasta la obtención del raster (tan solo faltan las fases de validación y de creación de los metadatos). Algunas de estas etapas no se han completado del todo y ofrecen resultados con errores más o menos importantes que deben corregirse, pero a partir de lo que se ha hecho se puede tener una idea de lo que falta por hacer, de los problemas que deben abordarse y más o menos del modo de resolverlos. Por otra parte, a partir del análisis de los resultados obtenidos y de los errores cometidos se puede hacer una estimación del grado de precisión que se puede conseguir y fijar unos límites de error. Las tareas pendientes en cada etapa del proceso son las siguientes: Georreferenciación de las imágenes raster Corrección de errores. Como he propuesto en el apartado de resultados de la georreferenciación, o como se considere oportuno. Validación de cada imagen del el archivo .aux, incluyendo para cada una de ellas una tabla con los errores residuales de los puntos control enlazados y del ECM de la transformación. Creación de los Metadatos para cada imagen del archivo aux. Digitalización. Corrección de errores de digitalización en las curvas maestras y los puntos de cota afectados. Corrección de errores debidos a la incorrecta georreferenciación en todos los elementos afectados. Digitalización con inclusión de datos de altitud de curvas de nivel secundarias, puntos de cota, red hidrográfica, red de carreteras y caminos, y de líneas estructurales “implícitas”: líneas de cresta, de fondo de valle y de barranco, etc. Generalización de vértices de las curvas de nivel para utilizarlos como datos de entrada de tipo punto acotado (Mass point) en la generación del TIN. Es una opción que debe estudiarse. En caso de llevarse a cabo, debe hacerse sobre una copia del fichero “curvas de nivel”.

28

Creación de los Metadatos para cada archivo. Generación del MDE TIN y RASTER. Generación del TIN a partir de los datos digitalizados y generación del Ester a partir del TIN. Validación de los dos modelos incluyendo para cada uno una tabla con los errores locales para cada punto control, el error medio y el ECM. Creación de los Metadatos para cada modelo.

Bibliografía consultada ƒ

ArcGIS Desktop Help. Manual de ayuda del programa ArcView.

ƒ

Felicísimo, A.M.: Curso de Modelos Digitales del Terreno http://www.etsimo.uniovi.es/~feli/CursoMDT/CursoMDT.html

ƒ

Moldes Teo, F.J.: Tecnología de los Sistemas de Información Geográfica.RA-MA (1995)

ƒ

Meaza, G. y col. Metodología y práctica de la Biogeografía. Ediciones del Serbal (2000)

Un libro recomendado “La diversidad de la vida”, escrito por Edward O. Wilson. Crítica (Grijalbo) 1994.

Algunos recursos en Internet Bookmarks

Ir al índice

ƒ Behavioural Environmental Biology (U. Manchester) index.html Material on-line del master de Biología de la Conservación que imparte el Departamento de Ciencias Biológicas de la Universidad Metropolitana de Manchester. Incluye tres interesantes módulos: Research Methods, Multivariate Statistics y GIS Conservation Unit. El referente a SIG es muy interesante, con mucha información y muchos enlaces (merece mucho la pena entrar). Los tres módulos pueden consultarse (para acceder a ellos hay que pinchar en “web based teaching” ƒ

ESRI Conservation Program conservation geography, activism and multicultural social change Página web del Programa de Conservación de ESRI

ƒ

ESRI Conservation Program Resources General Conservation Documento donde se informa sobre el Programa de Conservación de ESRI. Dentro hay un vínculo a Un documento donde se detalla claramente todo lo referente a este programa, incluyéndose información muy interesante sobre la aplicación de un SIG en la conservación de la Naturaleza; está en “ESRI Conservation Program Charter PDF”

29

ƒ

ESRI Página web de ESRI España. Con mucha información útil y muchos enlaces interesantes.

ƒ

Página personal de A.M. Felicísimo Una página muy interesante. Incluye material extraído de un curso que imparte en distintas universidades sobre Modelos Digitales del Terreno y que he consultado hasta la saciedad (está en Cursos”). También incluye distintos trabajos sobre el tema, realizados por otras personas (está en “Invitados”) y un pequeño glosario.

ƒ

Gabriel Ortiz Analista de Cantabria en GIS SIG (Sistemas de Información Geográfica) y cartografía. Programación en ArcGIS, Ma Es una página muy interesante (y muy bien diseñada). Con mucha información. Merece la pena visitarla. Su autor trabaja en el CIMA (Centro de investigaciones Medio Ambientales), en Cantabria y tiene un amplísimo conocimiento sobre SIG.

ƒ

Pagina de presentacion del proyecto Fauna Iberica Ofrece amplia y detallada información sobre el proyecto Fauna Ibérica que lleva a cabo el Museo Nacional de Ciencias Naturales. Pinchando en la casilla “organización” aparece el organigrama del proyecto. Pinchando luego en “Banco de Datos y SIG” y luego en el vínculo “Sistemas de Información Geográfica (SIG)”, se puede ver una pequeña explicación de cómo aplican el SIG en el proyecto.

ƒ

CTSP - The Conservation Technology Support Program Este programa concede (o concedía) becas para equipamientos, software y prácticas para ayudar a organizaciones conservacionistas a implantar un SIG. Independientemente de esto, la página tiene información interesante y enlaces con otras páginas de SIG aplicado a la conservación.

ƒ

GIS Lounge - Geographic Information Systems, Remote Sensing, Data, Analysis En esta página se ofrece información, noticias, recursos, etc. sobre SIG. (Al pinchar en el vínculo aparece un mensaje diciendo que la página no está disponible sin conexión, no hay que hacer caso; pinchando 2 ó 3 veces en “seguir desconectados” por fin se abre la página.

ƒ

Wildlife Ecology Software Es un servidor de softwares relacionados con la Ecología.

ƒ

Jenness Enterprises - GIS Consultation Ofrece servicios de consulta personalizados sobre SIG basado en ArcView y muchas aplicaciones gratuitamente (puede ser interesante).

30

Figuras Figura 1

Volver

31

Figura 2

Volver

Tomado del capítulo 6 “Métodos y Técnicas en Zoogeografía” del libro “Metodología y Práctica de la Biogeografía”. Ediciones del Serbal.

32

Volver Figura 3

33

Figura 4

Escala 1:2.500

Volver

Escala 1:2.500

34

Figura 5

Escala 1:25.000

Volver

Escala 1:10.000

35

Figura 6

Escala 1:2.500

Volver

36

Figura 7

Escala 1:10.000

Volver

37

Figura 8

Escala 1:10.000

Volver

38

Figura 9

Escala 1:10.000

Volver

raster imferior

Selaya 1

Escala 1:10.000

39

Figura 10

Escala 1:150.000

Volver

40

Figura 11

Escala 1:25.000

Volver

Figura 12

Escala 1:25.000

Volver

41

Figura 13

Escala 1:10.000

Volver

42

Figura 14

Escala 1:2.500

Volver

43

Figura 15

Escala 1:100.000

Volver

Escala 1:100.000

44

Figura 16

Escala 1:25.000

Volver

Volver

Imagen del raster mostrando el mismo área

Escala 1:25.000

45

Figura 17

Escala 1:10.000

Volver

Escala 1:10.000

46

Figura 18

Escala 1:10.000

Volver

Volver a Raster

Escala 1:10.000

47

Figura 19

Escala 1:10.000

Volver

Escala 1:10.000

48

Figura 20

Escala 1:125.000

Volver

1.537-1.700 1.374-1.537 1.211-1.374 1.048-1.211 885-1.048 722-885 559-722 396-559 233-396 70-233

49

Figura 21

Escala 1:25.000

50

Figura 22

Escala 1:10.000

Volver

51

Figura 23

Escala 1:5.000

Volver

52

Figura 24

Escala 1:5.000

Volver

1.211-1.374 1.048-1.211 885-1.048

53

Figura 25

Escala 1:10.000

Volver

54

Figura 26

Escala 1:10.000

Volver

55

Figura 27

Escala 1:10.000

Volver

56