O Quantum GIS combinado com as ferramentas do SAGA GIS gera um ganho significativo para as pessoas que trabalham com dados matriciais em ambiente de produção de dados. Nosso objetivo é utilizar estes aplicativos para manipulação dos recursos hídricos. Escrevi este tutorial para a Comunidade Quantum GIS Brasil a fim de realizar a instalação de todas as ferramentas necessárias para análises hidrológicas com software livre.

Acesse o tutorial completo através da Comunidade Quantum GIS Brasil: 

Quantum GIS 1.8: Instalação do SAGA GIS através da Ferramenta Sextante

Linux: Instalação do SAGA GIS no Quantum GIS através do Sextante

Procedimentos para instalação do SAGA GIS no Ubuntu 12.04

 1 - Acesse o Synaptic e instale o SAGA GIS fisicamente no computador (você deve marcar para instalação o saga, libsaga e o python-saga).

2 – Depois, com o Sextante instalado, basta acessar as configurações da ferramenta e marcar a opção Activate. Os algoritmos do SAGA GIS vão funcionar perfeitamente. Vou executar a mesma rotina para gerar um relevo sombreado com iluminação:

Terrain Analysis – Lighting – Analytical Hillshading

Resultado: http://goo.gl/amXg7

Utilizamos o SRTM TOPODATA como amostra para este tutorial. Em outras análises podemos efetuar testes com o ASTER GDEM.

Download do Tutorial

Escolha um dos formatos abaixo:

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

Clique nas opções abaixo para fazer o download do tutorial.

Download do Tutorial

Escolha um dos formatos abaixo:

Extração de Rede de Drenagem, Álgebra de Mapas e Delimitação de bacia Hidrográfica são os temas que serão abordados neste fantástico tutorial. Eu tentei construir esta dica de modo que o procedimento para extração de rede de drenagem possa ser reproduzido em outro SIG, pois as ferramentas são as mesmas.

A caixa de ferramentas Hydrology são algoritmos do Spatial Analyst para análise hidrológica no ArcGIS. Vamos utilizar essas ferramentas para extração de uma rede de drenagem a partir de um DEM SRTM.

Preparando os Dados

Download e Processamento do SRTM Topodata

No tutorial abaixo você pode obter uma grade vetorial contendo os links para download do DEM Topodata. O procedimento para obter este modelo de elevação processado pelo INPE está tudo muito bem documentado no tutorial abaixo:

Topodata – Modelo Digital de Elevação Processado pelo INPE

O DEM SRTM que vamos utilizar para extração da rede de drenagem precisa de um ajuste de projeção e reprojeção para o Sistema de Coordenadas Planas UTM de acordo com sua localização. O tipo de dado utilizado no tutorial é Unsigned Integer 16 bit (Números Inteiros Não Sinalizados de 16 bits) [1] [2].

Familiarize-se com alguns termos associados a Hidro

É importante tomar nota de alguns termos associados à análise hidrológica conhecendo os seus correspondentes em Português:

• Watershed = Bacia Hidrográfica
• Basin = Bacia Hidrográfica
• Drainage Network = Rede de Drenagem
• Stream Network = Rede de Drenagem
• Sub-Basin = Sub-Bacia
• Stream = Córrego, Riacho, Rio
• Upstream = Rio Acima (Em Geografia, “upstream” indica a fonte de um córrego ou rio contra a direção normal do fluxo de água)
• Downstream = Rio Abaixo

Tradução de algoritmos de Análise Hidrológica

Sink (Buraco, Depressão)
Cria um raster identificando todas as depressões, áreas de drenagem interna e vazios ou erros do SRTM ou outro raster de superfície.

Fill (Preencher)
Preenche pequenas imperfeições nos dados e remove todos os sinks do raster de superfície.

Flow Direction (Direção do Fluxo)
Gera um raster contendo a direção do fluxo de cada célula. No Flow Direction, cada pixel é potencialmente cercado por oito pixels vizinhos. A inclinação de cada uma destas oito direções pode ser calculada tomando-se a diferença em elevação indicada pelo valor de DEM para cada um destas oito localidades vizinhas e do valor no pixel a ser examinado.

Flow Accumulation (Acumulação de Fluxo)
O raster criado determina a acumulação de fluxo a partir do Flow Direction criado anteriormente.

Tutorial

Execute o ArcMap. A caixa de ferramentas Hydrology está localizada no menu Spatial Analyst – Hydrology. O pré-processamento que fizemos no Modelo de Elevação Topodata seguiu estes passos:

• O arquivo foi entregue com projeção/Datum arbitrária, sem referência espacial. Definimos a projeção para GCS_WGS1984;
• Reprojetamos o DEM com referência geográfica para WGS 1984 UTM Zone 22 S (esta é a projeção plana da região de Jaboticabal-SP, nossa área de estudo)
• Finalmente fizemos a conversão do tipo de dado de Float 32 bits para Unsigned Integer 16 bits.

Nosso DEM Topodata está pronto. Este PDF traz uma visão geral da região onde vamos executar nossas análises hidrólogicas no Spatial Analyst.

Flow Direction – Definindo as Direções dos Rios

Clique no menu Spatial Analyst – Hydrology – Flow Direction. Informe o DEM e o local de saída para o novo arquivo TIF:

Clique com o botão direito do mouse sobre o Flow Direction criado e selecione a opção Open Attribute Table. No campo Count da tabela, o Flow Direction deve conter apenas os valores 1, 2, 4, 8, 16, 32, 64 e 128 que representam as oito células adjacentes por onde a água deve fluir. A tabela de atributos do nosso raster possui muitos valores associados ao Flow Direction, o que significa que há vazios, depressões ou mesmo erro nos dados SRTM. Verifique o problema através desta imagem: | http://goo.gl/RfXeT

A ferramenta Sink cria um raster identificando possíveis erros nos dados SRTM.

Sink

Cria um raster identificando vazios e possíveis imperfeições no raster de superfície.

Clique no menu Spatial Analyst – Hydrology – Sinks

Informe o Flow Direction no campo de entrada e indique um arquivo TIF no campo de saída. Execute o processo.

Agora podemos identificar os problemas do raster deixando apenas o Sink habilitado. A ferramenta Fill deve ser utilizada para correção de todos os sinks.

Verifique o problema de sinks através desta imagem: | http://goo.gl/ZvPxc

Fill – Preenchimento de Vazios do DEM

A ferramenta Fill deve ser executada diretamente no DEM.

Clique no menu Spatial Analyst – Hydrology – Fill

Aguarde o término do processo. O raster Fill será adicionado no mapa.

Se o aspecto da imagem estiver totalmente cinza, acesse as Propriedades do Raster, guia Simbology e modifique o Stretch Type para Minimum/Maximum.

Visualização do DEM corrigido pela ferramenta Fill: | http://goo.gl/TAkHK

Com o novo DEM “Filled”, será preciso gerar um novo Flow Direction. Vamos organizar nosso ambiente de trabalho removendo do ArcMap o antigo DEM, o Flow Direction com problema e o próprio Sink. Só vamos manter no mapa o DEM preenchido pela ferramenta Fill.

Apenas o DEM corrigido pela ferramenta Fill deve permanecer no mapa: | http://goo.gl/Ew9a5

Execute novamente o algoritmo Flow Direction e aguarde o término do processo.

Este é o Flow Direction CORRETO. Tabela de valores cujo campo Cont registra valor o máximo 128:

Visualização do Flow Direction: | http://goo.gl/ETrrB

Flow Accumulation – Criando a Acumulação do Fluxo

Clique no menu Spatial Analyst – Hydrology – Flow Accumulation

Execute o processo Flow Accumulation. Nas opções, informe a entrada, a saída e modifique o tipo de dado para Inteiro (INTEGER).

O Flow Accumulation será criado, mas o raster pode parecer escuro. É preciso aproximar o mapa com a ferramenta Zoom para visualizar os locais onde ocorre o fluxo de água acumulada:

Visualização do Flow Accumulation: | http://goo.gl/BCaHd

Para ampliar a visualização, acesse as Propriedades do Raster. Na guia Simbology, marque a categoria Classified e selecione 8 classes. Modifique a rampa de cores para Azul Piscina – Lilás. Será necessário alterar a primeira cor da classificação para branco (clique sobre o retângulo). Clique no botão OK.

Nova Visualização do Flow Accumulation: | http://goo.gl/cvKdL

Os rios agora podem ser visualizados no mapa.

Álgebra de Mapas para Extração da Rede de Drenagem

Para criar uma rede de drenagem, é preciso estabelecer um limiar para descobrir quantos pixels adjacentes devem formar os rios. Em nosso caso, vamos estabelecer um limiar de 500 pixels.

Restaure a simbologia cinza do Flow Accumulation através das Propriedades do Raster, guia Symbology, categoria Stretch. Modifique o Stretch Type para Minimum/Maximum.

Formulação de Expressões Algébricas no SIG

A expressão matemática é criada com base em proposições que podem surgir no decorrer do projeto.  A função CON (CONDITIONAL) determina um teste lógico e retorna um resultado. Simples assim. Em relação à rede de drenagem, poderíamos estabelecer uma condição para criar o limiar, filtrando apenas os pixels significativos. A condicinal do limitar torna-se intuitiva quando é lida desta maneira:

Se o valor dos pixels do raster Flow Accumulation for superior a 500, deve-se substituí-los pelo valor 1.
Caso contrário, transformar em NODATA.

Com esta condição, todas as células do raster com valor superior a 500 serão consideradas durante a criação de um novo mapa raster. As demais células serão transformadas em NODATA. Esta é a questão lógica: o velho IF > ELSE. O próximo passo consiste em redigir essa informação na linguagem do computador (expressão matemática).

Clique no menu Spatial Analyst – Map Algebra – Raster Calculator

Digite a fórmula abaixo:

Con("03_FLOW_ACCUMULATION.tif" > 500,1)

Substitua “03_FLOW_ACCUMULATION.tif” pelo nome do seu raster que representa o Flow Accumulation. Informe um local de saída para a Rede de Drenagem.

