Recursos naturales renovables

 

Sobre la clasificación de sistemas arbolados usando información espectral multi–angular

 

On the classification of tree systems using spectral multi–angular information

 

Alejandro Cano–González, Fernando Paz–Pellat* , Martín Bolaños–González, Enrique Palacios–Vélez, Enrique Mejía–Sáenz, José L. Oropeza–Mota, Rene Valdez–Lazalde, José Chávez–Chan, Alfonso Zarco–Hidalgo

 

Colegio de Postgraduados, Campus Montecillo, Km. 36.5, Carretera México–Texcoco, 56230. Montecillo, Texcoco. Estado de México. *Autor responsable: (pellat@colpos.mx).

 

Recibido: Agosto, 2007.
Aprobado: Marzo, 2009.

 

Resumen

El uso de información espectral multi–angular se ha justificado como una estrategia para aumentar la precisión de los sistemas de clasificación de cultivos y de la vegetación natural, los cuales ahora sólo usan información espectral, con un esquema de análisis multivariado o de árboles de decisión, entre otros. En este trabajo se discuten los esquemas de caracterización de la información espectral multi–angular asociada a la vegetación, asi como su uso en los sistemas de clasificación. Para revisar la factibilidad de uso de datos espectrales multi–angulares, se diseñó un experimento tipo maqueta con sistemas arbolados, cinco especies forestales y se realizaron mediciones asociadas a la función de distribución bidireccional de las reflectancias. La información multi–angular fue modelada en forma compacta y usada para definir un parámetro general, la pendiente, que engloba toda la variación angular de las reflectancias. Los resultados muestran que, dejando fijo el fondo de la vegetación, es posible discriminar sistemas arbolados y que cuando el fondo de la vegetación varía hay confusión con coberturas aéreas bajas, la cual se reduce al incrementarse ésta.

Palabras clave: BRDF, clasificación, sistemas arbolados, sensores remotos.

 

Abstract

Multi–angular spectral information has been used to increase precision in classifying crops and natural vegetation. These classification systems now use only spectral information with a scheme of multivariate analysis or of decision trees, among others. In this paper, the characterization schemes of multi–angular spectral information associated with vegetation and its use in classification systems are discussed. To review the feasibility of using multi–angular spectral data, a maquette–type experiment was designed for tree systems with five forest species. Measurements taken were associated with the bidirectional reflectance distribution function (BRDF). The multi–angular information was modeled in a compact form and used to define a general parameter, the slope, which comprises all the angular variation of the reflectances. The results show that, with the vegetation background remaining fixed, it is possible to discriminate tree systems, and when the vegetation background varies and there is little canopy cover, confusion arises but decreases in the measure that canopy cover increases.

Key words: BRDF, classification, tree systems, remote sensors.

 

INTRODUCCIÓN

La clasificación de la vegetación, asignación de tipos homogéneos de vegetación a una clase u otorgar una etiqueta de clase a un píxel de una imagen satelital, usando sensores remotos ha sido un objetivo de largo plazo de la tecnología satelital. Para clasificar la vegetación, entendida en términos de la clase asociada a la cobertura (parcial o total) de las plantas sobre una parcela o píxel de una imagen satelital, se ha usado una gran variedad de sensores remotos en plataformas espaciales, donde predominan los que son gratuitos y con resolución espacial suficiente para análisis regionales o globales. Dentro de los sensores de uso gratuito, uno de los más usados para clasificar la vegetación es el AVHRR, con resolución espacial a nadir de 1.1 km, el cual ha estado disponible desde 1978. El sensor MODIS, con resolución a nadir de 250, 500 y 1000 m, en operación desde 2000, es usado intensivamente en el mapeo de la vegetación a escalas regionales y globales. Ambos sensores son de barrido ancho, por lo que los píxeles en una escena son vistos con diferentes ángulos cenitales de visión e iluminación solar, generando un problema de geometría sol–sensor. La información multi–angular asociada a la vegetación, geometría sol–sensor, ha sido la base para desarrollar sensores operacionales como el POL–DER, que tiene hasta catorce ángulos de observación y una resolución espacial de 7 km; el ATSR–2 con dos ángulos de visión y una resolución de 1 km; y, el MISR, con nueve ángulos de visión y resolución espacial de 275 m a 1.1 km. Así, la información espectral multi–angular es inherente a los sensores actuales de cubrimiento global y alta resolución temporal.

