Re: [QGISBrasil:4393] Problema na elaboração de mapa de declividade

[Leônidas Descovi Filho]

Olá Marcelo,

Bem vindo.

Não entendi claramente qual é a dúvida dentro do QGIS, poderia explicar melhor seu problema/dúvida?

Abs.

Em 22 de maio de 2013 16:23, <marcelo...@hotmail.com> escreveu:

Boa tarde a todos! Está é minha primeira publicação aqui, e  gostaria de parabenizar pela ideia de criação de um grupo destinado a difusão de conhecimentos sobre o QGIS.
Sou iniciante em SIG's e estou tendo dificuldades para elaboração de um mapa temático de declividade. Os primeiros passos realizados foram: Obtenção das imagens SRTM/TOPODATA, recorte da área de interesse, Raster > Análise > MDE modo: declividade (em porcentagem), com a imagem gerada, utilizei  o estilo mapa de cores, e criei os estilos manualmente seguindo a classificação do INCRA, para declividade(plano, suave ondulado, moderadamente ondulado...)Porém esse se trata de um arquivo para visualização, sendo necessários outros procedimentos para que seja possível calcular a área de cada uma das classes unicamente a partir do QGIS. Lendo o material do INCRA, ele fala em RECLASSIFICAR a imagem. A imagem com os valores de declividade foi gerada, mas não temos os
valores agrupados por classe dentro da imagem, o que temos é uma visualização dessas classes criada a partir do Mapa de Cores. Sendo necessario atribuir novos valores para cada um dos pixels, e é ai que entra a parte em que se tem que utilizar a calculadora raster, se alguém ja leu essa apostila do INCRA e entendeu o processo teria como me explicar melhor?
Em outro material que pesquisei tinha esse processo de forma mais facil utilizando a Reclassificação de dados raster com QGIS > SEXTANTE > GRASS, porém ao tentar utilizá-lo deu o erro :

Um erro ocorreu enquanto executava o seguinte código Python: ...(com várias linhas).

