Estoy tratando de combinar dos capas de trama (una son datos satelitales nocturnos, la otra son datos de densidad de población GPW). Estoy usando la calculadora ráster para combinar los dos (lo siguiente es para una operación de división, pero el mismo problema ocurre incluso para la suma o la resta), pero obtengo franjas de valores, que tienden a ser 0 o nodata. Adjunto una captura de pantalla para mostrar cómo se ve esto:
¿Por qué está pasando esto? Me he asegurado de que las capas ráster estén georreferenciadas con el mismo sistema (WGS84) e incluso he intentado alinear los rásteres, pero fue en vano. Lo siento si esta es una pregunta para novatos: soy bastante nuevo en el software.
EDITAR: http://beta.sedac.ciesin.columbia.edu/data/set/gpw-v4-population-density gpw es la versión 2015 de estos datos
los datos satelitales están aquí: http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html particularmente la capa de luces estables de 2013
fuente
gdal_info filename.tiffseparado en sus dos rásteres. Esta es una herramienta de línea de comandos que proporciona metadatos detallados y puede mostrar algo (como resoluciones no coincidentes o CRS) que son más difíciles de comparar solo en QGIS. (Ejecútelo desde una ventana de terminal u osgeo4w shell si está utilizando Windows)Respuestas:
Esto fue de alguna manera difícil debido a los conjuntos de datos originales que está utilizando. Si verifica la información de la capa pop_density, parece que todo está bien.
Para la capa de luces sucede algo. Obtiene una advertencia y también puede ver que el origen del eje y la compresión de esta capa son diferentes. De hecho, este tiff está cerca de 3 gb, mientras que el anterior tiene solo 250 mb.
Lo más importante aquí para poder realizar una calculadora ráster es que todos los rásteres deben estar en la misma proyección y las "celdas" deben superponerse exactamente. Si verifica los detalles, estos no coinciden con la segunda premisa.
Entonces, como dijiste que lo habías intentado, intenté alinear también pero no funcionó. Puede ser debido a la advertencia de los datos. Tienes que usar traduct firts y hacer algo de compresión en los datos. Eso funcionó para mí. Ráster -> Conversión -> Traducir (convertir formato).
Ahora el archivo tif debería estar alrededor de 1 gb. Ahora intente alinear su ráster utilizando uno de esos como "plantilla". Ráster -> Alinear ráster ...
Ahora debe tener 2 ráster listos para usar en la calculadora ráster. Sin embargo, dependiendo de su carnero, la calculadora ráster no querría hacer el trabajo. Si este es el caso, puede usar una herramienta fácil de grass r.mapcalculator . Como ves en la imagen es muy fácil de usar y de crear tu cálculo.
El resultado de dividir las luces / pop_density debería parecerse a esto.
No sé si esto es lo que estabas buscando, pero ahora es tu turno.
Una última nota: tenga cuidado con los datos sobre los bordes de la línea costera. Puede cubrir estas áreas con interpolar, pero esta es otra pregunta.
fuente