Los sistemas de clasificación de la vegetación basados en sensores remotos se han fundamentado en análisis de agrupaciones (Alvarez et al., 2003), análisis logístico de regresión (Seto y Kaufmann, 2005), conjuntos borrosos (Foody, 1996), conjuntos rugosos (Yun y Jianwen, 2006), teoría matemática de la evidencia (Cayuela et al., 2006), redes neuronales (Carpenter et al., 1979), árboles de decisión (Running et al., 1995; Friedl y Brodley, 1997) y reglas generadas por sistemas expertos (Daniels, 2006). La evolución temporal de índices de vegetación como el NDVI o las reflectancias o temperaturas radiativas es usada para mejorar las clasificaciones (Justice et al., 1985; Nemani y Running, 1997). Uno de los sistemas de clasificación de la vegetación con mayor uso es la del grupo del sensor MODIS/MISR (Knyazikhin et al., 1998a y b; Zhang et al., 2002a y b). Un punto crítico de este sistema es el uso de las relaciones del índice espectral NDVI con la RFAA (radiación foto–sinteticamente activa absorbida) y el IAF (índice de área foliar), en sus patrones temporales y espectrales. La validez de la relación entre el NDVI y la RFAA y el IAF ha sido discutida y, en ciertas condiciones, se ha establecido una relación más o menos lineal (Asrar et al., 1992; Myneni y Williams, 1994). Con un enfoque de clasificación basado en árboles de decisión, Nemani y Runing (1997) usaron la relación entre el NDVI y la temperatura de superficie, Ts, para discriminar diferentes clases de vegetación, usando balances de energía y agua como paradigma. El problema de este enfoque es que requiere un canal térmico de los sensores y los multi–angulares rara vez cuentan con este tipo de canal.

Los sistemas de clasificación de la vegetación basados en sensores remotos discutidos, utilizan el NDVI como parte importante de sus algoritmos de discriminación, lo cual los hace poco atractivos para ecosistemas con mucho follaje, por el problema de saturación de la banda del rojo (R), la cual está incluida en el NDVI junto con la banda del infrarrojo cercano (IRC). Paz et al. (2007) discuten el diseño de índices de vegetación y muestran que las hipótesis implícitas del NDVI son de las peores para capturar la dinámica del crecimiento de la vegetación.

Para poder definir los alcances y limitaciones de posibles sistemas de clasificación de la vegetación usando información espectral multi–angular, uno de los objetivos principales de esta investigación fue evaluar el potencial de discriminación de diferentes tipos de vegetación, diferentes tamaños, tipos de hojas y distribución espacial y angular, mediante la información obtenida por mediciones angulares de la reflectancia.

 

INFORMACIÓN ESPECTRAL MULTI–ANGULAR

Los objetos terrestres generalmente reflejan la energía solar en forma anisotrópica, de tal forma que esta señal angular puede ser usada para distinguir o clasificar los diferentes objetos en una escena con un satélite multi–angular (Diner et al., 1999; Asner, 2000). La caracterización de las mezclas suelo–vegetación ha sido estudiada (Kimes, 1983; Ranson et al., 1985) y es la base para que la información espectral multi–angular puede ser usada para una mejor discriminación y caracterización de los cultivos y ecosistemas.

Modelos de la BRDF

