Captulo 5

Formato raster. Algebra de mapas.

En el tema anterior se expusieron las ideas bsicas referentes a la codicacin y almacenamiento de informacin espacial en formato raster, as como los fundamentos de la visualizacin y consulta de dicha informacin. Ms all de estas utilidades bsicas, uno de los usos fundamentales de los Sistemas de Informacin Geogrca es la obtencin de nuevas capas de informacin a partir de otras previamente disponibles. Para ello se dispone de un conjunto de herramientas de clculo con matrices de datos que reciben el nombre genrico de lgebra de mapas. El lgebra de mapas incluye un amplio conjunto de operadores que se ejecutan sobre una o varias capas raster de entrada para producir una o varias capas raster de salida. Por operador se entiende un algoritmo que realiza una misma operacin en todas las celdillas de una capa raster. Estos operadores se denen mediante ecuaciones, por ejemplo el operador B = A 100 genera una nueva capa (B) asignando a cada celdilla el valor de la celdilla correspondiente en la capa (A) multiplicado por 100. (gura ??). Cada capa raster es una matriz de nmeros y la operacin se realiza para todos los nmeros de la matriz, por tanto para todas las celdillas de la capa raster. Sin embargo aunque ambas operan sobre matrices, el lgebra de mapas tiene poco ms que ver con el concepto matemtico de lgebra matricial.

5.1. Operadores en lgebra de mapas

Pueden denirse innitos operadores, aunque normalmente se clasican en funcin de las celdillas implicadas en el clculo en: Operadores locales Operadores de vecindad o focales Operadores de bloque Operadores de area 79

80 Operadores de area extendida Operadores globales

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

5.1.1. Operadores locales

Los operadores locales generan una nueva capa a partir de una o ms capas previamente existentes. Cada celdilla de la nueva capa recibe un valor que es funcin de los valores de esa mismo celdilla en las dems capas. Zmx,y = f (Z1x,y , Z2x,y , .., Znx,y ) la funcin representada por f () es una combinacin de operadores y funciones de distinto tipo: Aritmtico (suma, resta, multiplicacin, divisin, raiz cuadrada, potencia, ...) Lgico (AND, OR, XOR, NOT) Relacional (>, >=, <, <=, ==, !=) Trigonomtrico (sen, cos, tan, ...) Condicional (si condicin entonces instrucciones sino instrucciones) En cuanto a operadores de tipo aritmtico, puede tratarse de operaciones sencillas como la multiplicacin por 100 de una capa de altitudes en metros para obtener una capa de altitudes en centmetros o la suma de 12 capas de precipitacin mensual para obtener una capa de precipitacin anual. Un ejemplo sencillo de utilizacin del lgebra de mapas con un modelo real, sera la estimacin de la erosin potencial mediante el modelo de la USLE cuya ecuacin fundamental es: Ex,y = Rx,y Kx,y Lx,y Sx,y Cx,y Px,y (5.2) (5.1)

donde R es la erosividad de la lluvia, K la susceptibilidad del suelo a ser erosionado, L la longitud de la pendiente, S la pendiente, C un factor que depende del tipo de cultivo y P un factor que depende de las prcticas de cultivo. Si se dispone de una capa raster para cada uno de estos factores resulta muy simple obtener una estimacin de la erosin1 . Otro ejemplo algo ms complejo sera el clculo del parmetro I de Topmodel2 que se calcula mediante la ecuacin:
Sin embargo obtener cada una de estas capas requiere operaciones bastante ms complejas y muchas capas de informacin Topmodel es un modelo hidrolgico muy extendido en el que el parmetro I representa la capacidad de cada celdilla para producir escorrenta
2 1

5.1. OPERADORES EN LGEBRA DE MAPAS

81

I = ln(

a ) tan

(5.3)

