COMITÉ SATELITAL – MINAGRIbiblioteca.inia.cl/medios/platina/descarga/Taller... · proveniente de...

COMITÉ SATELITAL – MINAGRI (CIREN, ODEPA, INIA, SAG, CONAF, INFOR, INDAP, CNR, SRM AGRICULTURA RM.)

TALLER DE FUNDAMENTOS DE PERCEPCIÓN REMOTA DOCUMENTO DOCENTE: VISUALIZACIÓN DE IMAGENES

SATELITALES.

Roberto Castro R. Ingeniero Forestal. Ms. en Ingeniería y Gestión Ambiental (UPC. España).

Doctorado en Cartografía Teledetección y SIG (UAH. España). Coordinador Area Satelital y Ambiental CIREN.

LA PLATINA. SANTIAGO. 28 y 29 de MAYO 2009

Visualización de datos imágenes satelitales

Cuaderno Pedagógico del GDTA traducido al español por: Ricardo Honorato P.

Pontificia Universidad Católica de Chile

Adaptación

Roberto Castro R. Angélica Berríos S.

Este documento fue adaptado por los autores del cuaderno pedagógico elaborado por la Agrupación para el Desarrollo de la Teledetección Aeroespacial (GDTA) de

Francia. 1992.

3

PROLOGO

El GDTA

La Agrupación para el Desarrollo de la Teledetección Aeroespacial se ubicada en Toulouse, en el Centro Espacial del complejo científico de Rangueil. Es una Agrupación de Interés Económico (GIE), fundada en 1973 y que se disolvio en el año 2002, integrad por los siguientes organismos: El CNES (Centro Nacional de Estudios Espaciales), el IGN (Instituto Geográfico Nacional), el BRGM (Oficina de Investigaciones Geológicas y Mineras), el IFP (Instituto Francés del Petróleo), el BDPA SCET-AGRI (Oficina para el Desarrollo de la Producción Agrícola) y el IFREMER (Instituto Francés de Investigaciones para la Explotación del Mar). El GDTA tenía como finalidad favorecer el desarrollo de la teledetección y la utilización de imágenes espaciales en Francia y en el extranjero. Su principal actividad es la formación, tanto básica ,como continua y de profundización destinada a especialistas en cuyos campos de actividad profesional pueden ser útiles los métodos de teledetección espacial: cartografía, agricultura, urbanismo, oceanografía, geología, prospección minera, medio ambiente, planificación forestal, etc. Los que conocimos y fuimos parte de los profesionales capacitados por el GDTA agradecemos la entrega del conocimiento en una tecnología que aporta directamente al desarrollo de los paises. Roberto Castro Ríos.

Concepción gráfica y maquetaje

Dominique LASSELIN Florence NOIRET Bruno THOLLET

Comité de lectura

Jean-Pierre ANTIKIDIS, Gérard BEGNI, Jean-Paul DARTEYRE, Guy FLOUZAT, André FONTANEL, Pierre-Yves REVILLION, Bruno THOLLET, Gilbert WEILL

4

Indice ____________________________________________________________________________

PROLOGO________________________________________________________3

INDICE___________________________________________________________4

RESUMEN________________________________________________________6

INTRODUCCIÓN___________________________________________________7

DESARROLLO Y PREGUNTAS_______________________________________8

CAPÍTULO 1 ______________________________________________________9

PRINCIPIO DE VISUALIZACIÓN DE UNA IMAGEN DIGITAL ________________9

1.1. Repaso sobre la imagen digital de teledetección __________________________________________ 9

1.2 Principio de visualización ____________________________________________________________ 12

1.3 Esquema de un sistema de visualización ________________________________________________ 13

CAPÍTULO 2 _____________________________________________________15

VISUALIZACIÓN DE UN CANAL DE UNA IMAGEN_______________________15

2.1 Visualización en niveles de gris _______________________________________________________ 15 2.1.1 Visualización bruta_______________________________________________________________ 15 2.1.2 Definición e interpretación del histograma de un canal ___________________________________ 26 2.1.3 Modificaciones de dinámica________________________________________________________ 32 • Principio __________________________________________________________________________ 32 • Modificación de tipo lineal ____________________________________________________________ 32 • Equipoblación (o equidistribución) por niveles. ____________________________________________ 35 • Equipoblación ______________________________________________________________________ 38 • Otras modificaciones de dinámica ______________________________________________________ 39

2.2 Visualización de un canal en colores ___________________________________________________ 40

5

2.2.1 Principio _______________________________________________________________________ 40 2.2.2 Definición de una tabla a color______________________________________________________ 40 2.2.3 Modificación de dinámica: Comparación con la visualización en niveles de gris ______________ 41

CAPÍTULO 3 _____________________________________________________42

GENERALIZACIÓN EN DOS O TRES DIMENSIONES: VISUALIZACIÓN SIMULTÁNEA DE VARIOS CANALES O COMPOSICIÓN COLOREADA ______42

3.1 Principio __________________________________________________________________________ 42

3.2 Composición coloreada normal _______________________________________________________ 47

3.3. Interpretación de una composición coloreada___________________________________________ 49

CONCLUSIONES _________________________________________________50

"RESPUESTAS" __________________________________________________51

ANEXOS ____________________________________________________________________________ 61

Anexo 1______________________________________________________________________________ 62

Anexo 2______________________________________________________________________________ 63

Anexo 3______________________________________________________________________________ 64

Anexo 4______________________________________________________________________________ 65

Anexo 5______________________________________________________________________________ 66

Anexo 6______________________________________________________________________________ 67

Anexo 7______________________________________________________________________________ 68

Anexo 8______________________________________________________________________________ 69

Anexo 9______________________________________________________________________________ 70

BIBLIOGRAFÍA___________________________________________________71

6

Resumen

Este documento pedagógico supone conocidas algunas bases físicas elementales de la teledetección, las nociones de satélites y de sensores, así como las nociones básicas relacionadas con el color y su representación. Estas nociones constituyen las bases necesarias para la comprensión de la visualización de una imagen digital proveniente de sensores satelitales. Vamos a abordar sucesivamente: • El principio fundamental de visualización, • La visualización de un canal y las diferentes maneras de modificar globalmente la dinámica de una imagen, a saber: • El despliegue lineal • La equipoblación (con o sin niveles), el despliegue de tipo gaussiano, etc. • La visualización de varios canales o composición coloreada que no es más que una ampliación de la visualización de un sólo canal, basada enteramente en las nociones de reproducción trícroma. Se terminará este Trabajo Dirigido con las bases de la interpretación de los principales colores presentes en una composición coloreada y constatarán de nuevo que las nociones básicas relacionadas con el color son indispensables.

7

Introducción ____________________________________________________________________________ Introducción

La visualización de imágenes es el primer problema que se plantea al usuario de datos digitales de teledetección. Este querrá igualmente mejorar globalmente la dinámica de esta imagen, con miras a un análisis de calidad de los datos, de una localización del o de los centros de interés y de una primera interpretación visual de la información contenida. Vamos entonces a estudiar el principio fundamental de la visualización, así como las diferentes técnicas que permiten modificar la dinámica de la imagen de manera global, sin cambiar el orden relativo de los valores de los pixeles de la imagen original. Todas estas nociones son las mismas, sea el soporte de visualización un film fotográfico, papel fotográfico, un producto de impresión o una pantalla de un sistema de tratamiento de imágenes.

8

Desarrollo y Preguntas

9

Capítulo 1

Principio de visualización de una imagen digital

1.1. Repaso sobre la imagen digital de teledetección Se llama función imagen toda representación de un objeto o de una escena en un plano. Existen tres formas físicas de existencia de una imagen: • forma binaria (trazo), • en niveles de gris, • en colores Una imagen digital representa la expresión numérica1 de la función imagen en un sentido general. Una imagen digital de teledetección es una matriz geométrica en dos dimensiones. Esta imagen se obtiene en general, mediante un sensor electrónico2 ubicado en una plataforma3 que capta la energía proveniente de la superficie terrestre. Esta energía recibida por el sensor es transformada en señales eléctricas por los detectores4 Estas señales continuas (analógicas) son amplificadas y digitalizadas, es decir cuantificadas de manera discontinua5 .