Una función utilizada para caracterizar la variación de las reflectancias con la geometría sol–sensor es la BRDF (Nicodemus et al., 1977) o función de distribución de la reflectancia bidireccional, que caracteriza las reflectancias en función de los ángulos (cenitales y acimutales) de visión y de iluminación. La modelación de la BRDF ha tenido un enfoque empírico (Minnaert, 1941; Hapke, 1981; Walthall et al., 1985) o basado en la teoría de transferencia radiativa de medios turbios u homogéneos (Ross, 1981; Goel, 1988) o basado en óptica geométrica (Li y Strahler, 1985, 1992; Strahler y Jupp, 1991).

Los modelos de la BRDF más usados están basados en funciones matemáticas que representan las componentes esenciales del proceso de transferencia de radiación (kernels) que combinan modelos de medios turbios y de óptica geométrica (Wanner et al., 1995). Estos modelos semi–empíricos de kernels pueden ser lineales (Roujean et al; 1992; Lucht y Roujean, 2000) o no lineales (Staylor y Suttles, 1986; Rahman et al., 1993). Los modelos de la BRDF requieren al menos tres parámetros, lo que obliga a contar con información multi–angular con mínimas variaciones de las reflectancias de las superficies a caracterizar.

Pinty et al. (2002) y Gobron et al. (2002) usaron el modelo no lineal de la BRDF de Rahman et al. (1993) para proponer que k (parámetro de forma de la función de Minnaert modificada) de este modelo puede usarse para discriminar tipos de vegetación con fondos muy reflectivos (nieve o desierto). Así, la forma de la BRDF puede ser caracterizada en función de ser cóncava (k>1), convexa (k<1) o lineal (k=1). Este enfoque ha sido usado por Nolin (2004).

Bolaños et al. (2007) proponen una nueva familia de modelos de la BRDF que simplifican el problema de modelación, ya que sólo requieren un parámetro para su caracterización.

Clasificación de la vegetación con información espectral multi–angular

Se ha argumentado que la información espectral multi–angular mejora las clasificaciones de la vegetación (Bicheron et al., 1997; Asner, 2000). El uso de información multi–angular en una sola banda espectral en la clasificación de la vegetación, da resultados similares al uso de información de múltiples bandas espectrales con un solo ángulo de visión; aunque la discriminación de clases de vegetación está fuertemente influenciada por las diferencias a nivel espectral y, en menor medida, a nivel angular (Barnsley et al., 1997).

Zhang et al. (2002a y b) argumentan que la clasificación de la vegetación debe basarse en principios de consistencia entre las propiedades radiativas de las clases y sus firmas espectrales distintivas, manteniendo fija la resolución espacial (problema de mezclas). Estos autores utilizan espacios espectrales (R–IRC), y muestran que las clases de vegetación (radiativamente diferentes) están diferenciadas en función de su localización en el espacio espectral y tres métricas: inclinación, longitud e intercepción de los patrones cuasi–lineales en ese espacio. Asi, estos autores favorecen la hipótesis de firmas espectrales multi–angulares únicas para las clases de vegetación.

 

EXPERIMENTOS CON SISTEMAS ARBOLADOS

Para analizar si la información espectral multi–angular puede discriminar los diferentes tipos de vegetación, particularmente los arbolados, se diseñó un experimento tipo maqueta (superficie de suelo con vegetación de altura o edad pequeña, como una maqueta arquitectónica) en el Colegio de Postgraduados, Campus Montecillo, Estado de México, en septiembre de 2006. En el experimento se usaron cinco especies forestales, con plantas relativamente homogéneas (Cuadro 1). En cada especie forestal se midió BRDF usando tres configuraciones geométricas (Figura 1).

