Geotecnologias Luís Lopes

2012: ano de mudanças! Gostaria de opiniões...

2012-01-16T08:00:00.001-03:00

Olá pessoal,

2012, ano que promete fortes emoções! Como alguns sabem (pelo twitter @geoluislopes), já estou morando em minha nova cidade: Macapá, capital do Amapá (moro no hemisfério norte do planeta rs), ocasião esta, assumi o cargo efetivo de Especialista em Geoprocessamento do Estado.

Espero que tenham um pouco de paciência, pois aqui (acreditem!) não chegou a banda larga ainda. O Governador prometeu para este ano...vamos ver!

Fim de 2011, disponibilizei o Arquivo do blog que está separado por mês e dia da postagem, para facilitar a busca de algum artigo publicado aqui.

Confira o Arquivo do blog.

Tivemos também postagens referentes ao projeto aqui implementado: o Seja Parceiro. Esta iniciativa consiste dar espaços para divulgar artigos, tutoriais, materiais de pessoas/profissionais da área de geotecnologias.

Divulgue seu material através do Seja Parceiro, veja as regras.

Iniciamos um "movimento" em prol do uso de software livre em geoprocessamento. Trata-se da "Série Quabrando tabus: software livre também faz!", utilizando a hashtag nas redes sociais #SLGeoTbFaz. Vale destacar que esta ideia contou como apoio dos GeoBlogueiros Eliazer Kosciuk (@eliazerk), Sadeck (@sadeckgeo), Anderson Medeiros (@ClickGeo), Esdras Andrade (@geoparalinux) e José Carlos (@geojcarlos), com publicações em seus respectivos blogs. A Série também recebeu uma nota na 3ª Edição da Revista FOSSGIS Brasil (@fossgis).

Primeira postagem de 2012, anunciamos mais um colunista para o blog: Christian Nunes.

Essa postagem é para "ouvir" você, leitor!

O blog está correspondendo a sua demanda?
Quer ver mais um colunista de alguma área mais específica?
Percebe-se a evolução do blog?
O que acharam da iniciativa do ano passado "Quebrando tabus: #SLGeoTbFaz?

São algumas questões, citadas acima, que gostaríamos de saber para melhorar cada vez mais o blog.

Bom 2012 e aguardo seus comentários!

Apresentação: coluna de Christian Nunes

2012-01-09T08:00:00.023-03:00

Olá companheiros...

Estou iniciando hoje uma pequena contribuição, no blog do Luis Lopes, sobre a temática da cartografia e de como ela vem sendo entendida no espaço acadêmico e profissional. Inicialmente vou relatar minha experiência na área.

Meu primeiro contato com a cartografia na universidade foi no ano de 1998, quando entrei no curso de licenciatura e bacharelado em Geografia na Faculdade de Geografia e Cartografia da Universidade Federal do Pará, no campus do Guamá, em Belém, estado do Pará. Naqueles anos da graduação tive uma aproximação muito tímida com os estudos cartográficos, pois as disciplinas somente demonstravam a parte teórica da cartografia, não tínhamos um laboratório e nem equipamentos onde os alunos pudessem trabalhar com softwares de geoprocessamento e de sensoriamento remoto, desse modo, os estudos eram muito superficiais, sem uma carga-horária prática. Felizmente, hoje, o curso de geografia da UFPA já oferece um instrumental em geotecnologias que antes não dispunha.

A aproximação com a cartografia começou a se estreitar em minha vida no ano de 2004, depois que me formei e tive oportunidade de trabalhar como consultor ambiental da Secretaria de Meio Ambiente do município de Inhangapí, no estado do Pará. Naquele momento, que coincidiu com minha entrada em um curso de especialização em Gestão Ambiental, também na UFPA, aprofundei minha leituras sobre a cartografia e os principais instrumentos de gestão e ordenamento territorial (EIA/RIMA, Cadastro Técnico Multifinalitário, Unidades de Conservação, ZEE, etc) e comecei a entender, na prática, onde eu poderia utilizar os conceitos e categorias da geografia e da cartografia que tinha aprendido na graduação.

fonte: FFLCH, USP

No ano de 2005 tive que me ausentar do trabalho como consultor e ingressar no mestrado em Geografia na UFPA. Sendo que, quando retornei ao curso de geografia, já tínhamos um aparato tecnológico que oferecia o suficiente para eu aprofundar meus conhecimentos, agora de forma prática, na cartografia. Com ajuda de colegas da graduação que atuavam como monitores do Laboratório de Análise da Informação Geográfica (LAIG), com destaque para o amigo Luis Sadeck, pude conhecer ainda mais as importantes ferramentas disponíveis para a cartografia dos dias de hoje.

No ano de 2006 ingressei, por um concurso de cargo temporário, como analista intelectual no Sistema de Proteção da Amazônia – SIPAM, no Centro Regional de Belém, onde, entre outras funções, tive que elaborar diversos mapas e ministrar cursos de softwares de geoprocessamento em vários municípios da Amazônia. A experiência no SIPAM durou três anos e meio e me possibilitou um conhecimento significativo/prático do que podemos extrair e manipular da informação espacial em meio computacional ou não, além do que, a experiência no SIPAM também me aproximou de pessoas que se preocupam a disponibilizar/democratizar a informação e o aprendizado sobre a cartografia, como o colega Luís Lopes.

Em setembro de 2009 fui aprovado e convocado para exercer o cargo de Professor da graduação em Geografia da UFPA, onde ministro disciplinas ligadas à cartografia, geoprocessamento e sensoriamento remoto, além de coordenar o LAIG e disponibilizar a diversas pessoas – alunos da UFPA ou não, cursos gratuitos na área de cartografia. Também na UFPA venho desenvolvendo, ou participando, de projetos e atividades em que o uso da cartografia é imprescindível, enfocando não somente a parte técnica de manipulação das geotecnologias, mas também analisando os conceitos e categorias importantes para se entender onde a cartografia pode ser aplicada.

Minha coluna será direcionada a alunos de cursos de graduação ou profissionais que estão iniciando seus estudos em cartografia. Os textos terão como principal objetivo discutir sobre a importância do uso das ferramentas de cartografia – mais precisamente suas geotecnologias, nos mais variados setores da atividade humana, enfocando, principalmente a educação cartográfica e a ocupação humana do espaço geográfico. Para isso, faremos uma análise inicial sobre a importância da ciência cartográfica desde seus primórdios, buscando elucidar o leitor que está iniciando suas leituras de forma bem básica, mostrando os principais tipos de representação e da linguagem da cartografia para posteriormente demonstrar, na pratica estudos onde a cartografia e suas ferramentas são analisadas e podem subsidiar os trabalhos dos leitores. Espero contribuir na discussão e estimular ainda mais a curiosidade de quem deseja aprender.

Série "Quebrando tabus": gerar relevo sombreado usando o QGIS | #SLGeoTbFaz

2011-12-26T07:30:00.008-03:00

Olá pessoal,

Como foram de Natal? Espero que tenha acontecido tudo dentro do normal! Final de ano chegando, festividades com família e amigos, mas ainda dá tempo de aprender/conhecer mais uma ferramenta do QGIS, dando continuidade a Série "Quebrando tabus: software livre também faz!".

Para esta Série, adotamos a hashgat #SLGeoTbFaz e publiquei, aqui no blog, 3 tutorias da mesma, confira:

Vale ressaltar que esta Série conta com o apoio dos amigos "geoblogueiros" Eliazer Kosciuk (@eliazerk), Sadeck (@sadeckgeo), Anderson Medeiros (@ClickGeo), Esdras Andrade (@geoparalinux) e José Carlos (@geojcarlos), com publicações em seus respectivos blogs. A Série também recebeu uma nota na 3ª Edição da Revista FOSSGIS Brasil (@fossgis)

3ª Edição da Revista FOSSGIS Brasil

Recentemente, publiquei o tutorial mostrando como gerar relevo sombreado no ArcGIS. Esta função obtém a iluminação hipotética de uma superfície, de elevação, para determinação dos valores de iluminação para cada célula em um raster.

Veja como gerar relevo sombreado usando o ArcGIS

Para este tutorial usaremos uma SRTM da mesma região usada para o ArcGIS. Clique aqui para fazer o download e acompanhar o passo a passo. Adicione o raster ao QGIS. Siga o caminho Raster > MDE (modelo digital de elevação), como mostra a figura 1.

Figura 1. Caminho para gerar relevo sombreado.

A janela correspondente ao passo acima irá se abrir (Figura 2). Em Arquivo de entrada (MDE raster), insira a SRTM. Em Arquivo de saída dê um nome e um destino ao novo arquivo que será gerado. No drop list Modo selecione a opção "sombreamento". Mantendo as configurações padrões já apresentadas nesta janela o seu relevo sombreado é gerado, mas vamos ver algumas delas:

Em Azimute de insolação insira, em graus, a posição angular do sol (o padrão utilizado é 315°, medido em sentido horário a partir do norte, ou seja, 45° a noroeste);

Em Altitude de insolação insira a inclinação do sol, em graus, em relação ao horizonte, de 0° a 90° (o padrão utilizado é 45°);

Por fim, clique em carregar na tela quando concluído para que ao final do processo o arquivo seja adicionado automaticamente à lista de camadas e clique em OK.

Figura 2. Janela MDE (modelo digital de elevação)

Veja o resultado na figura 3. E então? Alguma diferença para um software proprietário, como o ArcGIS?

>> Software livre também faz! <<

Como é de praxe, vou disponibilizar este tutorial para download em formato PDF, mas agora será através do PagSocial. Clique no banner abaixo, faça o download e automaticamente haverá divulgação em sua rede social: twitter, facebook ou orkut. Experimente e valorize nosso trabalho:

Esta é a minha última postagem de 2011. Obrigado pela ótima audiência neste ano: uma média de 50 países acessam, todos os meses, este blog!

Aproveite e curta nossa página no facebook

Desejo a todos um Feliz Ano Novo, e um 2012 repleto de felicidades e realizações!!!

Material de apoio:
GISTutor

Seja parceiro | Autoria do Profº Christian Nunes

2011-12-19T07:30:00.073-03:00

Olá pessoal,

Essa é mais uma postagem do nosso programa Seja Parceiro. Através deste programa o leitor pode contribuir divulgando o seu material. Para tal, estabelecemos algumas regras para este tipo de publicação, confira:

Veja como enviar seu material para divulgação, leia as regras do Seja Parceiro

Pois bem, recentemente recebi um e-mail do Profº (e amigo) Christian Nunes, da Universidade Federal do Pará - UFPa. O Professor já é conhecido dos leitores, pois já foi divulgado aqui no blog, também pelo Seja Parceiro. Recomendo que leiam:

Seja Parceiro | O ensino da cartografia na era da (geo)informação, por Christian Nunes

Neste e-mail, o professor compartilhou cinco textos, que estão disponíveis abaixo:

Currículo do Profº Christian Nunes:

