# Reunião de orientação (25/06/2026) - João **Channel:** Ricardo Michael Pinheiro Silveira **Source:** https://www.youtube.com/watch?v=BcTPhuTv7Mk **Transcript page:** https://www.withtranscript.ai/video/BcTPhuTv7Mk ## Chapters - 0:09 — Visão geral do fluxo operacional - 3:56 — Geração dos atributos geomorfométricos - 17:32 — Quantificação das bacias por estatísticas - 26:02 — Análise multitemporal do uso do solo - 29:34 — Seleção de datas para turbidez - 35:09 — Processamento de imagens e turbidez ## Transcript **[0:09] Speaker A:** João, eu vou fazer aqui essa gravação para te dar um apoio na parte operacional do desenvolvimento do trabalho. Aí, claro, qualquer dúvida, tem um monte de etapas aqui que eh eu vou mostrar só o básico, daí questões assim mais detalhadas, você pode me perguntando que eu vou auxiliando, tá? Então vou dividir aqui em três etapas que são para cada um dos teus três objetivos específicos. Primeiro vai ser a parte aqui do levantamento de dados da geomorfometria. Depois vou falar um pouco sobre o levantamento dos dados de uso do solo. E o terceiro vai ser a parte do das imagens de satélite pro cálculo de turbidez, tá? Então a primeira questão que quero destacar é o seguinte. você vai ter um grande conjunto de de dados e você vai ter muitas tabelas. Então, o que eu vou passar aqui são instruções para você ter essa base de dados organizada, estruturada e já ter as tabelas prontas, porque também tendo essas planilhas prontas, você já vai ter uma série de aí de resultados **[1:10] Speaker A:** quantificados que posteriormente, depois que terminar tudo isso, fazer o levantamento, aí vai pra fase da interpretação e análise dos resultados, tá? Então, vou começar falando aqui da parte da geomorfometria e quais vão ser os atributos geomorfométricos que você vai usar para fazer essa avaliação do relevo em relação à suscetibilidade e à produção de sedimentos nos reservatórios. Aí tem uma questão que é interessante, João, que eh claro que os reservatórios em si que você vai trabalhar, eles não contemplam. Então aqui em cada uma dessas subbacias do Alto Iguaçu tem um reservatório, eles vão estar aqui, né? Tem a tem vai ter aqui então o Piraquara, aqui é o Passaúa, né? Que vão est aqui na vai est na bacia do Barigui. Então claro que a quantificação vai ser mais importante nessas bacias onde estão os reservatórios. De qualquer forma, a quantificação vai ser feita também para essas outras bacias aqui do entorno que não fazem parte do sistema de abastecimento, mas até para fins de comparação, tá? **[2:15] Speaker A:** Só para poder dizer, olha, então nas áreas que atualmente tem o a aqui as bacias para para abastecimento e comparando com aquelas que não possuem, tá bom? Então a base de dados vai começar, eu vou passar para você aqui essa base que vai ter a hidrografia, as curvas de nível, os pontos cotados e essa delimitação das subbacias. Aí tendo o dado de entrada, que é o modelo digital do terreno de 5 m processado, que eu vou mandar aqui para você, a partir dele é que vão ser gerados os atributos geomorfométricos. Então vou mostrar aqui um por um como que é o processamento deles. Altimetria é o próprio MDT, né, que já tá aqui representado. Claro que não vou falar agora de fins de eh representação gráfica de como vai ser a simbologia deles pros mapas finais. Isso daí vem depois. Tô falando aqui primeiro no processamento e como que eles vão ser quantificados, tá? **[3:12] Speaker A:** Então aqui a altimetria, a declividade para fazer o processamento, raster, análise declividade. A declividade a gente vai trabalhar, dado de entrada o MDT, a gente vai trabalhar ela com, ela pode estar em graus ou porcentagem, né? As unidades a gente vai trabalhar com elas em porcentagem. Então é importante marcar essa opção aqui para expressar a declividade como porcentagem. Marcando isso, executar. Eu vou deixar aqui tudo como temporário, mas aí no para você ter essa base de dados processada, aí então você coloca vai bem organizado tudo numa única pasta, né? Aí você vê qual que vai ser a nomenclatura, tudo certinho para facilitar. **[3:56] Speaker A:** Maravilha. Ó, temos aqui a declividade em porcentagem gerada. Próximo atributo, tô pegando aqui pela lista do teu projeto, vão ser os relevos sombreados. O relevo sombreado ele não vai ser quantificado comparado ali nas subbacias, ele tem uma função mais qualitativa, mais visual para representação, né, visualização das áreas e depois também paraa elaboração dos produtos cartográficos. Então, mas para gerar os relevos sombreados, análise, sombreamento. Aqui é interessante trabalhar. Você pode fazer vários relevos sombreados. Então, dado de entrada vai ser o MDT. Aí você vai, você vai fazer esses vários relevos sombreados. Primeiro, trocando os azimuts de insolação. É interessante fazer um com 315 e outro com 45 de azimuta de insulação. Então, no mínimo dois elevos sombreados. Mas pode testar, fazer outros também, aumentando o exagero vertical. Você vê que ele vai dar aumentar, né, fazer um exagero das sombras. E também ver como que fica na representação, marcando essa opção aqui do sombreamento multidirecional. **[5:04] Speaker A:** Tá? Para você ter alguns relevos sombreados aí que vão te favorecer a visualização da área. Próximo atributo, deixa eu ver aqui na lista, é a área de contribuição. Área de contribuição vai ser um atributo bem importante. Ele vai mostrar então para cada pixel, ele vai fazer a contagem de quanto de área que drena para aquele pixel. Como ele pegando então os valores de elevação de cada pixel, ele vem fazendo uma contagem de qual seria o caminho da água. A partir do momento que ele conta o número de pixels que drenam para aquele ponto, então como a gente tem a resolução do MDT que é de 5 m, a gente sabe que cada pixel tem 25 m² aqui de área de contribuição. Então ele vem fazendo essa contagem, sabendo que cada pixel tem 25. Ele vem fazendo a contagem 25 25. E dessa maneira que ele vai ter um cálculo de área de contribuição para toda todo recorte. **[5:59] Speaker A:** Aí aqui as ferramentas que a gente vai usar são do saga. do Saga Gis. Se o teu KGIS não tiver instalado o Saga, daí eu te mostro como instalar manualmente, tá bom? Aí a ferramenta que a gente vai usar é o flow accumulation. Aí tem uma ferramenta que é aqui, ó, flow accumulation step, que ele faz tudo num único passo. Como dado de entrada, a gente vai colocar o MDT. Aqui a gente pode colocar essa opção do sync removal mesmo, porque é o padrão que ele vai remover o que achar de erros assim dos erros mais e achar pequenos buracos, né? Tanto que ele faz aqui nesse nessa opção e aqui o principal que é qual que vai ser o algoritmo para fazer esse cálculo da área de contribuição. Então a gente pode trabalhar com esse multiple flow direction mesmo, porque ele vai gerar um ele a diferença para esse D8, o D8 ele faz, ele entende assim, cada pixel central ele vai forçar como se todo fluxo de **[6:59] Speaker A:** água fosse para um único pixel. vai ficar um negócio mais quadradão. Eu recomendo usar esse aqui, o múltiplo, porque ele vai considerar, ele faz uma avaliação ponderada. Ele vai dizer assim: "Olha, maior parte do fluxo pode ir para esse pixel, mas pode ser que uma outra parte do fluxo iria pro pixel que vai estar aqui ao lado." Então ele vai gerar uma mais próximo do como é a representação disso na realidade. O que vai interessar pra gente é só o flow acumulation, que esse específico nem precisa processar. Feito isso, só vir aqui em executar, tá? Aí teremos mais um atributo que daí vai ser então a área de contribuição. Próximo é o índice de posição topográfica. Aí pro índice de posição topográfica, assim como para outros atributos, a gente vai usar uma ferramenta que é a do grass chamada R.nbors, que ele vai fazer uma análise de vizinhança a partir, né, de um dado raster. **[7:56] Speaker A:** Então aqui o dado que a gente vai entrar também vai, isso é calculado a partir do MDT. E o que que é o índice de posição topográfica? Ele é o pixel central menos a média da elevação do entorno. Então o que que a gente vai usar essa ferramenta aqui? É para calcular qual que é a média da elevação do entorno que a gente vai pensando em suscetibilidade, a gente tá pensando aqui quais vão ser as áreas que estão mais altas no relevo, que predominaria então aqui a erosão e as áreas mais baixas vai predominar a deposição, seguindo essa lógica, né? E então dado de entrada vai ser o mdt. Como que a gente quer calcular a é a média? Eu vou marcar aqui average, que average é média. tamanho, que tá aqui na versão traduzida, tamanho do bairro, mas é o tamanho da análise de vizinhança que a gente vai fazer aqui. **[8:45] Speaker A:** É importante que esse valor que é atribuído aqui, ele é o número de pixels. Então, veja aqui, ele tá com valor três. Três pixels nesse modelo, sabendo que ele tem resolução espacial de 5 m, equivale a 15 m, que é 3 x 5, né? Então, o que que eu sugiro? O índice de posição topográfica, ele é multalar. Então é interessante trabalhar fazer pelo menos três índice posição topográfica. Dá para fazer um com 100 m, outro com fazer dois a princípio, um com 100 m, um com 1 km e depois a gente avalia fazer aí um com 2,5 k, né, de de raio para fazer essa análise de média de vizingança. Então vou dar um exemplo aqui. Se a gente quer fazer calcular uma disposição topográfica, um IPT, com raio de 100 m, então a gente vai ter que fazer o seguinte, a seguinte continha. 100 m dividido pela resolução do tamanho do pixel, que nesse caso é 5, isso equivale a 20 pixels. **[9:52] Speaker A:** Ou seja, 20 pixels aqui nessa área, nesse recorte, equivale a um raio de 100 m. Então, só que aqui ele só aceita número ímpar. Então, a gente coloca um a mais, vamos colocar 21, ou seja, 21 pixels que vale aqui e então esse raio em metros que a gente quer determinar. Próximo passo e último é marcar aqui, ó, usar vizinhança circular, porque usando vizinhança circular, ele vai fazer o cálculo a partir de um raio. Se não marcar essa opção, ele vai fazer um cálculo como se fosse uma janela, ou seja, ele vai considerar assim um pixel central e uma um quadrado de 21 por 21 pixels. Então é dessa maneira que a gente vai calcular, tá calculando aqui a média da elevação do entorno. Então marquei aqui. Aí é só vir em executar, né? Eu não vou fazer aqui porque ele vai levar um leva um tex aqui. **[10:47] Speaker A:** Ele pode demorar um pouquinho pro processamento quanto maior for o raio, tá? Mas depois que terminar aqui o processamento, o último passo para gerar o índice de posição topográfica é raster, calculadora has aí aqui na calculadora raster, como o IPT ele é o pixel central menos a média do entorno, o pixel central é o quê? É o próprio MDT. a gente vai fazer mdtos. Aí aqui o arquivo que foi gerado do R neighbors, a média da elevação, é o que eu colocaria aqui na expressão. Então, digamos que ele tivesse aqui, né? Vou até colocar aqui, digamos que fosse assim, elevação média sem M, tá? Digamos que esse daqu foi um arquivo que foi gerado, então vou fazer o MDT menos o arquivo gerado do R neighbors. Perfeito. Aí eu posso salvar aqui o arquivo, posso vir aqui colocar na na área de trabalho, né? E vir aqui e colocar, então, chamar ele de IPT 100 M, né? **[11:51] Speaker A:** Porque eu sei que é um índice poção topográfica com raio de 100 m e salvar. e ok, ele vai calcular, ele vai gerar um arquivo com valores positivos, que vão ser as áreas mais altas do relevo, valores negativos, que vão ser as áreas que estão mais baixas em relação ao entorno. Feito isso, a gente pode fazer mais um índice de posição topográfica. Então, venho aqui no R neigors de novo, vou calcular de novo a média da elevação. Average. Digamos que agora eu queira fazer um com raio de 1000 m. Então, 1000 m dividido por 5, que é o tamanho do pixel do modelo que a gente tá usando, dá 200. Ou seja, aqui na ferramenta eu vou colocar 2011, 2011 pixels. Usar vizinhança circular, executar. Depois de executar, has calculadora has mdtado. **[12:45] Speaker A:** Nome de saída. Vou salvar aqui como IPT1 km. Porque daí então esse vai ser um índice poição topográfica mais generalizado com um raio de 1 km. Ele vai mostrar só as principais crisitas e os principais vales, tá? Então só não se assuste. Quanto maior o tamanho do raio que a gente coloca aqui, mais ele demora. Até para esses processamentos aqui do L Neighbors, até sugiro fazer no momento que você não vai usar o computador para nada, tá? Então esse é do em disposição topográfica. Você pode testar vários aí fica, né? Só tô mostrando aqui o passo a passo. Próximo atributo, amplitude altimétrica. O que que é amplitude altimétrica? É a diferença do ponto mais alto do relevo subtraído pelo ponto mais baixo. Mais alto pelo mais baixo. Para isso, a gente também vai usar essa ferramenta aqui, o R neigors. E a amplitude altimétrica, ele também é multiescalar, que a gente pode determinar qual que é o amplitude altimétrica no raio de 1 km. **[13:46] Speaker A:** Então, camada de entrada MDT. Aí aqui que vai mudar, ó, na operação de vizinhança, o range, que essa opção aqui, ó, range, o range ele já faz, ele já subtrai o máximo pelo mínimo. Então, a gente vai selecionar range. Se a gente quiser fazer o qual que é a amplitude altimétrica no raio de então 1 km, é a mesma coisa. Colocar aqui range 2011 pixels, usar vizinhança circular. Marca essa opção. E aqui o resultado disso, pode até colocar como arquivo definitivo. Esse daqui já vai ser a amplitude altimétrica. Então você pode daí salvar ele no salvar aqui o nome e colocar, né? Amplitude 1 km. Aí você que tá mostrando qual que é a variação de relevo no raio de 1 km, tá? Próximo atributo é o índice topográfico de umidade. Esse índice ele leva em consideração. Aí a gente vai usar uma ferramenta do Saga. **[14:49] Speaker A:** Tem acho que tem várias ferramentas que fazem isso, mas eu acho que fica com melhor representação é esse daqui, ó. Saga Wetness Index. que a área de contribuição, aliás, o índice topográfico do unidade, ele faz, ele é uma divisão da área de contribuição dividido pela declividade. Ou seja, o objetivo dele é mostrar o seguinte, é qual que é a tendência de acumulação de água no terreno a partir só considerando aí características topográficas. Eu até tenho um processado aqui que eu vou mostrar só para ver a carinha dele, né? Deixa eu pegar aqui, pegar aqui o TW. Esse daqui é um atributo bem importante, né? Para isso, pegar. Esse daqui é um que é processado para todo o Paraná. Vai até demorar um pouquinho para para carregar, ó. Mas então ele vai dar, quanto maior for o valor dele, significa que maior é a área de contribuição e menor é a declividade. **[15:49] Speaker A:** Porque claro, se a gente tem um um ponto, né, que ele vai ter uma grande área de contribuição e ele é plano, a chance dele acumular água vai ser maior, porque vai ter muita água que drena para aquele ponto e pelo fato de ser plano, a água vai ter uma chance maior de se acumular ali. Então é isso que esse atributo aqui ele ele representa, tá? Por isso que ele há uma relação entre área de contribuição e declividade. E aí quanto maior for o valor dele, significa que a partir de características topográficas, maior vai ser a chance de acumulação, tá? Aí tem mais um atributo que é o teu último desse pacote de dados que é o índice de potência do escoamento. Esse a gente vai fazer manualmente. Daí não vou colocar aqui porque esse daí a gente processa numa reunião pessoalmente para para eu te mostrar como que ele é, mas é bem simples. **[16:40] Speaker A:** Ele ele também leva em consideração a área de contribuição e a declividade, só que ao invés de fazer uma divisão, ele multiplica a área de contribuição e adequidade. Porque daí o que que o índice de potência de escoamento ele faz? A lógica que ele vai mostrar é outra. Ele vai querer dizer assim, ó, quanto maior for a área de contribuição e quanto maior for a declividade, maior é o potencial erosivo daquele rio. Isso também vai ser bem importante pra produção de sedimentos, né? Então pense, a gente tem um rio que ele tá tem um valor muito alto deado de contribuição determinado ponto e ele tá muito inclinado, muito declivoso, putz, ali. Então, quando tiver uma chuva muito forte, o potencial dele de erosão vai ser muito grande. É isso que esse último atributo topográfico que você vai trabalhar mostra, tá? Então isso aqui é encerrando essa parte de geração dos atributos. **[17:32] Speaker A:** Depois de ter os atributos, a gente vai paraa quantificação aqui das bacias, claro, com ênfase nas bacias, que é vai ser o teu olhar, né, que tem aqui os reservatórios, mas também como dá para processar a mesma coisa para todas elas, daí a gente vai para poder comparar, né? Aí o que que a gente vai fazer então, né, na sequência? Qual que é a lógica? É, considerando aqui esses polígonos que são aqui as subbacias, extrair estatísticas. Por exemplo, a gente tem aqui a declividade, né? Vamos colocar a declividade aqui por baixo, só para entender aqui a lógica. Deixa eu desligar a hidrografia, ó. Qual que é a declividade média da bacia do Barigui? Qual que é a declividade média aqui da bacia do do Piraquara? Porque a gente tem um vessa avaliação de relevo pensando no como, né, já que principalmente aqui das áreas a montante dos reservatórios, tá? **[18:32] Speaker A:** Então aqui a gente vai fazer primeiro para toda as bacias e depois ao longo aí da sua pesquisa você pode até duplicar esse arquivo das subbacias e fazer um um recorte nos polígonos vindo assim só até onde começa os reservatórios, sabe? para ver assim qual que é a área, as características aí a montante antes da água chegar nos reservatórios. Então você poderia cortar aqui o polígono seguindo aqui, né, essa bacia só aqui paraas áreas, ó. Opa, que eu tô no valigui aqui, ó. Ou seja, aqui, ó, aqui começa o reservatório do Passaúna. Bem aqui nessa área, né? Então aqui para passar um, né, poderia fazer uma avaliação do polígono somente a partir daqui, daqui para cima, cortando o polígono pra gente ter uma estatística das características de produção, né, de sedimentos que vai ser da cadela hidrográfica nessa porção aqui, a montante. Mas fazendo para uma primeira caracterização para todos os polígonos, a ferramenta que a gente vai usar vai ser, a gente pode usar essa daqui, ó, V. **[19:51] Speaker A:** Hest stats, ele faz o quê? o seguinte, ele faz então uma estatísticas de dados, de dados matriciais, né, que é o que a gente tem aqui. Tem outras ferramentas que fazem isso, mas é, vou mostrar porque a gente vai usar essa. Então aqui, primeira coisa que ele pergunta, qual que é o arquivo de polígono? Nesse primeiro caso, primeira caracterização vai ser as subbacias inteiro, né? Aí aqui que a gente vai escolher qual que vai ser o atributo que a gente vai fazer essas estatísticas. A gente vai fazer para todos esses atributos. Aliás, a gente vai fazer paraa declividade, pra área de contribuição, pro índice de posição topográfica, amplitude altimétrica, o índice topográfico unidade e o índice de potência de escoamento. Talvez até paraa área de contribuição não seja tão necessário assim. Então, mas vou fazer aqui primeiro para declividade para você ver a lógica, né? **[20:41] Speaker A:** Então, aí aqui então selecionei o arquivo de polígono, selecionei o arquivo que eu quero calcular estatísticas e aqui nós temos quais são as estatísticas possíveis de serem calculadas. Então o que a gente vai marcar vai ser mínimo, máximo, range que é a amplitude, average que é a média, deixa desvio padrão, variância, deixa esse daqui, não precisa soma, não precisa primeiro quartil, mediana, terceiro quartil fazer e o percentil não precisa. Então a gente vai deixar marcado pros atributos vão ser esses aqui, tá? Volta, tem que vir aqui, ó, em voltar. Maravilha. Aí ele pode demorar um pouquinho, tá? Esse processamento aqui para fazer. Eu levo alguns minutinhos. Vou aí em executar. Eu já deixei um pronto aqui para mostrar, por isso que daí já tá feito, ó. Daí ele tem aqui, ó, vai gerar um arquivo chamado estatísticas raster. E olha que interessante, quando eu abro a tabela de atributos, ele é o arquivo das bacias, das subbacias. **[21:44] Speaker A:** Tá aqui as bacias, ó. A bacia do rio Palmital, do Passaúa, tá aqui, ó. pega do passarona que a gente estava olhando ali como exemplo. E só que ele acrescentou aqui novas colunas. Nesse caso são aqui as estatísticas da declividade. Então, por exemplo, pro Passaú, o valor mínimo aqui, ó, é de 847. O máximo aqui tá de o desvio padrão é 260, né? Bem, **[22:15] Speaker B:** aí quando você sair com essa porta ali, o Robson tá ali trabalhando, ele pediu para deixar fechar. **[22:20] Speaker A:** Beleza, obrigado. **[22:24] Speaker B:** Então, amplitude 260, né? Aqui a média 937, o desvio padrão 41, variância 1721. Então, estão aqui as informações estatísticas aqui. Eu fiz esse tem uns valores aqui que são da é do MDT, é da altimetria, não é daquividade, tá? Por isso que ele deu esses valores aqui bem diferentes. Então isso daqui é valor de elevação. E mas o que que é interessante? Então tenho aqui isso tabelado já para todas as bacias, todas as bacias aí que fazem ão fazem parte do desse contexto. E deixa eu até separar aqui, ó. **[23:08] Speaker B:** vai ter isso. Só que daí a vantagem a gente pode trabalhar com isso aqui no Excel para daí fazer ali uma tabela síntese. Então posso, a gente pode exportar is aqui no formato de planilha do Excel, clicar num arquivo, exportar, guardar elementos como e aqui eu vou salvar no formato XLSX. Posso ir lá e vou escolher aqui a pasta e salvar aqui. Então, subbacias. Declividade, porque ele vai ter um para cada atributo e salvar. Ele vai exportar isso como uma planilha que daí eu vou abrir aqui. Ah, não exportei. Feito isso, aí a gente pode repetir esse. Ah, daí é interessante salvar também esse arquivo, exportar, né? O arquivo, então, salvar o arquivo espacial, pode salvar em gel package ou shape file com o mesmo nome. Pode ser esse daqui, subbacias declividade. Porque depois se você quiser fazer uma representação, colocar assim, né, quantitativa, ele já vai est aqui pronto, ó. **[24:14] Speaker B:** Quero colocar a média da elevação aqui de cada subbacia. Vou pegar aqui, ó, o campo average, que é média. Opa, graduado. Average, classificar. colocar **[24:33] Speaker B:** sete classes separados aqui por intervalo igual aqui, ó, né? Por como que é possível isso? Porque eu tenho na tabela de atributos essas informações. Então, a gente consegue fazer classificada daí dessa maneira. Então, você vai fazer isso para esses para cada atributo. Ah, terminou para declividade, o próximo vai ser, vamos fazer aqui pro índice, posição topográfica. O arquivo aqui vai ser o arquivo das subbacias e vai selecionar aqui o seria aqui o arquivo processado do índice posição topográfica. Ó, digamos, tem aqui o do índice topográfico de umidade, né? Vou selecionar ali as estatísticas. Esse é mínimo, máximo, variância, aliás, variância vai, esse daqui não vai, soma não vai e percentil não vai. Volta e executar, né? Ou você pode vir aqui e salvar, vai gerar um arquivo vetorial que na tabela de atributos vai ter aí essas estatísticas. Daí você consegue converter isso pro Excel e no Excel consegue fica muito mais fácil para você trabalhar e organizar. **[25:37] Speaker B:** Então vejo que você vai ter muitas tabelas ali assim, né? Tu primeiro tudo que vai facilitar a tua comparação e você vai ver como vai ficar até muito mais fácil para você descrever, falar sobre os resultados disso daqui depois, tá? Então, dessa maneira, encerra aqui a parte de processamento e de geração dos dados gerais aqui dessa primeira fase da parte de análise do relevo geomorfométrica. **[26:02] Speaker B:** Pra segunda análise é a parte da análise multitemporal aí do do uso do solo para ver o que que mudou, né? Para isso, a gente vai usar os dados do MAP biomas. Deixa eu separar aqui, já deixei separadinho. Então, aqui no map biomas, colocar aqui map biom, a gente vai vir aqui em downloads. Aí aqui em downloads tem aqui o mapa das coleções. **[26:29] Speaker B:** Aqui ele tem a opção para fazer o download aqui direto, ó. Link para download direto. E aí você vai fazer o download. Ele tem o link aqui, ó, que é só você trocar o ano. Deixa eu dar um zoom aqui para você ver, tá? Tem aqui, ó, 2024, né, os anos que você vai usar seu projeto. Deixa eu até confirmar aqui. Então você vai fazer assim, vai baixar para vai colocar aqui 2000. Eh, o primeiro vai ser 1985, você colocar aqui 1985, enter, né? Ele vai pegar o arquivo, só que isso daqui é o arquivo para todo o Brasil. Aí depois só vai precisar cortar paraa tua área de estudo. Então vai gerar aqui o arquivo, a cobertura para todo ano 1985. Depois de 85 você vai fazer eh vamos fazer de 20 em 20 anos. Então vai ser 1985, depois 2005 e vai baixar aqui. **[27:28] Speaker B:** E o último, por enquanto, que tem disponível é 2024. Acho que não tem 2025 ainda. Deixa eu confirmar. É, não tem. Então você vai pegar aqui 2024, a gente vai analisar 3 anos. 85, 2005, 2024. Vai baixar aqui então pros 3 anos. Feito o download, ele vai ter aqui para todo o Brasil. Você vai ver aqui has vai recortar só pra área de estudo. Haster extração recortar pela máscara. Camada de entrada vai ser o arquivo do uso do solo do MEP biomas. Camada máscara vai ser esse arquivo das subbacias. E aí o arquivo que você vai colocar aqui na para salvar, né? Salvar como arquivo vai ser o uso do solo para 1985. Você fazer a mesma coisa para 2005 e a mesma coisa para 2024 porque ele vai ficar recortadinho aqui já pra tua área de estudo para esses 3 anos. Aí eu vou colocar, vou deixar disponível um arquivo lá na pasta que a gente usa ali no doc.fp.br, né? **[28:38] Speaker B:** Tem a tua pasta lá específica. Vou deixar um arquivo já com a simbologia do uso do solo, porque isso facilita bastante, né, para você colocar aqui já no colocar carregar automaticamente a legenda do uso do solo. Depois que tiver essa base de dados organizada, a gente vai fazer essa quantificação para cada uma dessas subbacias. O que que mudou no uso do solo? Então a gente vai ter também isso tabelado, tudo tabelado então para 85, 2005, 2024, quais que foram as variações das classes de uso do solo para, né, verificar. Aí a gente vai fazer isso tanto para as bacias quanto para os rios. que a gente vai daí vou passar depois um arquivo de polígono só das áreas de preservação permanente do entorno dos rios para ter essas tabelinhas aqui, tá? Então esse é o segundo objetivo específico do teu trabalho. **[29:34] Speaker B:** Aí o terceiro e último é referente à análise da turbidez ali dos reservatórios a partir do sanamento remoto, né? E aí que que é interessante até dei essa sugestão, né? Ele vai ter que pegar algumas datas específicas para fazer essa análise, né? E é interessante pegar em condições aí de chuva diversas ao longo desse período. Como a gente tem os dados do Sentinel a partir de 2016, então pra análise biturbidez a gente vai trabalhar com os dados a partir o recorte temporal a partir de 2016. Então, a primeira questão que é interessante é avaliar assim, ó, como que se dá a diferenciação de turbidez aí de sedimentos em suspensão nos reservatórios em condições de pouca chuva. Para isso a gente teve aí um período que foi 2020, que ficou foi um ano de seca que os reservatórios ficaram no nível baixíssimo. Então, pô, esse é um caso legal, pegar uma imagem aí de 2020 ali do inverno do auge, quando tava super seco, para fazer uma avaliação. **[30:42] Speaker B:** Agora, como que fica os reservatórios em condição de muita chuva, assim, de um foi um mês muito chuvoso? Eh, opa. Aí, então isso também é um, é interessante pegar uma imagem de satélite após essas chuvas, quando não tiver cobertura de nuvem ali, né, para avaliar qual que é o impacto disso, da pra gente poder comparar como que fica o sedimento em suspensão após um grande período de chuva, como que fica em condições meio normais assim, tem então é essa um pode ir pegando várias imagens desses contextos assim, sabendo a importância da chuva, que que eu deixo como recomendação, deixei aqui já como referência, eu tô aqui no site do Instituto Nacional de Meteorologia, aqui em dados meteorológicos, banco de dados meteorológicos, aqui você vai vir em prosseguir, vai colocar aqui o teu e-mail, prosseguir, porque daí eles mandam os dados pro e-mail, tá? Aí tem a opção aqui de você fazer o download, né? **[31:46] Speaker B:** Separar aqui, ó, com vírgula e pegar dados horários ou dados mensais. Então você quer saber, putz, qual que foi um mês aí bastante chuvoso ou um mês normal, dá para selecionar aqui dados mensais. Ele vai gerar uma planilha do Excel que vai ser mandada por e-mail para você com dados mensais de chuva. Eu fiz uma aqui pegando dados horários de estação automática, região sul. Aqui dá para selecionar quais vão ser as as informações que vão ter na tabela de da planilha. Dá para selecionar tudo, né? Não tem problema. E eu selecionei aqui, ó, Curitiba. Peguei da estação Curitiba, a estação central. Aí vim prosseguir. Ele vai mandar um link por e-mail para vocês. Tem que clicar daí lá no link para poder gerar a confirmação. Depois ele vai processar e vai mandar um novo e-mail com isso daqui, ó. Uma uma planilha do Excel. Eu coloquei ali aqui na data, deixa eu mostrar para você. **[32:47] Speaker B:** Cadê? Ele tem a opção aqui, ó. Data. Data de início, eu comecei aqui 1eo de janeiro de 2016 e data final coloquei a data de hoje, tá aqui, ó, 25/06. Então, não vai trabalhar porque é só a partir de 2016, né, que a gente vai analisar aí. Tá aqui. Então, tá vendo? Deixa eu dar um zoom. Ele tem aqui a data, ó. E como eu selecionei lá diário, ele tem os dados diários, 1o de janeiro, 2 de janeiro, 3 de janeiro. Aí ele tem aqui as variáveis. Então, por exemplo, esse daqui é precipitação total, que é o que interessa pra gente. Ele tem aqui os valores, quanto que choveu em cada um desses dias até hoje, até a data de hoje. Tanto que, ó, tá marcando aqui que pegar lá o último dia, ó. dia 21, 4 dias atrás, choveu 25 mm. Tá aqui, tá marcado, né? **[33:43] Speaker B:** Então isso daqui é um, pode ser um bom parâmetro. Putz, vamos filtrar aqui. Filtrar aqui. Vou pegar, por exemplo, qual que foi a maior chuva? Deixa eu colocar aqui do maior. **[34:01] Speaker B:** É quando tem nu é que deu erro. Não, não tem registro na estação, tá? o erro. Mas, ó, dos dias que tem registro, a maior chuva desde 2016 foi em 22 de fevereiro de 2019. Essa estação registrou 119 mm. Putz, esse daqui é um bem legal para comparar, né? Como que ficou os reservatórios após essa esse acumulado de chuva. Então você pode colocar daí na busca para procurar a imagem de satélite lá do Sentinel, procurar nos dias aqui a partir do dia 22 de de fevereiro de 2019, ver qual imagem de satélite não tem cobertura de novenos reservatórios para poder fazer o cálculo do NDTI, né? Então, ah, vamos pegar aqui algum outra data, né? Que teve uma chuva também significativa, entendeu? Então, só para pegar imagens assim que depois você vai poder para explicar nos materiais e métodos o por que você selecionou algumas dessas imagens de satélite, né, para comparar como que os reservatórios se comportam com distintos acumulados pluviométricos. **[35:09] Speaker B:** Selecionando mais ou menos os períodos que você vai querer fazer a análise da turbidez, a gente vai pro site do Copérnicos Browser. Aí, então aqui no site eu fiz o meu login aqui na busca, né? Vim aqui em search na aba para buscar as imagens para fazer o download. Aí eu selecionei aqui, ó, Sentinel 2, peguei com 10% de cobertura de nuvem. Mas isso daqui pode até controlar, tá? Porque o importante pra gente é que não tenha cobertura de nuvem nesses caras aqui, ó, nos reservatórios. Se tiver um monte de cobertura aqui na mancha, na área urbana de Curitiba, não tem problema. O que não pode ter paraa gente é você vai ter que verificar se não tem nuvem no reservatório. Vai ter que dar uma olhadinha aqui para eles, ó. Ó, no passador não tem nenhuma nuvem. Perfeito. Dá para usar. Então é isso que você vai ter que controlar aqui para baixar as imagens. **[36:03] Speaker B:** Ó, aqui no piringuava sem nuvem. Perfeito. Então, seleciona aí sentinel o nível dois, que ele já vem corrigido com correção geométrica e atmosférica. E aqui no período. Então daí que que você vai colocar, né? Qual que é a data que você vai querer fazer a busca? Então a gente viu aquela chuva lá de fevereiro de 2019, colocar aqui, ó, o ano 2019 a partir do dia, ó, fevereiro, a partir do dia 2, agora não lembro se foi, não, acho que foi dia 22 a chuva, né? Então eu coloco, se for 22, coloco que a partir do dia 22, sei lá, colocar aí um período de as duas semanas, três semanas. Quanto mais perto da data, melhor. Só que é mais difícil porque às vezes pode ser que choveu muito num dia e nos outros dias ficou nublado. Daí tem, né? Então você vai ter que daí fazer essa avaliação. **[36:54] Speaker B:** Colocou aqui a data próxima, sei lá, colocar três semanas, um mês aqui, né? e search quando a gente vem. Vamos até ver aqui fazer ver agosto. Colocar aqui 2019. **[37:17] Speaker B:** Vou colocar até março. Colocar aqui até 23 de março. Buscar. Vamos ver se ele vai achar. Ele vai fazer aqui a busca. **[37:33] Speaker B:** Ó, diz que não, com base nesses parâmetros não foi achado nenhuma imagem. Vamos ver quanto que eu coloquei de cobertura de novo aqui. 21%. Deixar 26. Vou daí. Putz, não tem nenhuma. Tem que aumentar, colocar aqui um período maior. Vou colocar até dia 30 de março. Vamos ver se ele acha, vai ter alguma imagem. **[37:58] Speaker B:** tem uma, ó, com base nesses parâmetros, né? A gente pode vir aqui daí em visualizar. Daí aqui em visualizar a gente consegue até dar ter uma visualização prévia aqui antes para saber se tem cobertura de nuvem ou não nessas áreas nas nossas áreas aqui de de interesse. Ó, por exemplo, aqui já tem uma nuvenzinha aqui no Iraí já não daria para usar essa, né? Então vai ter que ver daí, ó. se essas nuvens tivessem um pouquinho aqui do lado, daria para usar certinho, porque tá vendo que aqui não tem. Então, daí é nesse sentido que vai ter que fazer essa essa busca, né? Pegar algumas imagens aí representativas desse período que não tem eh muita cobertura de novos reservatórios. E o legal é que aqui antes de você baixar, você pode até ter uma visualização do, olha aqui, ó, ndwi que é o índice de água, né? **[38:56] Speaker B:** Ó como que ele fica. Tá vendo que ele tem influência das nuvens aqui? E daí já daria problema pra gente. Então você consegue fazer essa visualização por aqui antes de baixar. Agora digamos que essa imagem aqui não tivesse nuvem. Você viu aqui? Tá tudo certinho em todos os reservatórios. Perfeito. A gente pode vir aqui em insert **[39:16] Speaker B:** até mostra pra gente qual que é o quadrante da imagem. Se tivesse mais de uma imagem aqui, quando a gente clica, ele apareceria aqui a lista de todas essas imagens. Se a gente quiser visualizar em tela, vem aqui em visualiz. Se quiser, ah, viu? Essa imagem aqui, ela é boa, vai dar para usar. Aí para fazer o download, vem aqui, ó, download. Ele vai demorar uns minutinhos carregando e depois vai aparecer a opção para fazer o download. Aí eu já trouxe uma aqui como exemplo de 2020, que 2020 foi o período da seca, né? Então, aí eu já baixei, peguei uma imagem aqui de 2020, já baixei ela, fiz o download, tá aqui, ó. Aí vamos abrir ela aqui no no queis granule. O que a gente vai usar são só quatro. Eh, deixa eu até confirmar, são três bandas só que vão ser usadas. **[40:10] Speaker B:** Então, posso ver aqui imagem R10 para pegar a resolução de 10 m, né? Então, a gente vai usar a banda três, a banda quatro e a banda oito. Só essas aqui são suficientes pra gente trabalhar. Vou adicionar elas, ó. Então, aqui. Aí, primeiro, o que que a gente faz? Primeiro faz o NDW pra gente mapear o as áreas dos reservatórios, né? Então o NDWI, ah, como ele tá na área de trabalho, ele ficou com esse nome grande aqui, né? Mas aí a gente vai fazer aqui no calculadora has o nd, tá? Então abre parênteses que daí aquilo que posso mandar depois a continha para você, né? Mas ele vai ser o NDW. É a banda três, que agora não sei qual delas é aqui, mas vai ser a banda três, que a gente vai selecionar aqui, menos a banda 8, fecha parênteses, dividido por banda 3 mais a banda 8. **[41:27] Speaker B:** Beleza? Vem aqui nos três pontinhos para salvar esse arquivo. Esse daqui vai ser o NDWI. Aí depois de gerar o NDW, a gente vai só ter que ver qual que vai ser o valor de cóte, que ele vai gerar um arquivo com variando de -1 a 1. Normalmente acima valores acima de zero, eles vão mostrar bem as áreas aqui que vai ser o vão ser os reservatórios, vão ser esses corpos d'água. Mas tem que só olhar se valor zero representa isso mesmo. Ou se vai ser o valor 0,1. Às vezes pode ser o 0,2. Aí só tem esse controle. Tendo aí então o NDWI, a gente vai converter ele para polígono, né? Para ter ele como formato de polígono. Então, para ter ele como polígono, a gente pode usar essa ferramenta aqui primeiro, o Rmap Calk Simple, para converter em arquivo binário o que que é corpo d'água e o que que não é corpo d'água. **[42:25] Speaker B:** Então, digamos que aqui tá o arquivo do NDWI, né, gerado. E, digamos que você viu que o valor acima de zero ficou bem representado, colocar a maior do que zero. Ele vai gerar um arquivo binário, 0 e 1. Se você colocou 0.1, então tem que colocar 0 maior do que 0.1. ele vai colocar como um, né, vai representar que você entendeu como corpo d'água acima desse valor aqui do arquivo que for gerado do desse RMP Calk Simple. A gente vai vir aqui em haser conversão, haser para vetor, né, para gerar polígono. Ah, aqui daí, então convertendo esse arquivo, ele vai gerar, vai ter um arquivo que na tabela de atributos desse polígono, ele vai ter ali zero, que não é água, e um, que vão ser os corpos da água. Aí, desse um é que vai excluir tudo que for zero para deixar somente os polígonos da dos corpos d'água. **[43:25] Speaker B:** Então esse aqui é o caminhamento daí para para isso, né? Beleza? Depois que a gente tiver os polígonos aqui do que é a massa dos corpos d'água, a gente vai calcular, então vem aqui calculador rastro e calcular o índice de turbidez para então a gente fazer uma avaliação dos sedimentos em suspensão. Aí para fazer aqui dos sedimentos em suspensão, NDTI, ele é a banda 4 menos a banda 3 dividido pela, claro, vamos selecionar tudo por aqui, né? pela banda quatro mais a banda três. O arquivo a gente vai salvar aqui como pode até deixar ele como temporário que ele vai gerar o NDTI para toda a área, mas a gente quer ver o NDTI, o Durbidd só para as áreas dos polígonos aqui, que daí vem o último passo, que o último passo é a gente vem e aqui em has extraction, extrair e recortar pela máscara. **[44:28] Speaker B:** Como assim recortar pela máscara? Então, a camada de entrada vai ser o arquivo que a gente acabou de gerar do NDTI para toda a área e a camada de máscara vai ser o arquivo de polígono dos corpos d'água gerado no passo anterior. Agora sim, daí quando então a saída dá até para colocar aqui em atribuir um valor sem dado, colocar o valor -999 para ele não ficar com nenhum pixel preto nas bordas, tá? Então, atribuir um valor sem dado -9 - 99909 aqui, tá? Colocar esse valor. E o arquivo de saída vai ser o índice de turbidez para avaliar os sedimentos em suspensão recortados já para as áreas dos que têm os corpos d'água. Esse vai ser é a etapa final, porque daí com base nisso a gente vai ter o eles aqui calculados para datas distintas que a gente vai conseguir correlacionar com o valor de chuva. **[45:28] Speaker B:** Fica até mais fácil para você colocar isso nos teus resultados, né, para fazer uma avaliação deles. A gente vai ter o valor, qual que é a média que deu de turbesse para esse reservatório, qual foi a média de turbesse que deu para esse. Ao longo do dos anos analisados, das imagens analisadas, esse daqui teve uma turbidez maior do que esse? Não sei, tu a pesquisa que vai avaliar isso. Ah, pode ser que você verifique que, sei lá, o reservatório aqui, o no de Piraquara, ele é o mais suscetível quando tem volume de chuva um pouco maior. Ele é o que dá maior turbidez. Não sei, pode ser um resultado que você detecte. Só fazemos processamentos, tendo isso nas tabelas, é que você vai conseguir avaliar, tá? Então isso daqui é a parte operacional do teu trabalho, vai ser isso daqui, né? Pensando aí no em todas as etapas. **[46:21] Speaker B:** Aí qualquer dúvida me avise. Claro que eu tô dei só um panorama geral. Deixa eu parar aqui a apresentação **[46:30] Speaker B:** e aí nós vamos conversando na sequência. Então isso aqui é só para deixar registrado um pouco desse passo a passo.