Declividade com valores muito elevados.

[Lucas Mota]

Olá pessoal, estou usando o qgis em um projeto de um município na Bahia, porém os meus valores de declividade pós processados estão muito elevados. Efetuei o mesmo processamento anteriormente e o valores estavam dentro do esperado, porém em um novo processamento os valores mudaram estranhamente. Gostaria de saber se isso tem relação com alguma configuração, ou eu terei de reinstalar o programa. Estou usando a versão 3.10.6 e já fiz o teste com a versão 2.18.7 encontrando o mesmo problema. Será que pode ser o GDAL? 

[yjmenezes]

Veja as unidades da saida ( graus, radianos, % ).

[Lucas Mota]

Já tentei, tanto em graus, quanto em % os valores ficam fora da realidade.

[yjmenezes]

Que raster esta usando como entrada ?

Veja as propriedades do Raster de entrada. tamanho de pixel, sistema de coordenadas ...

[Lucas Mota]

Já me atentei a estas informações, inclusive baixei a imagem novamente no Earth Explorer. Certamente é alguma configuração que não consigo identificar.

[yjmenezes]

Poste a saida das propriedades das imagens.   Ou  a saida do  gdalinfo   suaimagem.

Quais as classes de declive que pretende usar.

[]s

julio

cartognu.org/covid19

[Lucas Mota]

Segue a sequencia do processamento da imagem. 

1 - O MDE de entrada tem as seguintes propriedades:

Nome s13_w039_1arc_v3

Caminho D:\PastaDeTrabalho\Projetos\ProjetoTeste\s13_w039_1arc_v3.tif

SRC EPSG:4326 - WGS 84 - Geográfico

Extensão -39.0001388888888911,-13.0001388888888911 : -37.9998611111111160,-11.9998611111111124

Unidade graus

Largura 3601

Altura 3601

tipo de dado Int16 - Inteiro de 16 bits com sinal

Descrição do driver GDAL GTiff

Metadados do driver GDAL GeoTIFF

Descrição do registro D:/PastaDeTrabalho/Projetos/ProjetoTeste/s13_w039_1arc_v3.tif

Compressão

Banda 1

STATISTICS_APPROXIMATE=YES

STATISTICS_MAXIMUM=430

STATISTICS_MEAN=100.87001333873

STATISTICS_MINIMUM=-6

STATISTICS_STDDEV=84.470634331073

STATISTICS_VALID_PERCENT=100

Mais informação

AREA_OR_POINT=Point

DTED_CompilationDate=0002

DTED_DataEdition=02

DTED_DigitizingSystem=SRTM

DTED_HorizontalAccuracy=0009

DTED_HorizontalDatum=WGS84

DTED_MaintenanceDate=0000

DTED_MaintenanceDescription=0000

DTED_MatchMergeDate=0000

DTED_MatchMergeVersion=A

DTED_NimaDesignator=DTED2

DTED_OriginLatitude=0130000S

DTED_OriginLongitude=0390000W

DTED_Producer=USCNIMA

DTED_RelHorizontalAccuracy=NA

DTED_RelVerticalAccuracy=0004

DTED_SecurityCode_DSI=U

DTED_SecurityCode_UHL=U

DTED_UniqueRef_DSI=N28 043

DTED_UniqueRef_UHL=N28 043

DTED_VerticalAccuracy_ACC=0006

DTED_VerticalAccuracy_UHL=0006

DTED_VerticalDatum=E96

Dimensões X: 3601 Y: 3601 Bandas: 1

Origem -39.0001,-11.9999

Tamanho do Pixel 0.0002777777777777777775,-0.0002777777777777777775

2 - Fiz o processamento seguindo a imagem abaixo:

3 - O log gerado foi o seguinte:

4- A imagem após o processamento tem as seguintes propriedades, observe que o SRC muda depois do processamento:

Nome Teste

Caminho D:\PastaDeTrabalho\Projetos\ProjetoTeste\Teste.tif