Então é isso pessoal, agradeço pela atenção de vocês . (Obs: estou trabalhando com os dados Topodata, para gerar um mapa de declividade da Serra da Meruoca-CE(maciço residual).

Versão do QGIS
1.7.0-Wroclaw

-- 
Att.
Marcelo Henrique Viana Soares
-------------------------------------------------------------------

Universidade Estadual Vale do Acaraú - UVA
Centro de Ciências Humanas - CCH
Curso de Geografia

Currículo Lattes: http://lattes.cnpq.br/5200625061225143
Facebook: facebook.com/marcelo.h.viana

e-mail: marcelo...@hotmail.com
Celular: (88) 9230-9742

-------------------------------------------------------------------

--
--
---
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.
 
 

--
Leonidas

[Edson Bastos]

Olá pessoal,

Marcelo, eu trabalho no INCRA aqui em Aracaju-SE e a gente utiliza esse material que foi produzido pelo pessoal de Salvador e Minas Gerais, inclusive o Helder Gramacho e o Israel Ely fazem parte desta comunidade no google grupos. Essa reclassificação é uma transformação da matriz resultante pelo agrupamento dos pixels que fazem parte de cada intervalo de classe.

Pra gente compreender melhor, esse caminho metodológico, teoricamente, é semelhante ao que fazemos no Excel quando utilizamos a função "SE".  Imaginemos uma coluna chamada declividade. Cada linha tem um valor específico, mas nós precisamos agrupá-los. Então a gente cria uma segunda coluna "B" e estabelece uma função com as seguintes expressões aninhadas:

"Para qualquer célula da coluna A, considere o seguinte:

 - Se for menor ou igual a 2 será "Classe 1";

- Se for maior que 2 e menor ou igual a 5, será "Classe 2";

- Se for maior que 5 e menor ou igual a 10, será Classe 3;

- E assim segue até chegar ao último intervalo.

Esta explicação eles destacam nas páginas 79 e 80, no texto em vermelho, na 5ª versão, disponível no site do incra, nesse link: http://www.incra.gov.br/index.php/servicos/publicacoes/manuais-e-procedimentos

Como você colocou essa reclassificação é feita para que a matriz resultante tenha o total de valores correspondente ao número de classes para aquela região, ao invés de dezenas ou centenas de valores individuais.

Eu ainda não tentei esse processo com o GRASS, não vou poder ajudar muito nesse aspecto. Mas estamos trabalhando aqui seguindo a apostila e tá rodando beleza.

Quando você processa seguindo esses passos, ocorre algum tipo de erro? Se achar interessante, poderia enviar uma captura de tela pra comunidade aqui dar uma olhada.

 De qualquer maneira, depois da integração do sextante ao qgis, as possibilidades aumentaram exponencialmente, e certamente existem outras opções que levam ao mesmo resultado. Se eu descobrir algo, postarei aqui, ok?

Abraço.

Edson Bastos.

[marcelo...@hotmail.com]

Tem outra coisa que não entendi Edson Bastos. Na apostila quando ele fala para identificar o valor mínimo e máximo da declividade(em metadados)o meu valor aparece:

Val Máx: 219.0361785889

Val Min: -9999.0000000000

Acho que tem alguma coisa errada. E quando fiz toda a fórmula minuciosamente(usei as sete classes) e apliquei, a imagem que foi gerada não apresenta os valores min e max. qual terá sido o problema?

[marcelo...@hotmail.com]

AÊ! Agora sim. Fiz o seguinte: como suspeitei do valor mín(-999), resolvi elaborar a declividae de novo : Raster > Análise > MDE modo: declividade (em porcentagem) ai sim apareceu um valor mín aceitavel de 0,003%, o valor max tava ok. Logo após comecei a operação utilizando a calculadora raster para a reclassificação das classes(utilizei as sete), logo em seguida apliquei as outras etapas contidas no manual do INCRA(5ªversão) e saiu quase beleza, o problema agora é que apareceram algumas manchas brancas(acho que são pixels que tem valores de altitude acima do que foi colocado na expressão na calculadora,por isso não foram classificados), segue em anexo a imagem para observação e esclarecimento da minha dúvida em relaçao a essas manchinhas brancas.

[Israel Ely]

Vc pode com o identificador de feições e clicar em cima dos campos em branco e verificar se tem valores atribuídos a eles. Caso contrário, se for um valor acima do limite da fórmula que colocamos na apostila, vc pode alterar o limite superior..

Atenciosamente

Israel Ely

INCRA

[marcelo...@hotmail.com]

Israel Ely, esses espaços em branco contem valores atribuidos, então no caso como faria para aumentar o limite superior na fórmula?

[marcelo...@hotmail.com]

Outra dúvida, no meu caso precisei fazer a operação para utilizar as 7 categorias, e utilizei a fórmula:

(declividade_1@1 <= 2) + ((declividade_1@1 > 2 AND
declividade_1@1 <= 5) * 2) + ((declividade_1@1 > 5 AND
declividade_1@1 <= 10) * 3) + ((declividade_1@1 > 10 AND
declividade_1@1 <= 15) * 4) + ((declividade_1@1 > 15 AND
declividade_1@1 <= 45) * 5) + ((declividade_1@1 > 45 AND
declividade_1@1 <= 70) * 6) + ((declividade_1@1 > 70 AND 
declividade_1@1 <= 70) * 7) + ((declividade_1@1 > 70) * 8)

Só por motivo de segurança ela está correta?Por que fiquei com dúvida nesta última linha.(Segui o modelo da apostila do INCRA,porém na apostila o exemplo só utiliza 5 categorias)

[Israel Oliveira]

Desculpe a demora na resposta, a fórmula contém um erro no final, deve ficar assim:

(declividade_1@1 <= 2) + ((declividade_1@1 > 2 AND 
declividade_1@1 <= 5) * 2) + ((declividade_1@1 > 5 AND 
declividade_1@1 <= 10) * 3) + ((declividade_1@1 > 10 AND 
declividade_1@1 <= 15) * 4) + ((declividade_1@1 > 15 AND 
declividade_1@1 <= 45) * 5) + ((declividade_1@1 > 45 AND

declividade_1@1 <= 70) * 6) + ((declividade_1@1 > 70) * 7)

A ultima categoria é a classe ESCARPADO onde se encaixa toda a declividade acima de 70% ( declividade > 70).

Verifique se no arquivo raster, se os pixel que aparecem em branco após a reclassificação, possuem valor em porcentagem (arquivo de declividade).

Atenciosamente,

Israel Ely

Israel Ely de Almeida Oliveira

--
--
---
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ê recebeu esta mensagem porque está inscrito em um tópico do grupo "qgisbrasil" dos Grupos do Google.
Para cancelar a inscrição neste tópico, acesse https://groups.google.com/d/topic/qgisbrasil/IoYiKtAgmCY/unsubscribe?hl=pt-BR.
Para cancelar a inscrição neste grupo e todos os seus tópicos, envie um e-mail para qgisbrasil+...@googlegroups.com.

[marcelo...@hotmail.com]

Refiz todo o procedimento, e apliquei a fórmula que você indicou e deu certo. Os espaços em brancos sumirar, foram preenchidos pela classe escarpado. Muitíssimo obrigado pela atenção e ajuda Israel Ely. Abraços.
(Só pra demonstrar esse foi o resultado final)

Mapa de Declividade - Serra da Meruoca-CE

[marcelo...@hotmail.com]

Ops! Outra dúvida que queria que você me esclarecesse, é de como citar esta apostila do INCRA(5ªversao).

[Israel Oliveira]

Vou consultar o alexandre barberi q e mais apto nos assuntos academicos.... Rsrs, te respondo amaha

[Helder Gramacho]

Olá Pessoal,

só hoje via  discussão e como ela já foi resolvida, fica a dica: o Anderson Medeiros publicou um tutorial que ensina como realizar a reclassificação por meio do GRASS, é bem mais rápido. Segue o link:

http://andersonmedeiros.com/reclassificacao-dados-raster-sextante-grass/

Hélder

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.
 
 

--

Hélder Gramacho dos Santos

Engenheiro Agrônomo

Mestrando Ciênc. Geod. Tec. Geoinformação-UFPE

agroh...@gmail.com

[Felipe Dantas]

Olá, amigo.

Com toda essa questão, me deparei com um problema (no meu caso) e não vi nada falando na apostila do INCRA e aqui.

É o seguinte, Peguei 4 imagens SRTM - EMBRAPA - (duas na zona 23 e duas na zona 24), fiz o mosaico das imagens e as reprojetei para o sistema UTM das suas respectivas zonas.Fiz todo processo de geração da declividade em porcentagem e agrupamento das classes via calculadora raster.

Tenho a shape de um município piauiense (pois sou do piauí :)) e gostaria de analisar a declividade do mesmo. Porém, o município abrange mais de uma zona UTM (23 e 24).