Los sensores están concebidos para poder captar separadamente las radiaciones de acuerdo con diferentes intervalos de longitudes de onda del espectro electromagnético, llamados canales. Se dice que son multibandas o multiespectrales y se habla entonces de imagen multiespectral. Una imagen por lo tanto, generalmente contiene varios canales. Las bandas espectrales se escogen en función de la respuesta específica de los objetos que componen la superficie terrestre. Así, a cada superficie elemental o pixel corresponden tantas mediciones digitalizadas o cuentas digitales, como canales existentes. La resolución en tierra corresponde al registro de un pixel. Equivale a 20 metros por ejemplo, en el caso de una imagen SPOT obtenida en la forma multiespectral. Las cuentas digitales en los diferentes canales son posteriormente transmitidas a una estación de recepción donde son registradas línea tras línea. Ver en la figura 1 la formación de un canal de una imagen digital de teledetección.

9 1 En la forma digital, los valores son cuantificados de manera discontinua, por oposición a la forma analógica donde los valores son

continuos. 2 o radiómetro. Es el instrumento de la toma de vistas. Ejemplo: el instrumento HRV del satélite SPOT. 3 avión o satélite. Ejemplo: el satélite SPOT. 4 en general fotodiodos. 5 Esta digitalización normalmente en un código lineal de 8 bites (256 cuentas numéricas posibles).

Capítulo 2: Visualización de un canal de una imagen

Cultivo

Agua

Bosques

Suelo desnudo

Paisaje

Representación espacial (1).

(1) Un número de pixeles representacada una de las zonas.Existen pixeles mixtos y mezclados, laproporción de esto está encorrespondencia con la resolución delsensor.

10

Números correspondientes a la señal eléctrica registrada por el sensor (2).

(2) En la frontera entre dos zonas los pixeles mixtos toman valores intermedios.

Figura 1: Formación numérica de un canal de una imagen.


A un pixel cualquiera de la imagen están asociados tres valores por canal: • su coordenada en línea Li • su coordenada en columna Cj • la cuenta digital correspondiente a la medición física realizada por el sensor sobre este pixel en un rango de longitudes de onda dado: CN(ij). En la figura 2 encontraran el ejemplo de una imagen del satélite SPOT, obtenida en el modo multiespectral (o XS). Esta imagen de la figura 2 posee tres canales denominados XS1, XS2 y XS1 que corresponden respectivamente a los rangos de longitud de ondas del verde, rojo y del infrarrojo cercano. Un pixel Pij tiene como coordenadas Li Cj y tres cuentas digitales asociadas XS1ij, XS2ij, XS1ij correspondiendo cada uno a una intensidad o valor radiométrico en el canal considerado. Observación: Una imagen digital de teledetección posee mayor o menor grado de deformación geométrica en relación a una carta, según sea su nivel de pre-tratamiento geométrico.

Figura 2. Ejemplo de estructura de una imagen multiespectral SPOT.

L = Línea, C = Columna, i = índice de línea, j = índice de columna, n = número de líneas, p = número de columnas.

11


12

1.2 Principio de visualización El principio de la visualización de un canal de una imagen digital es de asociar un nivel de gris o de un color a cada cuenta digital. Comenzaremos por abordar la visualización de un canal en niveles de gris, antes de generalizar a la visualización de este canal según una tabla de colores cualesquiera. Mejorar posteriormente el contraste de una visualización y modificar su dinámica1 se traduce de hecho, en introducir una función de correspondencia entre cuentas digitales y una gama de colores, con el fin de aumentar el contraste global de la imagen. Esto se hace en la práctica por intermedio de las tablas de codificación o LUT2. El valor posible de la cuenta digital de un canal de la imagen constituye el valor de entrada de esta tabla, el valor visualizado3 constituye el valor de salida. Las figuras 3 y 4 muestran un ejemplo de tabla de codificación presentada en forma digital y gráfica respectivamente. Sin embargo, es igualmente posible visualizar varios canales simultáneamente en síntesis aditiva de los colores. Es lo que se llama una composición coloreada. Esta constituye una generalización de la visualización de un canal en dos o tres dimensiones.

Valor de Entrada

Valor de Salida

Color o nivel de gris asociado

0 0 Negro 1 0 Negro 2 0 Negro 3 0 Negro

25 0 Negro 26 2 Gris muy fuerte 27 4 28 8

152 250 153 252 154 255 Blanco 155 255 Blanco 156 255 Blanco

255 255 Blanco

12 1 La dinámica representa la relación del nivel medio máximo al nivel medio mínimo de la potencia de una señal. 2 Look-Up Table. 3 Visto en la pantalla. Ref.Cf 1.3 Estructura de un sistema de visualización.


Figura 3: Ejemplo de código de tabla de visualización (LUT).

Figura 4: Ejemplo de tabla de visualización representada en forma gráfica. Tabla correspondiente a la figura 3.

1.3 Esquema de un sistema de visualización Un sistema de visualización de una imagen digital está constituido por: • una unidad central que representa el "corazón" del sistema, • una memoria asociada a esta unidad central, • una unidad de almacenamiento de información digital, donde está contenida la imagen digital. Esta

unidad puede ser un disco duro magnético, por ejemplo. • una carta gráfica que está constituida por dos elementos principales: • una memoria que permite almacenar los valores de las cuentas digitales de la

imagen visualizada,

13


14

• un convertidor numérico-analógico que transforma los valores de las cuentas

Esta señal es enviada a la pantalla de visualización,

cual recibe señal video.

e la unidad central recibe la orden de visualizar una imagen , va a:

transferir los datos a la memoria de la carta gráfica, rmente los datos numéricos

n señal eléctrica video enviada a la pantalla de visualización.

er en la figura 5 un ejemplo de esquema de un sistema de visualización.

ráfica asociada a esta pantalla. Un sistema de visualización de agen digital puede en ciertos casos utilizar esta carta gráfica y esta pantalla para mostrar las ágenes. Sin embargo se utiliza generalmente una segunda carta gráfica mejor adaptada y por lo tanto

na segunda pan

digitales almacenadas en señal eléctrica analógica vía una tabla de codificación.

• una pantalla de visualización o monitor. Esta pantalla está unida a la carta gráfica de lala Cuando el procesador d • leer sobre la unidad de almacenamiento los datos digitales que corresponden a esa imagen y pasarlos a la memoria central, • • pasando por las cartas de codificación, la carta gráfica convierte posterioe V • Observaciones importantes en relación con las cartas gráficas: • Según su tipo, una carta gráfica puede ser monócroma o a color. Además, está limitada en cuanto al número de colores o de niveles de intensidad restituidos: 16, 64, 128, 256, 2563. • Todo computador posee como estándar al menos una pantalla de control que permite manejarlo y por lo tanto, una carta gimimu talla para mostrarlas.

Pantalla de visualización de imagen

Computador

Pantalla de control

Señal de video

Unidad central

Teclado

Figura 5: Ejemplo de sistema de visualización de una imagen digital: microsistema de tratamiento de imágenes. El computador está integrado por la unidad central, el disco duro, el lector de disquete, la carta de imagen gráfica y la carta gráfica de la pantalla de control. La pantalla de visualización de imagen está separado.


Capítulo 2

Visualización de un canal de una imagen

2.1 Visualización en niveles de gris

