Estou tentando combinar duas camadas raster (uma é a de dados noturnos de satélite, a outra são os dados de densidade populacional GPW). Estou usando a calculadora raster para combinar os dois (o seguinte é para uma operação de divisão, mas o mesmo problema ocorre mesmo para adição ou subtração), mas obtém faixas de valores, que tendem a ser 0 ou nodata. Anexo uma captura de tela para mostrar como é isso:
Por que isso está acontecendo? Verifiquei se as camadas de varredura são georreferenciadas com o mesmo sistema (WGS84) e até tentei alinhar os rasters, mas sem sucesso. Desculpe se esta é uma pergunta de novato - eu sou relativamente novo no software.
EDIT: http://beta.sedac.ciesin.columbia.edu/data/set/gpw-v4-population-density gpw é a versão 2015 desses dados
Os dados de satélite estão aqui: http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html particularmente a camada de luzes estáveis de 2013
fonte
gdal_info filename.tiffseparadamente em suas duas rasters. Essa é uma ferramenta de linha de comando que fornece metadados detalhados e pode aparecer algo (como resoluções incompatíveis ou CRSes) que são mais difíceis de comparar apenas no QGIS. (Execute-o de uma janela do terminal, ou shell osgeo4w se você estiver usando o Windows)Respostas:
De alguma forma, isso foi difícil por causa dos conjuntos de dados originais que você está usando. Se você verificar as informações para a camada pop_density, parece que está tudo bem.
Para a camada de luzes, algo acontece. Você recebe um aviso e também pode ver que a origem do eixo e a compressão dessa camada são diferentes. De fato, esse preço é de quase 3 GB, enquanto o anterior tem apenas 250 MB.
A coisa mais importante aqui para poder executar uma calculadora raster é que todos os rasters devem estar na mesma projeção e "células" devem se sobrepor exatamente. Se você verificar os detalhes, eles não correspondem à segunda premissa.
Então, como você disse que tentou, tentei alinhar também, mas não funcionou. Pode ser por causa do aviso dos dados. É necessário usar as primeiras transações e fazer alguma compactação nos dados. Isso funcionou para mim. Raster -> Conversão -> Traduzir (converter formato).
Agora, o arquivo tif deve estar em torno de 1 gb. Agora tente alinhar sua varredura usando uma delas como um "modelo". Varredura -> Alinhar varredura ...
Agora você deve ter 2 raster prontos para serem usados no calculador raster. No entanto, dependendo da sua ram, a calculadora raster não gostaria de fazer o trabalho. Se for esse o caso, você pode usar uma ferramenta fácil de grass r.mapcalculator . Como você vê na imagem, é muito fácil de usar e criar seu cálculo.
O resultado da divisão de lights / pop_density deve parecer algo assim.
Não sei se é isso que você estava procurando, mas agora é a sua vez.
Uma última observação: tenha cuidado com os dados sobre as margens da linha costeira. Você pode cobrir essas áreas com interpolar, mas essa é outra questão.
fonte