Para cada especie y configuración geométrica se usaron dos suelos (claro y oscuro) como fondo y se realizaron mediciones durante tres o cuatro ángulos cenitales (y acimutales) de iluminación. En todos los casos se usó un equipo de BRDF diseñado ex profeso para medir ángulos de visión de 0 a 50°, con incrementos de 10°, tanto en la dirección de iluminación solar como en la contraria, obteniéndose 12 mediciones por efectos de visión (los ángulos de nadir se midieron dos veces: iluminación y sombreado). El equipo de BRDF mantiene relativamente constante el área de observación reduciendo la altura de mención (para el ángulo de visión del radiómetro) en función del ángulo de visión. El radiómetro usado fue el modelo FR Jr de ASD^MR, que es un radiómetro hiperespectral con rango de medición de 350 a 2500 nm y resolución de proceso de 1 nm. En todas las mediciones se tomaron fotografías de la cobertura de la vegetación con una cámara digital modelo Cybershot DSC–V1 de Sony^TM. Con estas fotografías se obtuvo la cobertura aérea (COB) de la vegetación, con un proceso de clasificación supervisada usando puntos de control de los objetos suelo o follaje en la fotografía digital.

Además se efectuaron mediciones para el caso de sólo suelo (COB = 0), para tener cuatro coberturas diferentes para cada especie analizada. El uso de suelo como fondo de la vegetación fue para ver los efectos de fondos con reflectancias contrastantes (propiedades ópticas del material debajo de las plantas, que define las interacciones múltiples de los fotones que chocan contra éste y rebotan para interceptar de nuevo al follaje de las plantas). En los sistemas naturales el fondo de los sistemas arbolados (bosques) varía estacionalmente (Miller et al., 1997) y tiene su propia función BRDF (Peltonieri et al., 2005).

Las reflectancias hiperespectrales se convirtieron a bandas espectrales del sensor TM (Landsat 5) usando las funciones de respuesta espectrales correspondientes, después de eliminar valores con problemas de vapor de agua atmosférico detectables fácilmente por el rompimiento del patrón espectral (caso de la banda de 1000 nm) o por el ruido asociado a la señal del radiómetro. En el Cuadro 2 se muestra las bandas espectrales usadas.

Modelo de la BRDF usado

El modelo de la BRDF utiliza un solo parámetro (Bolaños et al., 2007), basado en la siguiente convención para el ángulo de visión:

donde, la geometría de iluminación está definida por el ángulo cenital solar (θs) y el ángulo acimutal solar (øs) mientras que la geometría de visión está definida por el ángulo cenital de visión (θv) y el ángulo acimutal de visión (øν).

Dado que durante las mediciones de reflectancia hubo muchas variaciones en las componentes de la radiación solar (difusa y directa) y que el mayor contraste en las reflectancias angulares se da en la dirección contraria a la de la iluminación o θvpositivo (Barnsley et al., 1997), en este trabajo sólo se analiza el modelo de la BRDF para ese caso (Bolaños et al., 2007):

donde, Rn es una reflectancia normalizada, multiplicada por cos(X), y X es una variable angular de posición que reduce la complejidad de la BRDF al usar simetrías angulares.

 

ANÁLISIS DE LOS RESULTADOS

En la Figura 2 se muestra el ajuste del modelo de la BRDF para el caso del fresno, COB = 91.76%, θs = 14.4° y en la Figura 3 el caso del suelo desnudo claro, θs = 56.1°. En ambos casos la regresión lineal fue forzada a tener una ordenada al origen de 90°. Los ajustes del modelo analizado de la BRDF fueron adecuados (R²>0.98) para la gran mayoría de los casos, mostrando que la BRDF puede ser aproximada en forma simple al usar simetrías angulares. Resultados similares fueron encontrados para el caso de pastizales (Bolaños et al., 2007).

Dado que las pendientes b de la función BRDF de las bandas espectrales están en función del ángulo cenital (y acimutal) solar, una forma de compactar (generalizar) la función BRDF es analizar la relación entre estas pendientes, dejando el ángulo cenital solar implícito. En este estudio las relaciones entre las b de la banda del R y del IRC fueron consideradas, aunque las relaciones entre el resto de las bandas presentan patrones similares al de éstas.

