J'essaie de combiner deux couches raster (l'une est des données satellites nocturnes, l'autre est les données de densité de population GPW). J'utilise une calculatrice raster pour combiner les deux (ce qui suit est pour une opération de division, mais le même problème se produit même pour l'addition ou la soustraction) mais j'obtiens des bandes de valeurs, qui ont tendance à être 0 ou nodata. Je joins une capture d'écran pour montrer à quoi cela ressemble:
Pourquoi cela arrive-t-il? J'ai fait en sorte que les couches raster soient géoréférencées avec le même système (WGS84) et j'ai même essayé d'aligner les rasters, mais en vain. Désolé si c'est une question de débutant - je suis relativement nouveau dans le logiciel.
EDIT: http://beta.sedac.ciesin.columbia.edu/data/set/gpw-v4-population-density gpw est la version 2015 de ces données
les données satellites sont ici: http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html en particulier la couche de lumières stables de 2013
la source
gdal_info filename.tiffséparément sur vos deux rasters. Il s'agit d'un outil en ligne de commande qui fournit des métadonnées détaillées et peut afficher quelque chose (comme des résolutions ou des CRS incompatibles) qui sont plus difficiles à comparer dans QGIS seul. (Exécutez-le à partir d'une fenêtre de terminal ou d'un shell osgeo4w si vous utilisez Windows)Réponses:
Cela a été difficile à cause des jeux de données originaux que vous utilisez. Si vous vérifiez les informations pour la couche pop_density, il semble que tout va bien.
Pour la couche de lumière, quelque chose se passe. Vous obtenez un avertissement et vous pouvez également voir que l'origine de l'axe et la compression de ce calque sont différentes. En fait, ce tiff est proche de 3 Go alors que le précédent ne fait que 250 Mo.
La chose la plus importante ici pour pouvoir effectuer une calculatrice raster est que tous les rasters doivent être dans la même projection et les "cellules" doivent se chevaucher exactement. Si vous vérifiez les détails, ceux-ci ne correspondent pas à la deuxième prémisse.
Donc, comme vous avez dit que vous aviez essayé, j'ai également essayé d'aligner mais cela n'a pas fonctionné. Peut-être à cause de l'avertissement des données. Vous devez utiliser les premières traduct et effectuer une compression sur les données. Cela a fonctionné pour moi. Raster -> Conversion -> Traduire (format de conversion).
Maintenant, le fichier tif devrait être d'environ 1 Go. Essayez maintenant d'aligner votre raster en utilisant l'un d'eux comme "modèle". Raster -> Aligner le raster ...
Vous devez maintenant avoir 2 rasters prêts à être utilisés dans la calculatrice raster. Cependant, selon votre RAM, la calculatrice raster ne voudrait pas faire le travail. Si tel est le cas, vous pouvez utiliser un outil simple de grass r.mapcalculator . Comme vous le voyez dans l'image est très facile à utiliser et à créer votre calcul.
Le résultat de la division des lumières / pop_density devrait ressembler à ceci.
Je ne sais pas si c'est ce que vous cherchiez, mais maintenant, c'est votre tour.
Une dernière note: soyez prudent avec les données sur les bords de la ligne côtière. Vous pouvez couvrir ces zones avec interpolation mais c'est une autre question.
la source