Problemas com MDTs do aerolevantamento de SC
173 views
Skip to first unread message
Ana Franco
Oct 27, 2017, 2:19:01 PM10/27/17
to qgisbrasil
Boa tarde pessoal,
Tenho verificado que alguns MDTs do aerolevantamento SDS 2012 de Santa Catarina apresentam um intervalo de cotas muito discrepante. Por exemplo, a carta MDT_SH-22-X-B-II-3-SO-B, apresenta como valor mais baixo de cota -32767, e o mais alto 171.644 (esse parece ok). Liguei na Diretoria de Cartografia do Estado e o técnico que me atendeu não é da área, não sobe responder, passou um número de outra secretaria, que não atende. Gostaria de saber se mais alguém está usando esse material e já encontrou esse problema, e como posso solucioná-lo?
Agradeço desde já,
Abraços,
Ana
Tenho verificado que alguns MDTs do aerolevantamento SDS 2012 de Santa Catarina apresentam um intervalo de cotas muito discrepante. Por exemplo, a carta MDT_SH-22-X-B-II-3-SO-B, apresenta como valor mais baixo de cota -32767, e o mais alto 171.644 (esse parece ok). Liguei na Diretoria de Cartografia do Estado e o técnico que me atendeu não é da área, não sobe responder, passou um número de outra secretaria, que não atende. Gostaria de saber se mais alguém está usando esse material e já encontrou esse problema, e como posso solucioná-lo?
Agradeço desde já,
Abraços,
Ana
Jorge Santos
Oct 27, 2017, 2:39:11 PM10/27/17
to qgisb...@googlegroups.com
Esta não parece ser uma questão que envolve o QGIS. Em questões assim, os moderadores orientam escrever [OFF-TÓPIC] na frente do assunto.
Eu já usei esse MDE de Santa Catarina. O que eu usei tinha resolução espacial de 1 metro. Foi um voo da Engemap na época em que eu trabalhei na empresa.
Pode repassar o link para download dessa folha?
Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Analista de Geoprocessamento
----------------------------------------------
--
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+unsubscribe@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.
yjmenezes
Oct 27, 2017, 4:07:51 PM10/27/17
to qgisbrasil
Ana,
provavelmente o -32767 representa locais sem informacao. NoData.
veja com o GDAL.
gdalinfo suacarta.tif
[]s
julio
provavelmente o -32767 representa locais sem informacao. NoData.
veja com o GDAL.
gdalinfo suacarta.tif
[]s
julio
Ana Franco
Oct 27, 2017, 6:15:05 PM10/27/17
to qgisbrasil
Oi Jorge,
Obrigada pelo aviso sobre o OFF-TOPIC.
Te encaminho um link para a pasta do dropbox onde coloquei o arquivo problemático. Para o download desse material é necessário realizar um cadastro no site sigsc.sds.sc.gov.br e solicitar o arquivo zipado 1:25.000. Cada pacote desses vai ter entre 2,5 a 10 GB (dependendo se for MDT ou ortofoto).
Segue o link para o dropbox:
https://www.dropbox.com/sh/86tclb9063qex99/AACJXG73diDVs1miphsT-VHCa?dl=0
Você já conseguiu corrigir esse tipo de problema? Parece um problema de aquisição de dado que não foi corrigido, embora a SDS só tenha disponibilizado os arquivos depois de um tempo de tratamento dos dados...
Abraço e obrigada,
Ana
Obrigada pelo aviso sobre o OFF-TOPIC.
Te encaminho um link para a pasta do dropbox onde coloquei o arquivo problemático. Para o download desse material é necessário realizar um cadastro no site sigsc.sds.sc.gov.br e solicitar o arquivo zipado 1:25.000. Cada pacote desses vai ter entre 2,5 a 10 GB (dependendo se for MDT ou ortofoto).
Segue o link para o dropbox:
https://www.dropbox.com/sh/86tclb9063qex99/AACJXG73diDVs1miphsT-VHCa?dl=0
Você já conseguiu corrigir esse tipo de problema? Parece um problema de aquisição de dado que não foi corrigido, embora a SDS só tenha disponibilizado os arquivos depois de um tempo de tratamento dos dados...
Abraço e obrigada,
Ana
On Friday, October 27, 2017 at 4:39:11 PM UTC-2, Jorge Santos wrote:
Esta não parece ser uma questão que envolve o QGIS. Em questões assim, os moderadores orientam escrever [OFF-TÓPIC] na frente do assunto.Eu já usei esse MDE de Santa Catarina. O que eu usei tinha resolução espacial de 1 metro. Foi um voo da Engemap na época em que eu trabalhei na empresa.Pode repassar o link para download dessa folha?
Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Em 27 de outubro de 2017 16:19, Ana Franco <anavi...@gmail.com> escreveu:
Boa tarde pessoal,
Tenho verificado que alguns MDTs do aerolevantamento SDS 2012 de Santa Catarina apresentam um intervalo de cotas muito discrepante. Por exemplo, a carta MDT_SH-22-X-B-II-3-SO-B, apresenta como valor mais baixo de cota -32767, e o mais alto 171.644 (esse parece ok). Liguei na Diretoria de Cartografia do Estado e o técnico que me atendeu não é da área, não sobe responder, passou um número de outra secretaria, que não atende. Gostaria de saber se mais alguém está usando esse material e já encontrou esse problema, e como posso solucioná-lo?
Agradeço desde já,
Abraços,
Ana
--
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.
Ana Franco
Nov 15, 2017, 9:45:00 AM11/15/17
to qgisbrasil
Bom dia,
Com o GDAL não consegui, mas acabei conseguindo com o GRASS, agora há poucos dias.
Entretanto outros MDTs do mesmo aerolevantamento apresentaram problema para geração das curvas mesmo no GRASS (no ArcGIS também não roda). Descobri que o valores negativos (-32767) estavam concentrados nas bordas (onde não deveria haver dado) e num pequeno conjunto de pixels adjacentes no topo de um morro.
O que fiz para corrigir (gostaria de uma opinião sobre se o método de correção é adequado) foi transformar o tipo do pixel de float (32 bit) para unsigned (32 bit), para que ele excluísse os valores negativos. Aí gerei as curvas com equidistância de 2 metros no GDAL e no ArcGis (só no GRASS que não deu certo).
Gerei um shapefile de pontos aleatórios sobre a cena e extraí os valores de cota do MDT original e corrigido, para verificar a distorção que possa ter havido nos valores de cota durante a transformação.
Em um ponto cuja cota era 85,190918 o valor foi arredondado para 85, ok.
Mas em um ponto onde o valores era 131,96505 o valor foi arredondado para 131.
Praticamente um metro de distorção, que dependendo da declividade da encosta, desloca substancialmente a curva.
Alguma observação ou sugestão pessoal? Gostaria de ter certeza antes de usar os dados...
Abraços,
Ana
Com o GDAL não consegui, mas acabei conseguindo com o GRASS, agora há poucos dias.
Entretanto outros MDTs do mesmo aerolevantamento apresentaram problema para geração das curvas mesmo no GRASS (no ArcGIS também não roda). Descobri que o valores negativos (-32767) estavam concentrados nas bordas (onde não deveria haver dado) e num pequeno conjunto de pixels adjacentes no topo de um morro.
O que fiz para corrigir (gostaria de uma opinião sobre se o método de correção é adequado) foi transformar o tipo do pixel de float (32 bit) para unsigned (32 bit), para que ele excluísse os valores negativos. Aí gerei as curvas com equidistância de 2 metros no GDAL e no ArcGis (só no GRASS que não deu certo).
Gerei um shapefile de pontos aleatórios sobre a cena e extraí os valores de cota do MDT original e corrigido, para verificar a distorção que possa ter havido nos valores de cota durante a transformação.
Em um ponto cuja cota era 85,190918 o valor foi arredondado para 85, ok.
Mas em um ponto onde o valores era 131,96505 o valor foi arredondado para 131.
Praticamente um metro de distorção, que dependendo da declividade da encosta, desloca substancialmente a curva.
Alguma observação ou sugestão pessoal? Gostaria de ter certeza antes de usar os dados...
Abraços,
Ana
Jorge Santos
Nov 15, 2017, 10:33:59 AM11/15/17
to qgisb...@googlegroups.com
Geralmente a computação dos processos em Ponto Flutuante é mais onerosa para o computador, mas parece que o Júlio conseguiu gerar o MDE em Float32. Os valores das bordas não deveriam atrapalhar
Eu vou transformar o seu MDT em Inteiro Não Sinalizado de 16 Bits.
Segue o vídeo. Agora, vamos partir para a avaliação das curvas.
Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Analista de Geoprocessamento
----------------------------------------------
--
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+unsubscribe@googlegroups.com.
Jorge Santos
Nov 15, 2017, 11:11:08 AM11/15/17
to qgisb...@googlegroups.com
Esta é uma região relativamente plana. O intervalo máximo de elevação aqui é 243.
Nos links [1] [2], seguem as curvas de nível geradas no QGIS. No vídeo consta todo o procedimento realizado.
Envie seus pontos aleatórios para que eu possa avaliar a elevação.
Uma suavização para as curvas será bem-vinda. Pesquisarei à respeito.
Uma curiosidade: esses dados são de domínio público? Eu posso ter acesso ao link?
Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Analista de Geoprocessamento
----------------------------------------------
Em 15 de novembro de 2017 12:44, Ana Franco <anavi...@gmail.com> escreveu:
--
Jorge Santos
Nov 15, 2017, 11:16:34 AM11/15/17
to qgisb...@googlegroups.com
Consegui uma boa suavização para as curvas de nível com auxílio do plugin Generalizer.
Este é um exemplo de ferramenta que deveria ser incorporada ao programa. Link do vídeo: https://www.youtube.com/watch?v=j9fny1Fm2-E
Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Analista de Geoprocessamento
----------------------------------------------
Em 15 de novembro de 2017 12:44, Ana Franco <anavi...@gmail.com> escreveu:
--
Ana Franco
Nov 15, 2017, 1:56:02 PM11/15/17
to qgisbrasil
Oi Jorge, então, na verdade estamos falando de arquivos diferentes, pois aquele que o Julio conseguiu gerar as curvas com o GDAL (e eu não, mas com o GRASS eu gerei) é o MDT_SH-22-X-B-II-3-SO-B. Agora estou tentando resolver em outros dois, em que nem no GRASS foi possível gerar (nem ArcGIS). Estes são as cenas MDT_SH-22-X-B-II-3-NE-E e MDT_SH-22-X-B-II-3-NE-F. Pintei os valores negativos de vermelho pra identificar onde estavam localizados, e tirando um pequeno agrupamento dentro da cena, o resto estava nas bordas da imagem, onde deveria ser nodata.
Encaminho link para o dropbox onde disponibilizo as cenas originais, meu arquivo corrigido, as curvas geradas e os pontos de validação: https://www.dropbox.com/sh/640eg3qmexrhwph/AAAbiIpkuwWPmyR0pbR6Qdvka?dl=0
Esses arquivos são públicos sim, basta acessar o site e criar um login em http://sigsc.sds.sc.gov.br/ e baixar dados matriciais e vetoriais (restituição) que cobrem todo o estado de Santa Catarina.
Obrigada novamente,
Abraço,
Ana
Encaminho link para o dropbox onde disponibilizo as cenas originais, meu arquivo corrigido, as curvas geradas e os pontos de validação: https://www.dropbox.com/sh/640eg3qmexrhwph/AAAbiIpkuwWPmyR0pbR6Qdvka?dl=0
Esses arquivos são públicos sim, basta acessar o site e criar um login em http://sigsc.sds.sc.gov.br/ e baixar dados matriciais e vetoriais (restituição) que cobrem todo o estado de Santa Catarina.
Obrigada novamente,
Abraço,
Ana
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.
Ana Franco
Nov 15, 2017, 2:54:59 PM11/15/17
to qgisbrasil
Olá Jorge, a conversão que você indicou funcionou bem com essas cenas também!
Muito obrigada!!!
Abraço,
Ana
Muito obrigada!!!
Abraço,
Ana
Jorge Santos
Nov 15, 2017, 3:04:43 PM11/15/17
to qgisb...@googlegroups.com
Isso que eu iria dizer, Ana. Você não havia visto meus vídeos. Fiz as correções com eles.
Estou no portal baixando os dados originais. É lento pacas.
Jorge Santos
Analista de Geoprocessamento
----------------------------------------------
Analista de Geoprocessamento
----------------------------------------------
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+unsubscribe@googlegroups.com.
Reply all
Reply to author
Forward
0 new messages