A rede de drenagem será visualizada no mapa. Um aspecto interessante é sobre o limiar: ele deve ser utilizado com uma espécie de filtro para os rios. Dependendo da aplicação, o uso da hierarquia dos rios pode ser ampliada ou omitida. Os algoritmos Strahler e Shreve são mais específicos para classificar rios segundo a ordem. Os testes que realizei aqui demonstram que, se o limiar for aumentado, os riachos poderão ser omitidos. Segue um layout com a drenagem extraída de acordo com os limiares 500, 1000 e 2000. O maior limiar gera córregos mais genéricos.

Visualização dos Limiares 500, 1000 e 2000: | http://goo.gl/f4dlu

Stream to Feature – Drenagem para Shapefile

Clique no menu Spatial Analyst – Hydrology – Stream to Feature

Informe a rede de drenagem, o Flow Direction e o nome de saída para o arquivo shapefile de linhas.

A extração/vetorização da Drenagem, para a resolução espacial do DEM SRTM, ficou excepcional. Muito bom mesmo.

Alta resolução: http://goo.gl/4krQ9

Basin – Bacia Hidrográfica

Spatial Analyst – Hydrology – Basin

Último processo. Informe o Flow Direction e o arquivo TIF de saída.

A bacia será criada. Acesse as Propriedades do Raster e modifique a simbologia e classes como foi feito anteriormente.

Resultado Final

Bacias Hidrográficas e Rede de Drenagem (PDF)

Agora posso dar prosseguimento ao meu trabalho em Jaboticabal-SP, pois eu tenho a delimitação perfeita da bacia hidrográfica daquela região.

Que tal um desafio? É possível reproduzir essa mesma dica em outro SIG, pois agora você está um pouco mais familiarizado com os algoritmos. A ordem para gerar a bacia é esta:

DEM UTM - FILL DEM UTM -FLOW DIRECTION – FLOW ACCUMULATION – STREAM NETWORK – BASIN

Boa sorte! Jorge Santos

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

O PostGIS pode ser definido dos maiores aplicativos livres para análise de banco de dados espaciais da atualidade . Tive que instalá-lo para para cumprir um pequena tarefa e estou repassando o tutorial para que outra pessoa possa alcançar o mesmo objetivo.

Roteiro de instalação

Em nosso roteiro de instalação, vamos utilizar a versão 9.1 do PostgreSQL combinada com a versão 2.0 do PostGIS. O PostgreSQL 9.1 não possui versão 64 bits. A sequência das instalações é sistemática e deve ser seguida de acordo com a ordem abaixo:

• Instalação do PostgreSQL 9.1
• Instalação do PostGIS 2.0
• Conexão e Gerenciamento da Base de Dados com pgAdmin III.

Clique nas opções abaixo abaixo para baixar o tutorial.

Download do Tutorial

Escolha um dos formatos abaixo:

Outros tutoriais para o PostGIS podem ser acessados nos sites do Anderson Medeiros e Geoprocessamento.net.

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

O ERDAS 2013 está repleto de novidades. Hoje vamos aprender como utilizar as bases OpenStreetMap e Bing Maps.

Adicionando a Base OpenStreetMap

Tenho alguns remanescentes de vegetação no formato shapefile e quero validar o posicionamento dessas informações.  Estou com os polígonos abertos no ERDAS 2013.

Na guia Home, clique na categoria Basemap e selecione a base de mapas OpenStreetMap:

Aguarde o carregamento das imagens. As estradas da base OpenStreetMap serão exibidas no mapa:

Visualização em Tela Cheia: http://goo.gl/Vas9S

A base OSM deve apresentar serrilhados por causa do posicionamento on-the-fly das camadas. Mude o método de reamostragem dos pixels para suavizar o serrilhado. Na categoria Basemap existe um botão Options. Modifique a reamostragem para Cubic Convolution e o problema será resolvido:

Adicionando a Base Bing Maps

A cobertura Bing Maps inclui três mapas básicos:  aéreo, estrada e mapas híbridos. O uso gratuito do serviço Bing Maps requer autenticação e o acesso pode ser obtido através de um cadastro simples. Todo usuário que deseja ter acesso às imagens do Bing Maps vai precisar de uma chave obtida diretamente da Microsoft para ter acesso a base de mapas. Para mais informações, consulte o licenciamento do Bing Maps.

Cadastro no site do Bing Maps

Visite este link http://www.bingmapsportal.com/  e faça login com sua conta da Microsoft no campo Sign In. Vou utilizar minha conta do Microsoft Outlook.

Há dois tipos de licenciamento: Trial e Basic. Cada conta só pode utilizar duas chaves apenas. A licença Trial permite acessar a base Bing Maps por 90 dias, sendo que a chave não pode ser excluída nem modificada. Vou criar os dois tipos de licença para esta conta:

No cadastro, informe o nome da aplicação, o tipo de licença e a finalidade da aplicação. Eu cadastrei as aplicações como ERDAS 2013 GIS Web, indiquei as licenças Trial/Basic  para cada uma delas e  informei a finalidade como Educacional e Site Público. Digite o captcha e envie o formulário através do botão Submit. As chaves aparecem na parte de baixo e você ter acesso a elas quando quiser.

 Autenticação no ERDAS 2013

 Na guia Home, clique na categoria Basemap e selecione a base de mapas Bing Aerial. Digite a chave na janela de autenticação e clique no botão Test. Em seguida, clique no botão OK:

Aguarde o carregamento das imagens de satélite e faça bom uso da sensacional base de mapas Bing Maps.

