Krigagem no SAGA/QGIS

[Pedro Celestino Dantas]

Caros Colegas da comunidade,

Com muita dificuldade e muitos testes consegui fazer uma mapa de dispersão de Nitrato, com base nas informações de 130 poços de captação de Natal/RN, utilizando o módulo de Kigagem do SAGA. O método que escolhi foi a Krigagem Ordinária, que segundo a bibliografia "estima pontos não amostrados. Este ponto é resultado da combinação linear dos valores encontrados na vizinhança próxima"

Não tenho certeza se foi o mais correto (não é fácil para iniciante a geoestatística), no entanto o mapa ficou até bonito.

Ao trazer o Raster para o QGIS percebi que as informação dos níveis de Nitrato que deveriam conter nos pixels não veio, aparecerem valores do RGB.

Onde foi que eu Errei?
--
Pedro Junior
Geógrafo - CREA/RN: 2107854874
Mestre em Engenharia Sanitária e Ambiental
Setor de Monitoramento de Veículos - URBANA
Telefones de Contato
(84) 3232-8767 Trab. - 8817-2927(Oi) e 9671-1277(Tim) Cel.

Respeite o ciclista.

Uma bike a mais, um carro a menos.
..... __@
...._ \ >_
...(_)/ (_)____________________________________
░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

[Leônidas Descovi Filho]

Oi Pedro,

Legal tua estimativa por krigagem, antes de mais nada parabéns pela socialização. Como voce gerou o mapa? poderia nos ensinar/explicar os passos?

Abs

--
--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

--

Leônidas

[Pedro Celestino Dantas]

Caro Leônidas,

Com uma tabela de poços com níveis de nitrato por ano (anexo) adiciono os pontos através da ferramenta "Adicionar camada através de texto delimitado no QGIS.

Salvei em SHP.

No SAGA, talvez vc consiga pelo SEXTANTE, no momento não consigo instalar o plugin, abro o shp com os dados de Nitrato na caixa onde consta as pastas do computador (C, D, E).

Seleciono o shp na árvore (tree) - Ver Imagem1 Anexa

Modules > Geoestatistc > Kriging > Ordinary Kriging (VF Global) - Imagem2

Na caixa que se abre, onde tem a palavra "not set" em points selecione a camada Nitrato e em "not set" abaixo selecione a coluna 2007. Podendo fazer também para os anos de 2008 e 2009 - Imagem3

Dê 2 cliques em Ok nas duas caixas seguinte.Na terceira caixa selecione as coordenadas dos vértices do mapa. Para facilitar eu defini e salvei as coordenadas para as futuras utilizações. Para isso abra o arquivo anexo (Limite Natal.sprm) no Botão "Load"  e Okay - Imagem4

Proto a Krigagem está Pronta.

Dê dois clique no primeiro arquivos gerado que aparecem em "Tree" e na caixa que abre selecione "New".

Não sei porque geram dois ao invés de um.

Sobre o nome do arquivo clique com o botão direito e em "Save As Image" e escolha o formato. Eu gosto do formato Tiff. - Imagem5

Abra o Quntum e Edite seus mapas.

Finalizei o meu conforme anexo (Dispersão Nitrato 2008_2009.jpg). Não está ainda finalizado como eu quero, mas já está 80% concluído.

Segue também o modelo de Layout que fiz.

Espero que consiga fazer os passos.

Como falei no Email anterior não sei se este método de Krigagem é a mais correta, por isso pedi ajuda aos colegas mais experientes.

Espero ter ajudado

[Pedro Celestino Dantas]

Esqueci de enviar o shp dos limites dos bairros de Natal. Anexo

[Felipe Dantas]

Tu salvou a shape de nitrato em qual sistema de coordenadas, Pedro ?

--
--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

--

Felipe Ramos Dantas

Graduando em Tecnologia  em Geoprocessamento - IFPI

Contatos: tim (86) 9921-5742 /  claro (86) 9411-4372

Respeite o ciclista.

Uma bike a mais, um carro a menos.
..... __@
...._ \ >_
...(_)/ (_)___________________________

_________
░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

[Pedro Celestino Dantas]

Sistema de Coordenadas UTM, no Datum Sad69, Zona 25S

[Leônidas Descovi Filho]

Valeu Pedro ótimo tuto.. no sextante é possível sim, só que na minha tentativa acho que configurei algo incorretamente e o algoritmo ficou rodando um tempão e não pude concluir o processo pois estava de saída..

Usei algoritmo Ordinady Kriging (Global) que está aninhado no menu Analysis --> Sextante toolbox--> SAGA --> Kriging.

Abs

--
--
---
Você recebeu esta mensagem do Grupo "QGISBrasil"
Para postar neste grupo, envie um e-mail para qgisb...@googlegroups.com
Para cancelar a sua inscrição neste grupo, envie um e-mail paragisbrasi...@googlegroups.com
Para ver mais opções, visite este grupo em
http://groups.google.com/group/qgisbrasil?hl=pt?hl=pt-BR
 
---
Você está recebendo esta mensagem porque se inscreveu no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste grupo e parar de receber seus e-mails, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para obter mais opções, acesse https://groups.google.com/groups/opt_out.
 
 

--

Leônidas