{"id":7751,"date":"2019-01-10T04:04:29","date_gmt":"2019-01-10T03:04:29","guid":{"rendered":"http:\/\/www.sigterritoires.fr\/?p=7751"},"modified":"2023-04-29T06:48:19","modified_gmt":"2023-04-29T04:48:19","slug":"tutorial-de-clasificacion-de-imagenes-con-qgis-2-1-ajuste-atmosferico-de-imagenes-landsat-8","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/es\/tutorial-de-clasificacion-de-imagenes-con-qgis-2-1-ajuste-atmosferico-de-imagenes-landsat-8\/","title":{"rendered":"Tutorial de clasificaci\u00f3n de im\u00e1genes con QGis: 2.1 Ajuste atmosf\u00e9rico de im\u00e1genes Landsat 8"},"content":{"rendered":"\n

\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nHemos\nvisto anteriormente c\u00f3mo corregir las im\u00e1genes Landsat para considerar la\nreflectanciaTOA (tope de la atm\u00f3sfera). Aqu\u00ed hay un diagrama que nos ayudar\u00e1 a\nentender este fen\u00f3meno:  \n\n\n\n<\/p>\n\n\n\n\n\n\n\n

\"\"<\/figure>\n\n\n\n

La\nreflectancia que hemos calculado es por lo tanto, el % de luz reflejada con\nrespecto a la totalidad de la luz visible incidente . Pero como se ha visto en\nel diagrama , el sensor del sat\u00e9lite mide dos cosas al mismo tiempo: la luz\nreflejada por nuestros objetivos en la superficie de la tierra ,y la luz dispersada\npor las part\u00edculas suspendidas en el aire. <\/p>\n\n\n\n

Se puede\nrealizar la correcci\u00f3n atmosf\u00e9rica de las im\u00e1genes satelitales para eliminar la\nporci\u00f3n de luz debida a la difusi\u00f3n. <\/p>\n\n\n\n

Como ya\nhemos dicho, \u00e9sto s\u00f3lo tiene inter\u00e9s si usted trabaja en im\u00e1genes separadas en\nel tiempo El porcentaje de difusi\u00f3n es el mismo para todo el conjunto de una\nimagen, realizar esta correcci\u00f3n es perder el tiempo. <\/p>\n\n\n\n


\n Las instrucciones para realizar esa correcci\u00f3n en im\u00e1genes Landsat est\u00e1\ndisponible en <\/p>\n\n\n\n

Pero si usted es un principiante seguramente se le dificultar\u00e1 entender\nel procedimiento. Veremos aqu\u00ed  el procedimiento  usando las\nherramientas de QGis. <\/p>\n\n\n\n

M\u00e9todo DOS<\/strong> (Dark Object Subtraction)<\/strong>\n<\/p>\n\n\n\n


\nEl m\u00e9todo que vamos a emplear es el m\u00e1s com\u00fanmente empleado.
\nLa resta de objetos negros es una m\u00e9todo de correcci\u00f3n atmosf\u00e9rica empirico\nsimple para la imagen satelital , que asume que la reflectancia de los objetos\noscuros incluye una componente importante de la difusi\u00f3n atmosf\u00e9rica. <\/p>\n\n\n\n

El\nconcepto es simple: si se tiene un objeto completamente negro y, por lo tanto ,\nno refleja la luz, este objeto que deber\u00eda tener un valor 0 de reflectancia\ntendr\u00e1 un valor no nulo (diferente de 0) en el sensor del sat\u00e9lite. Este valor\nser\u00e1 la reflectancia debida a las part\u00edculas presentes en la atm\u00f3sfera. <\/p>\n\n\n\n

La resta\nde objetos oscuros investiga en cada banda el valor de p\u00edxel m\u00e1s oscuro. La\ndifusi\u00f3n es suprimida restando este valor a cada p\u00edxel de la banda . <\/p>\n\n\n\n

Lo\nprimero: la correcci\u00f3n que vamos a hacer solo aplica a la luz visible y el\ninfrarrojo cercano. Vamos a trabajar por lo tanto , para im\u00e1genes Landsat8, en\nlas bandas 2 (Azul), 3 (verde), 4 (Rojo) y 5 (infrarrojo cercano). <\/p>\n\n\n\n

Lo\nsegundo que hay que saber: el principio en uso es calcular la distribuci\u00f3n para\nla banda 4 (rojo) y usar una f\u00f3rmula de correspondencia a fin de determinar la\ncorrecci\u00f3n para las otras tres bandas. <\/p>\n\n\n\n

Lo\ntercero: vamos a trabajar con las bandas 2 , 3, 4 y 5 en reflectancia TOA y\ntambi\u00e9n en radiancia bruta (DN) para la banda 4 (rojo). En el manual del\nusuario citado previamente no es evidente  a primera vista . Veamos el\nprocedimiento paso a paso: <\/p>\n\n\n\n

1- se\ncalcula la reflectancia TOA para las cuatro bandas con la f\u00f3rmula <\/p>\n\n\n\n

TOA banda<\/strong> = (<\/strong> (DN * 0.00002)-0.1)\n\/ sen (<\/strong> elevaci\u00f3n<\/strong> solar<\/strong> )<\/strong> <\/p>\n\n\n\n

2-\nanalizando los valores brutos de radiancia (DN) de la banda roja, vamos a\ndeterminar una radiancia debida a la difusi\u00f3n atmosf\u00e9rica ( DNred ) <\/p>\n\n\n\n

3-\ntranformaremos esta difusi\u00f3n en reflectancia TOA <\/p>\n\n\n\n

CORrojo<\/strong> = (<\/strong> (<\/strong> DNrojo*\n0.00002) -0.1) \/ sen (<\/strong> elevaci\u00f3n<\/strong> solar<\/strong> )<\/strong> <\/p>\n\n\n\n

4- con un\n\u00e1baco se ingresa esta difusi\u00f3n de la banda roja y se obtiene la difusi\u00f3n para\nlas bandas 2,3 y 5 ( CORazul , CORverde , CORir). <\/p>\n\n\n\n

5- se\nrestan los valores de reflectancias obtenidos de TOA para obtener los valores\nde reflectancia de las cuatro bandas corregido por el efecto atmosf\u00e9rico.<\/p>\n\n\n\n

Rrojo<\/strong> =<\/strong> TOA (<\/strong> rojo) \u2013<\/strong>\nCORRrojo<\/strong> <\/p>\n\n\n\n

Razul<\/strong> =<\/strong> TOA (<\/strong> azul) \u2013<\/strong>\nCORazul<\/strong> <\/p>\n\n\n\n

Rverde<\/strong> =<\/strong> TOA (<\/strong> verde)\n\u2013<\/strong> CORverde<\/strong> <\/p>\n\n\n\n

Rir<\/strong> =<\/strong> TOA (<\/strong> ir<\/strong> )\n\u2013<\/strong> CORir<\/strong> <\/p>\n\n\n\n

Veamos en\ndetalle los pasos 2 y 4. <\/p>\n\n\n\n

C\u00e1lculo de la difusi\u00f3n de la banda roja.<\/strong> <\/p>\n\n\n\n

Hay varios m\u00e9todos para determinar la correcci\u00f3n a realizar en la banda\nroja. El m\u00e9todo recomendado es el conocido como Bin 5 – 0.008.<\/strong> Usted encontrar\u00e1 aqu\u00ed una descripci\u00f3n detallada del m\u00e9todo.\n<\/p>\n\n\n\n

En la pr\u00e1ctica, con QGis, usaremos el histograma de la banda para determinar el valor umbral de la difusi\u00f3n atmosf\u00e9rica. 1- Primero debe ejecutar la herramienta Estad\u00edsticas de capa raster<\/strong> en los datos en la banda 4 (rojo). \u00a0 <\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

Una vez ejecutado, tendr\u00e1 una l\u00ednea de resultados en el panel del Visualizar resultados\u00a0 <\/strong><\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nHaga\ndoble clic en esta l\u00ednea para abrir los resultados en su navegador:  \n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

Tenga en cuenta el valor m\u00ednimo aqu\u00ed 3916. <\/p>\n\n\n\n

Buscaremos el valor que se asignar\u00e1 a la correcci\u00f3n de la difusi\u00f3n\natmosf\u00e9rica en el histograma de la imagen. \n\nHaga clic en la capa de la imagen ->\nPropiedades -> Histogram\n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nIngrese\nel valor m\u00ednimo que se encuentra en las estad\u00edsticas (aqu\u00ed 3916) en el cuadro\nMin, luego haga clic en el cuadro Max para mover el valor m\u00ednimo de la l\u00ednea de\npuntos en el histograma. \n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

Puede usar el formato de l\u00ednea para encontrar la base del histograma,\npero en este ejemplo, se usar\u00e1 (y se recomienda) el formato de barra. Si el\nhistograma est\u00e1 en formato de l\u00ednea, haga clic en Preferencias \/ Acciones<\/strong>\ny desactive la casilla de verificaci\u00f3n Marcar como l\u00edneas<\/strong> . <\/p>\n\n\n\n

Aparecer\u00e1 un histograma en formato de barra. El histograma de QGIS no representa\/\nno muestra todos los valores como el histograma de ArcGIS, pero muestra muy\nbien la forma que puede ser apropiada para establecer un valor de dispersi\u00f3n\nque represente la base del histograma (en caso de un aumento relativamente\nabrupto de las frecuencias).   Los siguientes pasos muestran c\u00f3mo se puede\nusar el histograma de QGIS para establecer la dispersi\u00f3n. \n\nPase el cursor del mouse sobre el histograma\npara ver una lupa. Arrastre la lupa hasta el extremo inferior del histograma\npara ampliar la base del mismo. Si hace demasiado zoom o si desea volver a\nempezar por cualquier raz\u00f3n, haga clic en Preferencias \/ Acciones<\/strong> y Recalcular\nhistograma<\/strong> . El histograma se ampl\u00eda a continuaci\u00f3n.   \n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nObservar\u00e1\nque la barra con el valor m\u00e1s bajo sigue siendo bastante superior al valor de\nimagen bajo que ingres\u00f3 anteriormente en el \u00e1rea de Min (3619); el histograma\nde QGIS no muestra muchos valores en una cola estad\u00edstica (Landsat 8\ngeneralmente tiene largas colas estad\u00edsticas). La base del histograma es el\n\u00e1rea utilizada para establecer la dispersi\u00f3n. Haciendo zoom en las tres\nprimeras barras del histograma:  \n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nEn\neste ejemplo, los valores comprendidos, aproximadamente, entre 5570 y 5595\npueden acercarse a la base del histograma; \u00e9stos son valores muy bajos y no\nest\u00e1n demasiado en la cola estad\u00edstica (este proceso puede ser ligeramente\nsubjetivo). Para las necesidades de este documento, el DN 5569 ser\u00e1 utilizado\ncomo base para el histograma. El histograma se ampliar\u00e1 una vez m\u00e1s para este ejemplo.\n\n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nLa\nregla a utilizar es tomar como valor de dispersi\u00f3n el valor m\u00ednimo de la\nprimera barra que exceda una frecuencia de 5 (de ah\u00ed el nombre Bin 5) y que no\ntenga frecuencias inferiores a 5 a su derecha. En el siguiente ejemplo:   \n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

Las dos primeras barras exceden el valor 5, pero como la tercera es m\u00e1s\nbaja, se retendr\u00e1 el valor m\u00ednimo de la cuarta barra. <\/p>\n\n\n\n

En nuestro ejemplo, el valor de radiancia que se retendr\u00e1 para la\ncorrecci\u00f3n atmosf\u00e9rica de la banda roja ser\u00e1 5569.<\/strong> <\/p>\n\n\n\n

Podemos pasar a la tercera etapa, <\/p>\n\n\n\n

3- transformamos este resplandor en reflectancia TOA. <\/p>\n\n\n\n

COR(rojo) = ((DNrojo * 0.00002)\n-0.1) \/ sen (elevaci\u00f3n solar)<\/strong> <\/p>\n\n\n\n

  COR(rojo) = ((5569 *\n0,00002) -0,1) \/0,42631886<\/strong>   = 0.026694<\/strong> <\/p>\n\n\n\n

Este valor corresponde al Bin 5. Dado que el m\u00e9todo es el Bin5 – 0.008,\neliminamos 0.008 de este valor: 0.018694<\/strong> <\/p>\n\n\n\n

Es este valor el que debe restarse de los valores de TOA del rojo previamente\ncalculado para tener los valores corregidos de la difusi\u00f3n atmosf\u00e9rica. <\/p>\n\n\n\n

Correcci\u00f3n de otras bandas.<\/strong> <\/p>\n\n\n\n

A partir del valor de correcci\u00f3n para el rojo obtenemos los valores que\nse aplicar\u00e1n a cada una de las otras bandas: <\/p>\n\n\n\n

Ir a la siguiente p\u00e1gina \n\n\nEncontrar\u00e1 una ventana donde puede ingresar el\nvalor de correcci\u00f3n para el rojo   \n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\nIngrese\nsu valor (aqu\u00ed 0.018647<\/strong> ), autom\u00e1ticamente tendr\u00e1 los valores de las\notras correcciones   \n\n\n\n<\/p>\n\n\n\n

\"\"<\/figure>\n\n\n\n

As\u00ed obtenemos las correcciones de <\/p>\n\n\n\n

  • 0.06309 para la banda azul (2) <\/li>
  • 0.03442 para la banda verde (3) <\/li>
  • 0.00626 para la banda de\ninfrarrojo cercano (5) <\/li><\/ul>\n\n\n\n

    Para tener las bandas corregidas usamos la calculadora raster restando a\ncada banda corregida TOA la constante de correcci\u00f3n que hemos encontrado. <\/p>\n","protected":false},"excerpt":{"rendered":"

    Hemos visto anteriormente c\u00f3mo corregir las im\u00e1genes Landsat para considerar la reflectanciaTOA (tope de la atm\u00f3sfera). Aqu\u00ed hay un diagrama que nos ayudar\u00e1 a entender este fen\u00f3meno: <\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[842],"tags":[],"class_list":["post-7751","post","type-post","status-publish","format-standard","hentry","category-non-classe-es"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-211","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/7751","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=7751"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/7751\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=7751"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=7751"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=7751"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}