Para analizar las BRDF, en forma independiente del ángulo cenital solar, se usó la relación entre las pendientes del modelo BRDF de las bandas del R e IRC:

El uso de sólo las bandas del R e IRC es debido a que este par de bandas espectrales se presenta en prácticamente todos los sensores remotos a bordo de plataformas espaciales y en el alto contraste entre estas dos bandas para el caso de la vegetación (Tucker, 1979).

En la Figura 4 se muestra la relación (3) para el caso del fresno. Aunque los ajustes resultaron razonables (R²>0.9) para los casos analizados, el modelo lineal propuesto puede considerarse como una consecuencia del modelo de la BRDF generalizado (Bolaños et al., 2007).

En la Figura 5 se muestra la variación de las pendientes B, relación (3), en función de la cobertura aérea, para las especies analizadas y el caso del suelo claro. En la Figura 6 se muestra esta misma información para el caso del suelo oscuro.

En las Figuras 5 y 6 se usó un polinomio de segundo grado para ajustar la relación entre B y COB, sólo para mostrar patrones. El uso de la pendiente B, relación (3), está de acuerdo en que este parámetro captura las diferencias producto de cambios en la iluminación solar. El análisis de las Figuras 5 y 6 muestra que los patrones espectrales multi–angulares de los sistemas arbolados analizados son separables (se pueden clasificar correctamente), si queda fijo el tipo de suelo como fondo de la vegetación. Esta evidencia experimental soporta la idea de utilizar información espectral multi–angular para una mejor discriminación de los tipos de vegetación naturales, particularmente los arbolados.

El caso de que el fondo de la vegetación varíe (caso normal), dejando fija la cantidad de vegetación, hay zonas de traslape de información (confusión en términos de clasificaciones), que no permiten separar las señales usando la información espectral multi–angular (Figuras 7 y 8). Esto se da particularmente para el caso de coberturas aéreas baja y no se presenta si la cobertura aérea aumenta.

En las Figuras 7 y 8 se muestra que, independientemente del tipo de fondo (suelo) de la vegetación, los patrones espectrales angulares compactados de las especies analizadas muestran dos formas: cóncavos y convexos, definiendo así un camino de análisis de las señales espectrales multi–angulares que puede usarse en forma similar a la propuesta por Pinty et al. (2002) y Gobron et al. (2002).

Para tener una perspectiva del problema de clasificación de sistemas arbolados usando un enfoque clásico espectral, en la Figura 9 se muestra el espacio del R–IRC para las especies analizadas (sólo el caso de ángulos de visión contrarios al de iluminación), teniendo como fondo un suelo claro y usando todas las coberturas aéreas analizadas por especie. Se observa que hay importantes traslapes de información, por lo que la clasificación de especies es muy problemática.

En la Figura 10 se presenta la información de la Figura 9, pero utilizando reflectancias normalizadas, asociadas al modelo de la BRDF discutido. Se observa que los patrones espectrales son más separables usando algoritmos de discriminación basados en análisis angulares lineales, por lo que puede inferirse que con el formateo angular, normalización, de los espacios espectrales puede obtenerse una mejoría en las clasificaciones de la vegetación, cuando menos en los extremos.

 

CONCLUSIONES

Partiendo de una discusión del uso de la información espectral multi–angular para clasificar y caracterizar a la vegetación, se han analizado las diferentes propuestas teóricas que soportan esta estrategia, así como los modelos de la función de distribución bidireccional de las reflectancias (BRDF).

Usando un modelo simplificado de la BRDF, dependiente de un solo parámetro, se han analizado los resultados de un experimento de sistemas arbolados en maquetas, donde se realizaron mediciones de las reflectancias variando los ángulos cenitales de visión y de iluminación. El modelo ajustado de la BRDF resultó adecuado para todas las bandas del sensor TM (Landsat 5), determinadas por la convolución de las mediciones hiperespectrales a las de este sensor usando sus funciones de respuesta.