Doutorando do Programa de Pós-Graduação em Ecologia. Mestre em Geografia. Especialista em Gestão Ambiental. Bacharel e Licenciado em Geografia.
Pesquisador do Grupo Acadêmico a Produção do Território e Meio Ambiente na Amazônia - GAPTA. Atualmente é Professor Assistente na Faculdade de Geografia e Cartografia da Universidade Federal do Pará e Coordenador do Laboratório de Análise da Informação Geográfica - LAIG.

Estamos realizando uma pesquisa sobre o perfil de nossos leitores. É importante sua participação. É rápido! São apenas 10 perguntas. Clique aqui e participe!

Até a próxima!

Luís Lopes

Siga @geoluislopes
Skype: geo.luislopes

6 maneiras de acompanhar o blog GeoLuisLopes [reeditado]

2011-11-17T07:30:00.015-03:00

Olá pessoal,

Hoje a postagem será diferente do normal, mas não deixa de ser um "tutorial".

Sabemos hoje que a internet é uma das principais fontes de conhecimento. Diversos blogs/sites da nossa área, de geotecnologias, tem surgido e bastante novidade. Isto é notável!

Deviso a este boom de informações, nós, blogueiros, temos que buscar alternativas para que nossos leitores fiquem a par de todo o conteúdo divulgado aqui, nesta página. Outro fator crucial para escrever esse post foi um e-mail que recebi pelo formulário de contato de Geiny Viegas, perguntando como poderia acompanhar o blog.

Então mostrarei 6 formas que você pode ficar por dentro do que é postado. Então vamos lá!

fonte: mestreseo

1. Curtir nossa página no Facebook

O blog GeoLuisLopes possui uma fan page no Facebook. Curtindo a página você acompanha todas as novidades em seu próprio perfil da rede social.

Curta nossa página no Facebook

2. Assinar o nosso Feed

Feed nada mais é que facilitar a leitura de um ou mais blogs/sites concentrado em um único ponto. A exemplo disto é para quem possui conta Gmail, pois é possível ler os feeds assinados, através do Google Reader. A imagem abaixo indica que o blog tem um feed. No blog ele encontra-se no canto superior direito. Caso tenha problemas, clique aqui.

http://feeds.feedburner.com/GeotecnologiasLusLopes

3. Rede sociais

Constantemente, divulgo as novas e antigas postagens por meio do Twitter e, recentemente, no Google+. Converso também com leitores/seguidores por meio do Skype.

Perfil do Google+
Twitter Siga @geoluislopes
Skype: geo.luislopes

4. Cadastrar para receber atualizações em seu e-mail

Cadastre-se agora, inserindo seu e-mail no campo disponível abaixo, e receberá, no e-mail cadastrado, um link de verificação para confirmar a subscrição. Este sem dúvida é o mais interessante, pela comodidade.

Esta opção, também estará ao final de cada postagem.

5. Salvando em seus favoritos

Uma forma relativamente fácil e rápida de seguir. Independente do navegador usado há opção de salvar o conteúdo do blog de modo favorito, pois assim com um simples clique o página será aberta

6. Google Friend Connect

"Através deste recurso, um usuário usando uma conta Google, Yahoo! ou OpenID, pode interagir com outros visitantes do site ou de sua própria rede social, compartilhando recursos, arquivos ou aplicações OpenSocial" (wikipedia).

Na própria página do blog (fim da página) há um gadget que você pode se cadastrar (figura abaixo) e e receber as atualizações diretamente do seu perfil do blogger, caso escolha seguir publicamente o blog, ou interagir com o Google Reader, modo privado para seguir. Para fazer parte da rede Google Friend Connect, clique aqui.

Então isso pessoal,

Espero vocês agora!

Até mais...

Luís Lopes

Siga @geoluislopes

skype: geo.luislopes

Material de apoio para esta postagem:
http://www.escritosideologicos.com/ e
pt.wikipedia.org

Série "Quebrando tabus": visualizar orientação da declividade usando o QGIS

2011-10-28T07:30:00.013-03:00

Olá pessoal,

Antes de iniciar, gostaria de externar minha felicidade com a repercussão que está Série de postagens causou. Claro, que isto não seria possível sem a adesão dos colegas Eliazer Kosciuk (@eliazerk), Sadeck (@sadeckgeo), Esdras Andrade (@geoparalinux), Anderson Medeiros (@ClickGeo) e José Carlos (@geojcarlos) em postar em seus respectivos blogs as mesmas ferramentas existentes em outros softwares livres. Meus sinceros agradecimentos! Lembrando que a Série ganhou destaque, também, na 3ª Edição da Revista FOSSGIS Brasil

3ª Edição da Revista FOSSGIS Brasil

Mostraremos, neste tutorial, como visualizar a Orientação da declividade usando o Quantum GIS (QGIS), por meio da ferramenta Aspecto. Como a proposta de Série é mostrar que o software livre realiza funções com a mesma qualidade que um software proprietário, este mesmo procedimento foi feito no ArcGIS e publicado aqui no blog mesmo.

Visualizar orientação da declividade usando o ArcGIS

Para este tutorial usaremos uma SRTM, download aqui, adquirido no site da EMBRAPA em Relevo, devidamente corrigida.

Faça parte desse movimento! Mostre que software livre tem a mesma eficiência de um software proprietário. Quebre tabus! Não se esqueça de inserir a hashtag #SLGeoTbFaz para divulgar.

Este tutorial está disponível em PDF, clique aqui para download.

Gostou do tutorial? Deixem seus comentários e contribuições.

Luís Lopes

Siga @geoluislopes
skype: geo.luislopes

Material de apoio: GISTutor.com, manual QGIS 1.6.0

Prêmio TopBlog 2011 | Vote GeoLuisLopes

2011-10-07T14:11:00.003-03:00

Olá pessoal,

Como muitos já devem ter percebido, no canto superior direito do blog há um selo Top100, do prêmio TopBlog. Selo conquistado em 2010 devido a uma votação popular, na qual o blog estava concorrendo na categoria Tecnologia e Internet. Nesta ocasião ficamos entre os 30 melhores blogs do Brasil nesta categoria: Top30!

Em 2011, estamos concorrendo mais uma vez no Prêmio TopBlog.

“(...) Top Blog Prêmio é um sistema destinado a reconhecer e premiar, mediante a votação popular e acadêmica os Blogs Brasileiros mais populares (...)” - topblog.com.br

Se você gosta das postagens deste blog, ajudem-nos votando. É fácil!

Clique aqui e vote ou basta clicar no selo TopBlog, na coluna lateral do blog, e posteriormente em Votar.

Até mais,

Luís Lopes

Twitter: @geoluislopes

Skype: geo.luislopes

Lançada a 3a edição da Revista FOSSGIS Brasil

2011-09-30T08:31:00.002-03:00

Caros leitores,

É com grande satisfação que anunciamos o lançamento da terceira edição da Revista FOSSGIS Brasil. Nesta edição falaremos da situação das das Geotecnologias livres no Brasil, onde serão abordados os motivos para a sua utilização, desafios da implantação em instituições públicas, além de dicas para uma migração saudável. Há ainda um breve relato de como está o uso das Geotecnologias livres em Portugal.

Ainda no âmbito do governo falaremos sobre o i3Geo na gestão pública, um software livre brasileiro criado no Ministério do Meio Ambiente por Edmar Moretti, que é o entrevistado desta edição, onde fala um pouco mais sobre esta fantástica ferramenta.

Você ainda poderá saber um pouco mais sobre o ILWIS, um Desktop GIS que já foi proprietário e hoje é opensource; além do GISVM, uma máquina virtual GIS.

Não deixe de ler também os artigos sobre os projetos GeoMajas e MapServer, a matéria sobre como utilizar o padrão WPS com o OpenLayers e a nova seção, com o tema “Quebrando Tabus: Software livre também faz”.

Para fazer o download desta edição, basta acessar o seguinte link:

http://fossgisbrasil.com.br/download/

Espero que gostem desta edição e nos enviem suas opiniões e sugestões que são de fundamental importância para o nosso trabalho.

Atenciosamente,

Equipe FOSSGIS Brasil

Montando banco de dados espacial com a GDAL - Parte 4 (Obtendo os dados para o banco de dados (Shapefile do PRODES))

2011-09-20T16:13:00.000-03:00

Obtendo os dados para o banco de dados (shapefile do PRODES)

Os polígonos de desmatamento do PRODES são disponibilizados na forma de shapefile (SHP). Pode-se baixar os SHP por cena ou na forma de mosaico (http://www.dpi.inpe.br/prodesdigital/dadosn/).

O PRODES utiliza uma grade própria de cenas para fazer ao levantamento do desmatamento, sendo essa a unidade de trabalho do PRODES. Nessa postagem vamos trabalhar com os SHPs de PRODES por cena, sendo o ano base de 2010.

Podemos verificar no endereço das cenas de 2010 (http://www.dpi.inpe.br/prodesdigital/dadosn/2010/), que podemos baixar arquivos de imagem e SHPs.

Para organizar os dados, vamos criar uma pasta para o PRODES 2010 e baixar os SHPs nela. A seguir os passos no terminal configurado para usar o Linux.

- Criar o diretório e torna-lo como diretório corrente

mkdir prodes2010

cd prodes2010

- Obter a página do PRODES 2010 (o nome do arquivo é index.html)

wget http://www.dpi.inpe.br/prodesdigital/dadosn/2010/

- O arquivo index.html possui o nome dos SHPs, e vamos criar uma listagem desses arquivos com o endereço completo.

A listagem vai conter em cada linha o endereço completo de cada SHP, como exemplo, http://www.dpi.inpe.br/prodesdigital/dadosn/2010/PDigital2010_00165_shp.zip

Vamos criar um script shell, 'list_zip.sh', para listar os endereços completos, e coloca-lo no diretório onde está o arquivo 'index.html'.

Conteúdo do arquivo 'list_zip.sh':

for file in `grep _shp.zip index.html | cut -d= -f4 | cut -d'"' -f2`;

do echo "http://www.dpi.inpe.br/prodesdigital/dadosn/2010/$file";

done;

A sequência de processamento do script é:

1) Executa a sequência de comandos conectados com pipe '|', definida entre os caracteres '`'

`grep _shp.zip index.html | cut -d= -f4 | cut -d'"' -f2`

Obs.: A saída de cada comando passa a ser a entrada de outro comando (uso do pipe)

1.1) Lista as linhas do arquivo index.html que contém o texto '_shp.zip'

grep _shp.zip index.html

1.2) Corta o 4o campo (field) de cada linha, onde os campos são delimitados pelo caracter '='

cut -d= -f4

1.3) Corta o 2o campo (field) de cada linha, onde os campos são delimitados pelo caracter '"'

cut -d'"' -f2

2) Cada linha da saída da sequência de comandos do ítem (1) passa a ser o valor da variável

'file' dentro do comando 'for'

for file in `Lista de valores do ítem (1)`;

do "Comando onde utiliza o valor de file";

done;

Obs.: O comando de interação 'for' inicia com o comando 'do' e termina o laço de interação com o comando 'done'

2.1) Cada valor de 'file' no interior do comando 'for' é utilizado para escrever na saída padrão (o janela do terminal) o caminho completo do SHP.