2.1.1 Visualización bruta Las intensidades radiométricas de un canal bruto de una imagen digital de teledetección, se traducen en cuentas que toman valores comprendidos entre dos límites extremos, que son normalmente 0 y 2551. En este caso, hay 256 valores posibles. Se dice que los valores están codificados en 8 bites o un octeto2. Al momento de la visualización en intensidades de gris, la correspondencia establecida por las tablas de codificación entre cuentas digitales y los niveles de gris, se hace generalmente como sigue: el Negro es asignado al valor mínimo 0, el Blanco al valor máximo 255 y una gama de intensidades crecientes de gris, corresponden a los valores intermedios. Ver figura 6.

Negro Blanco

0

Escala de intensidades crecientes de

Cuentas digitales 255

Figura 6: Correspondencia entre cuentas digitales e intensidad de gris. Las cartas gráficas elaboradas permiten restituir en la pantalla los 256 niveles de gris. Ciertas cartas gráficas a color restituyen de hecho, 256 niveles en cada uno de los colores primarios Azul, Verde y Rojo, lo que permite obtener en el monitor a color, 256 niveles de gris por síntesis aditiva de los tres haces coloreados con señales idénticas.

15

15 1 Ciertos datos satelitales son sin embargo codificados originalmente de 0 a 1.023, es decir sobre 1.024 valores posibles. Se trata por ejemplo de datos del sensor AVHRR del satélite NOAA. 2 256 = 28

Sin embargo, el ojo humano no puede distinguir estas 256 intensidades de gris (o de un color cualesquiera en general). Puede distinguir alrededor de una quincena de niveles diferentes. Además, se puede estar limitado por el sistema de visualización de la imagen. Por ejemplo: ciertas cartas gráficas no restituyen más que 4, 16 o 64 intensidades de gris3. Los sistemas de restitución de imágenes digitales4 pueden presentar la misma limitación. Se puede entonces estar obligado a comprimir la información de partida, codificada en 256 valores. La figura 7 presenta como ejemplo la escala de gris obtenida en una impresora, limitada a 5 niveles de gris, incluyendo el Negro y el Blanco. En este caso, la correspondencia entre cuentas digitales y estas cinco intensidades se obtiene de la manera siguiente: Los valores comprendidos entre 0 y 51 son asignados al Negro aquellos entre 52 y 102 a un gris obscuro, etc... El intervalo [0,255] está compartido en cinco "tramos" de igual magnitud.

0 51 52 102 103 153 154 204 205 255

Figura 7: Ejemplo de correspondencia entre cuentas digitales en una escala de cinco intensidades de grises.

APLICACION: Tomemos como ejemplo, la visualización del canal XS1 de una parte5de una imagen obtenida por el satélite SPOT en la forma multiespectral sobre la ciudad de La Estrella en la Costa de Marfil con fecha 04/06/88. En adelante reutilizaremos siempre esta imagen para las diferentes aplicaciones. Observen en el anexo 1 la visualización bruta de este canal, en niveles de gris, según la correspondencia de la figura 6. PREGUNTA Nº 1: • Están restituídos todos los niveles en esta visualización bruta? Justifique su respuesta. • Qué sería necesario cambiar para mejorar esta visualización bruta? En este caso, qué sería necesario conocer a priori, en relación a los valores de la cuentas digitales?

−

25 3 Esto puede ser igualmente niveles de verde o de ámbar, que corresponden al color de la pantalla monocromática. 4 Impresores, restituidor fotográfico... 5 Las dimensiones de la parte son: 384 columnas por 512 líneas.

26

2.1.2 Definición e interpretación del histograma de un canal El histograma de valores de un canal de una imagen digital es una función que expresa el número de pixeles asociados a cada cuenta digital. Se puede deducir el número de pixeles que pertenecen a cada una de las intensidades de la escala de gris adoptada. El histograma puede representarse en forma gráfica o en forma digital: • Existen tres formas de representación gráfica de un histograma, que presentan generalmente un eje horizontal con las cuentas digitales de 0 a 255 y un eje vertical con el número de pixeles: • la forma "barras" donde el número de pixeles asociados a cada cuenta digital están representado por trazos finos verticales. Estos trazos no se juntan para dos cuentas digitales consecutivas. • la forma "columnas" donde el número de pixeles está representado por líneas verticales gruesas. Estas columnas se juntan para dos cuentas digitales consecutivas. • la forma "gráfica". El gráfico es una función continua que une cada uno de los extremos superiores de las "barras" descritas anteriormente. Estas tres formas de representación están representadas respectivamente en las figuras 8.1, 8.2 y 8.3 sobre el ejemplo ficticio de la imagen digital de la figura 1. • la forma digital de este histograma es un cuadro, que presenta de entrada los valores posibles de las cuentas digitales de la imagen, de 0 a 255, y como salida el número de pixeles asignados a cada una de ellas. Observación: la suma de los pixeles asociados a cada cuenta digital es evidentemente igual al número total de pixeles de la imagen, es decir al producto del número de líneas por el número de columnas de esta imagen. La representación digital del histograma de valores de la imagen convencional de la figura 1 se presenta en la figura 9. APLICACIÓN: Retomemos el canal XS1 de nuestra imagen de la ciudad de La Estrella y calculemos su histograma. La figura 10 les presenta el histograma de este canal bruto expresado bajo forma gráfica6 Esto permite una visualización de sus diversos modos7. El mismo histograma se expresa bajo forma digital en la figura 11.

26 6 Forma gráfica. 7 Un modo es un “pic” del histograma, que corresponde a un máximo de frecuencia, para valores característicos de la ocupación

27

Valores de cuentas digitales de la imagen presentada en la figura 1.

Cuentas digitales

Número de pixeles

Figura 8.1. Representación gráfica en forma de “barras” del histograma de valores de la imagen presentada en la figura 1.

28

Figura 8.2. Representación gráfica del mismo histograma, pero en forma de columna.

Cuentas digitales

Número de pixeles

Figura 8.3. Representación gráfica del mismo histograma, pero en forma gráfica.

Cuentas digitales

Número de pixeles

29

Figura 9. Histograma de valores de la imagen de la figura 1, representada en forma numérica.

30

Figura 10. Histograma de canal XS1 de la imagen de La Estrella representada en forma gráfica.

Cuentas digitales

Número de pixeles

PREGUNTA Nº 2 • Cómo interpretan ustedes el histograma presentado en las figuras 10 y 11? Tomando en cuenta la visualización bruta de la imagen8, a qué entidades simples corresponden los dos modos principales del histograma de este canal XS1? • Al analizar la repartición de las cuentas digitales presentes en este histograma y considerando su correspondencia con las intensidades de gris, la visualización del canal bruto del anexo 1, les parece normal, optimizada?

−

30 8 Anexo 1.

31

Figura 11. Histograma de canal XS1 de la imagen de La Estrella representada en forma numérica.

32

2.1.3 Modificaciones de dinámica1

• Principio Modificar la dinámica de una visualización, es modificar la correspondencia entre las cuentas digitales y la escala de colores escogida (en nuestro caso una escala de gris) de forma de aumentar el contraste global de esta visualización. Esto se traduce en la introducción de una tabla de codificación o LUT (Ref. 1.2 Principio de visualización). Se trata entonces de determinar una función f transformando las cuentas digitales, de manera de restituir la dinámica máxima, que va del negro al blanco en el caso de una visualización en niveles de gris. Falta ahora revisar las principales formas de cómo desplegar la dinámica de visualización de un canal de una imagen.

• Modificación de tipo lineal Efectuar una modificación de dinámica lineal se traduce en la introducción, como su nombre lo indica de una función lineal2entre los valores de entrada y los valores visualizados. Las figuras 3 y 4 representan tal función. Para determinar esta recta, conviene determinar dos límites mínimo y máximo, denominados Min y Max, entre las cuentas digitales de la imagen. Estos valores van a corresponder a los dos extremos del segmento de la recta, asociados respectivamente a los valores 0 y 255 en la tabla de codificación. Pueden definirse de manera interactiva considerando el histograma o de manera automática mediante un software, dejando por ejemplo 1% del número total de pixeles por debajo del valor Min y 1% por encima del valor Max. Los valores de salida f(x) de las cuentas digitales x comprendidas entre Min y Max serán desplegadas linealmente en el intervalo [0,255], de manera de corresponder a intensidades de gris repartidas del Negro al Blanco. La ecuación del segmento de la recta es la siguiente:

f (x) = 225 x - Min Max - Min

Observación importante: los pixeles que tienen valores x inferiores a Min serán asociados al valor f(x) = 03 y los que tienen valores x superiores a Max serán asociados a f(x) = 2554. PREGUNTA Nº 3 32 1 Todavía denominadas despliegues de dinámica, adaptaciones de dinámica, o stretching en inglés. Estas técnicas forman parte de los mejoramientos de la imagen. 2 Un segmento de recta. 3 Estarán saturados en Negro. 4 Estarán saturados en Blanco.

33

• Para mejorar lineal y globalmente el contraste del canal XS1 de la imagen de La Estrella, qué

valores Min y Max escogerían?. Observación: evidentemente no hay una respuesta exacta en lo absoluto. Estos valores son subjetivos y relativos al tema que se quiere resaltar. • Trazar encima la representación gráfica de la tabla de codificación correspondiente a la modificación

de dinámica lineal entre los límites Min y Max que acaban de determinar.

− Min =

Max =

Valor de salida o valor visual

255 200 150 100 50

Valor de entrada = valor de cuenta digital

0 50 100 150 200 250

Figura 12. Tabla de código correspondiente a la modificación de dinámica lineal aplicada al canal XS1 de la imagen de La Estrella. Representación gráfica.

34

• Los valores Min y Max escogidos corresponden a los límites absolutos de variación de las cuentas digitales? ¿Por qué?

−

El anexo 2 muestra la visualización del canal XS1 con modificación de dinámica lineal. Compárenla con la visualización bruta en el anexo 1. • Qué forma tendrá el histograma de valores de la imagen visualizada después de tal modificación lineal?. Comparen este histograma con el de la imagen bruta. Se reencuentran todos los valores comprendidos entre 0 y 255?

−

• Ayudándose del histograma, qué valores asociarían ahora a Min y Max para mejorar la visualización únicamente para la entidad "agua" de la imagen?.

− Min = Max =

En el anexo 3 encontrarán la visualización del canal XS1 que ha sufrido un despliegue lineal entre los valores Min y Max escogidos para mejorar la dinámica del agua. Compárenla con la imagen desplegada globalmente (anexo 2). • Misma pregunta para mejorar el contraste únicamente de la parte "no agua" de la imagen5.

− Min = Max =

En el anexo 4 encontrarán la visualización del canal XS1 que ha sufrido un despliegue lineal entre los valores Min y Max escogidos para mejorar la dinámica del campo terrestre. Compárenla con la imagen desplegada linealmente (anexo 2).

34 5 Campo terrestre.

35

e separar mejor visualmente los valores más representados y de aproximar los valores arginales.

• Equipoblación (o equidistribución) por niveles. Consideremos la curva de representación de un histograma cualquiera1. Sea H la integral normalizada de este histograma. En estadística, la función H se denomina función de repartición. Para una cuenta digital cualquiera x0, se tiene: H (x0) = Probabilidad {x < x0} Se denomina más comúnmente histograma acumulado del canal considerado. Es otra forma de representación del histograma de los valores de un canal. Esta función está representada en figura la 13 para nuestro ejemplo de aplicación2. Se puede leer por ejemplo que un 33.5% del número total de pixeles de este canal tienen una cuenta digital inferior a 38. Esto podría verificarse a partir del histograma del mismo canal representado bajo una forma digital (figura 11). Basta sumar el número total de pixeles que tienen cuentas digitales inferiores a 37 y dividir este número por el número de pixeles de esta imagen, o sea 512 * 384 = 196.6083.

Porcentaje de número total de pixeles

Cuentas digitales

Figura 13. Histograma acumulado del canal XS1 de la imagen de La Estrella representada en forma gráfica.

La equipoblación por niveles4. Es un despliegue de dinámica que tiende a asociar una cantidad igual de pixeles a cada nivel de gris (o nivel) de la imagen transformada. Esto tiene por consecuencia sobre la imagen, dm

35 1 El histograma del canal XS1 de nuestra imagen por ejemplo. 2 Canal XS1 de la imagen de La Estrella. 3 Producto del número de líneas por el número de columnas. 4 Como su nombre lo indica.

36

Para esto es suficiente con cortar el eje de las ordenadas de la representación de la función de repartición H(x) en un número de intervalos iguales correspondientes al número de niveles5 deseados6 sobre la imagen desplegada. La función H-1 nos da la correspondencia de cada uno de los intervalos en el eje de las cuentas digitales. Los intervalos de cuentas digitales que corresponden a cada una de las intensidades de gris son de longitud variable. Ejemplo: Con cuatro niveles, se divide el eje de las ordenadas de la representación gráfica del histograma acumulado en cuatro intervalos, representando cada uno un 25% del número total de los valores. Vía el histograma, a cada uno de estos intervalos corresponde sobre el eje de las cuentas digitales un intervalo de longitud variable. La figura 14 presenta tal división sobre el histograma acumulado del canal XS1 de nuestra imagen de referencia. El eje de las cuentas digitales se encuentra así dividido en cuatro intervalos conteniendo cada uno 25% de los pixeles: [0,28], [25,59], [60,70] y [71,255] La función de codificación asociada a tal modificación de dinámica es entonces una función llamada "en escaleras": todas las cuentas digitales de un mismo intervalo están asignadas al mismo valor y por lo tanto la misma intensidad de gris. Ver figura 15. Encontrarán en el anexo 5 la visualización del canal XS1 de la imagen de referencia, en equipoblación con cuatro niveles.

Cuentas digitales

Porcentaje de número total de

Figura 14. División sobre el histograma acumulado, el eje de las cuentas digitales se encuentra dividido en cuatro intervalos, cada uno correspondiente a 25% del número total de pixeles.

36 5 Número de niveles de gris. 6 O impuesto técnicamente.

37

Valor de salida o valor visualizado


Figura 15. Tabla de código correspondiente a una equipoblación sobre cuatro niveles aplicados a un canal XS1 de la imagen de La Estrella. Representación gráfica.

PREGUNTA Nº 4: Cuál es la forma de la representación gráfica del histograma de valores visualizados del canal XS1 de la imagen de La Estrella, que ha experimentado una equipoblación con cuatro niveles7?. - 37 -

−

37 7 Tabla de codificación de la figura 15.

V l d t d

38

• Equipoblación El método de equipoblación, tal como está definido en los sistemas de tratamiento de imágenes, puede ser considerado como un caso límite de equipoblación, en el cual el número de niveles sería fijado en 256. Una equipoblación se reduce de hecho a desplegar la dinámica de la imagen bruta según la curva del histograma acumulado. No hay reducción del número de niveles restituídos en la imagen desplegada, no hay reagrupación de valores. Contrariamente a la modificación de dinámica lineal donde los valores alcanzados por los pixeles son de intervalo constante, los espaciamientos entre los valores alcanzados en la imagen visualizada y desplegada llegan a ser de longitud variable. Observación: la diferencia visual entre los dos tipos de equipoblación se hace más pequeña mientras mayor es el número de rangos. En efecto nuestro ojo reagrupa visualmente los valores muy próximos que corresponden a niveles de gris casi idénticos, aglomerando estos valores en un mismo nivel visual8.



Figura 16. Tabla de código correspondiente a una equipoblación del canal XS3 de la imagen de La Estrella. Representación gráfica.

38 8 Nuestro ojo no es capaz de discriminar más de una veintena de niveles de gris. El número de niveles de gris perceptibles por el ojo

humano es por lo tanto inferior al número de niveles radiométricos que existen generalmente en una imagen.

39