La información del modelo de la BRDF fue compactada analizando las pendientes para las bandas del R e IRC, con el ángulo cenital solar implícito, para poder analizar los patrones de esta aproximación en función de la cobertura aérea de las cinco especies forestales. Los resultados mostraron que, dejando el suelo del fondo de la vegetación fijo, las señales espectrales multi–angulares son separables, por lo que la estrategia propuesta podría usarse para clasificar sistemas arbolados. En el caso de que el suelo varíe (sus propiedades ópticas), entonces se presenta confusión en la información espectral multi–angular, para el caso de coberturas aéreas bajas, reduciéndose este problema a medida que la cobertura aérea se incrementa.

Finalmente, es necesario enfatizar que los análisis realizados y conclusiones obtenidas están restringidos al diseño experimental utilizado, por lo que se requiere su revisión con datos de mediciones realizadas sobre sistemas naturales. Este último esquema está en progreso por el grupo de investigación usando un helicóptero de control remoto equipado con un radiómetro y cámara de video digital, controlando el resto de los factores que inciden en los valores de las reflectancias.

 

LITERATURA CITADA

Alvarez, R., R. Bonifaz, R. S. Lunetta, C. García, G. Gómez, R. Castro, A. Bernal, and A. L. Cabrera. 2003. Multitemporal landcover classification of Mexico using Landsat MSS imagery. Int. J. Remote Sensing 24: 2501–2514.        [ Links ]

Asner, G. P. 2000. Contributions of multi–view angle remote sensing to land–surface and biogeochemical research. Remote Sensing Rev. 18: 137–162.        [ Links ]

Asrar, G., R. Myneni, and B. Choudhury. 1992. Spatial heterogeneity in vegetation canopies and remote sensing of absorbed photosynthetically active radiation: A modeling study. Remote Sensing of Environ. 41: 85–103.        [ Links ]

Barnsley, M. J., D. Allison, and P. Lewis. 1997. On the information content of multiple view angle (MVA) images. Int. J. Remote Sensing 18: 1937–1960.        [ Links ]

Bicheron, P., M. Leroy, and F. M. Breon. 1997. Enhanced discrimination of boreal forest covers with direccional reflectances from the airborne polarization and directionality of earth reflectances (POLDER) instrument, J. Geophysics Res. 102: 29517–29531.        [ Links ]

Bolaños, M., F. Paz, E. Palacios, E. Mejía, y A. Huete. 2007. Modelación de los efectos de la geometría sol–sensor en la reflectancia de la vegetación. Agrociencia 41: 527–537.        [ Links ]

Carpenter, G. A., S. Gomal, S. Macomber, S. Martens, C. E. Woodlock, and J. Franklin. 1979. A neural network method for efficient vegetation mapping. Remote Sensing Environ. 70: 326–338.        [ Links ]

Cayuela, L., J. D. Golicher, J. Salas Rey, and J. M. Rey Benayas. 2006. Classification of a complex landscape using Dempster–Schafer theory of evidence. Int. J. Remote Sensing 27: 1951–1971.        [ Links ]

Daniels, A. E. 2006. Incorporating domain knowledge and spatial relationships into land cover classifications: a rule–based approach. Int. J. Remote Sensing 27: 2949–2975.        [ Links ]

Diner, D. J., G. P. Asner, R. Davies, Y. Knyazhikin, J. P. Muller, A. W. Nolin, B. Pinty, C. B. Schaaf, and J. Stroeve. 1999. New directions in earth observing: scientific applications of multiangle remote sensing. Bull. Metereol. Soc. 80: 2209–2228.        [ Links ]

Foody, G. M. 1996. Approaches for the production and evaluation of fuzzy land cover classification from remotely–sensed data. Int. J. Remote Sensing 17: 1317–1340.        [ Links ]

Friedl, M., and C. Brodley. 1997. Decision tree classification of land cover from remotely sensed data. Remote Sensing of Environ. 1: 399–409.        [ Links ]

