Delimitação de bacias - TauDEM

621 views
Skip to first unread message

Amanda Ap. Lima

unread,
Feb 24, 2016, 12:53:56 PM2/24/16
to qgisbrasil
Olá, boa tarde.

Gostaria de pedir ajuda, porque estou tentando delimitar uma bacia para estudo hidrológico e não estou obtendo sucesso.
Eu faço todos os passos com a utilização do STRM, que no caso não preciso utilizar a ferramenta mosaico, e depois de remover as depressões e os pixels espúrios, consigo gerar a direção do fluxo dos córregos, uso a ferramenta D8 Contributing Area para as informações de escoamento, e por fim a Stream Definition Thresheld pra destaque das nascentes.
Todo esse processo dá certo. O que acontece é que depois que gero o Shapefile do ponto de exutório, a bacia resultante da Segunda Versão-D8 Contributing Area não engloba toda a área de interesse. Ela gera uma sub-bacia pela metade da que realmente deve ser. O que pode estar acontecendo?

Obrigada desde já pela atenção!


Amanda Ap. Lima

unread,
Feb 24, 2016, 1:10:07 PM2/24/16
to qgisbrasil
Será que essa eliminação de depressão e pixels espúrios pode ter influenciado na delimitação da bacia?

Leônidas Descovi Filho

unread,
Feb 24, 2016, 4:47:19 PM2/24/16
to qgisbrasil
Oi Amanda

Qual o valor voce usou para o limite da bacia?

Se puder informe a área da bacia e os passos que voce usou. Posso tentar "repetir" o processo aqui e verificar onde pode estar o problema.

Abs

Leônidas



Em 24 de fevereiro de 2016 15:10, Amanda Ap. Lima <aalim...@gmail.com> escreveu:
Será que essa eliminação de depressão e pixels espúrios pode ter influenciado na delimitação da bacia?

--
Comunidade QGIS Brasil: http://qgisbrasil.org/
---
Você recebeu essa mensagem porque está inscrito no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.
Para postar nesse grupo, envie um e-mail para qgisb...@googlegroups.com.
Acesse esse grupo em https://groups.google.com/group/qgisbrasil.
Para mais opções, acesse https://groups.google.com/d/optout.

Amanda Ap. Lima

unread,
Feb 25, 2016, 5:56:38 AM2/25/16
to qgisbrasil
Bom dia, Leônidas!

Pra área da bacia, vou te mandar um kml*, tudo bem? Encontra-se em anexo.

Muito obrigada!
BACIA_WGS.kml

Vicente

unread,
Feb 25, 2016, 2:09:54 PM2/25/16
to qgisbrasil
Olá Amanda! Talvez você precise usar uma imagem SRTM maior do que essa que você está usando.
Muitas vezes, o rio principal da bacia delimitada nasce muito além da imagem usada e o TauDEM gera esse erro. Tente pegar uma imagem que englobe uma área maior do que essa que você pretende delimitar (mesmo que tenha que fazer o mosaico de duas ou mais imagens SRTM) e tente de novo.
Depois é só cortar a bacia do tamanho que você precisa e talvez realizar mais algumas operações simples de vetor e raster. Espero ter ajudado.
Boa sorte! 

Leônidas Descovi Filho