donde es la pendiente, que se calcula mediante un operador de vecindad, y a es el rea tributaria especca (el rea de la cuenca vertiente a la celdilla dividido por la anchura de la celdilla) que se calcula mediante un operador de vecindad extendida. La inclusin de operadores lgicos y condicionales permite elaborar operadores complejos en los que no slo los valores sino tambin las operaciones que denen la capa de salida dependan de los valores en las capas de entrada. El caso ms simple sera la reclasicacin o clculo del nuevo valor a partir del valor de la misma celdilla en otra capa (gura ??) en base a un conjunto de reglas sencillas de reclasicacin. Por ejemplo, muchas metodologas de estudios de impacto ambiental no trabajan directamente con los valores de pendiente sino con las pendientes reclasicadas. La operacin de reclasicacin se hara, suponiendo que la capa S contiene la pendiente en tanto por ciento y que queremos generar una capa de pendientes reclasicadas a la que llamaremos Sr , con rdenes del tipo: si S < 3 entonces Sr =1: Pendiente baja si 3 <= S AND S < 10 entonces Sr =2: Pendiente moderada si 10 <= S AND S < 25 entonces Sr =3: Pendiente media si 25 <= S AND S < 50 entonces Sr =4: Pendiente alta si S >= 50 entonces Sr =5: Pendiente muy alta Como puede verse, en este caso la reclasicacin implica pasar de una variable cuantitativa a una semicuantitativa de manera que pueden asignarse etiquetas de texto a los identicadores numricos utilizados. Todos los programas de SIG disponen de un mdulo que permite reclasicar capas de informacin siguiendo criterios similares. En GRASS, el mdulo encargado de reclasicar se denomina r.reclass, en la gura ?? aparece un ejemplo del tipo de codicacin que utiliza:

5.1.2. Operadores de vecindad o focales

Los operadores de vecindad adjudican a cada celdilla un valor que es funcin de los valores de un conjunto de celdillas contiguas, en una o varias capas. El conjunto de celdillas contiguas a la celdilla X ms ella misma constituye una vecindad. Generalmente se trabaja con vecindades de forma cuadrada y tamao variable, el tamao se dene como el nmero de celdillas que hay en el lado del cuadrado, siempre un nmero impar (3, 5, 7, etc.). Los ejemplos ms habituales de operador de vecindad son el ltrado de capas , los operadores estadsticos y los operadores direccionales.

82

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

Figura 5.1: Operador local aritmtico

Figura 5.2: Operador local lgico

5.1. OPERADORES EN LGEBRA DE MAPAS

83

Figura 5.3: Reclasicacin