echo "http://www.dpi.inpe.br/prodesdigital/dadosn/2010/$file"

Obs.: O valor de file é obtido pelo uso do caracter '$'

Você pode testar cada comando do script para estuda-los.

Para criar a lista contendo o endereço completo de cada SHP, execute o script shell, list_zip.sh, pelo comando 'sh', e canalize para um arquivo, ou seja, troque a saída padrão para um arquivo com uso do '>'.

sh list_zip.sh > list_zip.sh.txt

- O arquivo 'list_zip.txt' tem todos os SHPs que precisamos, assim, podemos baixar todos os arquivos usando o 'wget' , passando como argumento o arquivo com as listagem (uso da opção '-i')

wget -i list_zip.txt

A utilização dos comandos/programas via console LINUX é um recurso poderoso para obtermos dados de sites, para maiores informações vejam os detalhes do comando 'wget'.

Até a próxima,

Luiz Motta

motta.luiz@gmail.com

Série "Quebrando tabus" | Criar pontos aleatórios dentro de um polígono no QGIS

2011-08-23T05:30:00.003-03:00

Olá pessoal,

Dando continuidade a Série “Quebrando Tabus: SL também faz!”, postarei mais uma funcionalidade do Quantum GIS (QGIS). O tutorial dessa vez é ‘Criar pontos aleatórios dentro de um polígono’.

Funcionalidade esta, que já foi apresentada aqui no blog, porém, naquela ocasião, para o ArcGIS (confira aqui). Este tutorial obteve bastante audiência no blog, e por isso resolvi mostrar que o QGIS também é eficiente para o mesmo objetivo.

Antes começar o tutorial, gostaria de agradecer a receptividade da Série, especialmente ao Anderson Medeiros (@ClickGeo), Eliazer Kosciuk (@eliazerk), Esdras Andrade (@geoparalinux) e ao Sadeck (@sadeckgeo) por terem aderido a ideia, divulgado a Série em seus blogs e replicando os tutoriais para outros softwares.

Então vamos lá!

Para a criação de pontos aleatórios dentro de polígono, no QGIS, utilizarei um shapefile com limites de 5 municípios da região do Piemonte da Diamantina - Bahia (siga o tutorial, download do shape aqui).

Este tutorial está disponível em PDF, faça o download aqui.

Então pessoal, alguma diferença para o ArcGIS? O que estão achando da Série “Quebrando tabus: SL também faz!”?

Você domina ou trabalha em algum software livre em geoprocessamento? O que acha de colaborar nesta ideia?

Abraços e até mais!

Luís Lopes

Twitter: @geoluislopes

Skype: geo.luislopes

Série "Quebrando tabus" | Ferramenta Dissolver no QGIS

2011-08-19T07:30:00.080-03:00

Olá pessoal,

É evidente e notável o crescimento do uso de software livre em geoprocessamento. Tanto em instituições privadas como em públicas.

Em organizações/instituições públicas, órgãos do Governo, gastam-se centenas de milhares de reais em softwares proprietários, muitas vezes por indicações de compra de seus próprios técnicos que muitas vezes desconhecem a potencialidade de software livre de SIG, visto que são condicionados a trabalharem desde cedo com softwares proprietários.

Para aprofundar mais os conhecimentos ou ter os primeiros contatos em software livre para as geotecnologias, indico a leitura da revista FOSSGIS Brasil, lançada recentemente e que está em sua 2ª Edição.

Seguindo a linha de raciocínio, farei uma série de postagens, denominada “Quebrando tabus: SL também faz!”, mostrando que ferramentas presentes em software proprietário (no caso o ArcGIS, que é mais utilizado) também estão disponíveis em software livre (neste caso Quantum GIS, ou QGIS) com a mesma eficiência!

Tempos atrás publiquei o tutorial da ferramenta Dissolve: agregação baseado em um atributo, no ArcGIS, e agora mostrarei mesma ferramenta, dessa no QGIS: Dissolver.

Essa ferramenta irá agregar feições com base em um mesmo atributo. Neste caso usaremos como exemplo o shapefile de limite municipal do estado da Bahia (faça o download do shape e acompanhe o tutorial) e iremos gerar um novo shapefile de Mesorregião deste estado. Obs.: para isto o atributo Mesorregião deve existir na tabela para cada município (Figura 1).

Figura 1. Tabela de atributos com campo Mesorregião

Para dar início ao processo, siga o caminho Vetor > Ferramentas de geoprocessamento > Dissolver (Figura 2)

Figura 2. Caminho da ferramenta Dissolver no QGIS

A janela “Dissolver” abrir-se-á (Figura 3). Em “Entrar com camada vetorial” selecione o layer/camada que será usado no processo Dissolver (caso queira dissolver limites selecionados, marque a opção “Usar apenas feições selecionadas”). Em “Dissolver a partir do campo”, selecione o atributo que será usado para agregação, neste caso a Mesorregião. Clique em “Buscar” e dê um nome e um destino ao novo arquivo que será gerado. Após isto, clique em OK.

Figura 3. Janela Dissolver do QGIS

Ocorrerá o processamento e o QGIS perguntará se desejas adicionar automaticamente este novo arquivo gerado. Clique em “Yes” (Figura 4).

Figura 4. Adicionar o novo shape a lista de camadas

Observe o resultado final na figura 5, abaixo (Antes: municípios do estado da Bahia; Depois: agregação por mesorregião):

Figura 5. Resultado da ferramenta Dissolver no QGIS

E então? Gostaram da Série "Quebrando Tabus: SL também faz!"?

Deixem seus comentários!

Estamos concorrendo ao Prêmio TopBlog 2011: vote aqui e ajude-nos!

Abraços e até a próxima.

Luís Lopes

Twitter: @geoluislopes

Skype: geo.luislopes

GeoConnect People | Rede social do setor de geotecnologia

2011-08-15T19:42:00.002-03:00

Olá pessoal,

Essa postagem não se trata de mais um tutorial de algum software de SIG, mas para noticiar o lançamento (hoje, 15.08) da mais nova rede social voltado ao setor das geotecnologia: GeoConnect People

"Os conteúdos postados no GeoConnectPeople são ligados às tecnologias de geoprocessamento, Sistemas de Informação Geográfica (SIG), cartografia, agrimensura, fotogrametria, observação da Terra, geodésia, Sistemas Globais de Navegação por Satélites (GNSS), Serviços Baseados em Localização (LBS), Infraestruturas de Dados Espaciais (IDE) e suas diversas relações e aplicações." (MundoGEO)

Destaco que é de grande importância este novo canal. Aconselho a todos participarem, pois terá oportunidade de discutir, trocar informações e ideias com os grandes nomes do cenário das geotecnologias no Brasil.

Para fazer parte da rede social GeoConnect People, acesse aqui: http://geoconnectpeople.org/

Abraços e até a próxima!

Luís Lopes

Twitter: @geoluislopes

Skype: geo.luislopes

Seja Parceiro | Aplicando as Ferramentas Union e Field Calculator, por Renato Bossle

2011-08-05T08:00:00.008-03:00

Olá pessoal,

Como alguns já perceberam, este blog tem um espaço para as pessoas que queiram divulgar seus tutoriais, artigos, etc.

Poucos dias atrás recebi um e-mail de Renato Bossle, biólogo, da Prefeitura Municipal de São José dos Pinhais - PR, e mestre em Ecologia com ênfase na área de gestão de bacias hidrográficas e meio ambiente. O tutorial do Renato Bossle mostra a ferramenta Union, do ArcToolbox, e uma operação de soma usando o Field Calculator (software ArcGIS).

O tutorial elaborado por Renato Bossle está disponível neste link

Se você também tem/produziu um tutorial, artigo ou qualquer material relacionado a geotecnologia, veja como contribuir mandando-o para ser divulgado aqui no blog: leia as regras aqui.

Até mais,

Luís Lopes

geo.luislopes@gmail.com

twitter: @geoluislopes

skype: geo.luislopes

Montando banco de dados espacial com a GDAL - Parte 3 (Emulação do ambiente UNIX/SHELL no WINDOWS)

2011-07-18T08:00:00.005-03:00

>> Emulação do ambiente UNIX/SHELL no WINDOWS <<

O UNIX/SHELL, mas conhecido como SHELL, é um conjunto de programas que acessam o núcleo do sistema operacional UNIX, permitindo obter informações do sistema operacional. Num ambiente UNIX, como exemplo o LINUX, é comum o uso de um terminal, também chamado de console, que é uma janela onde chamamos os programas do SHELL. No caso do Windows acessamos os comandos do MS-DOS (Sistema Operacional da Microsoft) pelo programa "Prompt de Comando".

Nessa parte da postagem, vou mostrar os programas que foram utilizados para a montagem do banco de dados espaciais. Inicialmente vamos ver como podemos emular o SHELL no Windows com a obtenção dos programas pelos pacotes da OSGEO4W e pelo projeto GnuWin32.

1) Pacote da OSGEO4W

Na Parte 2 da postagem foi mostrado a opção de utilizar o instalador da OSGEO para Windows (OSGEO4W) para obter o pacote (conjunto de arquivos) da GDAL.

O OSGEO4W possui um pacote chamado MSYS que contém um conjunto mínimo de programas do SHELL para Windows, logo, deve-se marcar esse opção no instalador da OSGEO4W para termos alguns comandos que iremos utilizar nessa postagem (Figura 1).

Figura 1. Marque a opçção MSYS no instalador da OSGEO4W

O MSYS não possui todos os programas do SHELL que são utilizadas na postagem, assim, precisaremos instalar os pacotes que complementaram os comandos utilizados nessa postagem.

2) Pacote do projeto GnuWin32

Os pacotes do projeto GnuWin32 (http://gnuwin32.sourceforge.net/packages.html) são disponibilizados como arquivos ZIP.

Cada arquivo ZIP, corresponde a um programa do SHELL e deve ser descompactado de forma padronizada no diretório a ser instalado os pacotes do GnuWin32.

Os pacotes a serem utilizados são:

2.1) coreutils: comandos básicos do SHELL.

2.2) unzip: decompacta arquivos ZIP e

2.3) wget: baixa arquivos da Internet.

A seguir, figura 2, é ilustrado a árvore de diretórios e um exemplo do conteúdo zipado do programa "wget.exe".

Figura 2. Árvore de diretório (clique na figura para ampliar)

Na figura acima, o pacote de comandos do GnuWin32 foi colocado no diretório "C:\_gnu_windows", e nesse diretório é descompactado os arquivos de cada pacote.

Observa-se que no pacote temos os nomes dos diretórios que devem estar na árvore de diretórios da instalação do GnuWin32, logo, percebemos que não temos um SETUP, ou seja, um programa de instalação de pacotes, e sim, simplesmente os arquivos ZIPADOS que tem os programas que queremos ter no nosso ambiente de processamento e seus respectivos arquivos em diretórios padronizados (bin,lib,man, manifest, share, ...).

Para termos um console que emule o SHELL temos que configurar o ambiente de trabalho, assim, é necessário criar um arquivo BAT contendo as definições necessárias.

3) Configuração do ambiente de processamento no Windows

