1 / 16

Por Carmen Recondo

Corrección topográfica de la imagen para mejorar las clasificaciones en zonas montañosas. Corrección por el modelo de Minnaert. Por Carmen Recondo. SOLUCIÓN : En zonas montañosas , previamente a la clasificación , ha de realizarse la corrección topográfica de la imagen 

percy
Download Presentation

Por Carmen Recondo

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Corrección topográfica de la imagen para mejorar las clasificaciones en zonas montañosas. Corrección por el modelo de Minnaert. Por Carmen Recondo

  2. SOLUCIÓN: En zonas montañosas, previamente a la clasificación, ha de realizarse la corrección topográfica de la imagen Homogeneizar el valor de ND entre las zonas más y menos iluminadas, “aplanando” la superficie: ND(+ ) y ND( ) Una forma de hacerlo  El cociente  LIes la radiancia recibida por el sensor de una superficie inclinada y LHla que recibiría de la correspondiente superficie horizontal proyectada. Corrección:NDcorregido= NDobservado Casos: • Zona orientada al Sol LI > LH < 1 ND corregido< ND observado • Zona no orient.al Sol LI < LH > 1 ND corregido> ND observado • Zona horizontal  LI = LH = 1 ND corregido= ND observado

  3. El vector gradiente topográfico, El vector solar, El vector de observación, El ángulo i El ángulo e 2. Vectores y ángulos involucrados en la radiancia detectada por el sensor: SENSOR NORMAL SUPERFICIE DEL TERRENO e i HORIZONTAL

  4. a) Superficies lambertianas: dispersan la luz uniformemente en todas las direcciones La radiancia (L) recibida por el sensor se asume que es: L= Cte cosi e independiente de cos e 3. Tipos de superficie y modelos de radiancia. Fórmulas para  SENSOR NORMAL i •  es : • Casos: • Zona orientada al Sol  iI <V (cos i)I > cosV  < 1 • Zona no orient.al Sol  iI >V (cos i)I < cosV  > 1 • Zona horizontal  iI =V (cos i)I = cosV  = 1

  5. b) Superficies no lambertianas: más realistas, dispersan la luz de forma no uniforme Modelos de radiancia (L) diversos. El modelo de Minnaert: dondek () es la constante de Minnaert (entre 0 y 1=caso lambertiano),que puede ser obtenida de la propia imagen (una para cada banda o ) convirtiendo en lineal la ecuación anterior: donde L ND •  para cada banda sería:

  6. 4b) Imágenes derivadas del MDE: px , py y cos e • px =en sentido x (OE): • py=en sentido y (SN): • Cálculo de la pendiente pxy pyde cada punto (x,y) del MDE: Filtros

  7. Obtenidas pxy pyse calculacos een cada punto: para una superficie horizontal cos e=1

  8. 4c) Imagen obtenida del MDE y de los datos solares: cos i • El cálculo de cos i en cada punto requiere pxy py (distintas para cada punto) y las componentes del vector solar (ctes): Para una superficie horizontal: ; Se toman sólo los cos i0(i entre 0º-90º) Zona orientada al Sol (+ iluminada): cos i > cosV ya que i < V Zona no orient. al Sol (iluminada): cos i < cosV ya que i > V

  9. 4d) Obtención de la constante k de Minnaert para cada banda de la imagen donde L ND • Del modelo de Minnaert para L habíamos obtenido: • k () es entonces la pendiente del ajuste lineal y=b+kx, donde: Se obtiene una k (entre 0 y 1) para cada banda de la imagen yBANDA n = ln (NDBANDA n  cose) x = ln (cosi cose) En este ejemplo k crece con la superficie se hace más lambertiana con (k1= 0.16, k2 = 0.29, k3 = 0.37, k4 = 0.48, k5 = 0.67, k7 = 0.75) Banda 2 Banda 5

  10. 4e) Cálculo del cociente  para cada banda de la imagen • Obtenida k () para cada banda se calcula (x,y): Ejemplo para la banda 2: Para una superficie horizontal   =1 Zona orientada al Sol (+ iluminada)   < 1 Zona no orient. al Sol ( iluminada)   > 1

  11. 4f) Corrección topográfica de cada banda de la imagen • Obtenido  (x,y) para cada banda, cada una se corrige según: NDcorregido= NDobservado Se homogeneiza el valor de ND entre las zonas más y menos iluminadas, “aplanando” la superficie: ND(+ ) y ND( ) Ejemplo para la banda 5: Antes de la corrección:Después de la corrección:

  12. 4g) Comprobación visual de los resultados • Se observa finalmente que el problema inicial ha sido corregido, es decir, para cualquier composición en color, la topografía “desaparece” de la imagen Ejemplo: Imagen Landsat-TM del concejo de Mieres (Composición = color real; RGB=321) Antes de la corrección:Después de la corrección:

  13. 4g) Comprobación visual de los resultados • Se observa finalmente que el problema inicial ha sido corregido, es decir, para cualquier composición en color, la topografía “desaparece” de la imagen Ejemplo: Imagen Landsat-TM del concejo de Mieres (Composición = 3 primeras componentes principales) Antes de la corrección:Después de la corrección:

  14. ¿Modelo de Minnaert?

  15. ¿Modelo de Minnaert?

  16. ¿Modelo de Minnaert?

More Related