Gobron, N., B. Pinty, M. M. Verstraete, J. L. Widlowski, and D. J. Diner. 2002. Uniqueness of multiangular measurements–Part II: joint retrieval of vegetation structure and photosynthetic activity from MISR. IEEE Transactions on Geoscience and Remote Sensing 40: 1574–1592.        [ Links ]

Hapke, B. 1981. Bidirectional reflectance spectroscopy 1. Theory J. Geophysics Res. 86: 3039–3054.        [ Links ]

Justice, C. O., J. R. G. Townshend, B. N. Holben, and C. J. Tucker. 1985. Analysis of the phenology of global vegetation using meteorological satellite data. Int. J. emote Sensing 6: 1271–1381.        [ Links ]

Kimes, D. S. 1983. Dynamics of directional reflectance factor distributions for vegetation canopies. Applied Optics 22: 1364–1372.        [ Links ]

Knyazikhin, Y., J. V. Martonchik, R. B. Myneni, D. J. Diner, and S. W. Running. 1998a. Synergistic algorithm for estimating vegetation canopy leaf area index and fraction of absorbed photosynthetically active radiation from MODIS and MISR data. Theory J. Gephysics Res. 103: 32257–32275.        [ Links ]

Knyazikhin, Y., J. V. Martonchik, D. J. Diner, R. B. Myneni, M. Verstraete, B. Pinty, and N. Gobron. 1998b. Estimation of vegetation canopy leaf area index and fraction of absorbed photosyntheticallly active radiation from atmosphere–correct MISR data. Theory. J. Gephysics Res. 103: 32239–32256.        [ Links ]

Li, X, and A. H. Strahler. 1985. Geometrical–optical modeling of a conifer forest canopy. IEEE Transactions on Geoscience and Remote Sensing 23: 705–721.        [ Links ]

Li, X, and A. H. Strahler. 1992. Geometrical–optical bidirectional reflectance modeling of the discrete crown vegetation canopy: effect of crown shape and mutual shadowing IEEE Transactions on Geoscience and Remote Sensing 30: 276–292.        [ Links ]

Lucht, W., and J. L. Roujean. 2000. Considerations in the parametric modeling of BRDF and albedo from multiangle satellite sensor observations. Remote Sensing Rev. 18: 343–380.        [ Links ]

Miller, J. R., H. P. White, J. M. Chen, D. R. Peddle, G. McDermid, R. A. Fournier, P. Shepherd, I. Rubinstein, J. Freemantle, R. Soffer, and E. LeDrew. 1997. Seasonal change in understory reflectance of boreal forests and influence on canopy vegetation indices. J. Geophysics Res. 102: 29475–29482.        [ Links ]

Minnaert, M. 1941. The reciprocity principle in lunar photometry. Astrophysics J. 93: 403–410.         [ Links ]

Myneni, R. B., and D. L Williams. 1994. On the relationship between FAPAR and NDVI, Remote Sensing Environ. 49: 200–211.         [ Links ]

Nemani, R., and S. W. Running. 1997. Land cover characterization using multitemporal red, near–ir and thermal–ir data from NOAA/AVHRR. Ecol. Applications 7: 79–90.         [ Links ]

Nicodemus, F. E., J. C. Richmond, J. J. Hsia, I. W. Ginsberg, and T. Limperis. 1977. Geometrical considerations and nomenclature for reflectance. NBS Monogr. No. 160, National Bureau of Standards, U.S. Department of Commerce. 52 p.         [ Links ]

Nolin, A. W. 2004. Towards retrieval of forest cover density over snow from the multi–angle imaging spectroradiometer (MISR). Hydrological Processes 18: 3623–3636.         [ Links ]

Paz, F., E. Palacios, M. Bolaños, L. A. Palacios, M. Martínez, E. Mejía y A. Huete. 2007. Diseño de un índice espectral de la vegetación: NDVIcp. Agrociencia 41: 539–554.         [ Links ]