A configuração para processamento consistirá num arquivo BAT que terá todas as definições necessárias para executar os programas do SHELL e dos utilitários da GDAL, ou seja, será nele que iremos rodar scripts e "saborear" o poder do SHELL. Uma vez configurado o arquivo BAT, esse arquivo passará a ser o programa que iremos executar para emular o ambiente SHELL.

O conteúdo do BAT (clique aqui para baixar o BAT - OSGeo4W_PRODES.bat) consiste em:

@echo off

set OSGEO4W_ROOT=C:\OSGeo4W

PATH=%OSGEO4W_ROOT%\bin;%PATH%

for %%f in (%OSGEO4W_ROOT%\etc\ini\*.bat) do call %%f

REM *** // Gnu Windos and Scripts \\ ***

set http_proxy=172.21.12.1:3128

PATH=%path%;C:\_gnu_windows\package\bin;C:\_bd_prodes\sh

REM *** \\ Gnu Windos and Scripts // ***

@echo on

@cmd.exe

A seguir as observações referente as definições contidas no BAT:

3.1) O diretório onde está instalado os pacotes da OSGEO4W é "C:\OSGeo4W", pasta escolhida no momento da instalação da GDAL (veja na Parte 1 dessa postagem).

3.2) Foi definido a variável de ambiente http_proxy, com o valor do IP da maquina de PROXY na rede (valor = 172.21.12.1) e a porta do PROXY (valor = 3128). Essa variável só deve ser definida se o seu computador estiver numa rede corporativa, onde geralmente tem PROXY (veja com o administrador da rede os valores do seu PROXY), caso contrário, elimine ou faça um comentário, escrevendo REM no início da linha.

3.3) O pacote GnuWin32 foi instalado no diretório "C:\_gnu_windows".

3.4) Os Scripts criados para fazer o banco de dados espacial (clique aqui para baixar os Scripts - PRODES_sh.zip) estão no diretório "C:\_bd_prodes\sh", e podem ser chamado no terminal por estarem na variável PATH.

O BAT para configurar o console de processamento pode ser colocado em qualquer local, como sugestão, coloque no diretório "C:\_bd_prodes\bat".

Para facilitar o uso do terminal, você pode criar um atalho no seu Desktop apontado para esse arquivo (Figura 3).

Figura 3. Criar atalho em seu Desktop (clique na figura para ampliar)

Depois desse esforço todo está na hora de saborearmos o SHELL, assim, tome um cafezinho e volte animado para proximo ítem.

4) Saboreando o SHELL

Antes de colocar a mão na massa, gostaria de avisar que usei termo SHELL de forma simplista, para maior clareza deve-se consultar as referências sobre UNIX/LINUX. Nessa postagem preferi não aprofundar no UNIX/LINUX, e sim, mostrar uma estação de processamento GIS onde são manipulados dados via terminal.

Os programas do SHELL, também conhecidos como comandos, possui alguns conceitos, e para facilitar vamos explicar sucintamente o que é argumento de um comando.

Os comandos do SHELL podem conter argumentos, que são os parâmetros passados em sua chamada, fazendo uma analogia a uma função listar (x), x é o argumento da função, assim, a chamada para o comando 'listar' seria escrito no terminal como, "listar x".

Existem argumentos obrigatórios e opcionais, onde os argumentos opcionais são identificados com um prefixo de hífem ('-').

Para maior entendimento sobre o uso de um comando utilizamos o argumento de 'help', que mostra as informações sobre o comando, exemplo, "cp --help". Observa-se que o argumento 'help' possui no prefixo dois hífens ('--')

Agora vamos mostrar os scripts mais utilizados nessa postagem.

4.1) Baixando dados na Internet - wget

O 'wget' permite baixar conteúdos da Internet nos protocolos HTTP, HTTPS e FTP, e é por meio deles que vamos obter os dados do PRODES. O uso 'wget' em redes com PROXY necessita da definição da variável 'http_proxy' (ver item 3.2).

Para maiores informações sobre o 'wget' consulte o manual do comando utilizando o argumento padrão '--help', exemplo, 'wget --help'.

Vamos baixar os dados de PRODES do ano de 2010, para isso execute no seu terminal de processamento (OSGeo4W_PRODES.bat).

A seguir a sequência de comandos para 1) ir até o diretório de trabalho, 2) criar o diretório zip, 3) ir para diretório zip dentro do diretório atual e 4)fazer o download da página dos dados:

cd /home/lmotta/_work_2011/blog_lopes/banco_dados_prodes

mkdir zip

cd zip

wget http://www.dpi.inpe.br/prodesdigital/dadosn/2010

O 'wget' baixa o arquivo com o conteúdo da página. No Windows a implementação do 'wget' cria o arquivo index.html, no Linux/Ubuntu foi criado o arquivo '2010', logo, o comportamento do 'wget' pode variar conforme o sistema operacional (Macintosh, Distribuições LINUX,...).

Podemos dar uma rápida olhada no arquivo gerado utilizando o comando 'head', que permite ver as primeiras linhas de um arquivo.

head -n 20 2010

O comando acima exibe as 20 primeiras linhas do arquivo 2010. Observa-se o uso do argumento opcional '-n 20', para definir a quantidade de linhas a serem mostradas.

Podemos verificar que temos os endereços dos shapefiles que queremos baixar, como exemplo, HREF="PDigital2010_00265_shp.zip".

O endereço do arquivo 'PDigital2010_00265_shp.zip' ficará:

http://www.dpi.inpe.br/prodesdigital/dadosn/2010/PDigital2010_00265_shp.zip

O 'wget' pode baixar vários arquivos em um só comando, para isso devemos ter uma listagem dos endereços de cada shapefile. Exemplo do conteúdo do arquivo(exemplo list_shapefile.txt):

http://www.dpi.inpe.br/prodesdigital/dadosn/2010/PDigital2010_00165_shp.zip

http://www.dpi.inpe.br/prodesdigital/dadosn/2010/PDigital2010_00166_shp.zip

http://www.dpi.inpe.br/prodesdigital/dadosn/2010/PDigital2010_00167_shp.zip

...

Vamos executar um comando (composto de uma sequência de comandos) para termos os arquivos a serem utilizados nessa postagem.

4.1.1) Comando para ter a listagem dos endereços dos shapefiles

grep _shp.zip 2010 | cut -d= -f4 | cut -d'"' -f2 | sed 's/PDi/http:\/\wwww.dpi.inpe.br\/prodesdigital\/dadosn\/2010\/PDi/g' > list_shapefile.txt

Observações:

- Posteriormente será explicado os comandos utilizados (grep, cut e sed).

- Foi utilizado o recurso de 'pipe' (ver detalhes no próximo item).

- Foi utilizado o recurso de canalização, caracter '>', onde a saída passa a ser um arquivo (no caso foi o arquivo list_shapefile.txt), em vez da saída padrão (o terminal).

- 2010 é o nome do arquivo baixado pelo wget ao passar o endereço dos dados de 2010 do PRODES.

4.1.2) Comando para obter os shapefiles:

wget -i list_shapefile.txt

Obs.: o argumento '-i' define o uso do arquivo de leitura com os endereços.

Depois de baixar os arquivos zipados, crie um diretório, e descompacte os mesmos no diretório criado.

4.1.3) Comando para descompactar os shapefiles

cd ..

mkdir shp

cd shp