Visualização em Tela Cheia: http://goo.gl/WAA2Z

Algumas áreas possuem cobertura de imagens em Alta Resolução. Experimente aproximar a visualização coma ferramenta Zoom e novas imagens poderão surgir (Escala 1:10.000):

Visualização em Tela Cheia: http://goo.gl/GKOL6

Esta dica também funciona no ArcGIS 10.x pois a interface Bing Maps foi removida de todos os produtos ESRI. Um janela de autenticação deve surgir na sua tela durante a utilização do SIG. |Nos dois tipos de licenciamento eu pude verificar uma área no nível de Zoom nas imagens acima. Recomendo o uso da licença Basic. A Microsoft deseja ter um controle maior sobre a utilização dos seus dados e devemos reconhecer que neste ponto a empresa tem razão.

Compartilhe este post por ai!

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

Sobre o uDig

uDig é um aplicativo SIG Open Source orientado para a Web onde você pode descarregar suas camadas WMS, WFS, WPS e outros Web Services. O program possui uma Interface Gráfica de Usuário bastante amigável para trabalhar em ambiente Desktop.

Site Oficial: http://udig.refractions.net/

Sem sombra de dúvida, a maior fonte de tutoriais em Português-Br para o uDig  é o blog do Anderson Medeiros:

• Conheça Melhor o Udig
• Dando os Primeiros Passos no uDig
• Integração uDig com Base de Dados PostGIS
• Gerar Mapas Temáticos no uDig
• Geração de Estilos SLD com o uDig
• Definir Projeção de Shapefile no uDig

Download do uDig

Baixe o programa aqui: http://udig.refractions.net/download/

Conheça a Plataforma GXT

Escrito em Java com base na biblioteca GeoTools, GXT é uma plataforma para análise de dados espaciais com diversas ferramentas para Geoprocessamento, Estatística Espacial e  Análise Estatística de Algoritmos.  

Site Oficial: http://www.mangosystem.com:8080/gxt/docs/index.html

Integração uDig com a Plataforma GXT

O aplicativo uDig ganha musculatura com a utilização da plataforma GXT de tal forma que o resultado de suas análises espaciais pode ser comparado até mesmo com o ArcGIS. Com a instalação dessa sensacional ferramenta, é possível criar Slope (Encosta/Declive), Hill Shade (Relevo Sombreado), Buffer (Área de Influência) e outras centenas de algoritmos para análise espacial com base em vetor e raster.

Instalação da Plugin GXT no uDig

Faça o download do uDig e instale o programa. Clique no menu Help – Find and Install:

Marque na opção Search for new features to install e pressione o botão Next:

Clique no botão New Remote Site e digite os dados abaixo:

Name: GXT Toolbox for uDig
URL: http://1.234.82.19:8080/gxt/udig/

O servidor será adicionado. Pressione o botão Finish.

Após estabelecer a conexão, confirme as telas seguintes e o plugin será instalado. Para ativá-lo, clique no menu Window – Show View – Other. Marque o plugin GXT na janela Show View. Os comandos estarão disponíveis no menu direto da janela:

Agora podemos executar os algoritmos do plugin GXT sempre que for necessário. Mais ferramentas tornam o SIG uDig atrativo. Vou carregar um DEM e gerar um Slope para testar o processamento no uDig:

Visualização do uDig em  Tela Cheia: http://goo.gl/PXh5E

Eis o nosso Slope: http://goo.gl/0plE0. Este terreno é consideravelmente baixo, porém a simbologia em vermelho indica as partes mais inclinadas na imagem. O tempo de processamento está diretamente associado ao tamanho da imagem e a capacidade do equipamento disponível.

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

Clique nas opções abaixo para fazer o download do tutorial.

Download do Tutorial

Escolha um dos formatos abaixo:

Conversão de Bandas de 16 bits para 8 bits através da Matemática de Bandas

Nível de Usuário: Intermediário

Mudanças na resolução radiométrica de imagens de satélite podem ser realizadas no Quantum GIS através da Álgebra de Mapas.

Neste tutorial vamos conferir a fórmula utilizada para conversão de um raster de 16 bits para 8 bits. Este documento traz uma lista de todos os intervalos e suas possíveis aplicações.

Em PDI (Processamento Digital de Imagens), essa tarefa é conhecida como Rescale.

Pré-Requisitos

• Imagem de satélite com Resolução Radiométrica de 16 bits;
• Aplicativo SIG utilizado: Quantum GIS.

Conversão de Bandas