Peltonieri, J. I., S. Kaasalainen, J. Näränen, M. Rautiainen, P. Stenberg, H. Smolander, S. Smolander, and P. Voipio. 2005. BRDF measurement of understory vegetation in pine forests. Dwarf shrubs, lichen, and moss. Remote Sensing Environ. 94: 343–354.         [ Links ]

Pinty B., J. L. Widlowski, N. Gobron, M. M. Verstraete, and D. J. Diner. 2002. Uniqueness of multiangular measurements – Part I: an indicator of subpixel surface heterogeneity form MISR. IEEE Transactions on Geoscience and Remote Sensing 40: 1560–1573.         [ Links ]

Rahman H., B. Pinty, and M. M. Verstraete. 1993. Coupled surface–atmosphere reflectance (CSAR) model 2. Semiempirical surfaces model usable with NOAA advanced very high resolution radiometer data. J. Geophysics Res. 98: 20791–20801.         [ Links ]

Ranson, K. J., C. S. T. Daughtry, L. L. Biehl, and M. E. Bauer. 1985. Sun–view effects on reflectance factors of corn canopies. Remote Sensing Environ. 18: 147–161.         [ Links ]

Ross, J. K. 1981. The Radiation Regime and Architecture of Plant Stands. Dr. W. Junk Publishers, Norwell, Mass. 392 p.         [ Links ]

Roujeau, J. L., M. Leroy, and P. Y. Deschamps. 1992. A bidirectional reflectance model of the earth's surface for the correction of remote sensing data. J. Geophysics Res. 97: 20455–20468.         [ Links ]

Running, S. W., T. R. Loveland, L. L. Pierce, R. R. Nemani, and E. R. Hunt. 1995. A remote sensing based vegetation classification logic for global land cover analysis. Remote Sensing Environ. 51: 39–48.         [ Links ]

Seto, K. C., and R. K. Kaufmann. 2005. Using logit models to classify land cover and land–cover change from Landsat Thematic Mapper. Int. J. Remote Sensing 26: 563–577.         [ Links ]

Staylor, W. F., and J. T Suttles. 1986. Reflection and emission models for deserts derived from NIMBUS 7 ERB scanner measurements. J.0 Climate Application Meteorol. 25: 196–202.         [ Links ]

Strahler, A. H., and D. L. B. Jupp. 1991. Modeling bidirectional reflectance of forests and woodlands using Boolean models and geometric optics. Remote Sensing Environ. 34: 153–166.         [ Links ]

Tucker, C. J. 1979. Red and photographic infrared linear combination for monitoring vegetation. Remote Sensing of Environment 8: 127–150.         [ Links ]

Walthall, C. L., J. M. Norman, J. M. Welles, G. Campbell, and B. Blad. 1985. Simple equation to approximate the bidirectional reflectance from vegetation canopies and bare soil surfaces. Applied Optics 24: 383–387.         [ Links ]

Wanner W., X. Li, and A. H. Strahler. 1995. On the derivation of kernels for kernel–driven models of bidirectional reflectance. J. Geophysics Res. 100: 21077–21089.         [ Links ]

Yun, O., and M. Jianwen. 2006. Land cover classification based on tolerant rough set. Inter. J. Remote Sensing 27: 3041–3047.        [ Links ]

Zhang, Y., Y. Tian, R. B. Myneni, and Y. Knyazikhin. 2002a. Required consistency between biome definitions and signatures with the physics of remote sensing. I: empirical arguments. Remote Sensing Environ. 80: 418–434.        [ Links ]

Zhang, Y., N. Shabanov, Y. Knyazikhin, and R. B. Myneni. 2002b. Required consistency between biome definitions and signatures with the physics of remote sensing. II: theoretical arguments. Remote Sensing Environ. 80: 435–446.        [ Links ]