Encontrarán en el anexo 6 la visualización del canal XS1 de la imagen de referencia en equipoblación y en la figura 16 la tabla de codificación correspondiente. PREGUNTA Nº 5: Cuál es el aspecto de la representación gráfica del histograma del canal XS1 de la imagen de estudio visualizada con este tipo de equipoblación?

−

• Otras modificaciones de dinámica • el despliegue de tipo Gaussiano tiende a dar a cada tipo de histograma la forma de una curva de Gauss9 y restablece la media de los valores visualizados en 12710. Se escoge además la desviación tipo de las cuentas digitales de la imagen resultante. Esta transformación parece bien adaptada a la curva de sensibilidad del ojo. Este consideraría como armoniosa cualquier imagen que tiene este tipo de distribución. • Sin entrar en el detalle, uno mismo puede definir una infinidad de modificaciones de dinámica. Bastaría definir una función de correspondencia entre los valores digitales de la imagen bruta y los valores visualizados. La mayoría de los sistemas de visualización permiten al utilizador trazar por si mismo de manera simple e interactiva la representación gráfica de una tabla de codificación cualquiera. Por ejemplo: esta curva puede ser: • lineal por segmentos11

• en "escaleras"2

• de forma logarítmica3 o exponencial4, etc... • Las modificaciones de dinámica más utilizadas son sin embargo, la modificación de tipo lineal y la equipoblación. La primera aparece generalmente más armoniosa al ojo, pero no hace resaltar las diferencias entre cuentas digitales próximas y muy frecuentes en la imagen. Ejemplo: las zonas homogéneas como el agua, un bosque... si ellas ocupan una superficie relativamente importante en la imagen. Al contrario, la equipoblación mejorará la discriminación en las zonas homogéneas, si estas representan una proporción importante del número de pixeles de la imagen. Al revés, la discriminación de zonas

39 9 Forma de “campana”, correspondiente a una ley de distribución normal. 10 Centrada entre 0 y 255, correspondiendo en general al color más visible y al más suave para nuestro ojo, entre los colores de la tabla. 11 Con diferentes pendientes según los intervalos. 2 Constante por intervalo. 3 Función de pendiente decreciente. 4 Función de pendiente creciente.

40

heterogéneas será generalmente menos buena que con una modificación de tipo lineal, por ejemplo en el medio urbano. La equipoblación resultará generalmente preferible en el caso donde el número de niveles visualizables es limitado5. En este caso, se tratará de una equipoblación por niveles. En efecto más vale en este caso, asociar un mismo número de pixeles a cada uno de estos niveles. No es entonces posible aconsejar una u otra de las modificaciones. Todo depende de la imagen y del objetivo buscado. PREGUNTA Nº 6: Qué valores de la imagen serán preferidos en relación a otros con una modificación de dinámica de acuerdo a una función logarítmica?, exponencial?.

−

2.2 Visualización de un canal en colores

2.2.1 Principio El principio es prácticamente el mismo que el de la visualización en niveles de gris6, la única diferencia es que se asocia un color y no una intensidad de gris a cada cuenta digital.

2.2.2 Definición de una tabla a color Para obtener una tal visualización es primero necesario definir una paleta ordenada de colores llamada tabla a color7. La mayoría de los sistemas de visualización utilizan por defecto una tabla de colores denominada "arco iris" pues ella retoma los colores del espectro por orden creciente de longitud de onda:

Violeta, Azul, Cyan, Verde, Amarillo, Naranja, Rojo8

Una visualización de acuerdo a una tabla de color "arco iris" corresponde al término de equidensidad coloreada utilizada en fotografía. Puede entonces denominarse equidensidad coloreada digital.

40 5 A 8 o 16 niveles por ejemplo. 6 Ref. 1.2 Principio de visualización y 2.1 Visualización en niveles de gris. 7 O tabla de colores. 8 No están citados aquí más que algunos colores al interior de este espectro. El paso de un color a otro se hace evidentemente pasando por diferentes colores intermedios.

41

Pero uno mismo puede definir una tabla a color cualquiera escogiendo ciertos colores9 y ordenándolos de acuerdo a un orden cualquiera. Todo el problema será entonces la interpretación de la visualización de la imagen a través de dicha tabla: se requerirá generalmente de la escala de colores a mano como referencia para poder ordenar y por lo tanto analizar la información contenida en esta imagen.

2.2.3 Modificación de dinámica: Comparación con la visualización en niveles de gris

Todos los tipos de adaptación de dinámica descritos precedentemente se aplican de la misma manera que para la visualización en niveles de gris. El anexo 7 les presenta por ejemplo la visualización del canal XS1 de la imagen de La Estrella a través de una tabla de colores "arco iris"10, a la cual ha sido aplicada la misma modificación lineal de dinámica que en la imagen que aparece en el anexo 2. La escala de gris utilizada en la visualización en intensidades de gris es evidente en su definición y en su orden. Ella es interpretable instantáneamente por el ojo, no siendo necesaria una leyenda para entenderla. Por el contrario, la interpretación de un canal visualizado con una tabla a color, necesita una cierta experiencia y a menudo, la visualización de la escala de colores para poder ordenar ciertos fenómenos. Hay igualmente un paso obligatorio para los colores más obscuros o más claros de acuerdo a los matices. PREGUNTA Nº 7: Comparemos un canal de una imagen visualizado en niveles de gris y el mismo canal visualizado con una tabla a color11. Cuál es la ventaja de la segunda visualización en el caso en que las dos escalas tengan un mismo número de niveles.? −

41 9 256 como máximo, correspondiente a cada una de las cuentas digitales posibles de la imagen. 10 Definido por los tonos que van del azul al rojo. 11 Y con la misma adaptación dinámica.

42

Capítulo 3

Generalización en dos o tres dimensiones: visualización simultánea de varios canales o

composición coloreada

3.1 Principio La visualización simultánea de varios canales es posible en síntesis aditiva de los colores1. El principio de visualización de cada uno de los canales es similar al procedimiento de visualización de un canal en niveles de gris. La composición coloreada es una generalización de este principio en dos o tres dimensiones. Conviene en un primer lugar escoger dos o tres canales de una misma imagen multiespectral y de asociarles, uno de los tres colores primarios Azul, Verde y Rojo. El número máximo de canales visualizables simultáneamente es por supuesto de tres2. Una tabla de codificación está asociada a cada uno de los tres colores. El razonamiento que sigue está relacionado con la selección de los tres canales. Observación: En el caso de SPOT, los tres canales están asignados. Falta escoger la asignación de los colores a cada uno de ellos. Aplicación: Vamos a construir una composición coloreada (o CC) de la imagen de referencia de La Estrella, asociando el Azul al canal XS1, el Verde al canal XS2 y el Rojo al canal XS1. Vamos entonces a visualizar en síntesis aditiva: • el canal XS1 en niveles de Azul3

• el canal XS2 en niveles de Verde4

• el canal XS1 en niveles de Rojo5

PREGUNTA Nº 8 42 1 Ref. Cuaderno Pedagógico B3: Representación de colores. 2 Escoger dos es igualmente posible y simplifica el razonamiento que sigue. 3 Yendo del negro al azul intenso pasando por una gama ordenada de azules, como si se aplicara un filtro azul intenso sobre la visualización en niveles de gris de este canal. 4 idem que (47) con un filtro verde. 5 idem que (47) con un filtro rojo.

43

• La visualización bruta de los tres canales en síntesis aditiva sobre la pantalla será óptima?, por qué?. • Determinen para cada uno de los canales los valores Min y Max con miras a una modificación de dinámica lineal, sabiendo que las figuras 17 y 18 presentan los histogramas de los canales XS1 y XS2 bajo forma gráfica y la figuras 19 y 20 los representan bajo forma digital. Encontrarán en el anexo 8 el resultado de la composición coloreada que acabamos de elaborar, con modificación de dinámica lineal en los tres canales, restituída con cinco niveles posibles en cada uno de los tres colores primarios.

−

