Reprojeção ou Salvar como - Raster SRTM

235 views
Skip to first unread message

Leônidas Descovi Filho

unread,
Jan 4, 2017, 5:25:10 PM1/4/17
to qgisbrasil
Camaradas

Qual seria a melhor opção reprojetar raster em WGS 84 - EPSG: 4326 para projeção cônica de Albers EPSG: 102033? O Raster trata-se de um mosaico SRTM que compreende  áreas de partes de 4 fusos UTM: 20, 21, 22 e 23.

Tentei as opções de Reprojetar (RASTER) e Salvar como...(modo de saída: Dados Brutos e formato: Gtiff selecionando o SRC: EPSG: 102033 ). Os resultados foram satisfatórios apenas na segunda.

Haveriam perdas/ganhos entre estas diferentes formas de reprojeção? Será que posso assumir a forma que funcionou? Qual seria em tese a diferença na reprojeção para a seleção do SRC?

Na opção de reprojetar (menu raster) ocorreu o seguinte aviso:

Imagem inline 2Imagem inline 1

Ou seja, ela não funcionou.

Uso o Ubuntu 16.04 Xenial LTS e o Qgis 2.14.10 Essen LTR

Agradeço e aguardo dicas, contribuições e ajudas do camaradas.

Abs


Leônidas


Franklin Januário

unread,
Jan 5, 2017, 10:17:53 AM1/5/17
to qgisbrasil
Ao que parece não há informação dos parâmetros de transformação para o SRC alvo que você selecionou, você pode criar o mesmo manualmente se necessário, no segundo caso ele apenas deve ter "setado" o sistema de coordenadas sem nenhuma transformação.
Você chegou a testar se no segundo caso houve de fato alguma movimentação do raster? Provavelmente não.

Leônidas Descovi Filho

unread,
Jan 5, 2017, 1:43:19 PM1/5/17
to qgisbrasil
Franklin

