Correção atmosférica - Revision history

⚠️Samuel-rosa: /* Introdução */ Pequena correção ortográfica

2013-11-28T18:34:48Z

Introdução: Pequena correção ortográfica

← Older revision		Revision as of 18:34, 28 November 2013
Line 13:		Line 13:

	# Converter o valor de nível de cinza (DN - digital number = valor do pixel) para valor de radiância no topo da atmosfera (veja a fórmula abaixo ou use {{cmd\|i.landsat.toar}}). {{cmd\|i.landsat.toar}} suporta todas as versões do Landsat, desde MSS, até TM e ETM+.		# Converter o valor de nível de cinza (DN - digital number = valor do pixel) para valor de radiância no topo da atmosfera (veja a fórmula abaixo ou use {{cmd\|i.landsat.toar}}). {{cmd\|i.landsat.toar}} suporta todas as versões do Landsat, desde MSS, até TM e ETM+.
	# ~~Execute~~ a correção atmosférica (a qual converte o valor da radiância no topo da atmosfera para o valor de radiância na superfície): realizado através de {{cmd\|i.atcorr}}.		# Executar a correção atmosférica (a qual converte o valor da radiância no topo da atmosfera para o valor de radiância na superfície): realizado através de {{cmd\|i.atcorr}}.

	''Nota:'' Lembre de verificar se as áreas cobertas com água possuem valor de nível de reflectância > 0, uma vez que a reflectância é sempre positiva. Se negativa, você tem um pepino a resolver! Valores negativos de reflectância indicam que as entradas da equação de transformação não pertencem à imagem processada.		''Nota:'' Lembre de verificar se as áreas cobertas com água possuem valor de nível de reflectância > 0, uma vez que a reflectância é sempre positiva. Se negativa, você tem um pepino a resolver! Valores negativos de reflectância indicam que as entradas da equação de transformação não pertencem à imagem processada.

⚠️Samuel-rosa: Pequena correção ortográfica.

2013-11-28T18:33:05Z

Pequena correção ortográfica.

← Older revision		Revision as of 18:33, 28 November 2013
Line 5:		Line 5:
	'''Objetivo desse tutorial:'''		'''Objetivo desse tutorial:'''

	* ~~Transformção de um imagem Landsat ETM+ dos~~ valores de nível de cinza (DN) para valores de radiância		* Transformar os valores de nível de cinza (DN) de um imagem Landsat ETM+ para valores de radiância
	* Realizar a correção atmosférica de uma imagem Landsat utilizando o módulo {{cmd\|i.atcorr}} no GRASS GIS		* Realizar a correção atmosférica de uma imagem Landsat utilizando o módulo {{cmd\|i.atcorr}} no GRASS GIS

Neteler: i.landsat.toar in 6.4, too

2013-02-02T18:57:58Z

i.landsat.toar in 6.4, too

← Older revision		Revision as of 18:57, 2 February 2013
Line 12:		Line 12:
	Os principais passos para o processamento de imagens de satélite são:		Os principais passos para o processamento de imagens de satélite são:

	# Converter o valor de nível de cinza (DN - digital number = valor do pixel) para valor de radiância no topo da atmosfera (veja a fórmula abaixo ou use {{cmd\|i.landsat.toar~~\|version=65~~}}). {{cmd\|i.landsat.toar~~\|version=65~~}} suporta todas as versões do Landsat, desde MSS, até TM e ETM+.		# Converter o valor de nível de cinza (DN - digital number = valor do pixel) para valor de radiância no topo da atmosfera (veja a fórmula abaixo ou use {{cmd\|i.landsat.toar}}). {{cmd\|i.landsat.toar}} suporta todas as versões do Landsat, desde MSS, até TM e ETM+.
	# Execute a correção atmosférica (a qual converte o valor da radiância no topo da atmosfera para o valor de radiância na superfície): realizado através de {{cmd\|i.atcorr}}.		# Execute a correção atmosférica (a qual converte o valor da radiância no topo da atmosfera para o valor de radiância na superfície): realizado através de {{cmd\|i.atcorr}}.