Observación importante: Es naturalmente posible efectuar cualquier tipo de modificación de dinámica en cada uno de los canales que componen la composición coloreada. El anexo 9 les presenta la misma composición coloreada con equipoblación en los tres canales.

44


Cuenta digital

Número de pixeles


Cuenta digital

Número de pixeles

45Figura 19. Histograma de canal XS1 de la imagen de La Estrella presentada en forma

46

Figura 20. Histograma de canal XS2 de la imagen de La Estrella presentada en forma numérica.

47

3.2 Composición coloreada normal En numerosos sensores, se encuentran los tres canales correspondientes a los rangos de longitudes de ondas siguientes1. Verde, Rojo, Infrarrojo cercano • Para el instrumento HRV de SPOT, se trata respectivamente de los canales: XS1, XS2, XS1 • Para el Landsat MSS los canales: MSS4, MSS5, MSS7 • Para el Landsat Thematic Mapper; los canales TM2, TM3, TM4 El color “Infrarrojo cercano" no existe realmente, no es posible asociar a cada canal el color correspondiente a su rango de longitudes de onda. Una composición coloreada normal consiste entonces en "desplazar" estos rangos de longitudes de onda y a asociarles los colores primarios en el orden de clasificación correspondiente a su longitud de onda, a saber: Azul, Verde, Rojo Se trata de hecho de un simple deslizamiento de colores respetando el orden lógico de los canales:

Longitud de Onda Color asociado

Verde Rojo Infrarrojo Cercano

Azul Verde Rojo

Es este tipo de visualización que acabamos de elaborar2. Se la llama composición coloreada en falsos colores. Es análoga a otro tipo de imagen de teledetección: la fotografía aérea en infrarrojo color. El film fotográfico correspondiente denominado IRC, posee tres capas de gelatina sensibles a las partes Verde, Rojo e Infrarrojo cercano del espectro electromagnético3. Al momento del revelado, cada una 47 1 Clasificadas por orden creciente. 2 Figuras 8 y 9. 3 A diferencia de una película “colores naturales” utilizado por el fotógrafo aficionado y que presenta tres capas de gelatinas sensibles a los rangos Azul, Verde y Rojo del espectro visible.

48

de estas capas está respectivamente asignada a los colores primarios Azul, Verde y Rojo, obteniéndose por lo tanto una imagen análoga a una composición coloreada de falsos colores. Esto explica en parte la formación de la composición coloreada en falsos colores: las personas que ya han trabajado con tales fotografías, anteriores a las imágenes digitales, poseían ya las claves de la interpretación de una composición coloreada de este tipo. Los mismos elementos aparecen ahí en colores cercanos. En los dos casos, hay el mismo "deslizamiento" de los colores en relación al rango de longitudes de onda que les es asociado. Todas las otras asignaciones de colores son evidentemente posibles. En el caso de una imagen SPOT multiespectral, la elección está impuesta puesto que no hay más que tres canales, pero se puede al menos permutar los colores, asociando por ejemplo XS1, XS2, XS1 a los colores respectivos Rojo, Azul y Verde. Faltaría por cierto, interpretar visualmente los colores de tal composición coloreada... PREGUNTA Nº 9: • Cuántas composiciones coloreadas diferentes es posible elaborar con una imagen SPOT multiespectral? • Misma pregunta para una imagen Landsat Thematic Mapper4.

−

48 4 Seis canales sin el canal térmico.

49

3.3. Interpretación de una composición coloreada Al estar elaborada la composición coloreada por síntesis aditiva de los colores, es relativamente fácil de descomponer sus diferentes colores en "cantidades" relativas de colores elementales Azul, Verde, Rojo y de llegar así a una estimación relativa de la respuesta del elemento en los diferentes canales que componen esta visualización. Se llega así a crear una firma espectral y a emitir una hipótesis de la naturaleza del elemento estudiado. Se alcanza el procedimiento clásico de fotointerpretación consistente en definir las "claves" de identificación o de interpretación visual de los objetos. Ejemplo: en el caso de nuestra imagen de referencia, una parcela de color rojo vivo en la composición coloreada en falsos colores corresponde a una superficie que tiene una débil reflectancia en los canales XS1 y XS25 y una fuerte reflectancia en el canal XS16. Conociendo la firma espectral típica de la vegetación, está parcela debe por lo tanto, corresponder a una cobertura vegetal con actividad clorofiliana en fecha del 4/06/88. Si además, se conoce la región y el calendario agrícola, se puede deducir la(s) ocupación(es) del suelo posible(s). PREGUNTA Nº 10 • De la misma manera pueden ustedes analizar un color Negro o Azul Noche. • Misma pregunta para un color Blanco y Cyan muy claro.

−

49 5 A los cuales se ha asociado respectivamente los colores Azul y Verde. 6 Al cual se ha asociado el Rojo.

50

CONCLUSIONES

Uds. poseen ahora las nociones relacionadas con la visualización, la modificación de la dinámica de una imagen de teledetección, así como las bases de la interpretación de ésta. Les falta ahora aplicar estas técnicas en un sistema de tratamiento de imágenes. Los primeros tratamientos que Uds. acaban de descubrir, están enfocados a obtener el máximo de la discriminación visual de la información, sin perder el orden relativo entre los valores de los pixeles de la imagen original (modificación global). Existen técnicas matemáticas y físicas más sofisticadas para realizar este aumento del contraste global: este aumento del contrataste global: el análisis en componentes principales, el análisis factorial de las correspondencias Sin embargo, estas técnicas generan neo-canales donde la información contenida no está en correspondencia directa con la reflectancia en un rango dado de longitudes de onda. Siempre en esta categoría de tratamientos enfocados a discriminar la información, existen otros tipos de transformaciones, locales o espaciales esta vez, basadas en el tratamiento del vecindario de cada pixel. Hay una consideración de la homogeneidad o de la heterogeneidad local. Se trata de filtrados, de optimizaciones locales de densidad...

51

"RESPUESTAS"

52

− Respuesta Nº 1: * Todos los niveles no son restituidos en la visualización bruta del canal XS1 presentado en el anexo 1. No hay blanco y gris claro, salvo de manera muy puntual. Evidentemente, la dinámica global de esta imagen no está optimizada, esta visualización es mucho más obscura. Los pixeles de la imagen bruta no tienen por lo tanto valores repartidos en las 256 cuentas digitales posibles. Toda la dinámica posible no es por lo tanto utilizada. Se puede afirmar a la vista de esta imagen y de la escala de niveles de gris que le corresponden (Figura 6), que no hay (o muy pocas) cuentas digitales que tengan valores elevados en el canal bruto XS1. * Se podría mejorar esta visualización asociando el Negro al valor mínimo de las cuentas digitales de la imagen y el Blanco al valor máximo. * Pero para esto, las informaciones "estadísticas" que se acaban de deducir del análisis visual de esta visualización bruta, son insuficientes. Tendremos la necesidad de conocer la repartición exacta de los pixeles en los 256 valores posibles, es decir, el número asignado a cada valor. Esta función de repartición se llama el histograma del canal considerado y está definida en el punto 2.1.2. Esto nos permitiría conocer de manera precisa el intervalo de variación de las cuentas digitales.

− Respuesta Nº 2: * El histograma de este canal bruto nos muestra que los pixeles tienen valores comprendidos entre dos limites inferior y superior, y que al interior de este intervalo no hay equirepartición de valores en la imagen. Existen dos máximos locales principales, que corresponden a los valores 27 y 681; se dice que este histograma es bimodal2. Tomando en cuenta la visualización bruta de la imagen, el primer modo corresponde al agua3, el segundo al dominio terrestre. * Tomando en cuenta el histograma de los valores de este canal y de las interpretaciones precedentes, esta visualización es completamente normal, pero no está del todo optimizada. Todos los niveles de gris no son utilizados, la imagen no tiene entonces la dinámica óptima. Para mejorar esta visualización, sería necesario modificar la dinámica de la imagen de manera que los pixeles tengan valores visualizados repartidos en la totalidad del intervalo [0,255], y que de esta manera todos los niveles de gris estén repartidos del Negro al Blanco. 52 1 14.445 Y 4.478 pixeles son respectivamente asociados. 2 Posee dos “pics” desde el punto de vista gráfico. 3 Mar y laguna.