Gostaria de saber, no caso lá do pessoal que trabalha no INCRA, se um imóvel rural abranger mais de um zona, como é realizado o procedimento de geração da declividade e suas classes já que a mesma deve está em UTM.

Eu tava fazendo separadamente a declividade de cada zona e sua classes, porém ao colocar as duas raster no mapa, percebi uma grande diferenção de declividade. (anexo)

Valeu pessoal pela atenção, se alguém puder me ajudar !

[Eder Mileno S De Paula]

Olá,

Acho que a solução para seu problema é 'esticar' uma das zonas,
provavelmente a zona 24. Ou seja, torna-lá maior.
A melhor forma seria reprojetando a imagem para outra zona.
Esse alargamento de zona pode ser feito com restrições, tendo em vista
que quanto mais distante do meridiano central da zona mais distorções
irão ocorrer. O melhor seria reprojetar as duas imagens para um sistema
de coordenadas único, entretanto deve ser métrico para que se possa
fazer os cálculos no Qgis.

Espero ter contribuido.

Saudações Geográficas

Em 06-06-2013 21:12, Felipe Dantas escreveu:
> al pela atenção, se alguém puder me ajud

--
Eder Mileno S. De Paula
93 9157 2224

Respeite o ciclista.
Uma bike a mais, um carro a menos.
..... __@
...._ \ >_
...(_)/ (_)_______________________
░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░░\__

[Augusto Vieira]

Felipe,

Olhando a imagem percebi que vc utilizou cores diferentes para as faixas de declividade....sugiro, salvar o estilo da primeira e carregar na segunda, dessa forma nao terá a diferença visual apresentada,

Boa sorte

Augusto 

--

[Augusto Vieira]

Veja a faixa 2 - 5....uma esta em amarelo e a outra esta em laranja!

[Felipe Dantas]