Line 26:		Line 26:
	== Cálculo dos valores de radiância ==		== Cálculo dos valores de radiância ==

	A amostra de dados da Carolina do Norte - EUA possui, entre outras coisas, uma imagem Landsat ETM+ do dia 24 de maio de 2002. Cada pixel dessa imagem contém um valor de nível de cinza (''DN - digital number''). Para que seja possível utilizar a imagem de satélite para a realização de cálculos, ou comparar valores entre diferentes sensores, esses valores de nível de cinza precisam ser convertidos para valores de radiância ou reflectância. As fórmulas usadas para fazer essa conversão estão descritas aqui para imagens Landsat (ou utilize ~~[[GRASS_AddOns#i.landsat.toar~~\|i.landsat.toar]]).		A amostra de dados da Carolina do Norte - EUA possui, entre outras coisas, uma imagem Landsat ETM+ do dia 24 de maio de 2002. Cada pixel dessa imagem contém um valor de nível de cinza (''DN - digital number''). Para que seja possível utilizar a imagem de satélite para a realização de cálculos, ou comparar valores entre diferentes sensores, esses valores de nível de cinza precisam ser convertidos para valores de radiância ou reflectância. As fórmulas usadas para fazer essa conversão estão descritas aqui para imagens Landsat (ou utilize {{cmd\|i.landsat.toar}}).
	A conversão de valor de nível de cinza para radiância no topo da atmosfera é realizada da seguinte maneira:		A conversão de valor de nível de cinza para radiância no topo da atmosfera é realizada da seguinte maneira:

⚠️Samuel-rosa: /* Atmospheric correction */

2012-10-21T02:56:40Z

Atmospheric correction

← Older revision		Revision as of 02:56, 21 October 2012
Line 82:		Line 82:
	A hora de passagem local resultante 10:42:07 corresponde a 15:42 em GMT, que corresponde a 15,67 em horas GMT decimais (minuto decimal: 42 * 100 / 60). Esse valor é necessário para o arquivo de controle.		A hora de passagem local resultante 10:42:07 corresponde a 15:42 em GMT, que corresponde a 15,67 em horas GMT decimais (minuto decimal: 42 * 100 / 60). Esse valor é necessário para o arquivo de controle.

	== ~~Atmospheric correction~~ ==		== Correção atmosférica ==

	This radiance image can be used for the atmospheric correction with the '''6S algorithm'''. The algorithm will transform the top-of-the-atmosphere radiance values to top-of-canopy reflectance values using predefined information on the aerosol content and atmospheric composition of the time the image was taken. What follows describes the method to use this algorithm in GRASS GIS. Again, only the calculations for band 1 are shown, and the numbers to be changed for the other bands are indicated in red. The 'icnd_lsat1.txt' control file consists of the following parameters, and is written with a text editor:		This radiance image can be used for the atmospheric correction with the '''6S algorithm'''. The algorithm will transform the top-of-the-atmosphere radiance values to top-of-canopy reflectance values using predefined information on the aerosol content and atmospheric composition of the time the image was taken. What follows describes the method to use this algorithm in GRASS GIS. Again, only the calculations for band 1 are shown, and the numbers to be changed for the other bands are indicated in red. The 'icnd_lsat1.txt' control file consists of the following parameters, and is written with a text editor:

⚠️Samuel-rosa at 02:55, 21 October 2012

2012-10-21T02:55:23Z

← Older revision		Revision as of 02:55, 21 October 2012
Line 125:		Line 125:

	[[Category: Documentação]]		[[Category: Documentação]]
	[[Category: ~~Tutorial~~]]		[[Category: Tutoriais]]
	[[Category: Processamento de imagens]]		[[Category: Processamento de imagens]]
	[[Category: Languages/pt-br]]		[[Category: Languages/pt-br]]

⚠️Samuel-rosa at 02:50, 21 October 2012

2012-10-21T02:50:56Z

← Older revision		Revision as of 02:50, 21 October 2012
Line 124:		Line 124:
	* veja http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.atcorr/README		* veja http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.atcorr/README

	[[Category: ~~Documentation~~]]		[[Category: Documentação]]
	[[Category: Tutorial]]		[[Category: Tutorial]]
	[[Category: ~~Image processing~~]]		[[Category: Processamento de imagens]]
	[[Category:Languages/pt-br]]		[[Category: Languages/pt-br]]

⚠️Samuel-rosa: /* Cálculo dos valores de radiância */

2012-10-21T02:49:01Z

Cálculo dos valores de radiância

← Older revision		Revision as of 02:49, 21 October 2012
Line 65:		Line 65:
	Para processamento mais rápido, utilize um modelo digital de elevação com o valores inteiros:		Para processamento mais rápido, utilize um modelo digital de elevação com o valores inteiros:

	{{cmd\|r.mapcalc}} "elev_int = round(elevation)"		{{cmd\|r.mapcalc}} "elev_int = round(elevation)"

	Encontre o valor médio de elevação para iniciar o processamento (necessário no arquivo 'icnd' abaixo)		Encontre o valor médio de elevação para iniciar o processamento (necessário no arquivo 'icnd' abaixo)

⚠️Samuel-rosa: /* Estimating the overpass time from the sun position */

2012-10-21T02:48:30Z

Estimating the overpass time from the sun position

← Older revision		Revision as of 02:48, 21 October 2012
Line 71:		Line 71:
	{{cmd\|r.univar}} elev_int		{{cmd\|r.univar}} elev_int

	== ~~Estimating the overpass time from the sun position~~ ==		== Estimativa do horário de passagem a partir da posição do solo ==

	~~The satellite overpass time can be estimated rather precisely from the sun position reported in metadata using~~ {{cmd\|r.sunmask}}:		O horário de passagem do satélite pode ser estimado com bastante precisão a partir da posição do sol apresentada nos metadados utilizando {{cmd\|r.sunmask}}. Os [ftp://ftp.glcf.umiacs.umd.edu/glcf/Landsat/WRS2/p016/r035/p016r035_7x20020524.ETM-EarthSat-Orthorectified/p016r035_7x20020524.met metadados] desse exemplo são: SUN_AZIMUTH = 120,8810347; SUN_ELEVATION = 64,7730999.
	~~The~~ [ftp://ftp.glcf.umiacs.umd.edu/glcf/Landsat/WRS2/p016/r035/p016r035_7x20020524.ETM-EarthSat-Orthorectified/p016r035_7x20020524.met ~~metadata~~] ~~of this example report~~: SUN_AZIMUTH = 120.8810347, SUN_ELEVATION = 64.~~7730999~~

	~~# iteratively change~~ hour ~~and minutes to obtain a good result~~, timezone ~~needs to be correct of course~~:		Mude as horas (''hour'') e minutos (''min'') iterativamente para obter bons resultados; obviamente, o fuso horário (''timezone'') precisa estar correto:
	{{cmd\|r.sunmask}} -s elev=elevation out=dummy year=2002 month=5 day=24 hour=10 min=42 sec=7 timezone=-5		{{cmd\|r.sunmask}} -s elev=elevation out=dummy year=2002 month=5 day=24 hour=10 min=42 sec=7 timezone=-5

	~~Reports: sun azimuth~~: 121.342461~~, sun angle above horz.~~(~~refraction corrected~~): 65.~~396652~~		Relatório: azimute solar = 121,342461; ângulo solar acima do horizonte (corrigido para a refração) = 65,396652.

	~~The resulting overpass~~ local ~~time~~ 10:42:07 ~~corresponds to~~ 15:42 in GMT ~~which corresponds to~~ 15.67 ~~in decimal~~ GMT ~~hours~~ (decimal ~~minutes~~: 42 * 100 / 60). ~~This value is needed for the control file~~.		A hora de passagem local resultante 10:42:07 corresponde a 15:42 em GMT, que corresponde a 15,67 em horas GMT decimais (minuto decimal: 42 * 100 / 60). Esse valor é necessário para o arquivo de controle.

	== Atmospheric correction ==		== Atmospheric correction ==

⚠️Samuel-rosa: /* References */

2012-10-21T02:26:05Z

References

← Older revision		Revision as of 02:26, 21 October 2012
Line 113:		Line 113:
	More information on the use of {{cmd\|i.atcorr}} for other images can be found in the {{cmd\|i.atcorr}} manual.		More information on the use of {{cmd\|i.atcorr}} for other images can be found in the {{cmd\|i.atcorr}} manual.

	== ~~References~~ ==		== Referências ==

	* [http://6s.ltdri.org 6S Web site]		* [http://6s.ltdri.org 6S Web site]
Line 122:		Line 122:
	* Song, C., Woodcock, C.E., Seto, K.C., Lenney, M.P., Macomber, S.A., 2001. Classification and Change Detection Using Landsat TM Data: When and How to Correct Atmospheric Effects? Remote Sensing of Environment 75, 230-244. [http://www.unc.edu/courses/2008spring/geog/577/001/www/Song01_RSE.pdf PDF]		* Song, C., Woodcock, C.E., Seto, K.C., Lenney, M.P., Macomber, S.A., 2001. Classification and Change Detection Using Landsat TM Data: When and How to Correct Atmospheric Effects? Remote Sensing of Environment 75, 230-244. [http://www.unc.edu/courses/2008spring/geog/577/001/www/Song01_RSE.pdf PDF]

	~~How to add new sensors to~~ i.atcorr:		Como adicionar novos sensores ao i.atcorr:
	* ~~see~~ http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.atcorr/README		* veja http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.atcorr/README

	[[Category: Documentation]]		[[Category: Documentation]]

⚠️Samuel-rosa: /* Cálculo dos valores de radiância */

2012-10-21T02:24:20Z

Cálculo dos valores de radiância

← Older revision		Revision as of 02:24, 21 October 2012
Line 36:		Line 36:
	Lλ = ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN)) * (QCAL-QCALMIN) + LMINλ		Lλ = ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN)) * (QCAL-QCALMIN) + LMINλ

	~~Where~~:		Onde:

	* Lλ - ~~spectral Radiance at the~~ sensor~~'s aperture in~~ watts/(~~meter squared~~ * ster * μm), ~~the apparent radiance as seen by the satellite~~ sensor;		* Lλ - radiância espectral na abertura do sensor em watts/(metro quadrado * ster * μm), a radiação aparente tal como é vista pelo sensor do satélite;
	* QCAL - ~~the quantized calibrated~~ pixel ~~value in DN~~;		* QCAL - o valor quantizado calibrado de um pixel em nível de cinza;
	* LMINλ - ~~the spectral radiance that is scaled to~~ QCALMIN in watts/(~~meter squared~~ * ster * μm);		* LMINλ - a radiância espectral que é dimensionada para QCALMIN em watts/(metro quadrado * ster * μm);
	* LMAXλ - ~~the spectral radiance that is scaled to~~ QCALMAX in watts/(~~meter squared~~ * ster * μm);		* LMAXλ - a radiância espectral que é dimensionada para QCALMAX em watts/(metro quadrado * ster * μm);
	* QCALMIN - ~~the minimum quantized calibrated~~ pixel ~~value~~ (~~corresponding to~~ LMINλ) ~~in DN~~;		* QCALMIN - o valor mínimo quantizado calibrado de um pixel (correspondente à LMINλ) em nível de cinza;
	* QCALMAX - ~~the maximum quantized calibrated~~ pixel ~~value~~ (~~corresponding to~~ LMAXλ) ~~in DN=~~255.		* QCALMAX - o valor máximo quantizado calibrado de um pixel (correspondente à LMAXλ) em nível de cinza igual a 255.

	LMINλ e LMAXλ são os valores de radiância correspondentes aos valores mínimo e máximo de nível de cinza, e são apresentados no arquivo de metadados para cada imagem, ou na tabela 1. As informações sobre elevado ganho (''high gain'') ou baixo ganho (''low gain'') também são apresentadas no arquivo de metadados para cada imagem Landsat. O menor valor de nível de cinza (QCALMIN) é igual a 1 para imagens Landsat ETM+ (1), e o valor máximo de nível de cinza (QCALMAX) é igual a 255. QCAL é o valor de nível de cinza para cada pixel individual na imagem Landsat.		LMINλ e LMAXλ são os valores de radiância correspondentes aos valores mínimo e máximo de nível de cinza, e são apresentados no arquivo de metadados para cada imagem, ou na tabela 1. As informações sobre elevado ganho (''high gain'') ou baixo ganho (''low gain'') também são apresentadas no arquivo de metadados para cada imagem Landsat. O menor valor de nível de cinza (QCALMIN) é igual a 1 para imagens Landsat ETM+ (1), e o valor máximo de nível de cinza (QCALMAX) é igual a 255. QCAL é o valor de nível de cinza para cada pixel individual na imagem Landsat.
Line 51:		Line 51:
	{{cmd\|r.info}} lsat7_2002_xx		{{cmd\|r.info}} lsat7_2002_xx

	Na opção ‘comments’, os valores máximo e mínimo de radiância (LMAX ~~and~~ LMIN) para cada banda são fornecidos.		Na opção ‘comments’, os valores máximo e mínimo de radiância (LMAX e LMIN) para cada banda são fornecidos.

	Power users ~~podem usar este comando~~ para uma saída mais legível (exemplo para o canal 1):		Também é possível usar o comando abaixo (''Power users'') para uma saída mais legível (exemplo para o canal 1):

	CHAN=1		CHAN=1