Tutorial MDE ALOS-1 Huascar Vidal v1

Tutorial para obteno, correes e interpolao (opcional) do Modelo
Digital de Elevao (sensor PALSAR) do Satlite ALOS-1 de 2006 a 2011.
Por Huscar P. Vidal de Oliveira

Doutorando em Desenvolvimento e Meio Ambiente
(PRODEMA/UFC)
[email protected]
O satlite ALOS (Advanced Land Observing Satellite), chamado pelos japoneses de

Daichi, um projeto conjunto entre a Japan Aerospace Exploration Agency (JAXA) e a Japan
Resources Observation System Organization (JAROS). Seu desenvolvimento tem como um dos
principais objetivos proporcionar dados de observao da Terra que contribuam ao
desenvolvimento sustentvel, cartografia topogrfica, vigilncia do meio ambiente, desastres e
mudanas climticas em todo o mundo. Foi lanado em 24 de janeiro de 2006, no Centro
Espacial de Tanegashima (TNSC) e iniciou fornecimento de dados ao pblico em 24 de outubro
de 2006 e encerrado em abril de 2011. O satlite dispe de trs instrumentos de sensoriamento
remoto: dois instrumentos pticos, PRISM (Panchromatic Remote-sensing Instrument for Stereo
Mapping) e AVNIR-2 (Advanced Visible and Near-Infrared Radiometer type 2) e um radar
polarimtrico de abertura sinttica de banda larga PALSAR (Phased Array L-band Synthetic
Aperture Radar) para observao terrestre de dia e de noite.
A plataforma ALOS apresenta uma srie de melhorias ao seu antecessor JERS-1
(Japanese Earth Resources Satellite), que esteve em operao de 1992 a 1998.
O PALSAR um instrumento totalmente polarimtrico, operando em diferentes modos,
com polarizao simples FBS (Fine Beam Single polarization: HH ou VV), polarizao dupla FBD
(Fine Beam Dual polarization: HH + HV ou VV + VH) ou modo polarimtrico completo POL (HH +
HV + VH + VV). Ele tambm possui o modo ScanSAR, com polarizao nica (HH ou VV; 3/4/5-
beam).
Dentre as vantagens de um sensor imageador ativo de alta frequncia por microondas
como o PALSAR, reside na possibilidade de aquisio de dados topogrficos reais da superfcie
terrestre, uma vez que no sofrem interferncias de nuvens e facilidade em ultrapassar o dossel
vegetal. Como tambm para controle e monitoramento, visto que independe de iluminao
natural ou de emisso prpria do alvo.
Para a obteno do modelo, enumera-se os passos a seguir:
1) Cadastro gratuito no geoportal de descarregamento de imagens de satlite da Alaska

Satellite Facility (UAF/NASA), disponvel em: https://vertex.daac.asf.alaska.edu/
Recomendaes:
*No pode haver duplicidade de informaes no banco de dados da NASA, caso usurio j tenha
cadastro no geoportal Earth Explorer (https://earthexplorer.usgs.gov/) pode entrar direto com
seu login e senha, caso no lembre, pode recuperar seguindo os passos que pede na prpria
janela de Login e Senha;
*Para usurio iniciante, ateno no preenchimento obrigatrio das informaes pessoais que
pede no cadastro e principalmente ao definir uma senha forte com caracteres diversos;
*Ao final ser enviado um email automtico para autenticao do cadastro;
*Caso usurio no receba o email de autenticao implica em erros de cadastro ou duplicidade
de dados.
2) Delimitao da rea de interesse no ecl (mapa-mndi) do geoportal:

De forma a facilitar a navegao do usurio, o geoportal foi desenvolvido de forma bastante
dinmica com delimitao da rea de interesse atribuindo controle total do mouse nas funes
de Zoom e cliques sucessivos somatrios em vertex e ajustamento da poligonal.
3) Ao indicar a Plataforma ALOS PALSAR, o usurio ter que selecionar dentre as opes de
imagens disponveis, o ALOS PALSAR modo FBD, que o que contm o DEM.
Sugesto: Para facilitar o download dos arquivos zipados utilize o aplicativo DownThemAll!
(http://www.downthemall.net/) ou via [complementos/extenses/pesquisar_complementos]
do navegador FireFox. Instalado o programa, fica disponvel na janela do boto direito do mouse,
as funes de downloads sistemticos de arquivos.
Nos passos a seguir, so relacionados s correes de incongruncias resultantes desde

aquisio de microondas pelo sensor (radar), transmisso para a estao de recepo e gerao
dos dados. Essas imperfeies refletidas no raster, em seus pixels, precisam ser corrigidas, no
meramente por uma questo esttica, mas no intuito de corrigir todo o conjunto de dados do
clculo numrico de elevao. Sem o qual o modelo digital no iria condizer com a realidade do
territrio nele representado.
Para alguns procedimentos a seguir subentende-se alguns conhecimentos prvios dos
softwares ArcGIS e tambm QGIS, j que alguns passos podem ser realizados pelo mesmo.
4) Com base na rea litornea desse exemplo, identificou-se atravs da ferramenta Identify,
que todos os pixels de superfcie do mar apresentam pixel value (valor negativo) igual a -8. Essa
informao torna-se fator decisrio que ser utilizada mais adiante.
SUGESTO: Para reas de interesse sem referencial do mar, opo seria gerar shapefile de
pontos de controle, converter para formato KML e comparar com dados de altitude
representados nas imagens de satlite do programa Google Earth.
5) Com base na informao anterior, passaremos o raster por um filtro e converso do formato
GeoTIFF para IMG, para ampliao de range da imagem. No ArcToobox Data Management
Tools Raster Raster Dataset Copy Raster, segue-se os dados de configurao do
exemplo:
Observe que o NoData Value (optional) foi preenchido com valor inferior a -8 (que o valor
da superfcie do mar nesse exemplo e que ser corrigido para 0 nos passos seguintes).
Note o resultado da converso,

onde o formato IMG que nativo
no ArcGIS, amplia o range,
transparecendo valores antes
ocultos no formato GeoTIFF.
6) Para corrigir essa inconsistncia de valores de altitude, faz-se necessrio converter o DEM
que agora est no formato IMG para GRID, permitindo que esses valores estejam acessveis
numa tabela de atributos.
7) Remover valores negativos mediante aplicao da ferramenta Con. Nesse exemplo

praticado, o valor -8 corresponde ao valor 0, referente ao nvel do mar. Acesse a caixa de
ferramentas ArcToolbox, opo Spatial Analyst Tools Conditional.
Escreva a expresso VALUE >= -8, para que a ferramenta mantenha os pixels maiores e
iguais a -8, excluindo os demais pixels de valores incongruentes e salve em formato IMG.
8) Com a remoo dos pixels inconsistentes, utiliza-se a ferramenta Reclassify, que possui dupla
funo, primeiramente ir preencher os espaos vazios gerados anteriormente pela ferramenta
Con, atribuindo para os mesmos o valor 0 e por ltimo a reclassificao do range de valores de
altitude. No ArcToolbox Reclass Reclassify configure a ferramenta como segue:
Selecione o raster criando pela ferramenta Con, clique no boto Unique para reclassificao,
v para o final de tabela, substitua NoData da coluna New values pelo valor 0, nomeie e salve o
novo raster e execute a ferramenta.
9) Para finalizar as correes do MDE, antes de iniciar a interpolao, iremos utilizar a
ferramenta Fill. Sua funo pode ser melhor representada nos esquemas a seguir:
Pixels com anomilias de altitude:
com valores superiores
com valores inferiores
Essas duas ltimas ferramentas corretivas, atribuem novo valores aos pixels com anomalias,
com base nas informaes dos vizinhos mais prximos. Para tanto, no intuito de corrigir essas
grandes imperfeies, vamos ao ArcToolbox Spatial Analyst Tools Hydrology Fill.
Sugesto: Para os profissionais que pretendem extrair ao mximo todo o potencial de ALOS-1
sugere-se as seguintes fontes (links):
1. https://earth.esa.int/web/polsarpro/home;
2. http://www.dsr.inpe.br/sbsr2011/files/p0537.pdf;
3. http://repositorio.unesp.br/bitstream/handle/11449/86778/beneditti_ca_me_prud.p
df?sequence=1&isAllowed=y.
Foi demonstrada de forma sucinta as correes do MDE do ALOS Palsar com resoluo
espacial de 12,5m. Cabe ao usurio que nos acompanha avaliar seus objetivos, em relao
escala de trabalho, tempo e recursos, tendo em vista que as recomendaes a seguir limitam a
realizao da interpolao. Segue-se algumas recomendaes inicias:
*Nos passos demonstrativos seguintes, ser interpolado, um recorte da cena corrigida, de 12,5m
para 1m, o que implica na gerao de arquivo bem mais pesado que sua rea equivalente
anterior, sendo necessrio para tanto, recursos (processamento/armazenamento) e tempo
generosos do computador e do usurio;
SUGESTO: Caso o leitor pretenda trabalhar em escala municipal recomenda-se realizar recortes
individualizados do MDE com elaborao de shapefiles de articulaes em escala (cadastral)
1:10.000 (Sistema Cartogrfico Nacional), de sua rea de interesse. Para no haver
descontinuidade do MDE no resultado final, juno (mosaico) das partes, utilize do ArcToolbox
(Analysis Tools/Proximity) a ferramenta Buffer para expandir o limite das articulaes uma
distncia variavelmente pequena (exemplo 20m):
*No intuito de reduzir ainda mais o tempo de processamento do computador (durante

interpolao), o caro leitor poder restringir a rea, utilizando a ferramenta Clip (ArcToolbox/
Analysis Tools/Extract) na gerao de novo SHP (polgono) entre os shapes das articulaes e o
shape de delimitao de seu interesse. Em seguida, recorte o raster com SHP criado, com a
ferramenta Clip (ArcToolbox/Data Management Tools/Raster/Raster Processing).
Com o MDE recortado pelo shapefile (resultante entre articulao 1:10.000, por
desdobramento do mapeamento sistemtico nacional e o limite da sua rea de interesse),
passaremos para a interpolao.
10) Converta o MDE recortado para pontos atravs da ferramenta Raster to Point
(ArcToolbox/Conversion Tools/From Raster).
11) Finalmente, abra a ferramenta Topo to Raster (ArcToolbox/Spatial Analyst
Tools/Interpolation), configure de acordo com as marcaes na figura abaixo, nomeie/salve (sua
preferncia) e execute:
O resultado do exemplo do tutorial apresentou uma diferena vertical (Z) de 1m comparando
o mesmo ponto aleatrio (Coordenadas X e Y, datum WGS84 UTM) no Google Earth.
Segue abaixo o resultado (aspecto grfico) da interpolao:
MDE (original) 12,5m de Resoluo MDE (interpolado) 1,0m de Resoluo