O raster aberto no Quantum GIS é uma imagem de satélite Ikonos bruta (sem nenhum processamento), com três bandas e com tamanho de pixel de 1 metro (se a sua imagem estiver escura, acesse as Propriedades da Camada, guia Estilo. Modifique a opção Melhorar Contraste para Estender para MinMax:

Nas Propriedades da Camada, temos informações importantes sobre a imagem na guia Metadata. Para verificar a quantidade de níveis de cinza em cada banda, acesse as Propriedades da Camada, guia Estilo. No Campo Carregar Valores Min/Max da Banda, clique na opção Real e pressione o botão Carrega. Os valores Min/Max de cada banda serão exibidos. Guarde estes valores.

Fórmula

• Num Range de 16 bits, a escala é de 0-65535 níveis de cinza
• Num Range de 8 bits, a escala é de 0-255 níveis de cinza

Usaremos a seguinte equação para transformar cada banda de 16 bits para 8 bits:

( ( Raster - Valor Min Raster)*( Valor Max da Escala - Valor Min da Escala )/
( Valor Max Raster - Valor Min Raster ) ) + Valor Min da Escala

Execute o plugin RasterCalc (para instalar essa calculadora, basta clicar no menu Complementos – Buscar Complementos Python). Digite a fórmula abaixo para processar a primeira banda:

( ( [po_235348_rgb_0000000]@1  - 0)*( 255-0 )/( 2047 - 0 ) ) + 0

Selecione a opção Byte (Inteiro Não Sinalizado de 8 bits) e indique um local de saída para a banda 1. Clique no botão Calculate:

Segunda fórmula

( ( [po_235348_rgb_0000000]@2 - 0)*( 255-0 )/( 2047 - 0 ) ) + 0

Pressione no botão Calculate para processar a Banda 2 e modifique o nome do arquivo de saída. Há poucas alterações na fórmula em relação à expressão anterior.

Terceira Fórmula

( ( [po_235348_rgb_0000000]@3 - 0)*( 255-0 )/( 1982 - 0 ) ) + 0

Atenção: aqui o Valor Máximo da Banda 3 é 1982.  Clique no botão Fechar para sair da Calculadora Raster.

Camada Acumulada (Layer Stack)

Para verificar os resultados, faça um stack das bandas através do menu Raster – Miscelânea – Mosaico. Ajuste as opções abaixo:

• No campo Arquivos de Entrada, selecione as três bandas que foram criadas através do RasterCalc;
• No campo Arquivo de saída, selecione o diretório de saída do arquivo GeoTIFF;
• Marque o campo Camada Acumulada (Layer Stack);
• Marque o campo Opções de Criação e clique no botão Adicionar. digite as strings TFW  e YES nos campos apropriados.

Uma nova imagem em 8 bits será criada. Experimente abrir a imagem num aplicativo que não suporta dados em 16 bits como o Photoshop. O dado será visualizado perfeitamente:

Imagem colorida aberta no programa Adobe Photoshop

Agora você pode utilizar o Quantum GIS para realizar a Transformação Radiométrica. O processo pode ser executado no GDAL, GRASS, ENVI e muitos outros. A ESRI preparou um help para o ArcGIS que contém outras composições que podem ser realizadas através da Calculadora Raster. Experimente e compartilhe conosco o seu aprendizado.

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

Clique nas opções abaixo para fazer o download do tutorial.

Download do Tutorial

Escolha um dos formatos abaixo:

Informações Específicas

RapidEye é uma companhia que dispõe de uma constelação de cinco micro satélites que mapeiam sistematicamente a cobertura terreste.

O produto RapidEye possui Resolução Espacial de 5 metros ( Média-Alta Resolução).

Aquisição de Imagens RapidEye

As imagens RapidEye são ortorretificadas e podem ser adquiridas através da empresa AMS Kepler Engenharia de Sistemas ( conheça nossos parceiros). A empresa AMS Kepler disponibiliza dezenas de amostras de imagens RapidEye georreferenciadas, em 16 ou 08 bits (Resolução Radiométrica) que podem ser baixadas de acordo com a aplicação desejada (Área Rural, Área Urbana, Floresta, Ferrovia, Pastagem, Transportes, Mineração, Óleo e Gás, Áreas Litorâneas, etc).

Amostras RapidEye nas composições RGB True Color, False Color e Red Edge

As imagens RapidEye podem são comercializadas por telha (tile) ou por quilômetro quadrado. Entre em contato com a empresa AMS Kepler para obter maiores informações.

Pesquisa ao acervo de Imagens RapidEye

O Mecanismo de Busca EyeFind

Nível de Usuário: Intermediário

Para verificar a cobertura em acervo de imagens RapidEye sobre uma determinada área de interesse, utilize o mecanismo de busca conhecido como EyeFind.

http://eyefind.rapideye.com/

Submetendo uma Área de Interesse

Assumindo que você possui um conhecimento médio sobre SIG, desenhe sua área de interesse e exporte-a para KML/KMZ. O mecanismo de busca EyeFind suporta cinco tipos de formatos de área de interesse, os principais métodos são por Identificador de Telha (Tile ID) e  Upload de KMZ/KML.

Pesquisa RapidEye por Tile ID

Este método de pesquisa é bem eficiente se você possui as telhas do satélite RapidEye. 

• Grade de Cobertura RapidEye no Brasil
• Grade de Cobertura RapidEye Global

Carregue a grade no SIG e verifique quantas telhas RapidEye são necessárias para cobrir a sua área de interesse. Selecionamos o município de Aparecida do Taboado-MS: 

A cobertura deste município demanda pelo menos nove tiles RapidEye para ser totalmente imageada:

Anote os identificadores de cena ou exporte a Tabela de Atributos para o Excel e faça a transposição das células (Colar Especial – Transpor). Copie os dados do Excel para o Notepad e substitua os espaços pela vírgula. Os IDs devem ser separados por vírgula para que o mecanismo de busca do RapidEye aceite-os.

 Transposição de dados no Excel 2010

Substituição do Espaço pela Vírgula no Notepad

Pesquisa com EyeFind

Acesse o site EyeFind e marque as opções abaixo:

• Na guia Query, especifique o intervalo de tempo;
• O produto sempre será Level 3A (Ortorretificado);
• Em Area of Interest, selecione RapidEye IDs;
• Em Area of Interest, selecione Tile IDs.

Toda cobertura RapidEye de 01/01/2012 até 23/04/2013 será exibida no mapa. No lado esquerdo da tela são exibidas imagens em miniatura. Essas miniaturas das imagens recebem o nome de Quicklooks ou QuickViews.  No lado direito temos os Quickviews sobre o mapa. Repare que  EyeFind tentou cobrir a nossa área com as imagens mais recentes, porém algumas cenas estão com uma alta cobertura de nuvens:

As imagens também podem ser selecionadas de acordo com o período das chuvas e/ou da seca, mas isso depende da aplicação. A partir de agora você pode realizar a triagem das cenas, selecionando as imagens com um mínimo de nuvens possível.

Para verificar se existem imagens melhores, deve-se clicar duas vezes com o botão esquerdo do mouse sobre qualquer imagem coberta de nuvens e aguardar o carregamento da cena na relação à esquerda para desmarcá-la e selecionar outra cena melhor. Siga os passos abaixo para concluir esta etapa:

• Clique duas vezes com o botão esquerdo do mouse sobre a cena coberta de nuvens (1);
• Aguarde a exibição do Quickview no painel esquerdo (2);
• Desmarque o Quickview coberto de nuvens (3);
• Marque o próximo Quickview da lista com melhor visualização (4);
• Pressione o botão Update map (5).

As melhores cenas RapidEye neste intervalo de tempo são estas:

Consultas podem ser simples ou complexas. Podem envolver cobertura total, ângulo de elevação de cada cena, mínimo de nuvens, resolução temporal, etc. O critério de pesquisa quase sempre está associado a finalidade do projeto. Esta informação pode ser adquirida juntoao cliente final.

Download de Quickviews Georreferenciados

Para realizar o download dos Quickviews selecionados na pesquisa, clique no botão Download. Selecione as opções como mostra a imagem abaixo:

Após baixar todos os arquivos, decompacte o ZIP e carregue os Quickviews no SIG:

Os Quicklooks agora podem ser integrados ao GIS caso seja necessário realizar uma análise mais complexa. Eles não substituem a imagem, pois estão em baixa resolução e não possuem os parâmetros da imagem original. Ao terminar, encaminhe o shapefile da pesquisa para a empresa responsável pelo fornecimento das imagens. Nossa opção é pela a empresa AMS Kepler Engenharia de Sistemas.

Pesquisa RapidEye via KMZ/KML

O procedimento para pesquisa por KMZ/KML é muito mais simples. Clique no campo File Upload, localize o arquivo e clique no botão Query:

Este é o aspecto da área de interesse na consulta RapidEye. Repita as técnicas citadas acima para obter melhores imagens e envie o shapefile com o resultado para empresa fornecedora das imagens.

Resultado Final

Fim do tutorial. Este é o primeiro post sobre uma série que estou preparando sobre Processamento Digital de Imagens (PDI). Antes de entrar no tema principal, vamos conhecer os mecanismos de pesquisa de imagens de satélite. Continue acompanhado e recomende este post para seus amigos.

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

Clique nas opções abaixo para fazer o download do tutorial.

Download do Tutorial

Escolha um dos formatos abaixo:

ArcGIS 10: Configuração de Camadas Selecionáveis (Set Selectable Layers)

A seleção de feições é um dos recursos mais elementares durante a edição de objetos vetoriais no aplicativo SIG (Sistema de Informação Geográfica) ou GIS (Geographic Information System), pois permite que o operador destaque partes das entidades de ponto, linha ou polígono que serão editadas. Selecionar um dado é uma tarefa intuitiva e vale para qualquer aplicativo que você trabalha no computador. Se você seleciona um dado é porque em seguida seu próximo passo será recortar, copiar, colar, excluir ou acessar as propriedades desse dado. Este modo de atuação é satisfatório para qualquer ação realizada no computador.

A lógica do SIG consiste em gerenciar informação do espaço através de camadas sobrepostas (o espaço aqui mencionado pode ser considerado como o palco teatral, o lugar onde as grandes apresentações acontecem). Sobre uma mesma localidade do espaço é permitida a manipulação de entes que representam o terreno, a unidade da federação, o rio, a rodovia, etc. Geralmente estes elementos são logicamente distribuídos sobre o mapa de acordo com a importância de cada um. Para ampliar nosso conhecimento sobre as entidades que por nós serão editadas, vamos conhecer as principais ferramentas de seleção do SIG ArcGIS.

Ferramentas de Seleção do ArcGIS

Na instalação padrão, o ArcGIS fornece duas ferramentas básicas para manipulação de objetos através de seleções: Select Features (Selecionar Feições) e Clear Selected Features (Limpar Feições Selecionadas). Feições são as geometrias de ponto, linha ou polígono distribuídas no projeto do ArcGIS. As ferramentas podem ser encontradas na barra de ferramentas Tools:

A ferramenta de seleção Select Features possui algumas variedades:

• Seleção por Retângulo (Select By Retangle): é a ferramenta de seleção principal, utilizada em 95% dos casos;
• Seleção por Polígono (Select by Polygon): cria uma seleção irregular em forma de um polígono;
• Seleção por Laço (Select by Lasso): cria uma seleção em forma de um laço;
• Seleção por Círculo (Select by Circle): cria uma seleção em forma de um círculo;
• Seleção por Linha (Select by Line): cria uma seleção em forma de uma linha.

Ferramenta Zoom to Selected Features

Permite a aproximação das feições selecionadas. Útil para destacar na visualização um determinado trecho da seleção. A ferramenta Zoom to Selected Features (Aproximar Feições Selecionadas) é uma ferramenta de seleção do ArcGIS de grande valia. Esta ferramenta está presente na Tabela de Atributos do arquivo shapefile.

A instalação da ferramenta Zoom to Selected Features é importante para otimizar a edição no SIG.

Instalação da Ferramenta Zoom do Selected Features

Com o ArcMap aberto, clique no menu Customize – Toolbars – Customize…

Na janela Customize, clique na guia Commands e aguarde o carregamento das categorias. Na categoria Selection estão todas as ferramentas de seleção do programa. Deslize a barra de rolagem para baixo e localize a ferramenta Zoom to Selected Features:

A instalação é simples: clique e pressione com o botão esquerdo do mouse sobre a ferramenta e arraste-a para o menu superior (perceba que o cursor do mouse assume outra forma quando você pressiona e arrasta a ferramenta):

Arraste e solte a ferramenta Zoom to Selected Features ao lado das duas ferramentas de seleção padrão (Select Features e Clear Selected Features). Feche a janela Customize (Personalizar). Agora nós temos três ferramentas principais para seleção de feições:

Na tabela exibida mais acima, perceba que o estado do Maranhão estava selecionado. Ao pressionar a ferramenta Zoom to Selected Features, a visualização do mapa será direcionada para esta unidade da federação e para todos os arquivos shapefile que estiverem selecionados, pois estes elementos devem possuir informações selecionadas na tabela e na geometria:

Ampliar Visualização

Conheça a função Set Selectable Layers

Quando muitos objetos estão distribuídos no projeto do SIG – neste caso, o programa que estamos tratando é o ArcGIS – é comum desenhar um retângulo com a ferramenta Select Features sobre um determinado elemento. O problema é que o usuário acaba por selecionar outras entidades ao arrastar a ferramenta de seleção sobre o mapa e todos os objetos sobrepostos serão selecionados. Veja as situações abaixo onde pode ocorrer o problema de múltiplas seleções indesejadas:

Esta questão poderia facilmente ser solucionada se o operador desligasse todas as camadas, exceto a camada das rodovias. Somente esta camada das rodovias seria selecionada sobre o trecho desejado, e assim poderíamos ativar as camadas restantes:

Em projetos básicos, essa técnica daria certo, Infelizmente não funcionaria em projetos mais avançados com centenas de camadas sobrepostas, pois a seleção pode demandar a visualização do objeto que está por baixo. A forma correta para gerenciar seleções entre diversas camadas no ArcGIS é a utilização das Configurações de Camadas Selecionáveis (Set Selectable Layers).

Este recurso também está disponível na versão 9.3 do SIG ArcGIS, porém ele foi removido da versão 10 do programa. O procedimento para restaurar o Set Selectable Layers é a personalização de ferramentas, da mesma forma que executamos para adicionar a ferramenta Zoom to Selected Features:

• Clique no menu Customize – Toolbars – Customize…
• Clique na guia Commands e procure pela categoria Selection;
• Identifique a ferramenta Set Selectable Layers e arraste-a para o menu superior.

O procedimento é idêntico até aqui, mas tente arrastar essa ferramenta para o menu Selection. Este é o mesmo posicionamento presente nas versões anteriores do ArcGIS. Veja a imagem abaixo:

Ao terminar, encerre a janela Customize clicando no botão Close. 

Como funciona a Configuração de Camadas Selecionáveis (Set Selectable Layers)

• Clique no menu Selection – Set Selectable Layers;
• Clique no botão Clear All (Limpar Tudo);
• Marque apenas a camada que será selecionada.

A partir de agora você pode arrastar a seleção à vontade, somente a camada selecionada no menu anterior estará apta para receber a seleção, as outras servirão como orientação até que seja possível restaurá-las por intermédio do menu Set Selectable Layers.

Esta é uma ferramenta muito útil para o trabalho diário. Envie esta dica para os seus amigos!

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

Clique em uma das opções abaixo para fazer o download do tutorial.

Download do Tutorial

Escolha um dos formatos abaixo:

Buffer Wizard é um recurso opcional para criação de Buffer cujo objetivo é a geração ou ampliação de uma área de influência. Muitas pessoas desconhecem essa ferramenta porque costumam executar o geoprocesso Buffer que está localizado nas caixas de ferramentas do ArcToolBox.

Esta ferramenta também está disponível na versão 9.3 do SIG ArcGIS.

Instalação da Ferramenta Buffer Wizard

Abra o ArcMap e clique no menu Customize – Toolbars – Customize…

Na janela Customize, clique na guia Commands. Digite “Buffer Wizard” sem aspas. O aplicativo será exibido no campo Commands:

Clique e pressione com o botão esquerdo do mouse para arrastar o aplicativo para o menu superior:

Arraste o Buffer Wizard para uma barra de ferramentas existente e de fácil acesso. Vou fixar a ferramenta ao lado do botão Add Data, pois este é um local de fácil recordação da posição do Buffer Wizard:

O comando Buffer Wizard foi adicionado ao ArcMap. Clique no botão Close para fechar a janela Customize.

Tutorial: Criação de Buffer

Antes da criação do buffer, modifique a projeção do Dataframe para o Sistema de Coordenadas Planas UTM. Recomendo manter o shapefile na mesma projeção, pois vamos trabalhar com unidades métricas. O sistema de c0ordenadas da minha região é SAD 1969 UTM Zone 25 S:

Abaixo temos um exemplo de um traçado cuja extensão atravessa dois estados do território brasileiro.

[1257x783 Pixels]

Exercício para criação de um Buffer de 5 km sobre um traçado

Buffer sobre linhas

Uma determinada organização precisa adquirir imagens de satélite ou aerolevantamento para realizar um licenciamento ambiental que permite o mapeamento e análise dos impactos ambientais causados pela obra. Seu superior incumbiu você de criar uma área de abrangência de 10 km de largura em toda a extensão do traçado da linha.

O Buffer é o geoprocesso mais indicado para criação de feições do tipo corredores (rios, rodovias, ferrovias, linhas de transmissão de energia e assim por diante). Execute o Buffer Wizard e selecione o arquivo de linha na primeira tela. Clique em Avançar:

Buffer Distance: modifique a unidade de quilômetros para metros.
At a specified distance: digite o valor 5000. Para um buffer de 10 Km, 5 km devem ser definidos para cada lado da linha.
Based on distance from an attribute: permite especificar um valor para a distância com base num campo da tabela de atributos. Opção pouco utilizada.
As Multiple Buffer Rings: você pode criar múltiplos anéis sobre o buffer. Não vamos utilizar essa opção para esse trabalho.

Clique em Avançar. Novas opções devem ser definidas adiante.

Buffer Output Type: marque a opção Yes para dissolver as barreiras do buffer. Isso irá criar uma área contínua.
Create Buffer So They Are: Para linhas, essa opção não está disponível.
Where do You Want The Buffer to be Save: O assistente solicita um local para salvar o novo buffer.

Clique no botão Concluir. A área de influência será criada no ArcMap. Organize as camadas, subindo o traçado acima do Buffer. Este é o aspecto da área de abrangência:

[1306x824 Pixels]

Usando a ferramenta de medição, entendemos que o Buffer de 10 Km possui 5000 m para cada lado da faixa principal, como pode ser visto através da ilustração abaixo:

Buffer sobre Polígonos

As mesmas opções citadas acima são válidas, exceto:

Create Buffer So They Are: Para polígonos, selecione uma das opções para criar um Buffer Dentro (Inside), Fora (Outside), Dentro e Fora (Inside and Outside) e Envolvendo toda a Área (Outside an Include Inside).

Vamos optar pela opção que envolve toda a área, pois essa opção se adequa perfeitamente à nossa pesquisa. Veja algumas situações possíveis que podem ser criadas pelo Assistente de Buffer:

Além do Assistente para Criação de Buffer do ArcGIS 10, você pode utilizar o menu Geoprocessing – Buffer para criar suas áreas de influência. A escolha da melhor ferramenta fica por sua conta.

Resultado Final

Buffer sobre área, distâncias variadas

[1387x802 Pixels]

Bom trabalho!

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br

OSM2GIS é um aplicativo Web para conversão de dados OpenStreetMap para shapefile ou KML.

Acesse o site do projeto:  http://www.osm974.re/osm2gis/

Tutorial

1 – O processo consiste em navegar pelo mapa até encontrar sua região de interesse. Meu local selecionado é uma fração de uma localidade no estado do Rio de Janeiro:

2 – Clique no botão Give Me My Data! 

3 – Digite seu e-mail e marque a opção KML. O Sistema de Coordenadas padrão deve ser GCS_WGS1984, portanto, marque o código EPSG 4326 e clique no botão SEND:

4 – Verifique seu e-mail. Um link para download estará esperando na sua Caixa de Entrada. Clique nesse link para baixar os arquivos.

5 – Pronto, agora você tem o arruamento, locais de interesse e outras feições da área. Você vai precisar de um SIG para visualizar as ruas e locais. Vou executar o Quantum GIS:

[Tela Cheia - 1366x728]

6 – Podem ser adicionadas no QGIS imagens do Google Earth e a própria base OpenStreetMap . O procedimento para visualizar imagens do Google Earth, Bing Maps e outros serviços no Quantum GIS está descrito neste post.

[Tela Cheia - 1366x728]

7 – Não esqueça do arquivo KML: você pode abrir o arruamento selecionado no próprio Google Earth:

[Tela Cheia - 976x678]

Bom trabalho!

Dúvidas e sugestões podem ser enviadas para o e-mail contato@processamentodigital.com.br