Certo Eder, no meu caso, a precisão não é de grande importância, pois estou trabalhando o territorio de um municipio e trabalhando com uma pequena escala  acredito se eu der essa "esticada" não comprometerá o trabalho. Mas, se por acaso eu estivesse trabalhando com grande escala e a precisão do meu trabalho é fundamental, e se por exemplo eu tiver a área de uma fazendo e a mesma estivesse estre dois fusos, como seria o trabalho neste caso específico ?

Certo Augusto, realmente troquei a sequência das cores (propositalmente para que eu pudesse ter uma maior discrepância na imagem). Refiz aqui e coloquei na sequencia correta (anexo), mas ainda continua havendo uma grande diferença entre os dois mosaicos. Me dúvida é: será se fiz algo de errado, pois como áreas próximas podem ter valores bem diferentes de altitude ?; 

--

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.
..... __@
...._ \ >_

...(_)/ (_)___________________________

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

[Helder Gramacho]

Olá Felipe,

Nunca me deparei com uma situação como essa entretanto, cada fuso UTM pode ser prolongado em até 30' um sobre o outro justamente para evitar problemas em trabalhos de campo que ocorram nos limites de fusos, 30' dá aproximadamente 55,5km, dessa forma dificilmente ocorrerá de um imóvel rural se situar fora dessa faixa de prolongamento que acabei de citar.

No caso de um município, em que isso é mais provável de acontecer, o correto seria trabalhar com coordenadas geográficas ao invés de um sistema projetado.

att,

Hélder

--

--
---
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 Dantas]

Vlw, Helder. Não sabia desse detalhe que cada fuso pode ser prolongado em até 30'.

Para o exemplo que eu estava testando, eu fiz isso, "prolonguei" a zona 24 sobre a 23.

Obrigado a todos e bom final de semana !

[Gabriel Brito]

Boa tarde a tod@s.

Estou com um pequeno problema na elaboração do mapa de declividade, 

quando coloco a imagem sobre o imóvel uma pequena parte dele fica pra fora,

vi neste tópico que existe a possibilidade de esticar a imagem, como que faz isso??

Obrigado.

Gabriel Brito da Rocha

Graduando em Geografia

Estagiário do Instituto Nacional de Colonização e Reforma Agrária. 

[Harlan Rodrigo]

Esse MDT é o Topodata? Se for, baixe a região superior, faça o mosaico e depois crie a declividade. Acho a forma mais coerente de se fazer.

Se for um mdt mais específico, existem técnicas para inferir sobre esse área, especialmente se você tiver a malha de hidrografia...

---

Date: Wed, 2 Oct 2013 13:02:24 -0700
From: gabrie...@gmail.com
To: qgisb...@googlegroups.com
Subject: [QGISBrasil:5892] Re: Problema na elaboração de mapa de declividade

[Gabriel Brito]

Esse mdt que estou usando é o Topodata, tentei fazer o mosaico com as duas imagens mas o QGis da o erro que está na imagem.

[Harlan Rodrigo]

Um erro que eu percebi é que você marcou "camada acumulada". Essa opção serve para gerar uma imagem multiespectral. Para um mosaico como o seu, deixe essa opção desmarcada. Como não deu para perceber o caminho de seu arquivo, se for o caso, tire acentos e espaços dos nomes dos arquivos, das pastas e subpastas que estão dos dados, porque isso também dá erro para esses procedimentos.

Veja aí pra ver se funciona...

Até.

---

Date: Wed, 2 Oct 2013 13:46:51 -0700
From: gabrie...@gmail.com
To: qgisb...@googlegroups.com
Subject: [QGISBrasil:5894] Re: Problema na elaboração de mapa de declividade

Esse mdt que estou usando é o Topodata, tentei fazer o mosaico com as duas imagens mas o QGis da o erro que está na imagem.

[Gabriel Brito]

Harlan Rodrigo, agora deu certo =]

muito obrigado pela ajuda.

[Harlan Rodrigo]

Opa, estamos aí!...

---

Date: Wed, 2 Oct 2013 17:50:39 -0700
From: gabrie...@gmail.com
To: qgisb...@googlegroups.com
Subject: [QGISBrasil:5901] Re: Problema na elaboração de mapa de declividade

Harlan Rodrigo, agora deu certo =]

muito obrigado pela ajuda.