SRC EPSG:4326 - WGS 84 - Geográfico

Extensão -39.0001388888888911,-13.0001388888888911 : -37.9998611111111160,-11.9998611111111124

Unidade graus

Largura 3601

Altura 3601

tipo de dado Float32 - Ponto flutuante de 32 bits

Descrição do driver GDAL GTiff

Metadados do driver GDAL GeoTIFF

Descrição do registro D:/PastaDeTrabalho/Projetos/ProjetoTeste/Teste.tif

Compressão

Banda 1

STATISTICS_APPROXIMATE=YES

STATISTICS_MAXIMUM=8599828

STATISTICS_MEAN=997144.62699958

STATISTICS_MINIMUM=0

STATISTICS_STDDEV=976989.29361026

STATISTICS_VALID_PERCENT=96.67

Mais informação

AREA_OR_POINT=Area

Dimensões X: 3601 Y: 3601 Bandas: 1

Origem -39.0001,-11.9999

Tamanho do Pixel 0.0002777777777777777775,-0.0002777777777777777775

Observe também os valores na camada "Teste", estão muito elevados para declividade em %.

Em sexta-feira, 3 de julho de 2020 17:48:17 UTC-3, Lucas Mota escreveu:

[Márcio Midon]

Não teria de usar uma projeção plana?

--
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 ver essa discussão na Web, acesse https://groups.google.com/d/msgid/qgisbrasil/127b6e9c-fb74-4f87-9362-20b94ae230eao%40googlegroups.com.

--

Cordialmente,

Márcio Midon
Arquiteto e Urbanista UFPel
Especialista Geoprocessamento UFPR
Mestre Geomática UFSM
Tel 53 3305 7075 | Cel 53 98100 2020
SKYPE marcio_midon

[Lucas Mota]

Mesmo com projeção plana o resultado é esse.

[Márcio Midon]

Você reprojetou o raster antes de fazer o input?

Em 04/07/2020 14:19, Lucas Mota escreveu:

Mesmo com projeção plana o resultado é esse.  

[yjmenezes]

@Lucas,

O arquivo de entrada esta com arco de 1 segundo em sistema WGS84  Lon/Lat, e o ponto central seria lon=-38 e lat=-12. Na querida Salvador ?

Para esta latitude, 1 grau corresponde a 110946 metros.  Nao precisa reprojetar o raster !

0.000277777777 * 110946 daria o pixel em metros, coisa de 30 metros +-;

Use este valor como fator de escala, para converter de Graus para metros.

Configure o parametro e troque -s de 1.0 para 110946;

[]s

julio

cartognu.org/covid19    

BA ver grafico GR04

SUL esta feio, RS, SC,PR ver grafico GR07

[Lucas Mota]

O MDE pega toda a extensão da região metropolitana de Ssa, estou trabalhando com um município nessa região, haha. 

Mas então, me passou despercebido a unidade, alterando o parâmetro para 110946 eu volto a ter valores em metros, logo, após reprojetar eu resolvo meu problema. Agradeço a ajuda de vocês, finalmente vou poder dar sequência ao trabalho aqui.

Saudações 

Em sexta-feira, 3 de julho de 2020 17:48:17 UTC-3, Lucas Mota escreveu:

[yjmenezes]

@Lucas,

Vc nao precisa reprojetar, pode continuar com Lon/Lat.

Buffer, etc, *tudo funciona em Lon/Lat.  A dificuldade em usar estas unidades eh a mesma que temos em lidar com polegadas vs milimetros.

Grosso modo, considerando o raio da terra como 6370000 m,  um arco de 1 grau equivale a +- 111 Km.

[]s

julio

cartognu.org/covid19

[Lucas Mota]

Sim sim, entendi a essência. Vou estar a atento a isso e as essas conversões. Obrigado ymenezes.

Em sexta-feira, 3 de julho de 2020 17:48:17 UTC-3, Lucas Mota escreveu: