Mosaico Topodata

[Leônidas Descovi Filho]

Olá Qgiseiros

Não estou conseguindo montar um mosaico raster com imagens topodata.

Alguma dica?

Vejam a mensagem que recebo, em anexo.

Abs

Leônidas

[Jorge Santos]

Primeiro, você precisa Trocar a Projeção do SRTM Topodata.

1 - Coloque todos na mesma pasta e clique no menu Raster - Projeções - Trocar Projeção.

2 - Mande executar o processo na pasta inteira;

3 - Selecione a pasta para gravar os arquivos de saída; 

4 - Especifique o SRC Fonte para 4326 ou 4674;

5 - Especifique o SRC Alvo para 4326 ou 4674;

6 - Não marque a opção Adicioanr dados ao Projeto e execute o Processo.

7 - Crie um novo documento em branco, descarte as mudanças.

8 - Adicione todos os recortes Topodata Fixados;

9 - Faça seu mosaico. Se a quantidade de arquivos for muito grande, o algoritmo dá pau por causa da grande quantidade de parâmetros. Nesse caso, construa um mosaico virtual.

Verifique esses passos e me avise com está esse DEM.

Jorge Santos
Técnico em Geoprocessamento
http://www.processamentodigital.com.br/

Skype: jorgepsantos2002
Twitter: http://twitter.com/jpsantos2002
Facebook: http://www.facebook.com/ProcessamentoDigital

--
--
---
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 essa mensagem porque está inscrito no grupo quot;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 mais opções, acesse https://groups.google.com/d/optout.

[Leônidas Descovi Filho]

Salve Jorge Santos

Grande passo a passo.

Segue o mosaico DEM Topodata na tela do qgis - chugiak

Valeu pela força. Abs

OBS.: a projeção omitida é indicativo a ser considerado, nesse caso o alerta do Jorge. Vlw

Leônidas

[Leônidas Descovi Filho]

Observando com mais calma, percebi que as altitudes mínimas e máximas dos rasters separados (212 a 1284), diverge das altitudes do mosaico (de 272 a 1271).

Qual seria a explicação para essa diferença? método de mosaicagem? algum outro problema.

Outro detalhe, é sabido a existência de altitudes que superam 1800 metros dentro dessa área (mosaico). Haveriam problemas na krigagem realizada pelo INPE?

Alerto também que as imagens SRTM 90 m apresentam valores iguais a piores.

Alguém tem uma luz?

Leônidas

[Sidney]

Leônidas,
A Krigagem é um método de estimativa que usa médias, por isso é inerente à ela a suavização dos resultados.
Isso acontece com todos os métodos de interpolação, só a Triangulação honra os dados amostrais (mínimos e máximos), segundo o Prof. Landim.
Os dados da Missão Topografia por RADAR (Interferometria) do Ônibus Espacial - SRTM foram levantados com a resolução espacial original de 30 metros, e depois foram degradados para 90 metros (reamostragem = suavização) para serem distribuídos.
Daí se concluí que nenhum método, mesmo geoestatístico do mais avançado, mesmo usando vodoo ou magia wicca, pode-se restaurar a resolução original da Imagem.
Você pode consultar o Prof. Jorge Kazuo co-autor do livro Geoestatística: Conceitos e Aplicações se não quiser acreditar em mim.
Abraço.

[Jorge Santos]

A variação entre 1200 e 1800 é muito grande, algo não bate.

Bem estranho para mim. Faça um download do ASTER GDEM da mesma região e tire a prova real.

Quem pode ajudar a clarear isso são as estatísticas.

O mosaico foi criado com sucesso, porém o raster não calculou indicativos estatísticos como valores mínimo, máximo, média, desvio padrão, etc.

A única forma nativa do QGIS que eu conheço para computar as estatísticas é o recorte do raster.

WORKFLOW

1 - Crie um polígono bouding box a partir da extensão do mosaico (Vetor - Investigar - Polígono a partir da extensão da camada).

2 - Após criar o polígono, tente "Recortar" o mosaico com esse bounding box (Raster - Extração - Cortador).

3 - Use o bounding box como camada máscara para recortar o mosaico.

Me diga se os valores min/max foram modificados.

Jorge Santos
Técnico em Geoprocessamento
http://www.processamentodigital.com.br/

Skype: jorgepsantos2002
Twitter: http://twitter.com/jpsantos2002
Facebook: http://www.facebook.com/ProcessamentoDigital

[Leônidas Descovi Filho]

Buenas, testei aqui com um GDEM v2 que eu tenho da área e a altitude máxima dele marca 1243 metros. Confirmei novamente com a carta topográfica de São Joaquim, lá vejo cn e pontos cotados superiores aos 1700 metros. Fora o conhecimento empírico de campo, onde no morro da igreja, temos cotas ultrapassando os 1800 metros acima do nível médio dos mares.

Jorge seu primeiro link para o gdem não está funcionando. É Sidney a tal suavização estatística está matando com os dados reais.

Abs

Leônidas

[Jorge Santos]

É, mas de acordo com o SRTM e o ASTER (produtos diferentes), não há mapeamento de altitude de 1800 metros de altura na região da igreja, e você testou dados de duas fontes diferentes.

Se você tem a carta topográfica local, a melhor opção é ver se existem curvas no formato digital dessa região, ou mesmo gere essas curvas e interpole no QGIS para gerar o DEM.

Se você insistir em utilizar o SRTM/ASTER, considere a altitude de 1200 metros. Esta é a informação correta extraída da fonte de dados que você está utilizando.

Jorge Santos
Técnico em Geoprocessamento
http://www.processamentodigital.com.br/

Skype: jorgepsantos2002
Twitter: http://twitter.com/jpsantos2002
Facebook: http://www.facebook.com/ProcessamentoDigital

[Elias Mendes Costa]

Eu tive o mesmo problema do Leônidas quanto fiz o mosaico do MDE do estado de Rondônia no Qgis. Parece que na hora do mosaico ele gera valores errados. Peguei a mesma base de dados (MDE topodata) e fiz no ArcGis e deu valores bem diferentes (mais próximos do real). Assim concluí que parece ser problema na ferramenta do Qgis. Queria reportar o erro, mas fiquei meio perdido na hora de fazer isso na página do Qgis. Se alguém puder ajudar.