unread,
Feb 25, 2016, 3:51:40 PM2/25/16
to qgisbrasil
Oi Amanda

 A sua bacia tem menos que 50 km2, coloco em anexo um srtm 90x90 recortado (a folha inteira em wgs 84 pode ser obtida no link: http://www.relevobr.cnpm.embrapa.br/download/sp/sf-22-z-d.htm )

Veja se consegues rodar as rotinas usando mdt de base este anexo que te envio, em wgs 84 utm fuso 22s.

Vou tentar aqui também. Depois posto o resultado. Vou seguir o tutorial que voce enviou do Jorge.

Abs

Leônidas



SF-22-Z-D_wgs_corte.tif

Leônidas Descovi Filho

unread,
Feb 25, 2016, 3:53:18 PM2/25/16
to qgisbrasil
Outra opção seria usar srtm 30x30m, pois trata-se de uma microbacia hidrográfica.

Abs e boa sorte

Leônidas


Leônidas Descovi Filho

unread,
Feb 25, 2016, 7:05:29 PM2/25/16
to qgisbrasil
Oi Amanda

A sua pergunta inicial está diretamente relacionada no como voce definiu o ponto de exutório. Para facilitar colo aqui um print da minha tela. Voce definiu como ponto de exutório sendo próximo ao círculo vermelho 1 ou 2?
Imagem inline 1

Fiz os procedimentos do tutorial do Jorge, no qgis 2.12 Lyon, usando a imagem srtm 30 metros do EarthExplorer /USGS que coloco neste anexo e informo que aqui funcionaram todas as etapas. A unica ressalva reside na os  algoritmos que se repetem, em ingles e em portugues. Se realizares a mesma sequencia do tutorial do Jorge, ficaria assim:

Remover depressões (gera 1 tif)

Direções de fluxo D8 (gera 2 produtos tif - direcoes e grade declive)

Área de contribuição D8 (gera 1 tif)

Definição de limiar para fluxo (usei 100 para o limiar = threshold o palavrinha complicada sô!) a densidade da drenagem depende desse valor, quanto maior o limiar menor a densidade, e vice versa. Se alguem souber mais detalhes sobre o limiar, favor acrescente e colabore.

Depois destes processos tem-se a condição de visualizar toda rede de drenagem gerada sobre o raster. Nesse momento, voce define o ponto de exutório (de posse dessa visão sistêmica ΅todo" da rede de drenagem). Crie o shape de ponto na mesma projeção, edite e crie o ponto no exutório.

Passada a parte geral, chegou a parte do especifico: Bacia

1 - Marque o ponto de exutório (que sempre se localiza no encontro com outro curso dagua, lago ou mar ou ainda uma bacia piloto, bacia escola, que possui seu próprio ponto de exutório (atrelado a um experimento por ex.).

2 - Área de contribuição (desta bacia em pesquisa)

3 - Definição de limiar (para bacia)

4 - Alcance e fluxo da bacia hidrográfica (este gera vários arquivos tif, shp e txt referentes a bacia em pesquisa).

5 - Editar as propriedades das camadas raster e vetor (conforme o tuto).

6 - Poligonize a bacia (a drenagem vetorial já foi criada anteriormente no passo 4, é dela que vais estrair as nascentes).

Até esse ponto voce vai chegar a algo proximo a isso:

Imagem inline 2

Em vermelho o kml que vc enviou e em azul a bacia geoprocessada no taudem e poligonizada (menu raster).

Veja se essas dicas te ajudam e qualquer questão, poste aqui.

Abs e seja feliz


Leônidas


s24_w049_1arc_v3_cort_utm.tif

Amanda Ap. Lima

unread,
Feb 26, 2016, 12:03:20 PM2/26/16
to qgisb...@googlegroups.com
Leônidas,

Eu vou testar tudo isso nesse final de semana. 
Tenho trabalhado em outra coisa e agora vou voltar nisso.
Qualquer coisa volto a postar aqui.

Muitíssimo obrigada e vou ser muito feliz com certeza! Haha!
Tenha uma ótima tarde!



Você recebeu essa mensagem porque está inscrito em um tópico no grupo "qgisbrasil" dos Grupos do Google.
Para cancelar inscrição nesse tópico, acesse https://groups.google.com/d/topic/qgisbrasil/yCIxgdpc0DE/unsubscribe.
Para cancelar inscrição nesse grupo e todos os seus tópicos, envie um e-mail para qgisbrasil+...@googlegroups.com.

Leônidas Descovi Filho

unread,
Mar 1, 2016, 3:56:48 PM3/1/16
to qgisbrasil
Oi Amanda

Quando fizeres o teste, por favor, poste aqui as tuas impressões e resultados em geral. Seu feedback será muito importante para as pesquisas futuras dentro do grupo.

Abs e seja feliz

Leônidas


Amanda Ap. Lima

unread,
May 8, 2016, 8:35:46 PM5/8/16
to qgisbrasil
Leônidas, tudo bom? Há quanto tempo, não?! 

As pessoas só procuram resolver problemas somente quando se é necessário. Sou mais uma... --'
Testei novamente, e em outro arquivo e não obtive sucesso em nenhum. Não sei o que acontece.. 
Será que é na configuração? Eu passei por todos os passos do tutorial, sem mais e sem menos, fiz algumas vezes e nada. Você pode mandar um print das configurações do TauDEM dentro do QGIS? Estou fuçando pra ver o que está acontecendo.
No ponto do exutório, segundo o tutorial, pode ser colocado em um ponto localizado em uma barragem, por exemplo, por vias de pesquisa. E é o meu caso, que em estudo hidrológico precisamos colocar esse ponto em barragens. Acredito não ser esse o problema, pois onde se coloca o ponto, gera-se uma sub-bacia, ou estou errada? 
#Help I need somebody, help! 

Brincadeiras a parte, muito obrigada desde sempre! :)



Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Processamento Digital

unread,
May 8, 2016, 8:51:13 PM5/8/16
to qgisb...@googlegroups.com
O TauDEM foi desenvolvido para trabalhar com o ArcGIS. Os algoritmos do TauDEM emprestados para utilização no QGIS não trabalham com mais de um ponto de exutório para geração de várias sub-bacias.

Ainda não consegui explorar essa funcionalidade no QGIS. No ArcGIS é possível gerar várias sub-bacias simultaneamente. No QGIS não pesquisei ainda.

Você disse que a área de contribuição não corresponde com a realidade. Envie um print da sua área de estudo e o resultado do processo ou experimente gerar a bacia hidrográfica com os algoritmos do SAGA.

Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Skype: jorgepsantos2002
Twitter: http://twitter.com/jpsantos2002

Amanda Ap. Lima

unread,
May 9, 2016, 9:38:56 AM5/9/16
to qgisbrasil
Olá, Jorge.
Eu estou tentando agora pelo SAGA, ele gera os postos de nascente e exutório.
Meu grande problema está sendo depois do ponto, porque não gera bacia alguma, o que aparece é o ponto e um pixel preto (Isso no TauDEM). A correção seria no Threshold, alterar seu valor uma vez que a sub-bacia é pequena? 
No SAGA gerou as bacias, mas eu não tenho o controle de indicar o ponto de exutório, que no meu caso seria em uma barragem.

Abraço.

Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Leônidas Descovi Filho

unread,
May 9, 2016, 10:00:02 AM5/9/16
to qgisbrasil
Existe a possibilidade de uso do GRASS, o qual obtive sucesso em tentativas e usei em minha tese.
Voce necessita definir várias bacias ou sub-bacias simuotaneamente? acredito que essa ferramenta exita, mas ainda não necessitei, logo não procurei.

Veja se anima usar o grass, já compartilhei um passo a passo aqui no grupo, se não localizar eu reenvio aqui.

Abs e boa sorte no teu trabalho.

Leônidas



Processamento Digital

unread,
May 9, 2016, 10:10:21 AM5/9/16
to qgisb...@googlegroups.com
Mande seu ponto de exutório que eu gero aqui com base no Topodata.

Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Skype: jorgepsantos2002
Twitter: http://twitter.com/jpsantos2002

Em 9 de maio de 2016 10:38, Amanda Ap. Lima <aalim...@gmail.com> escreveu:

Amanda Ap. Lima

unread,
May 9, 2016, 10:14:45 AM5/9/16
to qgisbrasil
Oi, Leônidas!

Então, eu li sobre o GRASS,  só que ele assim como o SAGA, gera as bacias automaticamente. Eu insisto no TauDEM porque ele me dá a opção de colocar o ponto de exutório do local onde necessito. É pra estudo hidrológico de barragem para outorga d´água. Vou tentando aqui, se conseguir posto resultado. Vou tentar no ArcGIS também, ver se esse é o problema.

Muito obrigada!


Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Leônidas Descovi Filho

unread,
May 9, 2016, 11:07:42 AM5/9/16
to qgisbrasil
Não consegui entender qual é seu problema...

O passo a passo que postei aqui ajudou ou só confundiu?

Abs e boa semana!!1

Leônidas



Amanda Ap. Lima

unread,
May 9, 2016, 12:36:00 PM5/9/16
to qgisbrasil
Não! Super ajudou! É que em outra bacia ocorre o mesmo erro de antes.
Numa bacia menor do que aquela ele não gera o limite da bacia, o que ele gera é um pixel preto em cima do ponto criado. Mesmo erro que acontecia no outro arquivo.
Pois o que eu quero, e o que o plugin se propõe a fazer, é gerar uma bacia a partir do ponto em que eu determinar (a grosso modo), mas isso não acontece.
Tenso... 


Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Processamento Digital

unread,
May 9, 2016, 1:23:52 PM5/9/16
to qgisb...@googlegroups.com
Envie o ponto de exutório, Amanda. Verificarei aqui.

Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Skype: jorgepsantos2002
Twitter: http://twitter.com/jpsantos2002

Amanda Ap. Lima

unread,
May 9, 2016, 1:41:20 PM5/9/16
to qgisbrasil
Te enviei por e-mail.

Obrigada!


Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Amanda Ap. Lima

unread,
May 9, 2016, 2:32:35 PM5/9/16
to qgisbrasil
Caros,
No ArcGIS deu certo, a bacia foi gerada normalmente (ilustração em anexo) a partir de qualquer ponto na drenagem.
Não sei o que acontece no QGIS... :/


Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:
print.png

Vicente

unread,
May 10, 2016, 10:29:07 AM5/10/16
to qgisbrasil
Olá Amanda! Talvez você precise usar uma imagem SRTM maior do que essa que você está usando.
Muitas vezes, o rio principal da bacia delimitada nasce muito além da imagem usada e o TauDEM gera esse erro. Tente pegar uma imagem que englobe uma área maior do que essa que você pretende delimitar (mesmo que tenha que fazer o mosaico de duas ou mais imagens SRTM) e tente de novo.
Depois é só cortar a bacia do tamanho que você precisa e talvez realizar mais algumas operações simples de vetor e raster. Espero ter ajudado.
Boa sorte! 

Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Amanda Ap. Lima

unread,
May 12, 2016, 7:13:45 AM5/12/16
to qgisbrasil
Olá, Vicente!

Toda ajuda é válida. Tenho tido sucesso agora. Numa bacia pequena também nem tem muita necessidade de usar o arquivo SRTM para delimitação, vai a mão mesmo, mas valeu o teste.
Obrigada à todos!


Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Adriel Lucena

unread,
May 20, 2016, 1:48:45 PM5/20/16
to qgisbrasil
Olá, Amanda, pelo que percebe você conseguiu solucionar o seu problema. Contudo, só para esclarecer O SIG GRASS além de gerar dados hidrológicos automáticos da bacia hidrográfica, delimita uma bacia usando o ponto exutório definido a partir de um par de coordenadas, essa função é disponibilizada pelo módulo "r.water.outlet" do GRASS, para tanto você deve dispor do MDE da área de estudo, o arquivo "drainage direction" (obtido pelo módulo r.watershed) e o par de coordenadas que definem o exutório da sua bacia. Lembrando que não é necessário que seja feita correções nas áreas com depressões, o r.watershed usa o algorítmo de menor custo para contornar essas falhas (Sinks). Além disso o r.watershed gera subbacias na opção basin desse módulo que também poderia ser uma saida para seu problema, você exportaria as bacias em formato SHP para o QGIS (ou faria essa ação no GRASS) e selecionaria as subbacias do seu interesse (sua área de estudo)

att
Adriel Lucena



Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:

Amanda Ap. Lima

unread,
May 23, 2016, 7:34:07 AM5/23/16
to qgisbrasil
Oi, Adriel!

Obrigada. Em breve vou testar, quem sabe essa operação seja mais viável...
Valeu. Abraço!


Em quarta-feira, 24 de fevereiro de 2016 14:53:56 UTC-3, Amanda Ap. Lima escreveu:
Reply all
Reply to author
Forward
0 new messages