53

53

− Respuesta Nº 3: * La Figura 11 nos muestra que los limites absolutos de variación de las cuentas digitales de este canal son: Min = 22 Max = 215 El valor mínimo absoluto aparece significativo, la curva representativa del histograma4 presenta una fuerte pendiente a partir de la cuenta digital 21, siendo inmediatamente importantes los pixles asignados a estos valores bajos. Por el contrario, el valor máximo no es del todo representativo. Un sólo pixel le es asociado y este valor es aislado. De manera general, los valores comprendidos entre 96 y 215 no se distinguen de la representación gráfica debido a su escaso número comparado al de los valores fuertemente representados. Este un problema de escala gráfica, los 15 pixeles asociados al valor 103 por ejemplo representan apenas el 0.1% de los 14.445 pixeles asociados al valor 27. Escogeremos entonces como limites representativos: Min = 22 Max = 96 Observaciones:

* La representación gráfica del histograma es más fácil de interpretar. Expresa visualmente las proporciones relativas de los pixeles asignados a dos valores cualesquiera.

* El valor máximo escogido (96) no constituye la única respuesta posible. Se puede escoger un valor ligeramente diferente sin que la visualización se encuentre modificada de manera significativa.

* La representación gráfica de la tabla de codificación lineal aplicada al canal XS1 de la imagen de La Estrella, entre los valores Min = 22 y Max = 96 ,se presentan en la figura 21. Los pixeles que tienen un valor inferior a Min son asignados al valor 05; aquellos con un valor superior a Max son asociados a 255; los otros son distribuidos entre 0 y 255 de acuerdo a una recta. * Por las razones invocadas precedentemente, sólo el valor mínimo escogido 22 corresponde a un limite absoluto. El valor máximo escogido es netamente inferior al valor máximo absoluto. * Las figuras 22 y 23 presentan el histograma de los valores visualizados del canal XS1 de la imagen de Abdijan en el cual ha sido aplicada una modificación de dinámica lineal entre los valores 22 y 96.

4 Figura 10. 5 Esto no concierne a nuestro caso de estudio donde el valor Min corresponde al valor mínimo absoluto de los valores de la imagen.

54

54



Figura 21. Tabla de código correspondiente a una modificación de dinámica lineal aplicada al canal XS3 de la imagen de La Estrella entre los valores Min = 22 y Máx = 96. Representación gráfica.

Siendo lineal la modificación de contraste, no se cambia la forma general del histograma; al nivel gráfico, hay simplemente un estiramiento según el eje de las cuentas digitales en relación al histograma de la imagen bruta6, la disposición de la curva guarda exactamente la misma forma. Por el contrario, como acaban de constatarlo en la figura 22 ,hay una "ruptura" del histograma, todos los valores comprendidos entre 0 y 255 no están presentes. En efecto, cuando se despliegan (Max-Min) valores en el intervalo [0,255], todos los valores no pueden entonces ser considerados por las cuentas digitales transformadas: hay "hoyos". En nuestro caso particular por ejemplo ,se despliegan: (96-22) = 75 valores, hay entonces alrededor de una cuenta digital por cada tres o cuatro que está "ocupada" (1/4 < 75/256 < 1/3). La imagen que ha sufrido una modificación de dinámica no presenta más niveles de gris que la imagen bruta, pero estos están mejor repartidos del Negro al Blanco.

6 Cambio de escala según el eje de las X.

55

Cuentas digitales

Número de pixeles

Figura 22. Histograma de canal XS3 de la imagen de La Estrella que ha sufrido una modificación de dinámica lineal entre los valores Min = 22 y Máx = 96. Representación gráfica en forma de campana de los valores visualizados.

* Para mejorar el contraste únicamente de la entidad "agua" de la imagen, valores posibles para los límites de variación son: Min = 22 Max = 37 • Para mejorar el contraste únicamente de las partes "terrestres" de la imagen, valores posibles para

los limites de variación son: Min = 43 Max = 96

− Respuesta Nº 4: La representación gráfica del histograma de valores visualizados del canal XS1, que ha sufrido una equipoblación con cuatro niveles, será un conjunto de cuatro columnas o barras7, de tamaño sensiblemente equivalente y repartidos en un espaciamiento constante en el intervalo [0,255].

55 7 Representando cerca del mismo número de pixeles.

56

Figura 23. Histograma de canal XS3 de la imagen de La Estrella que ha sufrido una modificación de dinámica lineal entre los valores Min = 22 y Máx = 96. Representación gráfica en forma numérica.

57

Encontraran este histograma en la figura 24.

− Respuesta Nº 5: Econtrarán en la figura 25 la representación gráfica del histograma del canal XS1 visualizado en equipoblación, sin reducción del número de niveles visualizados. Podrán constatar que la escala de dilatación según el eje de cuentas digitales es variable y no es lineal como en la figura 22, que representa el histograma del mismo canal visualizado con una modificación de dinámica lineal. Los pixeles que tienen valores próximos y poco representados en la imagen bruta ven sus valores transformados muy próximos. Por el contrario, las cuentas digitales más frecuentes se encuentran discriminadas por esta transformación.

− Respuesta Nº 6: Teniendo en cuenta las formas respectivas de las funciones logarítmicas8 y exponenciales9, la primera transformación privilegiará la discriminación visual de los valores radiométricos débiles, y la segunda la de los fuertes.

− Respuesta Nº 7: Al no poder el ojo humano distinguir más que una quincena de niveles de gris diferentes en una misma paleta, la visualización de una imagen de una tabla de colores será interesante en el caso donde se desee visualizar nuestra imagen con un numero de niveles más importante. En efecto, nuestro ojo distingue evidentemente muchos más colores que niveles de gris, alrededor de 64 a 128 matices, bastante más fácilmente de acuerdo al entrenamiento personal y con una gama relativamente ordenada10. Por el contrario, si se utiliza un número de niveles más restringido, la interpretación de la imagen en niveles de gris será globalmente mejor. En efecto la escala de gris es evidente en cuanto a su orden y directamente interpretable por el ojo. Una escala de colores necesita un cierto aprendizaje visual.

− Respuesta Nº 8: * Lo mismo que para la visualización bruta de un canal, la visualización bruta de tres canales en síntesis aditiva no es óptima.

57 8 Curva de pendiente decreciente. 9 Curva de pendiente creciente. 10 Muchos menos si los ordenes lógicos no son respetados.

58

Cuentas digitales

Número de pixeles

Figura 24. Histograma de canal XS3 de la imagen de La Estrella que ha sufrido una equipoblación con cuatro niveles. Representación gráfica en forma de columnas de los valores

Cuentas digitales

Número de pixeles

Figura 25. Histograma de canal XS3 de la imagen de La Estrella que ha sufrido una equipoblación. Representación gráfica en forma de campana de los valores visualizados.

59

En efecto la dinámica máxima no es utilizada en ninguno de los tres colores, pues las cuentas digitales brutas de cada canal están repartidas más que en una fracción del intervalo [0,255]. Además, la media de estos valores brutos puede ser bien diferente de un canal al otro y la composición coloreada resultante tendrá entonces una dominante en el color asociado al canal donde la media de las cuentas digitales es la más fuerte. Para nuestro ejemplo de aplicación, la dominante será Roja, pues el canal XS1 presenta la media más elevada. Para mejorar la visualización de la composición coloreada precedente, es necesario naturalmente conocer antes que todo la repartición de los valores de los pixeles en cada canal, es decir el histograma de cada uno de ellos. * Los valores Min y Max de adaptación de dinámica lineal dados aquí no son más que indicativos, no tienen un carácter absoluto, de la misma for-a que para los limites Min y Max del canal XS1 que han sido determinados en el punto 2.1.3. Para XS1 se puede dar: Min = 39 y Max = 97 Para XS2: Min = 28 y Max = 93 Para XS1, se vuelve a considerar: Min = 22 y Max = 96