Filtrados El ltrado (gura ?? de capas acta moviendo una ventana, generalmente de 3x3 celdas, por toda la capa que recorre la capa entera. Esta ventana adjudica al celdilla central la media ponderada de los valores en las 9 celdillas que abarca la ventana. A partir de diferentes coecientes de ponderacin se consiguen diferentes resultados. La siguiente ecuacin calculara la media aritmtica. Zx1,y1 + Zx1,y + Zx1,y+1 + Zx,y1 + + Zx,y + Zx,y+1 + Zx+1,y1 + Zx+1,y + Zx+1,y+1 9 (5.4)

Zmx,y =

Utilizando diferentes coecientes de ponderacin podemos obtener ltros de distinto tipo. Este tipo de operadores se utiliza especialmente en el anlisis de imgenes de satlite. Operadores estadsticos Calculan un estadstico (media, desviacin tpica, mnimo, mximo, etc.) a partir de los valores de todas las celdillas que forman la vecindad y lo adjudican a la celdilla central en la capa de salida. Por ejemplo se puede establecer un ndice de diversidad del paisaje que sera el nmero de usos del suelo diferentes en las 25 (5x5) celdillas ms prximas a una celdilla dada. Operadores direccionales Permiten estimar un conjunto de parmetros relacionados con la ubicacin de los diferentes valores dentro de la vecindad. Su utilidad primordial es el anlisis de Modelos Digitales de Terreno (pendiente, orientacin, curvatura, etc.). Se estudiarn ms en detalle en el tema 6.

84

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

Figura 5.4: Operador de vecindad, ltro de la media

1 4 7 7 7 7

2 5 5 5 5 5

3 3 3 4 7 8

4 3 6 5 8 9

5 4 5 6 7 8

4 3 4 7 7 9

N N N N N N

N 3.7 4.8 5.5 6.1 N

N 3.8 4.3 5.3 6.2 N

N 4 4.3 5.7 6.9 N

N 4.2 4.8 6.1 7.3 N

N N N N N N

Cuadro 5.1: Ejemplo de operador de vecindad para el clculo de la media

5.1.3. Operadores de bloque

Actan de forma similar a los operadores de vecindad, pero en lugar de aplicar la ventana a cada una de las celdillas desplazandola celdilla a celdilla, se aplica a bloques completos del mismo tamao de la ventana y asignando el mismo resultado a todas las celdillas incluidas dentro del bloque. Las guras ?? y ?? muestran la diferencia entre un operador de vecindad y otro de bloque. En ambas se aplica un operador de media al mismo conjunto de datos. Los operadores de bloque permiten llevar a cabo cambios de escala en las capas raster, tal como aparece en la gura ?? el resultado puede considerarse como un aumento del tamao de las celdillas. Otra diferencia importante entre los operadores de vecindad y los de bloque es que mientras que los primeros exigen denir una vecindad de un nmero impar de celdillas, con los operadores de bloque pueden denirse ventanas de cualquier tamao.

5.1. OPERADORES EN LGEBRA DE MAPAS

1 4 7 7 7 7 2 5 5 5 5 5 3 3 3 4 7 8 4 3 6 5 8 9 5 4 5 6 7 8 4 3 4 7 7 9 3.7 3.7 3.7 6.1 6.1 6.1 3.7 3.7 3.7 6.1 6.1 6.1 3.7 3.7 3.7 6.1 6.1 6.1 4.2 4.2 4.2 7.3 7.3 7.3 4.2 4.2 4.2 7.3 7.3 7.3 4.2 4.2 4.2 7.3 7.3 7.3

85

Cuadro 5.2: Ejemplo de operador de bloque para el clculo de la media

5.1.4. Operadores de vecindad extendida, operadores extendidos

Son aquellos que afectan a zonas relativamente extensas que cumplen determinado criterio pero cuya localizacin precisa no se conoce a priori. Por tanto el operador (programa) debe determinar en primer lugar cual es el rea (conjunto de celdillas) que cumple dichas caractersticas. Entre los casos ms habituales estn: 1. Identicacin de polgonos homogeneos en mapas de variable cualitativa. Se parte de una variable cualitativa, como pueden ser los usos del suelo, y se obtiene un mapa en el que a cada una de las agrupaciones homogeneas con el mismo valor se le asigna un identicador nico (gura ??). De este modo se convierten en un mapa de entidades (polgonos). Por ejemplo, en un estudio de viabilidad de reintroduccin de una especie que habita en el bosque, no slo es importante saber las hectareas de bosque que tiene el territorio, sino tambin cual es el tamao de cada uno de los polgonos de bosque y cuales son los que cumplen los requerimientos de espacio de la especie a reintroducir. 2. Areas situadas a una distancia, inferior a un valor umbral, de un objeto denido por una o varias celdillas. Se genera as una zona tampn (buffer). El resultado sera una capa en la que se codicara de un modo el objeto de otro modo el rea tampn y de un tercero el rea exterior al tampon (gura ??). Siguiendo con el ejemplo anterior, la especie a reintroducir puede requerir la presencia de una supercie de agua a una distancia determinada. Este operador permitir discriminar cuales son las celdillas que cumplen esta condicin. 3. Lineas de ujo y cuencas de drenaje (gura ??). A partir del operador de direccin de drenaje visto anteriormente, puede construirse otro que de modo recursivo genere la linea de ujo que seguira un volumen de agua depositado sobre el territorio. La unin de todas las lineas de ujo que coinciden a un punto constituye la cuenca de drenaje de ese punto. Dada un celdilla, incluye todas las que drenan a este. Se trata de un operador fundamental en el desarrollo de modelos hidrolgicos. El operador linea de ujo permite determinar que celdillas se vern afectadas por un posible vertido de una industria sabiendo su ubicacin. El operador cuenca de drenaje permite evaluar el riesgo de inundacin y el clculo del parmetro I del modelo Topmodel que se ha visto con anterioridad.

86

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

Figura 5.5: Operador extendido. Transformacin de un mapa de variable cualitativa en un mapa raster de objetos

5.1. OPERADORES EN LGEBRA DE MAPAS

87

Figura 5.6: Operador extendido. Clculo de distancias

88

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

Figura 5.7: Operador extendido. Clculo de lineas de ujo y cuencas

Figura 5.8: Operador extendido. Clculo de cuencas visuales

5.1. OPERADORES EN LGEBRA DE MAPAS

89

Figura 5.9: Operador de rea: Altitud media por reas

4.

Cuencas visuales (gura ??). Determina el area visible desde un determinado punto e igualmente el area desde la que dicho punto es visible. Ejemplos de la utilidad de este operador pueden ser el estudio de impacto visual de un vertedero (cuanto mayor sea el rea desde la que es visible mayor ser su impacto) o el diseo de una red de torretas de vigilancia de incendios forestales (se tratara de maximizar la cuenca visual de cada una de ellas)

Tanto las lineas de ujo y cuencas de drenaje como las cuencas visuales se aplican de forma casi exclusiva a los Modelos Digitales de Terreno por lo que se tratarn con ms detalle en el tema correspondiente en el que, adems, se profundizar acerca de las aplicaciones de los SIG en el campo de la Hidrologa.

5.1.5. Operadores de rea o zonales

Son aquellos que calculan algn parmetro (supercie, permetro, ndices de forma, distancias, estadsticos) para una zona previamente conocida. Puede tratarse de diferentes niveles de una variable cualitativa (supercies con diferente litologa por ejemplo) o digitalizada e introducida por el usuario. Uno de los casos ms habituales es la obtencin del valor medio de una variable cuantitativa para diferentes valores de una variable cualitativa (gura ??). Por ejemplo obtener la altitud media para cada tipo de uso de suelo a partir de un mapa de elevaciones y de otro de usos del suelo. Otro caso interesante es el clculo de la supercie ocupada por cada uno de las entidades en un mapa de polgonos o por cada uno de los valores de una variable cualitativa (en ambos casos el operador es el mismo). La tilidad de este tipo de operadores es que permiten superar los anlisis celdilla a celdilla caractersticos de los operadores locales o de vecindad para llevar a cabo anlisis sobre agrupaciones homogeneas de celdillas. Permitirn por ejemplo calcular la pendiente media de una parcela de terreno o determinar la extensin de las diversas manchas de bosque existentes en un rea para determinar cuales tiene un tamao adecuado para la reintroduccin de especies.

90

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

5.1.6. Operadores globales

Son aquellos que afectan a toda la capa raster, se basan en el concepto de distancia. Operador distancia euclidiana, permite calcular para todas las celdillas su distancia a una serie de entidades en una capa raster. Operador distancia ponderada, permite introducir el concepto de friccin que es el coste de atravesar cada celdilla. Si por ejemplo tratamos de determinar el camino ms corto a pie hacia un yacimiento en mitad de una sierra, debera tenerse en cuenta la mayor o menor facilidad para atravesar las diferentes celdillas (debido a la pendiente, vegetacin, presencia de caminos, etc.) Operador entidad ms cercana, complementario de los anteriores ya que indica, para cada celdilla, cual es la entidad ms cercana. Resulta de utilidad para resolver problemas de distribucin de servicios (hospitales, colegios, etc.) La utilizacin del operador extendido linea de ujo sobre capas que representen distancias o distancias ponderadas permite determinar la distancia ms corta entre dos puntos teniendo en cuenta los diferentes factores que inuyen en el coste de atravesar una celdilla.

5.2. Modelizacin cartogrca en formato raster

El uso del lgebra de mapas supone la superacin de la fase inicial de presentacin y consulta de datos en un SIG e iniciar un uso ms avanzado (anlisis, modelizacin de procesos, toma de decisiones, etc.) que nos permita la resolucin de problemas espaciales concretos. Para ello, es necesario no quedarse en una concepcin del lgebra de mapas como un conjunto de operadores que se utilizan de forma separada, y entender el lgebra de mapas como un lenguaje con el que es posible expresar modelos ms o menos complejos que permiten resolver problemas.

5.2.1. El lgebra de mapas como lenguaje

Los diferentes programas de SIG han adoptado diferentes estrategias para la implemetacin del lgebra de mapas. Los primeros programas de SIG incorporaban mdulos para los operadores ms comunes (pendiente, orientacin, aritmtica). Con el tiempo se fueron desarrollando mdulos de propsito general que permitan establecer operadores locales simples, operaciones matemticas entre mapas o con un slo mapa; mientras que los operadores de vecindad o rea ms comunes seguan realizndose en mdulos aparte debido a su dicultad. Las ltimas versiones de los SIG ms avanzados como GRASS disponen de un mdulo que acta como intrprete de un lenguaje de programacin que permite realizar operaciones locales, de vecindad, de vecindad extendida o de rea. Los operadores de vecindad extendida ms complejos (cuencas de drenaje o cuencas visuales) siguen realizndose en mdulos aparte ya que son difciles de programar y siguen unos esquemas muy poco exibles que siempre se van a ejecutar del mismo modo.

5.2. MODELIZACIN CARTOGRFICA EN FORMATO RASTER

91

En el Departamento de Geografa de la Universidad e Utrech se ha desarrollado PCRaster (pcraster.geo.uu.nl). Se trata de un SIG relativamente barato que es exclusivamente un interprete de un lenguaje de lgebra de mapas muy potente que incluye operadores de vecindad extendida muy potentes y fciles de programar. Se utiliza fundamentalmente en estudios de tipo hidrolgico y geomorfolgico. En general pueden distinguirse dos posibles estrategias no exclusivas: Incluir un lenguaje de lgebra de mapas que permita incorporar los diferentes operadores como funciones con las que elaborar modelos complejos mediante una formulacin similar a la de la aritmtica. Mdulos (pequeos programas) que ejecuten cada uno de los operadores sobre una serie de capas de entrada generando diversas capas de salida. Estos operadores pueden activarse pulsando un botn en una Interfaz Grca de Usuario o bien escribiendo la orden en un terminal de texto. El conjunto de un mdulo junto con las capas de entrada y salida y los parmetros requeridos constituye una orden dada al sistema. La resolucin de un problema espacial mediante modelizacin cartogrca exigir la ejecucin ordenada de una serie de rdenes, en denitiva la implementacin de un determinado algoritmo para resolver el problema. De este modo el lgebra de mapas se convierte en un lenguaje de programacin Por ejemplo, en el captulo de lgebra de mapas del tutorial de GRASS aparece un conjunto de rdenes que permiten emular un operador de bloque (que no existe en GRASS), se trata en denitiva de un pequeo programa que constituye un operador de bloque.

5.2.2. El lgebra de mapas como sistema

Siguiendo un enfoque sistmico de lo que es un SIG, el lgebra de mapas es un conjunto de operadores, y cada operador constituye un componente del sistema de informacin que transforma una serie de entradas de informacin en una serie de salidas. Las entradas y salidas no son slo capas espaciales sino que tambin pueden ser tablas o visualizaciones en pantalla o impresora (gura ??). Este tipo de grcos se denominan diagramas de ujo y son de uso comn tanto en dinmica de sistemas como en programacin de ordenadores. Se han propuesto diversos mtodos de elaboracin de estos diagramas, una buena opcin es representar los mapas mediante rectngulos y los mdulos mediante valos. Las relaciones entre unos elementos y otros se representa mendiante echas. Pusto que los operadores de lgebra de mapas toman una o varias capas de entrada y producen una capa de salida, el anlisis SIG puede concebirse como una especie de mecano lgico mediante el cual diferentes piezas (operadores) se ensamblan formando anlisis complejos que, en denitiva, consituyen nuevos macrooperadores. Podemos pensar tambin en el conjunto de operadores como en componentes electrnicos o piezas de un mecano que pueden enlazarse para congurar el sistema de acuerdo a nuestras necesidades, es decir con el objeto de resolver un problema. En la parte inferior de la gura ?? puede verse como estos mdulos se combinan para solucionar un problema de lgebra de mapas planteado, se parte de 3 mapas y el resultado es un mapa y una tabla. La resolucin de problemas no va a involucrar slo a operadores-mdulos de SIG, sino tambin a otros programas presentes en

92

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

el sisema operativo e incluso herramientas desarrolladas por el usuario. En algunos casos estas herramientas se desarrollan a partir de la combinacin de algunos operadores que da lugar a un nuevo macro-operador (gura ??).

5.2.3. Conceptualizacin y formulacin de un problema espacial mediante modelizacin cartogrca

La aplicacin de estas ideas a la resolucin de problemas en un entorno SIG requiere, en primer lugar, la descomposicin del problema planteado en una serie de subproblemas y estos a su vez en otros, hasta que queden convertidos en una serie de operadores, interconectados en un diagrama de ujo, que se transeren y transforman informacin espacial. El primer paso para desarrollar un modelo que permita resolver un problema espacial es determinar como debe actuar el modelo, para ello hay que empezar teniendo claro cual es el resultado nal que se quiere obtener (generalmente uno o varios mapas y algunos informes) y aplicar el siguiente procedimiento recursivo: 1. 2. 3. Determinar que operador/es producira el mapa o los mapas resultado Determinar los mapas de entrada a dichos operadores Disponemos de esos mapas? SI. Terminado NO. Considerar estos mapas de entrada como nuevos mapas resultado y volver, para cada uno de ellos, al paso 1 En la gura ?? puede verse un ejemplo de aplicacin de esta metodologa, el problema estaba planteado as: Se requiere la seleccin de reas con una elevacin inferior a 200 m y pendiente inferior a 5 %. Estas reas deben a su vez estar a menos de 100 m de un curso uvial y no estar en cultivo. Para resolverlo se dispone de una capa vectorial con la red de drenaje (rios_vect), una capa raster de usos del suelo (usos) y una capa raster de elevaciones (elevaciones). Para resolverlo3 el primer paso sera determinar con que operacin podemos conseguir el mapa resultado y cuales son los mapas de entrada a dicha operacin. En este caso necesitamos un mapa de elevaciones (que ya tenemos), un mapa de pendientes (que podemos obtener a partir del de elevaciones) y un mapa de distancia a la red de drenaje. El mdulo que evaluar en que celdillas se cumplen las condiciones impuestas es r.mapcalc; y el mdulo que obtiene las pendientes a partir de las elevaciones ser r.slope.aspect. Para obtener un mapa de distancias a la red de drenaje tendremos que rasterizar esta (mdulo v.to.rast) y utilizar el mdulo r.buffer para establecer las distancias.
A continuacin se detalla la resolucin con GRASS, la resolucin con cualquier otro SIG ser igual pero cambiando los nombres de los mdulos
3

5.3. ESTUDIO DE VOLMENES Y MODELOS DINMICOS

93

El desarrollo de un proyecto SIG consistira de esta manera en la divisin sucesiva de un problema en subtareas cada vez ms simples hasta el momento en que cada una de estas subtareas pudiera expresarse como un operador de lgebra de mapas (gura ?? parte baja); este es el concepto fundamental del trabajo con un SIG. Una vez conceptualizado el modelo es necesario formularlo, es decir convertir el esquema en un conjunto de rdenes cuya ejecucin sucesiva resolver el problema, es decir en un programa. Por tanto el lgebra de mapas se convierte en un lenguaje de programacin y el trabajo con un SIG en el desarrollo de algoritmos. Existen programas de SIG que permiten convertir estos esquemas directamente en programas. Una vez que se ha explicitado un determinado trabajo en SIG como una serie de ordenes escritas en el lenguaje formal del lgebra de mapas lo que tenemos es un programa con todas las de la ley. Por tanto podemos escribirlo en un chero y utilizarlo con otra base de datos diferente. Evidentemente para hacer esto necesitamos un SIG basado en comandos, es decir basado en un lenguaje. Sin embargo un conjunto de rdenes no constituye un programa, es necesario en muchos casos disponer de otros elementos como estructuras condicionales o iterativas que evalen condicones. Por tanto es necesario incluir las ordenes del SIG en un lenguaje formal de suciente potencia. Generalmente se utilizan los denominados lenguajes de scripts. Finalmente cabe destacar cuatro ideas bsicas sobre el lgebra de mapas: Conocer el conjunto de operadores de lgebra de mapas no signica saber modelizar. El lgebra de mapas se aprende con la prctica y viendo ejemplos El lgebra de mapas crece ya que pueden crearse nuevos operadores a partir de la combinacin de otros operadores formado programas Los problemas complejos pueden dividirse siempre en subproblemas simples

5.3. Estudio de volmenes y modelos dinmicos

Debido a la fuerte inuencia que la tradicin cartogrca ha tenido sobre el desarrollo de los SIG, tanto la representacin de la realidad en formato raster como los operadores de lgebra de mapas tienen un carcter fundamentalmente bidimensional representando las diferentes capas diferentes aspectos del rea de trabajo. La incorporacin de la tercera dimensin (ms all de la representacin 2.5D que implican los modelos digitales de terreno) implicara la divisin de la misma en capas y cada una de estas capas se representa mediante una capa raster diferente4 . La implementacin de operadores tridimensionales no resulta demasiado complicada por lo que se reere a los operadores locales, bastara con repetir la misma operacin en todas las capas. Operadores ms complejos requeriran la repeticin de los mismos en diferentes capas y posteriormente un operador local entre las capas:
Ahora cada celdilla representara un volumen de dimensiones resx , resy y resz . A estos volmenens se les suele denominar voxels
4

94

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

Por ejemplo un operador de vecindad para calcular la media implicara el clculo de un operador de vecindad de media para cada capa y posteriormente calcular para cada celdilla la media entre los resultados de la capa, la capa superior y la capa inferior (gura ??). Los modelos llevados a cabo con tcnicas de modelizacin cartogrca son fundamentalmente estticos, es decir no incorporan la dimensin temporal. La implementacin de modelos dinmicos en un SIG es un problema an no del todo resuelto que requerira disponer de una capa raster para cada intervalo temporal. Esto supone la necesidad de manejar un nmero enorme de capas, especialmente por que en los modelos dinmicos se suelen considerar mltiples variables, lo que va normalmente ms all de las capacidades de un SIG y de los tamaos normales de los discos duros. Una solucin es la eliminacin de capas intermedia.

5.4. Errores
Toda operacin de lgebra de mapas implica la estimacin de una variable a partir de otras, por tanto implica la existencia de errores. Los errores pueden proceder de las fuentes u originarse en el proceso de digitalizacin y posterior procesado de los datos. Adems los errores pueden propagarse y aumentar al llevar a cabo operaciones con los datos. Sera deseable que las fuentes de los datos informaran acerca de los errores esperables en los datos. La digitalizacin supone errores en la localizacin de los objetos digitalizados. Un error de digitalizacin de un milmetro en un mapa 1:50.000 supone un error de 50 metros que se aade al propio error de partida del mapa. Algunos de los operadores previamente tratados tienen capacidad para propagar y multiplicar la magnitud de los errores. Si una pendiente tiene un error de + 5 % y elevamos al cuadrado un valor de 20, el resultado estar entre 225 y 625 con lo que la indeterminacin aumenta de forma considerable.

5.5. Bibliografa
De Meers, M.N. (2002) GIS Modeling in Raster John WIley & sons, Chichester 203 pp. Tomlin, C.D. (1990) Geographic Information Systems and Cartographic Modeling Prentice Hall,249 pp. Tomlin, C.D. (1991) Cartographic modeling en D.J. Maguire; M.F. Goodchild and D.W. Rhind (Eds.); 1991 Geographical Information Systems: Principles and Applications, John Wiley & sons, Chichester (www.wiley.co.uk/wileychi/gis/resources.html)

5.5. BIBLIOGRAFA

95

Figura 5.10: Ejemplos del concepto de operador

96

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.

Figura 5.11: Desarrollo de nuevos operadores con algebra de mapas

5.5. BIBLIOGRAFA

97

Figura 5.12: Esquema de resolucin de problemas con algebra de mapas.1

Figura 5.13: Emulacin de un operador tridimensional de media con 4 llamadas a operadores convencionales

98

CAPTULO 5. FORMATO RASTER. ALGEBRA DE MAPAS.