ls -1 ../zip/*.zip | sed 's/\.\./unzip -j \.\./g' > unzip_shp.sh

sh unzip_shp.sh

Observação: Primeiro criamos o diretório 'shp' dentro do diretorio 'banco_dados_prodes' e descompactamos todos os ZIP's.

A seguir são mostrados os comandos de forma mais didática utilizados nessa postagem.

4.2) Listar o conteúdo do diretório - 'ls'

O 'ls' lista as informações dos arquivos, para saber como formatar a saída e os tipos de informações, veja as opções do comando digitando 'ls --help'.

Geralmente, queremos aproveitar a lista (saída do comando ls) para extrair informações, para isso, existe uma forma de aproveitar a saída do comando como entrada para outro comando. Isso se chame 'pipe' , logo, vou mostrar a seguir um exemplo do comando 'ls' com o comando de contador de palavras ('wc').

4.3) Contar palavras -'wc'

Naturalmente, devemos dar uma olhada sobre o comando, 'wc --help', aqui vou direto a prática! Como contar os arquivos contidos no diretório, mas para ser mais realista contar o numero de shapefiles no diretório corrente do terminal:

ls -1 *.shp | wc -l

Para maior compreensão, execute cada comando, experimente, 'ls ', 'ls -l', 'ls -1', 'ls -1 *.shp', a listagem desses exemplos já mostram a finalidade de cada argumento.

Vamos fazer o papel do SHELL, ou seja, interpretar os comandos. O 'ls -1 *.shp' lista os arquivos shapefiles (por isso o uso do filtro *.shp) linha a linha e apenas o seu nome, comportamento da listagem definido pelo argumento '-1'.

Com a saída da listagem podemos passar como entrada para o comando 'wc -l', onde esse comando vai contar as linhas contidas na entrada, assim, temos o total de shapefiles num diretório.

Isso está ficando interessante! No Windows eu teria que selecionar apenas os arquivos com extensão 'shp', e olhar lá em baixo o total de arquivos selecionados, fica evidente que o uso do 'ls -1 *.shp | wc -l' é bem mais confiável do que fazer a operação na janela de gerenciador de arquivos do Windows.

O comando 'wc' pode contar outras informações, para saber mais, vc. já deve saber como! Help nele!

4.4) Concatenar e exibir arquivos - cat

O comando 'cat' é usado para mostrar o conteúdo de um ou mais arquivos, se vc. quer mostrar o conteúdo de mais de um arquivo, vc. está "concatenando" os arquivos, ou seja, juntando todos e mostrando de uma só vez.

O mais comum é ver o conteúdo de um arquivo, como exemplo, podemos ver o sistema de referência (SR) de um shapefile exibindo o valor do arquivo PRJ.

cat PDigital2010_00165_pol.prj

4.5) Pesquisando textos - grep

O grep é um poderoso comando, onde podemos fazer uma série buscas de textos em arquivos, ele merece ser pesquisado!

Como exercício vamos verificar se todos os shapefiles possuem o SR geográfico. No item anterior (exemplo de 'cat') você provavelmente viu o a palavra 'GEOGCS' no início da listagem, logo, o shapefile está no SR geográfico.

Para sabermos se todos os shapefiles estão no SR geográfico, vc. primeiro conta quantos arquivos PRJ existem e depois verifique se a contagem é a mesma com o comando:

grep -i "GEOGCS" *.prj | wc -l

No comando 'grep -i -c "GEOGCS" *.prj' fizemos a busca da palavra GEOGCS em todos os arquivos PRJ, e especificamos com o argumento 'i' para ignorar se o termo de busca está maiúscula ou minúscula (GEOGCS ou geogcs), a saída do comando mostra em cada linha o nome do arquivo e o texto com a palavra GEOGCS, como cada linha só aparece se existir a palavra GEOGCS, basta contarmos o número de linha ' | wc -l'.

O resultado é bom, mas não me agradou!

Tenho que usar dois comandos, um para contar o número de arquivos PRJ e o outro para contar se os arquivos PRJ possuem GEOGCS.

E se os resultados forem diferentes?

Vamos mostrar mais comandos e depois mostraremos como melhorar a pesquisa para obter os SR's.

4.6) Cortar texto - cut

O 'cut' permite extrair textos, delimitados pelo número da coluna ou delimitadores, em cada linha de um arquivo de texto.

Aproveitamos e vamos explicar o comando do item 4.1.1.

grep _shp.zip 2010 | cut -d= -f4 | cut -d'"' -f2

O grep filtra do arquivo '2010' (lembre-se no Windows é index.html) as linhas que possuem o texto '_shp.zip'. Essa saída passa para o comando 'cut -d= -f4'. O comando 'cut' usa como delimitador o sinal de igual '=-d' e extrai a informação que esta no 4º campo, ou seja, ele interpreta as linhas como se fossem uma tabela limitada por sinal de '=' e obtém os valores da 4ª coluna.

Como exercício, vamos usar os arquivos shapefiles (o usuário deve estar no diretório 'shp'), e obter uma lista com apenas as datas dos shapefiles (*.shp)

ls -l *.shp | cut -b45-

O 'cut' vai extrair da listagem (comando 'ls -l *.shp') os textos que começam a partir da coluna 45 (opção '-b45-').

4.7) Interpretador de texto - sed

O 'sed' oferece uma série de recursos para trabalhar com textos, o mais comum é a substituição de palavras no texto.

Como exemplo, vamos criar um arquivo texto que contém os comandos para mover arquivos dos shapefiles.

ls -1 *.shp | sed 's/\.shp/\.\* new_dir/g' | sed 's/PDi/mv PDi/g' > mv_shp.sh

O primeiro comando lista apenas os nomes dos arquivos com extensão 'shp', essa saída passa ser entrada para o próximo comando, sed 's/\.shp/\.\* new_dir/g'.

Vamos explicar o primeiro uso do 'sed', o argumento 's/\.shp/\.\* new_dir/g' é interpretado da seguinte forma:

o 's/' define o início da palavra a ser substituída, no caso '\.shp';
o caracter de escape '\' é para evitar a interpretação do '.', ou seja, trocar a palavra '.shp';
o próximo caracter '/' inicia a palavra que vai ser inserida '\.\* new_dir' ('.* new_dir); e
o caracter /g finaliza a substituição, resultando na substituição de '.shp' por '.* new_dir'.

O segundo uso do 'sed' substitui 'PDi' por 'mv PDi' e a canalização '>' cria o arquivo mv_shp.sh.

É padrão usarmos a extensão '.sh' para representar script do SHELL. Para mover os arquivos de cada shapefile executamos o script via interpretador do SHELL (comnando 'sh').

sh mv_shp.sh

4.8) Ordenar - sort

O 'sort' ordena as linhas de um texto. Seu uso dispensa exemplo, basta passar o nome de um arquivo 'sort nome_arquivo' ou receber a saída de um comando 'ls -1 | sort'

4.9) Listar apenas valores únicos - uniq

O 'uniq' mostra apenas os valores únicos das linhas de um texto.Um exemplo é verificar as datas que foram gerados os shapefiles de um diretório.

ls -l *.shp | cut -b37-46 | uniq

O primeiro comando lista os arquivos shapefiles com as suas informações, e o próximo comando retira apenas as datas de cada arquivo (colunas 37 à 46), e o 'uniq' mostra os únicos valores. Como resultado do comando temos as datas dos arquivos.

4.10) Descompactar - unzip

O unzip, trabalha com arquivos ZIP, permitindo, testa-los, mostrar os conteúdos e descompactar. Seu uso é bastante intuitivo, como exemplo, para descompactar o comando é 'unzip nome_arquivo'.

4.11) Finalmente vamos fazer a pesquisa sobre os tipos de sistema de referência.

O arquivo PRJ contém uma linha com o SR do shapefile no formato WKT (http://www.geoapi.org/3.0/javadoc/org/opengis/referencing/doc-files/WKT.html).

Se usarmos o comando 'cat' para concatenar todos os PRJ, podemos ter um problema, nem sempre o arquivo PRJ possui o carácter 'terminador de linha', assim, ao concatenar todos os PRJ, estaremos juntando todos os SR dos shapefiles numa mesma linha.

A construção de script SHELL pode ser feita numa única linha, e é o que estamos fazendo até agora, dessa forma, criamos o SHELL e executamos direto na linha do terminal.

O SHELL pode utilizar os recursos de programação, como definição de variáveis, interação (loops), condicional, funções e outros.

Nesse item vamos criar um script SHELL como um fonte de um programa e rodá-lo posteriormente. Isso já fizemos, criamos uma listagem de comandos para depois executar via interpretador SHELL (comando 'sh').

A ideia aqui é mostrar que existe esse recurso no seu PC, e para tirar proveito, o usuário avalia se o investimento em entendimento do SHELL pode ser uma boa opção para "domar seus dados espaciais".

Os arquivos PRJ não possui terminador de linha, assim, temos que colocar o terminador de linha para cada SR (conteúdo do PRJ) listado, e depois aplicar o comando 'uniq', assim, teremos a listagem dos SR existente nos shapefile, sendo mais interessante do que ficar comparando totais (exemplo do item 4.5).

Para termos a listagem dos tipos de SR's, vamos criar um arquivo de script SHELL denominado 'uniq_srs.sh'. O arquivo tem que estar no formato texto padrão, para isso use o 'Bloco de Notas' ou outro de sua preferência. A execução do script é feita com o comando 'sh uniq_srs.sh'.

A seguir a listagem do Script:

#!/bin/sh

# Return: List of Reference Systems types used by a set shapefiles

# Remove tempfile if exist

if [ -e "tempfile" ]

then

rm "tempfile"

# Add all reference systems inside tempfile

for file in `ls -1 *.prj`

echo "`cat $file`" >> tempfile

done

uniq tempfile

rm tempfile

Descrevendo o script:

A linha que inicia com o caracter '# ' não é interpretada, assim, pode ser usada para comentar o script.
A 1a linha define para o UNIX que se trata de um script SHELL, assim, se esse arquivo tem permissão de execução, pode ser chamado diretamente. Coisas do mundo UNIX/LINUX!
O script verifica se existe o arquivo 'tempfile', se existir, apaga o mesmo.
É criado uma lista contendo o nome dos arquivos PRJ, `ls -1 *.prj`. Observa-se o carácter de crase (`) no início e no final do comando ls -1 *.prj, força o comando a ser interpretado 1o, para depois fazer o loop.
O loop é feito com a instrução, for file in [Lista de arquivos], onde para cada arquivo PRJ é associado a variável file, e essa é utilizada dentro do loop. O loop está envolvido pelas palavras 'do' e 'done'.
Para cada valor de 'file' é executado o comando, echo "`cat $file`" >> tempfile, onde, o comando echo exibe o conteúdo de um texto (definido entre o caracter ") e coloca o final de linha. A saída é adicionada ao arquivo tempfile pelo dupla canalização '>>', logo, cada saída do comando echo é posta como uma nova linha no aquivo. O `cat $file`, é executado primeiro, ou seja, mostra o SR do arquivo, sendo colocado no texto (definido entre o caracter "), o sinal de $ é para ler o valor da variável.
Depois do loop, é utilizado o comando 'uniq' com o arquivo tempfile, mostrando apenas os tipo de SR's, depois removemos o tempfile.

Espero que tenha tido uma visão do uso de scripts para investigar grande quantidade de arquivos. Seu uso, passa a ser uma necessidade em ambientes servidores, mas acredito que vale a pena conhecer um pouco desse recurso tão escondido pelas "janelas" do seu Windows.

Até a próxima etapa e aguardo as dúvidas e comentários!

Luiz Motta

motta.luiz@gmail.com

Seja Parceiro | O ensino da cartografia na era da (geo)informação

2011-07-15T10:48:00.004-03:00

Olá pessoal,

Como é sabido, este blog, Geotecnologias Luís Lopes (ou GeoLuísLopes), abre espaço para publicação de materiais (artigos, textos, etc...) de qualquer pessoa que queira dar sua contribuição no ramo das geotecnologias, através do Seja Parceiro (clique aqui e contribua também).

O artigo de hoje foi enviado pelo amigo Prof. Christian Nunes (FGC/UFPA). O texto intitulado "O ensino da cartografia na era da (geo)informação" destaca a importância dos professores da área se aperfeiçoarem com as novas geotecnologias e usufruir das ferramentas disponíveis. Destaca também a interatividade professor/aluno em sala de aula de modo que a facilitar o entendimento da cartografia.

Segundo Christian Nunes:

"Nos dias de hoje, a popularização da informática e o acesso a textos e a produtos cartográficos em formato digital, além da flexibilidade/facilidade no manuseio das chamadas geotecnologias (hardwares, softwares e técnicas direcionadas para a geração de geoinformação, ou seja da informação espacial/geográfica), facilitou as atividades produtivas de diversos indivíduos, não somente para os professores de geografia ou de outras ciências humanas, mas de profissionais de várias ciências que integram às suas atividades docentes o ensino do/e pelo mapa, além de outros sujeitos que utilizam essas geotecnologias em sua localização cotidiana, como o Sistema de Posicionamento Global – GPS, ou a visualização das imagens do Google Earth, ou em outros sites de empresas privadas e governamentais."

Prof Christian destaca ainda o uso do software educacional Marble (conheça este software) que é um globo virtual e Atlas Mundial (Figura 1).

Figura 1. Software Educacional Marble

Para ler o artigo completo do Prof Christian Nunes, clique aqui e boa leitura!

Quer divulgar seu texto? Veja como você pode contribuir: Seja Parceiro!

Até mais!

Luís Lopes

Twitter: @geoluislopes

Skype: geo.luislopes

Definir e/ou mudar sistema de projeção e datum em lote | ArcGIS

2011-07-08T09:28:00.000-03:00

Olá pessoal ,

Depois de uma rápida votação pelo twitter (@geoluislopes) para escolher qual seria o próximo tutorial, o vencedor foi o ArcGIS.

Pois bem, hoje mostrarei como definir e/ou reprojetar seus dados (vetorial ou raster). Normalmente usamos as ferramentas Define Projection, para definir projeção, Project, para reprojetar um vetor, ou Project Raster, para reprojetar um raster.

Mas a boa destas ferramentas que vou apresentar é que você poderá processar os dados em lote. Isso mesmo! Em LOTE.

Quem nunca se deparou com situações em que diversos (dezenas ou centenas) dados precisam ser definidas as projeções ou reprojetados? As ferramentas que realizam estes procedimentos são Batch Define Coordinate System, para definir uma projeção, e Batch Project, para reprojetar.

Primeiro vamos conhecer a ferramenta para definir projeção em lote (Batch Define Coordinate System).

Para isto, ative seu ArcToolbox e siga o caminho Samples > Data Management > Projections > Batch Define Coordinate System (Figura 1).

Figura 1. Ferramenta para definir projeção em lote

A interface desta ferramenta será ativada (Figura 2). Em “Input Datasets” insira deus dados¹. Em “Coordinate System” defina a projeção desejada².

¹ Este campo não trabalhará com camadas como entrada, por isso realize a busca no diretório

² Esta definição deverá ser a mesma em que o dado foi construído.

Figura 2. Interface para definir projeção em lote

Para utilizar a ferramenta de reprojeção de dados em lote (Batch Project), ainda no ArcToolbox, siga em Data Management Tools > Projections and Transformations > Feature > Batch Project (Figura 3)

Figura 3. Ferramenta para reprojetar em lote

A interface desta ferramenta será ativada (Figura 4). Em "Input Features Class or Dataset” insira os dados que utilizados no processo. Em “Output Workspace” escolha o diretório de saída, ou seja, onde os arquivos serão armazenados. E, em “Output Coordinate System” escolha a nova projeção que os dados receberão.

Figura 4. Interface para reprojetar em lote

Então pessoal, é ou não ótimas ferramentas? Agiliza e muito nosso trabalho!

Aguardo os comentários!!!

Até mais,

Luís Lopes

twitter: @geoluislopes

e-mail: geo.luislopes@gmail.com

Revista FOSSGIS Brasil – Edição N° 2!

2011-06-30T10:08:00.001-03:00

Chegou o grande dia do lançamento da Revista FOSSGIS Brasil – Edição N° 2!

Após a grande repercussão da primeira edição da Revista, lançada em março com o temaRedescobrindo as Geotecnologias, iremos tratar nesta edição sobre um assunto muito importante, e que tem tudo a ver com as Geotecnologias livres: a Educação!

Capa da 2ª Edição da Revista FOSSGIS Brasil

É sabido que a educação no Brasil enfrenta grandes dificuldades, porém muitas instituições de ensino têm buscado alternativas para poder proporcionar um ensino de melhor qualidade, e com isso auxiliar na formação dos cidadãos.

Desta forma, diversas instituições têm visto no software livre uma opção, uma vez que este tipo de programa, além de não ter custo com licenciamento, auxilia os alunos no entendimento da palavra cooperação.

Seguindo esta vertente, iremos tratar nesta edição do uso dos sistemas de informações geográficas (SIG) no ensino, mostrando que é possível sim, proporcionar um ensino de qualidade utilizando o softwares livres, assim como é feito no curso superior de Tecnologia em Geoprocessamento do IFPB.

Também tivemos o prazer de realizar uma entrevista com Jorge Xavier da Silva, professor emérito da UFRJ e mentor do primeiro SIG nacional, o SAGA-UFRJ.

Ainda nesta edição você poderá conferir os artigos sobre o WKT Raster, uma alternativa doPostGIS para manipulação de informações Raster, Quantum GIS, Kosmo GIS, a continuação da coluna Mobile GIS, entre outros.

Para fazer o download acesse o site oficial do projeto (http://fossgisbrasil.com.br/) ou clique aqui ou ainda na imagem acima!

Resultado da promoção de Junho: GeoLuisLopes - Parêntese Editora

2011-06-17T21:56:00.003-03:00

Olá pessoal,

Já temos os sortudos que levaram os livros Calibração Absoluta de Sensores Orbitais e Novas Tecnologias para o Monitoramento

Os ganhadores foram:
@uedbn (Calibração Absoluta de Sensores Orbitais) e
@geodiversidade (Novas Tecnologias para o Monitoramento)

Aqui o link do resultado: http://sorteie.me/1ORZ3z

obs.: Pelo regulamento, somente residentes no Brasil poderiam participar da promoção. Diante disto, o livro "Calibração Absoluta de Sensores Orbitais", que teve como sorteado @uedbn, residente em Portugal, ficará disponível para outro sorteio!

Parabéns @geodiversidade! Entraremos em contato por DM. Para que a entrega seja efetivada, mande-me um e-mail (geo.luislopes@gmail.com) informando seu Nome, Endereço completo com CEP. Você receberá em seu e-mail o código de rastreamento da mercadoria.

O sorteio foi realizado através do site Sortei.me, especializado em sorteios via Twitter.

Vale lembrar que essa promoção é uma parceria entre a Parêntese Editora e o Blog Geotecnologias Luís Lopes.

Próximo mês, serão 3 livros sorteados e 3 ganhadores!!

Obrigado a todos que participaram e aguardem o novo sorteio!

--
Luís Lopes

Montando banco de dados espacial com a GDAL - Parte 2 (Pacote GDAL).

2011-06-09T09:36:00.018-03:00

Na parte 1 mostrei todo o conjunto de tecnologias(scripts, ambientes, bibliotecas e utilitários), ausentando os detalhes de cada uma delas, assim, privilegiei os usuários que já utilizam essas tecnologias, logo, esses usuários já poderiam experimentar os scripts disponibilizados para construir uma banco de dados espaciais.

As demais partes dessa postagem mostrará cada tecnologia, logo, podem ser aproveitada de forma independente pelo internauta.

A seguir a listagens das tecnologias a serem mostradas

Parte 2 - O Pacote GDAL

Parte 3 - Emulação do ambiente UNIX/SHELL no WINDOWS .

Parte 4 - Obtendo os dados para o banco de dados (Shapefile do PRODES)

Parte 5 - Banco de Dados Postgresql/Postgis

Parte 6 - Ingestão dos dados do PRODES

Parte 7 - Fazendo análise em banco de dados espaciais

O Pacote GDAL

A GDAL (www.gdal.org) é uma biblioteca para tratamento de dados matriciais, porém seu pacote, incorporou a biblioteca OGR, responsável para tratamento de dados vetorias. O pacote GDAL é constituído por bibliotecas e utilitários. Os utilitários para tratamento de dados vetoriais serão abordados nesse postagem.

1) Instalação do pacote da GDAL no WINDOWS:

A GDAL pode ser instalada de mais de uma forma, a mais comum é baixar o FWTools, o 'FW' faz referência ao autor da GDAL, Frank Wardeman. O endereço da FWTools para download é http://trac.osgeo.org/gdal/wiki/DownloadingGdalBinaries, onde podemos baixar o programa de instalação, http://download.osgeo.org/gdal/win32/1.6/gdalwin32exe160.zip.

Para os internautas que desejam seguir as outras partes da postagem, deve-se instalar a GDAL com o uso do programa de instalação OSGEO4W.

O programa OSGEO4W é um programa para instalar os projeto apoiados pela OSGEO(http://www.osgeo.org/) em ambiente WINDOWS, com um detalhe, todos os programas e bibliotecas ficam instaladas no diretório da OSGEO, evitando que a haja conflito com as bibliotecas dos seus aplicativos e sistemas, problema conhecido como "O inferno da DLL"(http://en.wikipedia.org/wiki/DLL_hell).

A página do projeto da OSGEO4W é http://trac.osgeo.org/osgeo4w/ , e o programa para instalação pode ser acessado diretamente no endereço http://download.osgeo.org/osgeo4w/osgeo4w-setup.exe

Na postagem irei mostrar como instalar a GDAL pela OSGEO4W.

1.1) Instalando a GDAL via OSGEO4W.

A OSGEO4W trabalha com dois diretórios, a) diretório de repositórios e b) diretório onde serão instalados os pacotes. Os pacotes podem ser bibliotecas, utilitários(acessados via console "Prompt do DOS") ou até mesmo programas completos, como exemplo, o QGIS.

Os diretórios (repositórios e de instalação) podem ser o mesmo, eu prefiro separá-los, assim, criei o dois diretórios, OSGeo4W_rep (repositórios) e OSGeo4W(instalação).

1.2) Selecionando os pacotes pela OSGEO4W

Execute o "osgeo4w-setup.exe" (coloquei no diretório OSGeo4W_rep) e sigas as instruções:

a) Escolha a opção "Advanced Install"

b) Install from Internet

c) Escolha como Root directory: C:\OSGeo4W (opção padrão). É recomendável deixar a opção para todos os usuários (opção padrão)

d) Local Package Directory: OSGeo4W_rep (opção padrão é o diretório onde está sendo executado o "osgeo4w-setup.exe").

e) Se estiver usando um Proxy, caso onde está numa rede corporativa, pode tentar usar as definições do Internet Explorer (opção Use IE5 Settings) ou definir manualmente seu Proxy (opção HTTP/FTP Proxy)

f) Marque a opção "Commandline_Utilities -> gdal: The GDAL/OGR libray and commandline tools" e continue o assistente de instalação.

Obs(s):

1) O programa de instalação marca automaticamente as dependências de cada item

selecionado, ou seja, para executar o utilitário(The GDAL/OGR libray and commandline tools) é necessário instalar os outros utilitários e as bibliotecas que são utilizadas por esse utilitário.

2) Os softwares livres utilizam as mesmas bibliotecas, assim, temos de forma organizada

as bibliotecas que são utilizadas pelos softwares.

* esse é uma das características do LINUX (meu ambiente preferido para um GIS Workstation)

g) Conclua a instalação selecionado as duas opções (Criar e adicionar o ícone).

2) Utilizando os comandos da GDAL /OGR

O pacote da OSGEO4W cria um console("Prompt do DOS") definido pelo arquivo OSGeo4W.bat que configura a ambiente para todos os pacotes por ele instalado. Alguns pacotes são executados na"linha de comando", ou seja, sem uma interface gráfica, esse é o caso dos utilitários da GDAL.

Cada sistema operacional possui seus comandos e a combinação desses com os comandos da GDAL proporciona uma poderosa ferramenta de manipulação e análise dos dados espaciais.

Na etapa 3 irei mostrar alguns comandos do sistema operacional de minha preferência, o UNIX/LINUX, e como obter esses comandos no WINDOWS.

Lembrando que temos duas tipos de utilitários, a GDAL (imagem) e o OGR (vetores), nessa postagem iremos abordar os utilitários para vetores.

A descrição dos comandos da OGR e os formatos suportados são descrito na página http://www.gdal.org/ogr/index.html.

Como exemplo, vamos utilizar o shapefile PDigital2010_00165_shp.zip do PRODES do ano de 2010, http://www.dpi.inpe.br/prodesdigital/dadosn/2010/.

Descompacte o shapefile, execute o console da OSGEO e digite o comando:

cd [c:\NOME DO DIRETÓRIO]

2.1) ogrinfo:

Mostra informações dos vetores (http://www.gdal.org/ogrinfo.html).

2.1.1) Ver os formatos suportados:

No console da OSGEO, digite "ogrinfo --formats"

2.1.2) Ver as informações de um shapefile (sumário)

No console da OSGEO, digite:

ogrinfo -so PDigital2010_00165_pol.shp PDigital2010_00165_pol

Observações:

2.1.2.1) A OGR possui o conceito de datasource para a fonte de dados espaciais, nessa fonte, pode ter um ou mais layers (camadas). O formato shapefile é um exemplo de um datasource contendo apenas um layer. Aplicando o conceito da OGR para mostrar as informações de um shapefile usamos dois argumentos, nome da fonte ( PDigital2010_00165_pol.shp) e nome da camada PDigital2010_00165_pol.

2.1.2.2) O utilitário ogrinfo, mostra todas as informações da camada, inclusive as coordenadas de cada feição no formato WKT(http://en.wikipedia.org/wiki/Well-known_text).

Para mostrar o sumário da camada foi utilizado a opção '-so'.

2.1.3) Ver as informações das duas primeiras feições de um shapefile

No console da OSGEO, digite:

ogrinfo -q -geom=NO -where "fid < 2" PDigital2010_00165_pol.shp PDigital2010_00165_pol

As opções utilizadas na camada foram:

a) -q: Não mostrar o sumário

b) -geom=NO: Não mostrar o campo geométrico

c) -where "fid < 2" : Mostrar as feições que possuim o campo fid (Feature Index) menor que 2.

Obs(s).:

- A expressão da clausula where é um string(sequência de texto) definido entre aspas duplas(").

- A 1a feição da camada inicia com fid igual a 0(zero)

2.2) ogr2ogr:

Cria novas camadas a partir de uma fonte : (http://www.gdal.org/ogr2ogr.html)

Pode ser utilizada para várias finalidades, como, uma simples conversão de formatos de arquivos vetoriais, como exemplo, inserir um shapefile num banco de dados Postgres/Postgis.

Como exemplo, iremos gerar um KML a partir do shapefile PDigital2010_00165_pol.shp.

As informações do KML serão filtradas, sendo:

1) Feições do tipo 'desmatamento' (sprclasse = 'desmatamento').

2) Os campos a serem mostrados no KML serão satelite e area.

O KML utiliza o sistema de referência EPSG: 4326 ().

No console da OSGEO, digite:

ogr2ogr -where "sprclasse = 'desmatamento'" -select "satelite,area" -t_srs "EPSG:4326" -f KML PDigital2010_00165_pol.kml PDigital2010_00165_pol.shp

A aplicação dos utilitários da GDAL ganha uma importância quando temos que fazer processamentos em vários arquivos, para isso, temos que usar os recursos dos comando do sistema operacional, assunto esse a ser tratado em outras etapas.

Até a próxima parte.

Luiz Motta

motta.luiz@gmail.com

Promoção do mês de Junho do GeoLuisLopes/Parêntese Editora

2011-06-08T14:09:00.002-03:00

Olá pessoal,

Como anunciado anteriormente, todo mês teremos um sorteio aqui no Blog, fruto de uma parceria firmada com a Parêntese Editora. A novidade desta vez é que serão 2 sorteados (um para cada livro). Os livros sorteados são:

1º sorteado - Calibração Absoluta de Sensores Orbitais
2º sorteado - Novas Tecnologias para o Monitoramento

Para participar é bem simples! Basta ser usuário do Twitter, seguidor de @geoluislopes e @parentese e escrever um tweet (copiar e colar) ou dar RT (retweet) na mensagem abaixo:

GEOsorteio | Calibração Absoluta de Sensores Orbitais e Novas Tecn. p Monitoramento. RT e siga @geoluislopes @parentese http://kingo.to/EOo

O sorteio será realizado no dia 17 de junho de 2011, às 21h (horário de Brasília).

Regulamento:

O participante da promoção deve residir no Brasil.
A participação é voluntária e gratuita.
O participante deve ser seguidor de @geoluislopes e @parentese e dar RT (retweet) na mensagem com o link do kingo.to da promoção.
O vencedor será escolhido por meio de sorteio realizado pelo http://sorteie.me/. O resultado será anunciado no blog e no Twitter.
Divulgado o resultado, o vencedor será contactado via DM (Mensagem Direta) no Twitter para solicitar os dados do mesmo. Caso não responda em até 48hs, perderá o prêmio e será realizado outro sorteio. O envio será de responsabilidade da Parêntese Editora.
A data e hora do sorteio podem ser alteradas caso seja necessário.
Perfis exclusivos para divulgação de promoções serão desclassificados.

Sobre o RT:

- Todos os retuítes, manuais ou automáticos, são contados. Desde que contenham o link kingo.to e a mensagem original em sua composição.

Boa sorte a todos!!

Obter elevação do terreno usando o QGIS

2011-05-26T20:38:00.001-03:00

Olá pessoal,

Recentemente lancei uma pequena enquete em meu perfil no twitter (@geoluislopes) para consultar sobre qual software faria o próximo tutorial: ArcGIS ou Quantum GIS? Então, a pedidos, farei mais um sobre o QGIS.

O tutorial de hoje é bastante interessante! Trata-se de obter dados de elevação. Recentemente fiz um tutorial semelhante utilizando o ArcGIS (calcular altitude de pontos através de um DEM ou TIN) e devido ao sucesso deste material, resolvi fazer também para o QGIS. Antes de começar o procedimento, é válido ressaltar que este plugin obtém a elevação usando o serviço do Google Maps.

O tutorial está disponível em PDF, clique aqui.

Deixem seus comentários e contribuições.

Abraços e até a próxima!

Luís Lopes

Twitter: @geoluislopes

E-mail: geo.luislopes@gmail.com

Livros de Geo e S.R. em promoção para os participantes do XV SBSR

2011-05-26T00:31:00.003-03:00

Olá pessoal,

Dando continuidade entre a parceria entre Geotecnologias Luís Lopes e a Parêntese Editora, gostaria de anunciar uma nova promoção da Parêntese:

"Dentro da política de apoio educacional a comunidade de Sensoriamento Remoto e Geoprocessamento no Brasil, a Parêntese Editora oferece aos participantes do XV SBSR uma oportunidade fantástica para aquisição de livros com 25% de descontos, frete gratuito (envio via registro módico) e pagamento em até 3 vezes sem juros no cartão de crédito.

Para participar, basta comprar os livros participantes da promoção através do site da editora e informar o código SBSR no fechamento do pedido. Promoção válida de 25/05/11 a 30/05/11 ou enquanto durarem os estoques específicos para a promoção.

Se o livro que você deseja comprar não estiver cadastrado no site, envie um e-mail para comercial@parentese.com.br e providenciaremos o cadastro da obra para aquisição dentro da promoção (se a obra estiver disponível no distribuidor local e/ou sujeito a prazo de entrega). Promoção válida exclusivamente para compras através do site www.parentese.com.br

Em junho começa a captação de obras para publicação em 2012, sem custos aos autores e com pagamento de direitos autorais. Se você tem material acadêmico para produção, entre em contato com o Editor de Projetos Acadêmicos através do e-mail academicos@parentese.com.br e apresente seu projeto/ideia. Assim o Editor entrará em contato e poderá ir ao seu encontro, em qualquer lugar do Brasil, para negociar a publicação da obra."

Carla Aldeia

Marketing Editorial

Parêntese Editora

www.parentese.com.br

________________________________________

Conhece alguém que participou do XV SBSR? Divulgue esta promoção!

Abraços a todos,

Luís Lopes

Twitter: @geoluislopes

E-mail: geo.luislopes@gmail.com

Resultado do 1º sorteio - Parêntese Editora e GeoLuisLopes

2011-05-20T20:06:00.001-03:00

Olá pessoal,

Já temos o sortudo que levou os livros Sensoriamento Remoto no Ambiente e Sensoriamento Remoto no Estudo da Vegetação.

O ganhador foi: Zenaldo Rodrigues @Zenaldovr

Aqui o link do resultado: http://sorteie.me/1M6Sii

Parabéns Zenaldo! Entraremos em contato por DM. Para que a entrega seja efetivada, mande-me um e-mail (geo.luislopes@gmail.com) informando seu Nome, Endereço completo com CEP. Você receberá em seu e-mail o código de rastreamento da mercadoria.

O sorteio foi realizado através do site Sortei.me, especializado em sorteios via Twitter.

Vale lembrar que essa promoção é uma parceria entre a Parêntese Editora e o Blog Geotecnologias Luis Lopes.

Próximo mês, também serão 2 livros sorteados, diferentes destes sorteados acima (serão anunciados no início do mês), e a novidade que serão 2 ganhadores!!

Obrigado a todos que participaram e aguardem o novo sorteio!

--
Luís Lopes

Conheça a nova parceria do Blog e participe!

2011-05-04T00:02:00.007-03:00

Olá pessoal, finalmente chegou 4 de maio! Chegou o fim do suspense!

É com grande satisfação que o Blog Geotecnologias Luís Lopes em parceria com a Parêntese Editora lança uma nova promoção, onde, todo mês, até dezembro, serão sorteados livros de Geoprocessamento e Sensoriamento Remoto.

De início, neste mês de maio, serão sorteados os livros Sensoriamento Remoto no Ambiente, de John Jensen (ver descrição do livro aqui), e Sensoriamento Remoto no Estudo da Vegetação, de Flávio Jorge Panzoni (ver descrição aqui), e terá UM único ganhador.

Para participar é bem simples! Basta ser residente no Brasil, usuário do Twitter, seguidor de @geoluislopes e @parentese e dar RT (retweet) ou escrever um tweet (copiar e colar) na mensagem abaixo:

Siga @geoluislopes e @parentese RT e concorra aos livros Sens. Remoto no Ambiente e Sens. Remoto no Estudo da Vegetação http://kingo.to/Avu

O sorteio será realizado no dia 20 de maio de 2011, às 20h.

Regulamento:

A participação é voluntária e gratuita.
O participante deve ser seguidor de @geoluislopes e @parentese e dar RT (retweet) na mensagem com o link do kingo.to da promoção.
O vencedor será escolhido por meio de sorteio realizado pelo http://sorteie.me/. E o resultado será anunciado no blog e no Twitter.
Divulgado o resultado, o vencedor será contactado via DM (Mensagem Direta) no Twitter para solicitar os dados do mesmo. Caso não responda em até 48hs, perderá o prêmio e será realizado outro sorteio. O envio será de responsabilidade da Parêntese Editora.
A data e hora do sorteio podem ser alteradas caso seja necessário.
Perfis exclusivos para divulgação de promoções serão desclassificados.

Sobre o RT:

- Todos os retuítes, manuais ou automáticos, são contados. Desde que contenham o link kingo.to e a mensagem original em sua composição.

Boa sorte a todos!!

Montando banco de dados espacial com a GDAL - Parte 1

2011-04-30T13:02:00.003-03:00

Olá pessoal,

A postagem deste mês mostra uma atividade de geoprocessamento pouco utilizada pelos usuários, a execução de tarefas automatizadas com o uso de scripts.

Os softwares de Geoprocessamento, apesar de oferecerem mecanismos de automatização de tarefas, possuem algumas limitações, e ainda, tende a ter uma linguagem própria, o que restringe seu uso ao próprio software.

No estudo de caso vai ser mostrado os scripts, no caso o SHELL do LINUX, que chamam os utilitários da GDAL/OGR que manipulam dados vetoriais.

Para usuários do WINDOWS é mostrado como baixar as utilitários para termos os mesmos comandos do LINUX.

Somando-se os recursos do Shell e os utilitários da GDAL/OGR , vamos criar um banco de dados espacial, com todos os shapefiles do PRODES o que seria uma atividade "humanamente" muito difícil.

Os scripts, além de automatizar as tarefas, tornam o trabalho mais controlado, garantido a padronização, e por fim, o script serve de documentação da atividades.

Por ser um conteúdo extenso, o mesmo foi dividido em partes. Esta é somente a PARTE - 1.

Este estudo de caso foi a montagem do banco de dados espacial, em PostgreSQL/PostGIS, dos dados do PRODES obtidos por shapefile de cenas em cada ano.

O resultado foi uma tabela com os desflorestamentos dos anos de 2000 a 2010, com os campos:

ogc_fid: Campo serial com o identificador da feição;
the-geom: campo geométrico do tipo polígono com o sistema de referência geográfico com o Datum SAD69;
ano: campo texto com o ano do desmatamento;
pathrow: campo texto com a órbita-ponto;
anoi: campo inteiro com o ano do desmatamento.

As fases de obtenção dos dados (shapefiles), verificação dos campos e valores e ingestão dos dados foram feitas por comandos/utilitários do UNIX e shell scripts chamando os utilitários da GDAL/OGR.

A seguir, a configuração do ambiente para executar os comandos e scripts.

1. GDAL/OGR (www.gdal.org):

A OGR é a biblioteca para manipulação de dados espaciais e utiliza os padrões da OGC. Sua instalação pode ser feita baixando os arquivos diretamente no site ou utilizando o instalador da OSGEO4W (http://trac.osgeo.org/osgeo4w/).

A OSGEO4W é um utilitário que baixa os projetos de software livre apoiados pela OSGEO e também diversas bibliotecas.No estudo de caso, utilizei os utilitários da GDAL/OGR instalados pela OSGEO4W.

2. Ambiente Shell:

Para emular o ambiente SHELL no Windows foram utilizados:

2.1 OSGEO4W possui o pacote MSYS que possui alguns utilitários do SHELL

2.1.1 Diretório: C:\OSGeo4W\apps\msys\bin

2.1.2 Comandos a serem utilizados:

2.1.2.1 sh: Interpretador do SHELL

2.1.2.2 grep e egrep: Procura palavras num fluxo de dados (Ex.: arquivos)

2.2 Para completar os utilitários do SHELL foi instalado alguns módulos do GnuWin32 (http://gnuwin32.sourceforge.net/packages.html)

2.2.1 Diretório: C:\_gnu_windows\package\bin

2.2.2 Em caso de proxy definir a variável 'http_proxy' set http_proxy=[endereço]:[porta]
2.2.3 Módulos: coreutils (comandos básicos), unzip (decompacta arquivos ZIP) e wget (baixa arquivos da Internet).

3. Os pacotes devem ser executados no console (prompt do DOS) com as definições de caminho e variáveis.

Definição da console está no arquivo 'OSGeo4W_PRODES.bat'

4. Foram criados Shell Script's (*.sh) para processamentos em lote.

Para executar os scripts utiliza-se o interpretador Shell, comando 'sh'

Exemplo: sh unzip_all.sh

Preparação para os dados espaciais do PRODES

Para seguir os próximos passos faça o download dos scripts, clique aqui.

Foram obtidos os shapefiles do PRODES por cena dos anos 2000 até 2010: http://www.dpi.inpe.br/prodesdigital/dadosn

Para cada ano foi criado um diretório de trabalho

G:\_bd_prodes: prodes2000, prodes2001,..., prodes2010

Listagem dos Scripts:

linkzip_all_need_link.sh

Utilização: sh linkzip_all_need_link.sh [index.html] [endereço dos ZIP´s]
index.html: Arquivo obtido com: wget [endereço dos ZIP´s]
Saída: Arquivo texto com a listagem dos endereços dos ZIP´s

unzip_all.sh

Utilização: sh unzip_all.sh
Saída: Descompacta todos os arquivos do diretório corrente

exist_CRS-in-shapefile_all.sh

Utilização: sh exist_CRS-in-shapefile_all.sh [String de Projeção]
Saída: Listagem na tela dos shapefiles com a ocorrência da projeção

exist_type-geom-in-shapefile_all.sh

Utilização: sh exist_type-geom-in-shapefile_all.sh [Tipo de geometria]
Saída: Listagem na tela dos shapefiles com a ocorrência da geometria

exist_field-in-shapefile_all.sh

Utilização: sh exist_field-in-shapefile_all.sh [Campo]
Saída: Listagem na tela dos shapefiles com a ocorrência do campo

exist_class-in-shapefile_all.sh

Utilização: sh exist_class-in-shapefile_all.sh [Campo] [Valor]
Saída: Listagem na tela dos shapefiles com a ocorrência do valor

ogr_get_values-of-field-in-shapefile.sh

Utilização: sh ogr_get_values-of-field-in-shapefile.sh [shapefile] [campo]
Saída: Listagem na tela com os valores do campo

ogr_get_values-of-field-in-shapefile_all.sh

Utilização: sh ogr_get_values-of-field-in-shapefile_all.sh [campo]
Saída: Listagem na tela com os valores do campo p/ todos os shapefiles

ogr_get_values-of-field-in-shapefile_by_list.sh

Utilização: sh ogr_get_values-of-field-in-shapefile_by_list.sh [arquivo] [campo]
Saída: Listagem na tela com os valores do campo p/ todos os shapefiles listados no arquivo

ogr_ingest_desflorestamento_postgis_all_p1.sh

Utilização: sh ogr_ingest_desflorestamento_postgis_all_p1.sh [ano]
Saída: Ingestão dos shapefiles (classe desflorestamento) no postgres/postgis
Existe variações do script (*_pN.sh) conforme o tipo de shapefile
Deve-se configurar a connecção do POSTGRES dentro so script.

Definição da tabela dos dados do PRODES

Os shapefiles do PRODES possuem várias classes e cobrem toda a a área da cena.

Todos os shapefiles foram inseridos em uma única tabela no banco de dados espacial, aplicando-se o filtro da classe de 'desflorestamento' nos shapefiles e criando apenas 4 campos: ogc_fid(index da tabela), ano (texto), anoi (inteiro), pathrow (definição da cena – tipo texto) e o campo Geométrico (polígono).

Ocorreram variações do filtro em função das mudanças do levantamento de cada ano, como exemplos:

2010: ano = '$1' and class_name = 'DESFLORESTAMENTO'

* $1 = 2010

2005: sprclasse = 'desmatamento'

* O ano e o pathrow foram obtidos pelo nome do arquivo
Exemplo: Pdigital2005_00157_359_pol → ano = 2005 e pathrow = 00157

Ingestão dos dados – Exemplo 2010

Baixar os ZIP's

1 - Obtenção dos 'Index.html' (Figura 1) do endereço http://www.dpi.inpe.br/prodesdigital/dadosn/2010:

Figura 1. Obtenção dos 'Index html'

Comando no console:

wget http://www.dpi.inpe.br/prodesdigital/dadosn/2010/

2 - Criação da listagem dos endereços dos ZIP´s (Figura 2):

Figura 2. Criação do arquivo dos endereços dos ZIP's

Comando no console:

sh linkzip_all_need_link.sh index.html "http://www.dpi.inpe.br/prodesdigital/dadosn/2010/" | grep _shp.zip > list_shape.txt

3 - Baixar arquivos zipados (Figura 3):

Figura 3. Baixar arquivos zipados.

Comando no console:

wget --input-file=list_shape.txt

4 - Descompactar os 97 arquivos baixados (Figura 4):

Figura 4. Descompactando os arquivos

Comando no console:

sh unzip_all.sh

5 - Alguns comandos do Shell e da OGR

5.1 - Ver se todos os shapefiles possuem projeção (contar os arquivos PRJ) (Figura 5):

Figura 5. Verificando os arquivos .PRJ

Comando no console:

ls -l *.prj | wc -l

5.2 - Ver as informações resumidas de um shapefile (utilitário da OGR) (Figura 6):

Figura 6. Resumo das informações dos shapefiles

Comando no console:

ogrinfo -so PDigital2010_00165_pol.shp PDigital2010_00165_pol

6 - Verificar se os Shapefiles tem a mesma projeção (Figura 7):

Figura 7. Verificando a mesma projeção nos shapefiles

* Foi feito o filtro na saída (uso do GREP), para listar apenas quem não tem projeção escolhida (': 0'), como não retornou nenhum elemento, todos possuem a projeção pesquisada.

Comando no console:

sh exist_CRS-in-shapefile_all.sh GCS_South_American_169 | grep ": 0"

7 - Verifica se os shapefiles são polígonos (o mesmo raciocínio anterior) (Figura 8):

Figura 8. Verificando a mesma geometria

Comando no console:

sh exist_type-geom-in-shapefile_all.sh Polygon | grep ": 0"

8 - Verifica se os shapefiles possuem o campo pathrow (Figura 9):

Figura 9. Verificação do campo pathrow nos shapefiles

Comando no console:

sh exist_field-in-shapefile_all.sh pathrow | grep ": 0"

9 - Verifica se os shapefiles possuem campo ano com valor 2010 (Figura 10):

Figura 10. Verificação do campo 2010 nos shapefiles

Comando no console:

sh exist_class-in-shapefile_all.sh ano 2010 | grep ": 0"

10 - Como todos os shapes foram conferidos (projeção, geometria, campos,...) é feita a ingestão no PostgreSQL/PostGIS (Figura 11):

Figura 11. Ingestão dos shapefiles no PostgreSQL/PostGIS

* O Script cria a tabela no Postgres se não existir.

Comando no console:

sh ogr_ingest_desflorestamento_postgis_all_p1.sh 2010

11 - Após inserir os primeiros dados, utiliza-se o POSTGRES para criar o campo 'anoi' do tipo inteiro e depois faz a atualização dos valores desse campo , fazendo a conversão dos valores do campo 'ano' (tipo texto) para inteiro com a instrução SQL :

UPDATE public.prodes

SET anoi = ano::integer

WHERE anoi IS NULL;

* Obs.: O campo 'ano' nos shapefiles é do tipo texto

A tabela do PRODES passa a ter apenas os desmatamento do ano corrente.

A utilização de SHELL script seria uma postagem a parte, mas a ideia é poder mostrar quando scripts são necessários, logo, para quem tem necessidade de tarefas repetitivas em dados vetoriais, pode aproveitar os scripts utilizados nesse artigo e fazer as suas adaptações.

O artigo procurou mostrar uma operação real, em toda sua fase, utilizando os scripts para o ano de 2010 do PRODES. As próximas etapas mostrarão as investigações dos dados, ou seja, perceber as diferenças de cada ano do PRODES, e assim, as mudanças dos scripts para poder inserir os dados no banco. Também teremos mais comentários sobre os comandos dos utilitários da GDAL/OGR, assim, vamos populando o banco de dados começando do ano mais recente até o ano de 2000.

A interação de vocês é fundamental, foram várias técnicas ilustradas, assim, sintam-se a vontade para perguntarem, tirarem suas dúvidas ou curiosidades.

O que acharam? Deixem seus comentários...

Até a segunda parte!

Luiz Motta
--
luiz.motta@gmail.com