Source: https://www.youscribe.com/catalogue/documents/autres/efecto-de-diferentes-metodos-de-clasificacion-de-imagenes-satelitales-2012893
Timestamp: 2018-08-19 22:58:40+00:00

Document:
EFECTO DE DIFERENTES MÉTODOS DE CLASIFICACIÓN DE IMÁGENES SATELITALES SOBRE ÍNDICES DE PAISAJE (Effect of different satellite image classification methods on landscape indices) - M. - Autres
EFECTO DE DIFERENTES MÉTODOS DE CLASIFICACIÓN DE IMÁGENES SATELITALES SOBRE ÍNDICES DE PAISAJE (Effect of different satellite image classification methods on landscape indices)
El estudio de patrones y procesos ecológicos en la escala de paisaje se basa frecuentemente en el uso de índices calculados a partir de clasificaciones de cobertura. Las clasificaciones dependen de la resolución espacial y espectral de las imágenes satelitales y pueden realizarse a traves de clasificadores “por píxeles” o “por regiones”. El primero solo considera la reflectancia de cada píxel para asignarlo a una clase determinada, y el segundo adicionalmente considera la vecindad del píxel mediante un proceso denominado segmentación, que agrupa píxeles vecinos en regiones homogéneas, previo a la clasificación. El objetivo del presente trabajo es comparar el método de clasificación por regiones con diferentes alternativas del método de clasificación por píxel sobre los índices de caracterización del paisaje en la región pampeana Argentina. Se realizaron clasificaciones no supervisadas utilizando técnicas por píxeles y por regiones, a partir de Indice de Vegetación Normalizado de cuatro imágenes LANDSAT 5 TM.
Se realizaron cuatro clasificaciones por píxeles (dos tamaños de píxel, 30x30m y 134x134m, y el filtro post-clasificación por mayoría, con o sin filtro). En la clasificación por regiones, se utilizó el algoritmo de “crecimiento de regiones” para el proceso de segmentación. Cada clasificación se disectó con una grilla de 8.3x8.3 km. En cada unidad se cuantificaron cinco índices de paisaje (número de parches, tamaño medio de parches, borde total, borde medio de parches, índice medio de forma) y se compararon entre clasificaciones. Para todos los índices las mayores diferencias se encontraron entre la clasificación por regiones y la clasificación por píxel de 30m sin filtro. Los valores de los índices para estas dos clasificaciones fueron, respectivamente: número de parches: 294 y 5688 parches
tamaño medio de parches: 24 y 1.2ha
borde total: 733769 y 2003763m
borde medio total: 2504 y 354m
índice medio de forma: 1.69 y 1.3.
La clasificación por regiones produce índices que corresponden mejor con la realidad debido a la incorporación de otros criterios de agrupación como el patrón de forma y la distancia a unidades vecinas. Por ejemplo el tamaño medio de parche calculado por este métod (promedio =24ha) fue mas cercano al tamaño real de los parches agrícolas de esta región, el cual varía entre 20 y 100ha.
The study of landscape-scale ecological patterns and processes is frequently based on the use of indices calculated from land cover classifications. Classifications, in turn, depend on satellite-imagery spatial and spectral resolution and can be created by “pixel-based” or “object (region)- based” classifiers. The former only considers each pixel´s reflectance, while the latter additionally considers pixel neighborhoods, through a process called “segmentation”, which groups neighbouring pixels in homogeneous regions. The objective of this work was to compare the object-based method with different alternatives of a pixel–based method on landscape indices, in the Pampas region of Argentina. Five regional, non-supervised classifications were made with pixel- and object-based methods, using Normalized Difference Vegetation Index (NDVI) values from four LANDSAT 5TM images: 4 pixel-based, varying pixel size (30x30m and 134x134m), with or without majority post-classification filters, and one objectbased. In the latter, we used the “region growing” algorithm for the segmetation process. Each classification was dissected by an 8.3km grid and for each 69km2-unit, we calculated the following indices: number of patches, mean patch size, total edge, mean patch edge, and mean shape index. These were compared between classifications. We found greatest differences between the object-based classification and the non-filtered 30m-pixel classification. The index values for these two methods were, respectively: number of patches: 294 and 5688
mean patch size: 24 and 1.2ha
total edge: 733769 and 2003463m
mean patch edge: 2504 and 354m
mean shape index: 1.69 and 1.3.
The object-based method produced index values that correspond better to reality, due to the incorporation of other classfication criteria, such as real-world shape patterns and neighbouring pixel´s spectral values. For example, the mean patch size calculated from the object-based method (mean = 24ha), was closest to true agricultural patch size in this region, which varies between 20 and 100ha.
Source : Revista de Teledetección 1988-8740 (2009) Num. 32
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-13
Efecto de diferentes métodos de clasificación
de imágenes satelitales sobre índices de paisaje
Effect of different satellite image classification
methods on landscape indices
M. C. Conde, S. B. Perelman y A. Cerezo
mconde@agro.uba.ar
Departamento de Métodos cuantitativos y sistema de información. Facultad de Agronomía-UBA,
Av. San Martín 4453,Ciudad de Buenos Aires, Argentina
Recibido el 26 de mayo de 2009, aceptado el 23 de noviembre de 2009
The study of landscape-scale ecological patternsEl estudio de patrones y procesos ecológicos en
and processes is frequently based on the use ofla escala de paisaje se basa frecuentemente en el
indices calculated from land cover classifica-uso de índices calculados a partir de clasifica-
tions. Classifications, in turn, depend on sate-ciones de cobertura. Las clasificaciones depen-
llite-imagery spatial and spectral resolution andden de la resolución espacial y espectral de las
can be created by “pixel-based” or “object (re-imágenes satelitales y pueden realizarse a traves
gion)-based” classifiers. The former only consi-de clasificadores “por píxeles” o “por regiones”.
ders each pixel´s reflectance, while the latterEl primero solo considera la reflectancia de cada
additionally considers pixel neighborhoods,píxel para asignarlo a una clase determinada, y el
through a process called “segmentation”, whichsegundo adicionalmente considera la vecindad
groups neighbouring pixels in homogeneous re-del píxel mediante un proceso denominado seg-
gions. The objective of this work was to com-mentación, que agrupa píxeles vecinos en regio-
pare the object-based method with differentnes homogéneas, previo a la clasificación. El
alternatives of a pixel–based method on lands-objetivo del presente trabajo es comparar el mé-
cape indices, in the Pampas region of Argentina.todo de clasificación por regiones con diferen-
Five regional, non-supervised classificationstes alternativas del método de clasificación por
were made with pixel- and object-based me-píxel sobre los índices de caracterización del pai-
thods, using Normalized Difference Vegetationsaje en la región pampeana Argentina.
Index (NDVI) values from four LANDSAT 5TMSe realizaron clasificaciones no supervisadas
images: 4 pixel-based, varying pixel sizeutilizando técnicas por píxeles y por regiones, a
(30x30m and 134x134m), with or without majo-partir de Indice de Vegetación Normalizado de
rity post-classification filters, and one object-cuatro imágenes LANDSAT 5 TM. Se realiza-
based. In the latter, we used the “regionron cuatro clasificaciones por píxeles (dos ta-
growing” algorithm for the segmetation process.maños de píxel, 30x30m y 134x134m, y el filtro
Each classification was dissected by an 8.3kmpost-clasificación por mayoría, con o sin filtro).
grid and for each 69km2-unit, we calculated theEn la clasificación por regiones, se utilizó el al-
following indices: number of patches, meangoritmo de “crecimiento de regiones” para el
patch size, total edge, mean patch edge, andproceso de segmentación. Cada clasificación se
mean shape index. These were compared bet-disectó con una grilla de 8.3x8.3 km. En cada
ween classifications. We found greatest diffe-unidad se cuantificaron cinco índices de paisaje
rences between the object-based classification(número de parches, tamaño medio de parches,
and the non-filtered 30m-pixel classification.borde total, borde medio de parches, índice
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-13 5M. C. Conde , S. B. Perelman y A. Cerezo
medio de forma) y se compararon entre clasifi- The index values for these two methods were,
caciones. Para todos los índices las mayores di- respectively: number of patches: 294 and 5688;
ferencias se encontraron entre la clasificación por mean patch size: 24 and 1.2ha; total edge:
regiones y la clasificación por píxel de 30m sin 733769 and 2003463m; mean patch edge: 2504
filtro. Los valores de los índices para estas dos and 354m; mean shape index: 1.69 and 1.3.
clasificaciones fueron, respectivamente: número The object-based method produced index va-
de parches: 294 y 5688 parches; tamaño medio lues that correspond better to reality, due to the
de parches: 24 y 1.2ha; borde total: 733769 y incorporation of other classfication criteria, such
2003763m; borde medio total: 2504 y 354m; ín- as real-world shape patterns and neighbouring
dice medio de forma: 1.69 y 1.3. pixel´s spectral values. For example, the mean
La clasificación por regiones produce índices patch size calculated from the object-based me-
que corresponden mejor con la realidad debido a thod (mean = 24ha), was closest to true agricul-
la incorporación de otros criterios de agrupación tural patch size in this region, which varies
como el patrón de forma y la distancia a unidades between 20 and 100ha.
vecinas. Por ejemplo el tamaño medio de parche
calculado por este método (promedio =24ha) fue
mas cercano al tamaño real de los parches agrí-
colas de esta región, el cual varía entre 20 y
PALABRAS CLAVE: índices de paisaje, clasi- KEYWORDS: landscape indices, satellite clas-
ficación satelital, segmentación, resolución es- sification, segmentation, spatial resolution
INTRODUCCIÓN comprender y gestionar su funcionamiento (Burel &
Baudry 2002). En la última década, el uso de imá-
La caracterización de un paisaje se relaciona di- genes satelitales y sistemas de información geogra-
rectamente con el conocimiento de su composición fica (SIG) se ha convertido en una herramienta
y de la configuración espacial de los elementos que fundamental para la cuantificación de estos índices
lo integran, aspectos relativamente independientes (Chuvieco 1999, Ares et al. 2001, Fuller 2001, Li &
que en combinación afectan los procesos ecológicos Wu 2004). Algunos autores encontraron diferente
y la supervivencia de las especies (Turner et al. sensibilidad entre distintos índices de paisaje a la ca-
2001, McGarigal & Marks 1995). La composición lidad de la información espacial y espectral sobre la
se refiere a las características asociadas con la pre- cual se realiza el análisis (Saura 2001, Saura 2004,
sencia y cantidad de cada tipo de parche dentro del Frohn & Hao 2006). Los índices pueden tener erro-
paisaje, sin considerar de manera explícita su distri- res derivados de la resolución espacial (tamaño del
bución espacial, mientras que la configuración des- grano de la información), de la extensión del paisaje
cribe la distribución física o el carácter espacial de en el estudio (área abarcada) y/o de la selección de
los parches dentro del paisaje (Fahrig 2003, Forman categorias en la clasificación, también denominada
& Godron 1986). Los parches son las unidades mí- resolución temática (Turner, Gardner & O`Neil
nimas de observación de un paisaje de acuerdo con 2001). Esto hace suponer que los índices de paisaje
la escala de detección y observación, y se caracteri- podrían ser sensibles a la metodología de clasifica-
zan por constituir una superficie homogénea (Lang ción y en particular al tamaño de las unidades que
& Blaschke 2009). sirven de base en la clasificación, es decir, que se
La cuantificación del patrón del paisaje a través de trate de píxeles aislados o de segmentos (polígonos
índices numéricos es muy utilizada en ecología para o regiones de píxeles colindantes previamente agru-
caracterizar y comparar los estados de fragmenta- pados).
ción y cobertura de la vegetación y sus cambios en La clasificación de las imágenes de la superficie te-
el tiempo. En los paisajes agrícolas, la cuantifica- rrestre se puede realizar a través de la fotointerpre-
ción del impacto que produce la expansión y la in- tación o de análisis cuantitativos basados en
tensificación de la agricultura, es fundamental para computación intensiva (Quattrochi & Pelletier
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-136Efecto de diferentes métodos de clasificación de imágenes satelitales sobre índices de paisaje
1991). Dentro de este último grupo los métodos tra- de alteración entre paisajes, o del mismo paisaje a
dicionalmente más utilizados de análisis de la infor- través del tiempo, es muy importante la consisten-
mación satelital son los de clasificación por píxel cia y exactitud de los índices. Además se requiere
(Blaschke & Strobl 2001, Guershman et al. 2003), que los índices presenten gran sensibilidad a los
los cuales sólo consideran la reflectancia de cada cambios en la composición y configuración del pai-
píxel para asignarlo a una clase determinada. Con saje, de manera de poder detectar señales tempranas
esta metodología suele ocurrir que dentro de un de degradación del sistema. El objetivo del presente
grupo de píxeles clasificados en una dada categoría trabajo es comparar el método de clasificación por
se encuentran píxeles aislados asignados a otras cla- regiones con diferentes alternativas de método de
ses. Este fenómeno es conocido como “sal y pi- clasificación por píxel sobre los índices de caracte-
mienta” y se debe a que en el criterio de asignación rización del paisaje en la región pampeana Argen-
no se tiene en cuenta a la vecindad entre píxeles. tina.
Para solucionar este problema se aplican filtros post
clasificación los cuales uniformizan las clases, re- MATERIALES Y METODOS
asignando los píxeles aislados o poco representati-
vos a las clases vecinas más representativas Área de estudio
(Chuvieco 2002). Uno de los problemas que pre-
senta esta técnica es la decisión sobre los paráme- Se seleccionó una zona en la Pampa Interior, una de
tros de filtrado ya que puede obtenerse el efecto las subregiones de la Llanura Pampeana en Argen-
deseado en algunas zonas de la imagen, pero ser ex- tina (León 1991). Esta zona se encuentra en la es-
cesivo o insuficiente en otras. Por otro lado, en los cena Landsat 5 TM Path/Row 227/85 y su superficie
últimos años se desarrollaron las técnicas de clasifi- es de 248000 ha (49kmx49km) (Figura 1). La re-
cación por regiones, que en lugar de clasificar píxe- gión posee un clima templado-húmedo, las precipi-
les aislados, utilizan como unidad mínima de taciones promedio anuales oscilan entre los 1200
clasificación regiones, poniendo mayor énfasis en la mm y los 600mm. Los suelos están dominados prin-
vecindad del píxel, al agrupar píxeles vecinos en re- cipalmente por molisoles y en menor proporción por
giones más homogéneas (proceso conocido como alfisoles, y oxisoles. El uso de la tierra en esta re-
“segmentación”, Meinel & Neubert 2004). Estas gión se caracteriza por una mezcla de cultivos anua-
técnicas utilizan la información espectral de cada les invernales (trigo) y estivales (soja, girasol, maíz)
píxel y la información espacial que abarca la rela- que corresponden al 33% del área, pasturas peren-
ción entre los píxeles y sus vecinos. Estos clasifica- nes y pastizales naturales o seminaturales que co-
dores procuran simular el comportamiento de un rresponden a un 56 % del área (Hall et al. 1992,
fotointérprete, al reconocer áreas homogéneas de la Baldi et al. 2006). Asimismo se caracteriza por po-
imagen satelital basadas tanto en las propiedades es- seer cuerpos de áreas permanentes y transitorios en
pectrales como en las espaciales (Camara et al. un 12 % de su superficie.
1996, Dragut & Blaschle 2006, Desclée et al. 2006). En las últimas décadas, la expansión e intensifica-
Este tipo de clasificaciones requiere un nivel de pro- ción agrícola tuvieron un fuerte impacto sobre esta
cesamiento que resulta más complejo ya que pre- heterogeneidad, disminuyendo la cobertura y frag-
senta una mayor demanda computacional que los mentando las pasturas naturales e implantadas, tal
clasificadores por píxel (Quattrochi & Pelletier como ocurrió en otros paisajes agrícolas del mundo
1991). Asimismo podemos definir los métodos de (Tscharntke et al. 2005, Baldi et al. 2006). La su-
clasificación como supervisados y no supervisados, perficie de cultivos anuales en la pampa interior ha
dependiendo si se utiliza verdad terrestre conocida aumentado en un 86% en los últimos 20 años.
para entrenar al clasificador. El uso de clasificado- (Fuente: Secretaría de Agricultura Ganadería Pesca
res supervisados es mucho mas confiable y exacto, y Alimentación, SAGPyA)
pero requiere documentación auxiliar, verdad te-
rrestre y/o fotografía aérea para entrenar al clasifi- Procesamiento de imágenes satelitales
cador. Tanto las clasificaciones por píxeles como las
clasificaciones con previa segmentación por regio- Se utilizaron cuatro imágenes LANDSAT 5 TM fa-
nes pueden utilizar métodos supervisados o no su- cilitadas por la Comision Nacional de Actividades
pervisados. Espaciales (CONAE), libre de nubes, de distintas fe-
En los estudios ambientales comparativos del grado chas de la misma estación de crecimiento (01/11/03,
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-13 7M. C. Conde , S. B. Perelman y A. Cerezo
Figura 1. Mapa de la Región Pampeana en Argentina. Se muestra el área de estudio ubicada en la subregión Pampa In-
terior. En el extremo superior izquierdo se ubica esta región en Sud América
19/12/03, 04/01/04, 08/03/04), de manera de captu- genes, debido a que este es un buen estimador de la
rar las diferencias fenológicas entre distinto tipos de presencia y condición de la vegetación presente (Pet-
cobertura. Las imágenes seleccionadas fueron co- torelli et al. 2005). Se realizaron clasificaciones no
rregidas geométricamente, radiométricamente, y at- supervisadas y se establecieron 11 clases funda-
mosféricamente. Las correcciones geométricas se mentadas en el conocimiento previo de los tipos de
realizaron en base una imágen ortorectificada dispo- cobertura que prevalecen en la zona (como por ejem-
nible en el sitio http://glcf.umiacs.umd.edu/data/land- plo soja, maíz, girasol, trigo, avena, pastizal natural,
sat/ (Universidad de Maryland). Cada imagen se pasturas implantadas, etc), sin embargo y para los
modificó según la referencia, utilizando 50 puntos fines de este trabajo la identidad de cada parche no
de control de tierra. El error asociado a la correc- fue de interés, dado que nos interesaba evaluar las
ción geométrica fue inferior al tamaño de un píxel. diferencias en los índices de paisaje entre clasifica-
El método de ajuste fue el polinomio de grado 1, el ciones y no entre clases. Por lo tanto, los índices
remuestreo fue el de “vecino más cercano” y la pro- fueron calculados para todo los parches sin discri-
yección utilizada fue UTM 20S-WGS84. Los valo- minar entre clases. Para ambos tipos de clasifica-
res de número digital (ND) fueron llevados a valores ción se utilizó el programa de distribución gratuita
de radiancia y posteriormente a reflectancia. Para el SPRING 5.1.3, desarrollado por el INPE (Instituto
calculo de radiancia se utilizaron los valores de Nacional de Pesquisas Espaciais) (Camara et al
“gains” y “biases” publicados por Chander & Mar- 1996).
kham (2003). Se realizaron correcciones atmosféri-
cas del efecto Rayleigh con el objetivo de llegar a Clasificación por regiones
valores aproximados de reflectancia en superficie.
Tanto la segmentación como la clasificación pos-
Clasificaciones terior se realizaron con las cuatro imágenes de IVN
con una resolución de 30m. La segmentación fue
Las clasificaciones se realizaron en base al Indice realizada con el algoritmo “crecimiento de regio-
de Vegetación Normalizado (IVN) de las cuatro imá- nes”, el cual se caracteriza por generar regiones ce-
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-138Efecto de diferentes métodos de clasificación de imágenes satelitales sobre índices de paisaje
rradas (Bins et a. 1996). Este algoritmo utiliza dos Por un lado se realizaron clasificaciones con la re-
parámetros: similitud y área. La similitud es la dis- solución espacial original de 30x 30m y por otro
tancia euclidiana entre los valores de ambos píxeles lado se realizó un redimensionamiento del tamaño
y el área es la mínima superficie que puede tener una del píxel a 134 x134m de manera que cada píxel cla-
región, la cual se define mediante el número de pí- sificado se asemejara al área mínima de los seg-
xeles. Se utilizaron valores de similitud de 20 y área mentos generados para la clasificación por regiones
mínima de 20 (aproximadamente 1.8ha) (Figura 2). (1.8ha). Se utilizó el filtro post clasificación de
La clasificación de los segmentos fue realizada con SPRING, que es un filtro de “mayoría” en donde el
un clasificador no supervisado disponible en píxel central se asigna a la categoría que cuente con
SPRING denominado Histograma, el cual permite la mitad más uno de los píxeles vecinos.
definir un número fijo de clases. En esta técnica se
calcula un histograma a partir de todos los segmen- Caracterización del paisaje
tos, luego los picos y valles del histograma se utili-
zan para localizar las clases (Fu & Mui, 1981). Con las clasificaciones obtenidas se evaluaron los
patrones de paisaje. Se disectó la región en una gri-
Clasificación por píxel lla de 8.3 x 8.3 km, de manera de generar 36 unida-
des de paisajes. En cada una de ella se cuantificaron
La clasificación no supervisada por píxel fue reali- los índices de paisaje, ampliamente utilizados en es-
zada con el algoritmo K-Medias. Este método de tudios de paisaje: número de parches, tamaño
análisis de conglomerados apunta al agrupamiento medio de los parches, borde total, borde medio por
de n observaciones en un número fijo de clases, parche, índice medio de forma de los parches. La
agrupando a cada observación en la clase con la descripción de los índices se presenta en la Tabla 1.
media más cercana. Se realizaron cuatro clasifica- (Reynolds & Wu 1999,O'Neil et al. 1997, MacGari-
ciones por píxeles combinando tamaño de píxel y gal & Marks 1995, Tischendorf 2001). Para el cál-
filtros post-clasificación para corregir el efecto sal y culo de los índices de paisaje se utilizó la extensión
pimienta. Patch Analyst del ArcGIS 9.0
Figura 2. Resultado de la segmentación en una subregión del área de estudio
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-13 9M. C. Conde , S. B. Perelman y A. Cerezo
Índice de paisaje Descripción Abreviación Fórmula Unidades
Es la cantidad de parches que hay en un
Número de parches NP parches
Tamaño medio de Es las suma de las áreas de todos los par- 2TMP m o ha
parches ches , divido el nro. de parches
Es la suma de la longitud de todos los bor-
Borde total des de los segmentos que envuelven a un BT m o km
Es la suma de la longitud de todos los bor-Borde medio de
des de los segmentos que envuelven a un BMP m o km
parche, divido el nro. de parches
Caracteriza el desvío de la forma actual de
Índice medio de un parche de una forma optimizada que es
forma un circulo. MSI incrementea cuando las
formas de los parches son mas irregulares
Tabla 1. Descripción de los cinco índices de paisaje considerados
RESULTADOS Y DISCUSIÓN parches de escaso tamaño, que resultaron más abun-
dantes en las clasificaciones por píxel con mayor re-
Para todos los índices las mayores diferencias se solución y sin agregación. Sin embargo, sí hubo
encontraron entre la clasificación por regiones y la diferencia en el índice medio de forma entre la cla-
clasificación por píxel de 30x30m sin filtro (Figura sificación por regiones y aquellas basadas en píxeles
3). Los valores de los índices para estas dos clasifi- (Figura 3.E).
caciones fueron, respectivamente: número de par- La clasificación por regiones produce índices que
ches (NP): 294 y 5688 parches; tamaño medio de corresponden mejor con la realidad. Por ejemplo, el
parches (TMP): 24 y 1.2ha; borde total (BT): tamaño medio de parche calculado por este método
733769 y 2003763.92m; borde medio total (BMT): (promedio= 24ha) fue mas cercano a la superficie
2504 y 354m; índice medio de forma (IMF): 1.69 y real de los parches agrícolas de esta región, la cual
1.3. varía entre 20 y 100ha (Ghersa et. al 2002).
Por otro lado, los índices de paisaje también varia-
ron entre las clasificaciones por píxeles. Al aumen- CONCLUSIONES
tar el tamaño de píxel, el NP disminuyó, mientras
que el TMP aumentó (Figura 3.A y 3.B). Asimismo La clasificación por regiones produce patrones
el borde total disminuyó y el borde promedio por menos fragmentados que los obtenidos con las cla-
parche aumentó ante la disminución de la resolución sificaciones por píxeles, aún aplicando los filtros de
espacial (Figura 3.C y 3.D), resultados que coinci- mayoría. La clasificación por regiones mejora el
den con aquellos obtenidos por Saura (2004a, comportamiento de los índices y la consecuente ca-
2004b) y Frohn & Hao (2006). En el caso del ín- racterización y comparación de los estados de frag-
dice medio de forma y en concordancia con Frohn & mentación y cobertura de la vegetación y sus
Hao (2006), éste no varío significativamente entre cambios en el tiempo. Esta mejora se debe a la in-
clasificaciones por píxel. Consecuentemente, coin- corporación de otros criterios de agrupación como
cidimos con Saura (2004b) en que los índices NP, el patrón de forma y la distancia a unidades vecinas.
TNP, BP y BMP son muy sensibles a la presencia de
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-1310Efecto de diferentes métodos de clasificación de imágenes satelitales sobre índices de paisaje
Figura 3: El valor medio y desvío estándar de los índices de paisaje para las distintas clasificaciones. (Reg: Regiones,
134s/f: 134x134m m sin filtro, 134 c/f: 134x134m con filtro, 30s/f: 30x30m con filtro, 30 s/f: 30x30 m con filtro). Figura 3.A:
Numero de parches (NP), Figura 3.B: Tamaño medio de parches (TMP), Figura 1.C: Borde Total (BT), Figura 3.D: Borde
medio de parche (BMP), Figura 3.E: Índice de forma (IMF)
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-13 11M. C. Conde , S. B. Perelman y A. Cerezo
ject-based method. Remote Sensing of Envi-AGRADECIMIENTOS
ronment 102:1-11.
DE WIT A. J. W. & CLEVERS J. G. P. W. 2004. Ef-Agradezco a la CONAE y Universidad de Mary-
ficiency and accuracy of per-field classifica-land por proveer en forma gratuita las imágenes sa-
tion for operational crop mapping.telitales. Constanza Caride por sus valiosos
International Journal Remote Sensing 25:comentarios.
4091-4112.
DRAGUT, L., & T. BLASCHKE. 2006. AutomatedREFERENCIAS
classification of landform elements using ob-
ject-based image analysis. GeomorphologyARES, J., BERTILLER, M., H. DEL VALLE.
81:330-344.2001. Functional and structural landscape in-
FAHRIG, L. 2003. Efects of habitat fragmentationdicators of intensification, resilience and re-
on biodiversity. Annual Review of Ecology,sistance in agroecosystems in southern
Evolution, and Systematics 34:487-515.Argentina based on remotely sensed data.
FROHN, R. C., & Y. HAO. 2006. Landscape metricLandscape Ecology 16(3): 221-234.
performance in analyzing two decades of de-BALDI, G., GUERSCHMAN, J.P., PARUELO,
forestation in the Amazon Basin of Rondonia,J.M..2006. Characterizing fragmentation in
Brazil. Remote Sensing of Environmenttemperate South America grasslands. Agricul-
100:237-251.ture, Ecosystems and Environment 116:197-
FORMAN, R. T. T. & M. GODRON. 1986. Lands-208.
cape Ecology. John Wiley, New York. 619 pg.BLASCHKE, T. & STROBL, J., 2001,What’s
FU, K. S. & J. K. MUI. 1981. A survey on imagewrong with pixels? Some recent development
segmentation. Pattern Recognition 13:3-16.interfacing remote sensing and GIS. Geo-
FULLER, D. O. 2001. Forest fragmentation in Lou-BIT/GIS, 6, 12-17.
doun County, Virginia, USA evaluated withBINS, L. S., FONSECA, L. M. G., ERTHALl, G. J.
multitemporal Landsat imagery. Landscape& li, F. M. 1996. Satellite imagery segmenta-
Ecology 16(7): 627-642tion: a region growing approach. In: Simpósio
GHERSA, C. M., D. O. FERRARO, M. OMACINI,Brasileiro de Sensoriamento Remoto, 8, Sal-
M. A. MARTÍNEZ-GHERSA, S. PEREL-vador, 1996. Anais. São José dos Campos,
MAN, E. H. SATORRE, & A. SORIANO.INPE..
2002. Farm and landscape level variables asBUREL, F. & J. BAUDRY. 2002. Ecología de pai-
indicators of sustainable land-use in the Ar-saje. Conceptos, métodos y aplicaciones. Bar-
gentine Inland-Pampa. Agriculture, Ecos-celona: Mundi-Prensa.
ystems & Environment 93:279-293.CAMARA G, SOUZA RCM, FREITAS UM, GA-
GUERSCHMAN, J.P., PARUELO, J.M., DIRRIDO J 1996. SPRING: Integrating remote
BELLA, C., GIALLORENZI, M.C., PACIN,sensing and GIS by object-oriented data mo-
F. 2003b. Land cover classification in the Ar-delling" Computers & Graphics, 20: (3) 395-
gentine Pampas using multi-temporal Landsat403.
TM data. International Journal of Remote Sen-CHANDER G. & MARKHAM B. 2003. Revised
sing 24 (17), pp. 3381-3402.landsat-5 tm radiometric calibration procedu-
HALL, A. J., C. M. REBELLA, C. M. GHERSA, yres and postcalibration dynamic ranges. IEEE
J. P. CULOT. 1992. Field-crop systems of theTransactions on Geoscience and Remote Sen-
Pampas. Páginas 413-440 en Ecosystems of thesing. Vol 41, Nro. 11, 2674-2677
World: Field Crop Ecosystems. Elsevier, Ams-CHUVIECO, E. 1999. Measuring changes in lands-
terdam.cape pattern from satellite images: short-term
LANG, S. & BLASCHKE T. 2009. Análise da pai-effects of fire on spatial diversity. International
sagem com SIG. Sao Paulo. Oficina de Tex-Journal of Remote Sensing 20:2331-2346.
tos.pp. 114.CHUVIECO 2002. Teledetección ambiental. Barce-
LEÓN, R. J. C. 1991. Geographic limits of the re-lona. Ariel, pag 410.
gion, Geomorphology and geology, RegionalDESCLÉE B., P. BOGAERT, & P. DEFOURNY.
subdivisions, Floristic aspects, Description of2006. Forest change detection by statistical ob-
the vegetation. Pages 369-387 in R. T. Cou-
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-1312Efecto de diferentes métodos de clasificación de imágenes satelitales sobre índices de paisaje
pland, editor. Natural Grasslands: Introduc- Landscape Ecology. Ed. M. G. Turner & R. H.
tion and Western Hemisphere. Elsevier, Ams- Gardner, New York: Springer-Verlag. pp 51-
terdam. 76.
LI, H., & J. WU. 2004. Use and misuse of landscape REYNOLDS, J.F., WU. J. 1999. Do landscape
indices. Landscape Ecology 19:389-399. structural and functional units exist? En: Inte-
MACGARIGAL, K. & B.J. MARKS. 1995. grating hydrology, ecosystem dynamics and
FRAGSTATS- Spatial pattern analysis pro- biogeochemistry in complex landscapes (Eds.
gram for quantifying landscape structure. U.S. J.D. Tenhunene & P. Kabat). John Willey &
Forest Service Technical Report PNW 351. Sons, UK.
MEINEL, G. & NEUBERT, M. 2004. A Compari- SAURA, S., MARTINEZ-MILLÁN J. 2001. Sensi-
sion of segmentation programs for high reso- tivity of landscape pattern metrics to map spa-
lution remote sensing data. Proceedings XXth tial extent. Photogrammetric Engineering and
ISPRS Congress, Istanbul, 14 -23. July 2004, Remote Sensing 67: 1027-1036
International Archives of Photogrammetry, SAURA, S. 2004a. Effects of remote sensor spatial
Remote Sensing and Spatial Information resolution and data aggregation on selected
Sciences, Vol. XXXV-B4, p. 1097-1102 (ISSN fragmentation indices. Landscape Ecology
1682-1750). 19:197-209.
O'NEILL, R.V., HUNSAKER, C.T., JONES, K.B., SAURA, S. 2004b. Comparación de índices de frag-
RITTERS, K.H., WICKAM, J.D., mentación estimados en imágenes de satélites
SCHWARTZ, P.M., GOODMAN, I.A., con distinta resolución espacial. Revista de Te-
JACKSON, B.L. Y W.S. BAILLARGEON. ledetección 21:89-93
1997. Monitoring environmental quality of the TSCHARNTKE T., A. M. KLEIN., A. KRUESS, I.
landscape scale. Bioscience 47: 513-519. STEFFAN-DEWENTER, C. THIES,. 2005.
PETTORELLI, N., J. O. VIK, A. MYSTERUD, J.- Landscape perspectives on agricultural inten-
M. GAILLARD, C. J. TUCKER, & N. C. sification and biodiversity-ecosystem service
STENSETH. 2005. Using the satellite-derived management. Ecology Letters 8:857-874.
NDVI to assess ecological responses to envi- TISCHENDORF, L. 2001. Can landscape indices
ronmental change. Trends in Ecology & Evo- predict ecological processes consistently?
lution 20:503-510. Landscape Ecology 16: 235-254.
QUATTROCHI D.A., PELLETIER R.E. 1991. Re- TURNER, M. G., GARDNER R. H &. O`NEILL,
mote sensign for Análisis of Landscape: An In- R.V. 2001. Landscape Ecology in theory and
troduction. En Quantitative Methods in practice. Springer .
Revista de Teledetección. ISSN: 1988-8740. 2009. 32: 5-13 13

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución