Suite

Somme des valeurs dans un raster (25 m) à l'aide d'un autre raster (1 km) ArcGIS Desktop 10.2


Je résume certaines données de population à l'aide d'une superposition raster de 1 km. J'ai converti les données de population, couche de polygones, en un raster de 25 mètres qui contient la valeur de population/m^2. Cette valeur est ensuite multipliée par 625 m^2 (la surface de chaque cellule), puis j'ai pensé que je pouvais faire un résumé zonal en utilisant ma couche raster de 1 km (qui est bien accrochée à la couche raster de 25 mètres), mais maintenant je suis coincé. Si je pouvais simplement additionner les valeurs de chaque cellule de 25 mètres qui se produit dans une cellule de 1 km, il me suffirait de multiplier par 625 pour obtenir la population totale approximative dans cette cellule… Je pourrais utiliser un rééchantillonnage majoritaire, mais cela devrait donner une meilleure estimation.

J'ai reçu d'excellents conseils sur le codage d'une solution, mais j'aimerais le comprendre en suivant cette ligne de pensée si possible.


Si vous voulez rester complètement en raster, ce serait bien si Resample avait une méthode Sum, mais tant que Spatial Analyst est disponible, Zonal Statistics devrait pouvoir résoudre votre problème.

Si vous utilisez un autre raster pour vos zones, chaque zone (qui sera chaque cellule dans ce cas) a besoin d'une valeur unique. Toutes les valeurs de cellule de 1 km sont-elles les mêmes ? Si c'est le cas, vous devrez modifier le raster et/ou ajouter une table d'attributs raster. Alternativement, plutôt que d'utiliser un raster pour l'entrée de zone, vous pouvez simplement utiliser votre filet de pêche (en supposant qu'il s'agisse de polygones). Voir le lien d'aide Statistiques zonales ci-dessus - il y a plusieurs considérations pour le raster ou le vecteur utilisé pour la définition de zone (taille de cellule, accrochage, nombre de canaux, présence d'une table d'attributs raster, etc.).

La sortie résultante doit être un raster de cellule de 1 km avec chaque cellule la somme de vos valeurs de cellule raster de 25 m. Vous pouvez ensuite exécuter ce raster en sortie via Raster Calculator pour votre opération [625 x chaque cellule] pour un raster récapitulatif final.

Si vous souhaitez travailler dans des tableaux, il existe également l'outil Statistiques zonales en tant que tableau qui crée un tableau plutôt qu'une sortie raster. Notez que la création d'une table d'attributs raster et celle-ci résument déjà des cellules identiques en tant qu'enregistrements uniques dans la table avec un nombre de fréquences, et non des enregistrements de cellules individuelles. Si vous voulez une table avec chaque cellule en tant qu'enregistrement (en supposant qu'elles ne soient pas déjà toutes des valeurs uniques), vous pouvez vous référer à Comment obtenir les valeurs de chaque cellule dans la table attributaire raster?


Vous avez déjà vos polygones de superposition, transformez maintenant votre raster de données en polygones en utilisant raster en polygone, sélectionnez « pas de simplification » et vos polygones correspondront exactement à vos cellules. Il existe des limites à la taille d'un raster pouvant être traité car le fichier de formes intermédiaire doit être inférieur à 2 Go pour tous les fichiers ; En raison de ces limitations de fichiers, vous pouvez effectuer la poligonisation par lots en utilisant l'extraction par masque et en sélectionnant des moitiés ou des quarts de vos polygones superposés comme masque. Cette approche accélérera également les choses car chaque lot peut être effectué simultanément avec plusieurs instances d'ArcCatalog.

Avant d'effectuer l'union, il serait préférable d'utiliser Intégrer pour minimiser/atténuer les légères disparités entre les couches de polygones, ce qui entraînerait des polygones en ruban. Assurez-vous de lire les mises en garde sur cet outil et de conserver une sauvegarde ; avec des paramètres trop enthousiastes, il peut détruire de bonnes données.

Union chaque partie à son tour avec les caractéristiques de superposition ; Union est limité à deux classes de fonctionnalités avec des niveaux de licence de base ou standard, ce qui devrait convenir à cette opération. Après l'union, toutes les valeurs des données superposées et du raster seront présentes, pour résumer cette utilisation des statistiques récapitulatives avec un champ de cas des attributs importants de la superposition et de la somme des attributs raster.

Remarque : ce processus fonctionnera tout aussi bien si vous l'exécutez sur des points générés à partir du raster, sauf si vous avez besoin d'une zone de chevauchement de polygones, auquel cas l'intégration serait plus préjudiciable que bénéfique.


Voir la vidéo: How to calculate area in ArcGIS from attribute table (Octobre 2021).