Nos parametros pode não estar havendo transformação, isso no caso do salvar como. Mas na geometria ocorre o básico com os Sistemas de Referência de Coordenadas (SRC), já que cada um está em um distinto sistema (Geográfico no caso do EPSG: 4326 e Plano no caso do EPSG: 102033. Um é inclinado e outro reto.


Caso tenhas interesse em reproduzir o erro e ou testar, segue link para imagem:

Imagem Original disponível:




Informo aqui os Metadados de cada imagem: a Original EPSG: 4326  e a Salva como EPSG: 102033:

Imagem Original EPSG: 4326:

Driver

GDAL provider

GTiff
GeoTIFF

Descrição do registro

/home/leonidas/Documentos/pos-doc/raster/teste_srtm/srtm_26_18.tif

AREA_OR_POINT=Area

Banda 1

STATISTICS_MAXIMUM=1386

STATISTICS_MEAN=577.27731011603

STATISTICS_MINIMUM=-1

STATISTICS_STDDEV=306.88562939715

Dimensões

X: 6001 Y: 6001 Bandas: 1

Origem

-55.0004,-24.9996

Tamanho do Pixel

0.000833333,-0.000833333

Sem valor de dados

-32768

Tipo de Dados

Int16 - Inteiro de 16 bits com sinal

Visão geral das pirâmides

SRC da Camada

+proj=longlat +datum=WGS84 +no_defs

Extensão da camada (fonte de projeção original da camada)

-55.0004165456048213,-30.0004166908911429 : -49.9995832122714887,-24.9995833575578104

Banda

Banda 1

No da Banda

1

Sem Informações

Nenhuma informação coletada até o momento



Imagem Salva como EPSG: 102033:

Driver

GDAL provider

GTiff
GeoTIFF

Descrição do registro

/home/leonidas/Documentos/pos-doc/raster/teste_srtm/SALVAR_COMO_90srtm_26_18.tif

AREA_OR_POINT=Area

Banda 1

STATISTICS_MAXIMUM=1409

STATISTICS_MEAN=575.68936012085

STATISTICS_MINIMUM=-1

STATISTICS_STDDEV=308.43965175361

Dimensões

X: 5528 Y: 6731 Bandas: 1

Origem

459421,807751

Tamanho do Pixel

90.0034,-90.0016

Sem valor de dados

-32768

Tipo de Dados

Int16 - Inteiro de 16 bits com sinal

Visão geral das pirâmides

SRC da Camada

+proj=aea +lat_1=-5 +lat_2=-42 +lat_0=-32 +lon_0=-60 +x_0=0 +y_0=0 +ellps=aust_SA +units=m +no_defs

Extensão da camada (fonte de projeção original da camada)

459421.3516898183734156,201950.5030881778802723 : 956960.0834587743738666,807751.2620764087187126

Banda

Banda 1

No da Banda

1

Sem Informações

Nenhuma informação coletada até o momento



Desde já grato pela ajuda.
Abs

Leônidas



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

Franklin Januário

unread,
Jan 6, 2017, 4:55:49 AM1/6/17
to qgisbrasil
Bom dia,

Fiz o procedimento por ambos os métodos, e reprojetou normalmente sem nenhuma mensagem de erro:

Versão do QGIS

2.14.6-Essen

Código da versão do QGIS

8e19f72

Compilado sobre Qt

4.8.5

Rodando sobre Qt

4.8.5

Compilado sobre GDAL/OGR

2.1.0

Rodando sobre GDAL/OGR

2.1.1

Compilado sobre GEOS

3.5.0-CAPI-1.9.0

Rodando sobre GEOS

3.5.0-CAPI-1.9.0 r4084

Versão do cliente PostgreSQL

9.2.4

Versão SpatiaLite

4.3.0

Versão do QWT

5.2.3

Versão PROJ.4

492

Versão QScintilla2

2.7.2





Deve ser algo em sua instalação.


Em quarta-feira, 4 de janeiro de 2017 20:25:10 UTC-2, Leônidas Descovi Filho escreveu:

Leônidas Descovi Filho

unread,
Jan 6, 2017, 9:02:48 AM1/6/17
to qgisbrasil
Grato Franklin

Os metadados são de imagens distintas: a bruta EPSG 4326 e a salva como Cônica de Albers EPSG 102033.

Creio que não tenha erro algum.

Abs


Leônidas



Leônidas Descovi Filho

unread,
Jan 6, 2017, 9:12:05 AM1/6/17
to qgisbrasil
Grato Franklin

Irei testar a reprojeção em outra versão do qgis.

Qual SO voce usou no teste?

Abs

PS. se os dois geraram o mesmo resultado, creio poder assumir o que funcionou aqui como correto!?


Leônidas


Franklin Januário

unread,
Jan 6, 2017, 1:20:41 PM1/6/17
to qgisbrasil
Boa tarde,

Windows 10 64 bits.

Att.

Leônidas




Leônidas



Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+...@googlegroups.com.

Leônidas Descovi Filho

unread,
Jan 6, 2017, 1:48:17 PM1/6/17
to qgisbrasil
Franklin

Teu resultado é muito interessante.

A partir dele posso assumir que os dois métodos apresentam o mesmo resultado.

Fabuloso!

Grato!

Abs

Leônidas



Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para qgisbrasil+unsubscribe@googlegroups.com.

Augusto Vieira

unread,
Jan 6, 2017, 1:55:22 PM1/6/17
to qgisb...@googlegroups.com
Leônidas

Confere se seu GDAL esta na versão 2.1.0 pode ser que tenha que atualizá-lo para 2.1.1.

Abraços

Leônidas Descovi Filho

unread,
Jan 6, 2017, 2:27:07 PM1/6/17
to qgisbrasil
Augusto

Como uso ubuntu, a versão que possuo é a mais atual para este sistema, ou seja, verão 1.11.3 de Set de 2015 equivalente a 2.0.1 do win. http://trac.osgeo.org/gdal/wiki/DownloadSource

A versão mais atual é a 1.11.5 para o ubuntu e 2.1.2 para win.

Pelo que percebi, necessito atualizar repositório.

Abs


Leônidas


Augusto Vieira

unread,
Jan 9, 2017, 6:43:15 AM1/9/17
to qgisb...@googlegroups.com
Olá Leonidas,

Desculpe a demora no retorno.... como sou novato no mundo Linux fiquei na dúvida sobre seu comentário, mas também rodo o Qgis 2.18 no Ubuntu 16.04 LTS e no meu Synaptic tenho a versão gdal-bin 2.1.0 instalada. Fiz o teste e também deu o mesmo problema no meu micro. Portanto, ou atualizamos para 2.1.2 ou deve ser outra biblioteca.

Imagem inline 1

Precisa conferir se a linha do comando esta montado certa, porque já aconteceu de dar esse problema para exportação de .tiff para .sqlite e ter que ajustar utilizando a ferramenta do lápis.

Valeu



Leônidas Descovi Filho

unread,
Jan 9, 2017, 8:22:18 AM1/9/17
to qgisbrasil
Augusto

Grato.

E o mais importante, acho que estamos detectando a possível fonte deste problema, e que está recorrendo em duas versões do QGIS, 2.14 e a 2.18.

Abs

Só corrigindo: a vesão que uso do QGIS é 2.14.10 Essen, ou seja, a versão de longa duração.

Leônidas


Reply all
Reply to author
Forward
0 new messages