.
CURSO INTRODUCTORIO DE ANÁLISIS DE IMÁGENES SATELITALES Y
SISTEMAS DE INFORMACIÓN GEOGRAFICA

2E: Operaciones booleanas elementales (parte del Sistema de Análisis Geográfico)

Vamos a realizar ahora un ejercicio simple en el cual buscaremos áreas hipotéticamente favorables para un cierto cultivo (por ejemplo alguna especie de flor exótica). Para esto supongamos que esta planta requiere una altura entre los 1000 y 2000 metros sobre el nivel del mar. Adicionalmente supongamos que las pendientes deben estar entre los 0 y 10% y que una plantación rentable no debe de bajar de 4 hectáreas.
Como ejemplo usaremos el DEM de una zona alrededor de la confluencia del río Ene con el río Cutivireni que contiene algunos de los mejores atractivos naturales del Perú.

Imagen B-1: DEM de zona Ene-Cutivireni

En base a este DEM podemos calcular las áreas con alturas entre 1000 y 2000 metros sobre el nivel del mar. Para esto se efectúa una reclasificación con los parámetros siguientes:

La imagen resultante se presenta a continuación. A este tipo de imágenes se les denomina "imagen booleana" pues sólo tiene valores de unos o ceros [0 = color negro = altitudes menores que 1000 o mayores que 2000 metros] y [1 = color rojo = altitudes entre 1000 y 2000 metros].

Imagen B-2: Areas (en rojo) entre 1000 y 2000 msnm

A continuación efectuamos una operación similar a la de la sección anterior para determinar todas las áreas con pendientes menores a 10 porciento, hallando la siguiente imagen booleana:

Imagen B-3: Areas (en amarillo) con pendientes menores a 10%

El siguiente paso es realizar una "operación booleana" sobre las dos imágenes anteriores, que consiste en hallar todas las areas que simultáneamente cumplen con las dos condiciones iniciales que hemos planteado para nuestro ejercicio. Esto equivale a hallar la intersección de dos conjuntos: el conjunto de las áreas con altitudes entre 1000 y 2000 metros y el conjunto de las áreas con pendientes menores a 10 porciento.
Esto lo realiza fácilmente un SIG ya que lo que hace es multiplicar los valores de los píxeles de una imagen por los de la otra. Como los píxeles tienen sólo valores de unos o ceros, la imagen booleana resultante tendrá píxeles con valor uno que simultánteamente cumplen con las dos condiciones planteadas y los demás píxeles tendrán valor cero, de acuerdo a la siguiente tabla:

  Imagen B-2 Imagen B-3 Resultado
Ninguna de las dos condiciones 0 0 0
Alturas entre 1000 y 2000 msnm 1 0 0
Pendientes menores a 10 0 1 0
Ambas condiciones 1 1 1

La siguiente imagen booleana es el resultado de la operación anterior:

Imagen B-4: Areas (en rojo) que cumplen dos condiciones simultáneamente

El área total que cumple con las dos restricciones planteadas es de 12,478 Hectáreas:

El problema con este resultado es que existen demasiadas áreas pequeñas que no son últiles para nuestro objetivo. Es necesario entonces separar las áreas en grupos individualizados como se ve en la siguiente imagen:

Imagen B-5: Areas de la imagen B-4 en grupos separados

El siguiente paso es determinar el área de cada uno de los 3,346 grupos formados y usar reclasificación para eliminar los grupos menores a 4 Hectáreas. La ventana de diálogo de reclass en el SIG IDRISI que usaremos es la siguiente:

La imagen resultante con los grupos individuales que cumplen con todas las condiciones planteadas inicialmente se muestra a continuación.

Imagen B-6: areas separadas con más de 4 Has

Haga click aquí para ver la lista de las áreas de los grupos mostrados en la imagen B-6. Los códigos de la 2ª columna son los que aparecen en el SIG al picar sobre los polígonos respectivos. El número total de áreas que cumplen las 3 condiciones de nuestro problema es de 449, cifra muy inferior a la que corresponde a la imagen B-5, lo cual reduce considerablemente la búsqueda de posibles áreas para sembrar nuestra exótica planta.