− Respuesta Nº 9: * Con una imagen SPOT multiespectral que posee tres canales: XS1, XS2 y XS1, es posible elaborar seis composiciones coloreadas diferentes. En efecto, si tomamos los tres colores primarios Azul, Verde y Rojo, tenemos tres canales posibles para la asociación con el primer color, nos quedan dos para el segundo y uno para el tercero. De donde 3x2x1 = 6 posibilidades. En términos de análisis combinatorio, este número representa la cantidad de permutaciones posibles de tres elementos11. * Para una imagen Thematic Mapper con seis canales12, el principio es el mismo. Hay seis posibilidades para la selección del primero, cinco para el segundo y cuatro para el tercero. Hay entonces 6x5x4 = 120 posibilidades para la selección de canales en un orden determinado, por lo tanto 120 composiciones coloreadas posibles. En el análisis combinatorio, este número corresponde al número de arreglos posibles de tres canales entre seis13.

59 11 Valor factorial de tres que se escribe 3¡ = 6. 12 TM1 a TM5, y TM7, el canal TM6 tiene una resolución diferente a la de los otros canales (120 m en lugar de 30) y no puede ser combinado directamente con los otros canales. 13 Siendo 6¡/(6¡-3¡) = 120.

60

− Respuesta Nº 10: * El color Azul-noche corresponde a una superficie que tiene valores relativamente débiles en el canal XS1, pero valores todavía más débiles en los otros dos canales asociados a los colores Verde y Rojo. Conociendo las propiedades espectrales de algunos elementos típicos14, se puede deducir que se trata probablemente del agua (mar o laguna) que absorbe mucho en los rangos de longitud de onda del visible y del infrarrojo cercano. * Una parcela de color Blanco o Cyan claro será una superficie que tiene reflectancias elevadas en los tres rangos de longitudes de onda15. Se tratará entonces de una superficie brillante. Un suelo desnudo claro, seco y liso, de arena a lo largo del litoral, techos muy brillantes en zona urbana...

60 14 Agua, suelo desnudo, vegetación activa. 15 Recuerden la formación del Blanco en síntesis aditiva de los colores...

61

ANEXOS

Anexo 1 Visualización bruta en niveles de gris del canal XS1.SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 2 Visualización en niveles de gris del canal SX1. Modificación de dinámica lineal entre valores : Min = 22 y Max = 96. SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 3 Visualización en niveles de gris del canal XS1. Modificación de dinámica lineal entre los valores: Min = 22 y Max = 37. SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 4 Visualización en niveles de gris del canal SX1. Modificación de dinámica lineal entre valores : Min = 43 y Max = 96. SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 5 Visualización en niveles de gris del canal SX1. Modificación de dinámica de tipo de equipoblación de cuatro niveles.SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 6 Visualización en niveles de gris del canal SX1. Modificación de dinámica de tipo de equipoblación. SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 7 Visualización de tabla de colores del canal XS1 de la banda 1. Modificación de dinámica de tipo lineal entre los valores: Min = 22 y Max = 96.SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 8 Modificación cuadrática del histograma la composición infrarrojo color de los canales 3, 2 y 1 (RGB). SPOT1. Sector La Estrella, VI Región. 27/04/98.

Anexo 9 Equalización lineal de la composición infrarrojo color del SPOT (3, 2 y 1, RGB).SPOT1. Sector La Estrella, VI Región. 27/04/98.

62

Anexo 1

Visualización bruta en niveles de gris del canal XS1. SPOT1. Sector La Estrella, VI Región. 27/04/98.

63

Anexo 2

Visualización en niveles de gris del canal SX1. Modificación de dinámica lineal entre valores : Min = 22 y Max = 96

SPOT1. Sector La Estrella, VI Región. 27/04/98.

64

Anexo 3

Visualización en niveles de gris del canal XS1. Modificación de dinámica lineal entre los valores: Min = 22 y Max = 37.


65

Anexo 4

Visualización en niveles de gris del canal SX1. Modificación de dinámica lineal entre valores : Min = 43 y Max = 96


66

Anexo 5 Visualización en niveles de gris del canal SX1. Modificación de dinámica de tipo de equipoblación

de cuatro niveles. SPOT1. Sector La Estrella, VI Región. 27/04/98.

67

Anexo 6 Visualización en niveles de gris del canal SX1. Modificación de dinámica de tipo de equipoblación.


68

Anexo 7

Visualización de tabla de colores del canal XS1 de la banda 1. Modificación de dinámica de tipo lineal entre los valores: Min = 22 y Max = 96.


69

Anexo 8

Modificación cuadrática del histograma la composición infrarrojo color de los canales 3, 2 y 1 (RGB).


70

Anexo 9

Equalización lineal de la composición infrarrojo color del SPOT (3, 2 y 1, RGB). SPOT1. Sector La Estrella, VI Región. 27/04/98.

71

Bibliografía • AEs Image, 1989. Manuel d'utilisation du logiciel Didactim, version 3.

• BLONDEL F.M., 1990. Manuel d'utilisation du logiciel de traitement d'images TTTUS, version

2.02. INRP (Institut National de la Recherche Pédagogique).

• CASTRO, R. 1999. BASES TEÓRICAS. Sistema para el Seguimiento y Análisis de Tierras

mediante Teledetección (TeleSAT). Proyecto: GCP/RLA/126/JPN, FAO.

• CAP SESA REGIONS, 1990. Manuel d'utilisation du logiciel Multiscope, version 2.0.

• FISCHER A., 1991. Bases physiques de la télédétection. Cahier pédagogique B 1, version 1.0.

GDTA.

• FOIN P., 1985. Cours de télédétection, fascicule 2 : Bases sur la connaissance du milieu naturel

et humain. Ecole Nationale des Sciences Géographiques (ENSG), Institut Géographique

National (IGN), St-Mandé.

• FOIN P., 1988. Cours de télédétection, fascicule 1 : Bases physiques. Ecole Nationale des

Sciences Géographiques (ENSG), Institut Géographique National (IGN), St-Mandé.

• LASSELIN D., 1991. Représentation des couleurs. Cahier pédagogique B3, version 1.2.

GDTA.

• LEO O. et DI~ER J.L., 1985. Télédétection, Techniques et applications cartographiques.

Editions FORHOM, Bmeau pour le Développement de la Production Agricole (BDPA).

• LEVESQUES A., 1987. Proposition d'utilisation des satellites Landsat 5 (TM) et SPOT 1

pottrl'inventaite et 1'aménagement foiestiers. Rapport de troisiéme année d'étude á 1'ENTTEF,

IAURIF.

• MALON J.F., 1984. Principaux traitements des images satellitaires et associées. BRGM

(Bureau des Recherches Géologiques et Miniéres).

• MS2i (Matra Sepimage Imagerie Informatique), 1988. Manuel d'utilisation du logiciel

Géopéricolor, version 2.2. Guide technique des systémes Péricolor 1500/2001.

• Open universiteit, 1989. Remote sensing course. Heerlen.

• Société FranGaise de Photo~tammétrie et Télédétection. 1988. Glossaire des tetmes ofFciels de

la télédétection aérospatiale. Commission enterministérielle de terminologie de la télédétection

aérospatiale (COMITAS). Commissariat général de la langue française. Bulletin nº 112/1988-4.

COMITÉ SATELITAL – MINAGRIbiblioteca.inia.cl/medios/platina/descarga/Taller... · proveniente de...

Documents

Transcript of COMITÉ SATELITAL – MINAGRIbiblioteca.inia.cl/medios/platina/descarga/Taller... · proveniente de...