V. F. Rodríguez-Galiano et al.Diferencia NVI Ratio NDVI Regresión NDVIU 1 U 2 U 3 U 1 U 2 U 3 U 1 U 2 U 3–0,159 –0,308 –0,458 –0,441 –0,927 –1,413 –0,441 –0,927 –1,413–0,140 –0,289 –0,438 –0,531 –1,018 –1,504 –0,531 –1,018 –1,504Tabla 4. Valores umbrales aplicados en las técnicas <strong>de</strong> realce <strong>de</strong> cambios.namiento por clase, con un número <strong>de</strong> píxel quevarió entre 9 y 12 (180 a 240 píxeles por clase), ypor tanto suficiente para calcular las matrices <strong>de</strong>varianza y covarianza por el clasificador. Para laevaluación <strong>de</strong> la exactitud <strong>de</strong> las clasificacionesse utilizaron un total <strong>de</strong> 145 píxeles por clase, in<strong>de</strong>pendientes<strong>de</strong> las áreas <strong>de</strong> entrenamiento.RESULTADOSY CONCLUSIONESLa forma i<strong>de</strong>al <strong>de</strong> evaluar la exactitud <strong>de</strong> unacorrección atmosférica es comparar medidas insitu <strong>de</strong> propieda<strong>de</strong>s atmosféricas y <strong>de</strong> reflectancia<strong>de</strong> superficie en el tiempo <strong>de</strong> adquisición <strong>de</strong>las imágenes con la estimación <strong>de</strong> los parámetrosresultantes <strong>de</strong> la corrección atmosférica(Moran et al., 1992). Por <strong>de</strong>sgracia, en estudiosmultitemporales en los que se utilizan imágenes<strong>de</strong> archivo esto no es siempre posible, por lo quehay que recurrir a otras formas <strong>de</strong> validación. Eneste estudio, se utilizó el error cuadrático medio(RMS) para medir la bondad <strong>de</strong> las diferentestécnicas <strong>de</strong> normalización aplicadas.El error cuadrático medio se calculó a partir<strong>de</strong> un conjunto <strong>de</strong> 15 muestras <strong>de</strong> reflectividadRMS (PIF)RMS (RCS)Banda 1 0,502 0,425Banda 2 0,457 0,667Banda 3 0,362 0,505Banda 4 0,590 1,035Banda 5 0,518 0,679Banda 7 0,023 0,029Tabla 5. Error cuadrático medio <strong>de</strong> los métodos PIF y RCS.constante, que no se habían usado en el cálculo<strong>de</strong> los coeficientes <strong>de</strong> calibración, comparandolas intensida<strong>de</strong>s <strong>de</strong> estos objetivos <strong>de</strong> referenciacon las <strong>de</strong> las imágenes calibradas.La Tabla 5 muestra los errores cuadráticos medioscalculados para cada método <strong>de</strong> normalizaciónen cada banda. Ambos métodos tienen el mayorerror en el infrarrojo cercano (banda 4), lo cualpue<strong>de</strong> ser <strong>de</strong>bido a diferencias en la presencia <strong>de</strong>vegetación entre las imágenes. El método PIF, conexcepción <strong>de</strong> la banda 1, presenta un error cuadráticomedio inferior al RCS, por lo que fue el elegidopara normalizar las imágenes utilizadas en elestudio <strong>de</strong> los cambios <strong>de</strong> usos <strong>de</strong>l suelo.Se ha <strong>de</strong> indicar, no obstante, que estos métodos<strong>de</strong> normalización absoluta son aproximados,pues admiten la linealidad <strong>de</strong> las correcciones yconsi<strong>de</strong>ran las condiciones atmosféricas comoconstantes en toda el área <strong>de</strong> estudio.Con carácter general, se pue<strong>de</strong> concluir que,los resultados <strong>de</strong> los métodos <strong>de</strong> <strong>de</strong>tección <strong>de</strong>cambio aplicados son similares. Los métodos <strong>de</strong>diferencia, cociente y regresión <strong>de</strong> imágenes hansido capaces <strong>de</strong> diferenciar las áreas <strong>de</strong> cambioa partir <strong>de</strong> los índices <strong>de</strong> vegetación normalizados(NDVI), con índices <strong>de</strong> exactitud global <strong>de</strong>l88,67%, 81,33% y 80,67%, e índices Kappa <strong>de</strong>l0,72, 0,56 y 0,50 respectivamente. La diferencia<strong>de</strong> imágenes y el cociente produjeron una estimaciónsimilar <strong>de</strong> las áreas <strong>de</strong> cambio, con estimaciones<strong>de</strong> cambio iguales a 1.686 y 1.697 km 2respectivamente (ver Tabla 6).Diferencia Cociente RegresiónE. global 88,67 81,33 80,67Kappa 0,72 0,56 0,50Cambio (km 2 ) 24 13 66Tabla 6. Evaluación <strong>de</strong> la exactitud y superficie <strong>de</strong> cambio<strong>de</strong> las técnicas <strong>de</strong> realce <strong>de</strong> cambios.12 Revista <strong>de</strong> Tele<strong>de</strong>tección. ISSN: 1988-8740. 2010. 34: 5-15
Análisis <strong>de</strong> cambios <strong>de</strong> usos <strong>de</strong>l suelo en la «Vega <strong>de</strong> Granada»: correcciones radiométricas y evaluación <strong>de</strong>l cambioLa clasificación supervisada <strong>de</strong> las imágenesLandsat normalizadas correspondientes a la Vega<strong>de</strong> Granada para los años 1998 y 2004 dio comoresultado unos mapas clasificados con una exactitudaceptable, con índices kappa iguales a 0,86y 0,92 respectivamente. La exactitud <strong>de</strong>l mapa <strong>de</strong>cambios final es igual al producto <strong>de</strong> las exactitu<strong>de</strong>sindividuales <strong>de</strong> cada mapa clasificado (Yuanet al., 2005), por lo que se obtuvo finalmente unmapa <strong>de</strong> cambios con un kappa igual a 0,79.N0 4.000 8.000MetrosLeyenda17-7-1998ClaseBosqueChoperasEmbalsesMatorralOlivarRegadíoSuelo <strong>de</strong>snudoSuelo urbanoFigura 3. Mapa clasificado <strong>de</strong> la Vega <strong>de</strong> Granada para el verano <strong>de</strong> 1998.N0 4.000 8.000MetrosLeyenda1-7-2004ClaseBosqueChoperasEmbalsesMatorralOlivarRegadíoSuelo <strong>de</strong>snudoSuelo urbanoFigura 4. Mapa clasificado <strong>de</strong> la Vega <strong>de</strong> Granada para el verano <strong>de</strong> 2004.Revista <strong>de</strong> Tele<strong>de</strong>tección. ISSN: 1988-8740. 2010. 34: 5-15 13