Espacial. Espacial.
Se voc quiser corrigir alguma diferena que for apresentado na sua rea, no intuito de
adequar a preciso, sugiro refazer o tutorial alterando o valor necessrio a partir do 5 Passo.
Espero que o tutorial tenha ajudado! Qualquer dvida ou sugesto, estou ao dispor. At a
prxima.
BIBLIOGRAFIA
BENEDITTI, C. A. Imagens multitemporais ALOS/PALSAR baseada em orientao a objeto na

discriminao da cobertura da terra. Programa de Ps-Graduao em Cincias Cartogrficas,
FCT/UNESP. Dissertao. Presidente Prudente, SP. Disponvel em:
http://repositorio.unesp.br/handle/11449/86778. Acesso em: 17 de dezembro de 2016.
ROSENQVIST, A.; SHIMADA, M. e WATANABE, M. ALOS PALSAR: Technical outline and mission
concepts. 4th International Symposium on Retrieval of Bio- and Geophysical Parameters from
SAR Data for Land ApplicationsInnsbruck, Austria, November 16-19, 2004. Disponvel em:
http://www.eorc.jaxa.jp/ALOS/kyoto/ref/ALOS_BioGeo-04.pdf. Acesso em: 16 de dezembro
de 2016.
SNCHEZ, G. A.; SHIMABUKURO, Y. E. e VALERIANO, D. de M. Processamento digital de imagens

de dupla-polarizao ALOS/PALSAR para a deteco de reas desmatadas na floresta
Amaznica. Anais XV Simpsio Brasileiro de Sensoriamento Remoto - SBSR, Curitiba, PR,
Brasil, 30 de abril a 05 de maio de 2011, INPE p.7525. Disponvel em:
http://www.dsr.inpe.br/sbsr2011/files/p0537.pdf. Acesso em: 27 de dezembro de 2016.

Tutorial MDE ALOS-1 Huascar Vidal v1

Enviado por

Direitos autorais:

Formatos disponíveis

Tutorial MDE ALOS-1 Huascar Vidal v1

Enviado por

Dados do documento

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Tutorial MDE ALOS-1 Huascar Vidal v1

Enviado por

Direitos autorais:

Formatos disponíveis

Tutorial para obteno, correes e interpolao (opcional) do Modelo

Digital de Elevao (sensor PALSAR) do Satlite ALOS-1 de 2006 a 2011.

Por Huscar P. Vidal de Oliveira

O satlite ALOS (Advanced Land Observing Satellite), chamado pelos japoneses de

1) Cadastro gratuito no geoportal de descarregamento de imagens de satlite da Alaska

2) Delimitao da rea de interesse no ecl (mapa-mndi) do geoportal:

Nos passos a seguir, so relacionados s correes de incongruncias resultantes desde

Note o resultado da converso,

7) Remover valores negativos mediante aplicao da ferramenta Con. Nesse exemplo

*No intuito de reduzir ainda mais o tempo de processamento do computador (durante

MDE (original) 12,5m de Resoluo MDE (interpolado) 1,0m de Resoluo

BENEDITTI, C. A. Imagens multitemporais ALOS/PALSAR baseada em orientao a objeto na

SNCHEZ, G. A.; SHIMABUKURO, Y. E. e VALERIANO, D. de M. Processamento digital de imagens

Você também pode gostar