Suite

Convertir le raster en polygone en fonction des plages de valeurs des cellules à l'aide d'ArcGIS for Desktop ?


J'ai un raster à valeur en virgule flottante que j'ai besoin de convertir en polygone comme suit

Placez les cellules avec VALUE < 0,5 dans (un petit/le plus petit nombre de) polygones

Cellules de la plage [0,5, 1] ​​en polygones

etc…

J'ai regardé l'outil Raster vers polygone, mais il n'y a pas de place pour spécifier des valeurs de plage. Existe-t-il un outil dans ArcGIS Desktop/Toolbox pour faire cela ?


j'utiliserais soit le Con ou alors Reclasser outils pour créer un nouveau raster avec la plage de valeurs qui vous intéresse, puis effectuez la conversion raster en vecteur sur ce raster nouvellement créé. Deux étapes et le tour est joué !


Le raster sera converti en une classe d'entités surfaciques à l'aide de l'outil ArcGIS Raster vers polygone, puis en entités linéaires à l'aide de l'outil Entité vers ligne. L'outil Raster vers polygone ne peut convertir que des rasters entiers en polygones. Si le raster en entrée est un raster à virgule flottante, vous devez utiliser le paramètre Map Algebra Expression pour le convertir en raster d'entiers.

Classe d'entités ligne en sortie.

Les répertoires manquants dans ce chemin seront créés s'ils n'existent pas.

Détermine si les polygones en sortie seront lissés en formes plus simples ou conformes aux bords des cellules du raster en entrée.

Vrai - Les polygones seront lissés en des formes plus simples. C'est la valeur par défaut.

False - Les polygones se conformeront aux bords des cellules du raster en entrée.

Champ utilisé pour affecter les valeurs des cellules du raster en entrée aux lignes du jeu de données en sortie. Il peut s'agir d'un nombre entier ou d'un champ de chaîne.

Nouveau système de coordonnées vers lequel projeter le raster.

Le raster ne peut être projeté dans un nouveau système de coordonnées que si la projection d'origine est définie. Une erreur sera générée si vous spécifiez un nouveau système de coordonnées sans définir le système de coordonnées d'origine.

L'outil ArcGIS Project Raster est utilisé pour effectuer la projection. La documentation de cet outil recommande de spécifier également une taille de cellule pour le nouveau système de coordonnées.

J'ai remarqué que pour certains systèmes de coordonnées, l'outil ArcGIS 9.2 Project Raster semble couper le raster projeté dans une mesure arbitraire trop petite. Par exemple, lors de la projection d'une image globale de chlorophylle MODIS Aqua de 4 km en coordonnées géographiques vers Lambert_Azimuthal_Equal_Area avec un méridien central de -60 et une latitude d'origine de -63, l'image résultante est tronquée pour ne montrer qu'un quart de la planète. Ce problème ne se produit pas lorsque Project Raster est appelé de manière interactive à partir de l'interface utilisateur d'ArcGIS, il se produit uniquement lorsque l'outil est appelé par programmation (la méthode ProjectRaster_management du géoprocesseur). Ainsi, vous ne le verrez peut-être pas lorsque vous utilisez Project Raster vous-même, mais cela peut arriver lorsque vous utilisez des outils MGET qui appellent Project Raster dans le cadre de leurs opérations de géotraitement.

Si vous rencontrez ce problème, vous pouvez le contourner comme ceci :

Tout d'abord, exécutez cet outil sans spécifier de nouveau système de coordonnées, pour obtenir le raster dans le système de coordonnées d'origine.

Dans ArcCatalog, utilisez l'outil Projeter un raster pour projeter le raster dans le nouveau système de coordonnées. Vérifiez que tout le raster est présent, qu'il n'a pas été découpé dans une mesure trop petite.

Dans ArcCatalog, recherchez l'étendue du raster projeté en cliquant dessus avec le bouton droit dans l'arborescence du catalogue, en sélectionnant Propriétés et en faisant défiler jusqu'à Etendue.

Maintenant, avant d'exécuter l'outil MGET qui projette le raster, définissez le paramètre d'environnement Etendue sur les valeurs que vous avez recherchées. Si vous appelez l'outil MGET de manière interactive à partir d'ArcCatalog ou d'ArcMap, cliquez sur le bouton Environnements dans la boîte de dialogue de l'outil, ouvrez Paramètres généraux, modifiez la liste déroulante Étendue sur "Comme spécifié ci-dessous" et saisissez les valeurs que vous avez recherchées. Si vous l'invoquez à partir d'un modèle de géotraitement, cliquez avec le bouton droit sur l'outil dans le modèle, sélectionnez Créer une variable, À partir de l'environnement, Paramètres généraux, Étendue. Cela placera l'étendue en tant que variable dans votre modèle, attachée à l'outil MGET. Ouvrez la variable Etendue, remplacez-la par "Comme spécifié ci-dessous" et saisissez les valeurs que vous avez recherchées. Si vous appelez l'outil MGET par programme, vous devez définir la propriété Extent du géoprocesseur sur les valeurs que vous avez recherchées. Veuillez consulter la documentation ArcGIS pour plus d'informations à ce sujet et les paramètres d'environnement en général.

Exécutez l'outil MGET. L'étendue du raster doit maintenant être de la bonne taille.

Une méthode de transformation utilisée pour convertir entre le système de coordonnées d'origine et le nouveau système de coordonnées.

Ce paramètre est une nouvelle option introduite par ArcGIS 9.2. Vous devez disposer d'ArcGIS 9.2 pour utiliser ce paramètre.

Ce paramètre n'est nécessaire que lorsque vous spécifiez que le raster doit être projeté dans un nouveau système de coordonnées et que ce nouveau système utilise un système de référence différent du système de coordonnées d'origine, ou qu'il existe une autre différence entre les deux systèmes de coordonnées qui nécessite une transformation. Pour déterminer si une transformation est nécessaire, je recommande la procédure suivante :

Tout d'abord, exécutez cet outil sans spécifier de nouveau système de coordonnées, pour obtenir le raster dans le système de coordonnées d'origine.

Ensuite, utilisez l'outil ArcGIS 9.2 Project Raster sur le raster pour le projeter dans le système de coordonnées souhaité. Si une transformation géographique est nécessaire, cet outil vous en demandera une. Notez le nom exact de la transformation que vous avez utilisée.

Enfin, si une transformation était nécessaire, saisissez le nom exact dans cet outil, réexécutez-le et vérifiez que le raster a été projeté comme vous le souhaitez.

Algorithme de rééchantillonnage à utiliser pour projeter le raster d'origine dans un nouveau système de coordonnées. L'outil ArcGIS Project Raster est utilisé pour effectuer la projection et accepte les valeurs suivantes :

NEAREST - interpolation du voisin le plus proche

BILINAIRE - interpolation bilinéaire

Vous devez spécifier l'un de ces algorithmes pour projeter vers un nouveau système de coordonnées. Une erreur sera générée si vous spécifiez un nouveau système de coordonnées sans sélectionner d'algorithme.

La taille de cellule du système de coordonnées projeté. Bien que ce paramètre soit facultatif, pour obtenir les meilleurs résultats, la documentation ArcGIS vous recommande de toujours le spécifier lors de la projection dans un nouveau système de coordonnées.

Les coordonnées x et y (dans l'espace de sortie) utilisées pour l'alignement des pixels.

Ce paramètre est une nouvelle option introduite par ArcGIS 9.2. Vous devez disposer d'ArcGIS 9.2 pour utiliser ce paramètre. Il est ignoré si vous ne spécifiez pas que le raster doit être projeté dans un nouveau système de coordonnées.

Classe d'entités, raster ou autre jeu de données géographiques existant ayant l'étendue à laquelle le raster doit être découpé.

AVERTISSEMENT : si vous utilisez cet outil dans un modèle de géotraitement ArcGIS et que vous sélectionnez un jeu de données en cliquant sur l'icône de dossier et en accédant au jeu de données, votre sélection peut mystérieusement disparaître de cette zone de texte après la fermeture de l'outil. Il s'agit d'un bogue dans ArcGIS. Pour contourner ce problème, faites glisser et déposez l'ensemble de données souhaité dans le modèle. Cela créera une couche dans le modèle pour cet ensemble de données. Sélectionnez ensuite cette couche dans cet outil en cliquant sur la liste déroulante plutôt que de cliquer sur l'icône du dossier. Le calque sélectionné ne doit pas disparaître lorsque vous fermez l'outil.

Rectangle auquel le raster doit être découpé.

Si un système de coordonnées projetées a été spécifié, le découpage est effectué après la projection et les coordonnées du rectangle doivent être spécifiées dans le nouveau système de coordonnées. Si aucun système de coordonnées projetées n'a été spécifié, les coordonnées doivent être spécifiées dans le système de coordonnées d'origine.

L'outil ArcGIS Clip est utilisé pour effectuer le clip. Le rectangle de découpage doit être transmis à cet outil sous la forme d'une chaîne de quatre nombres séparés par des espaces. L'interface utilisateur d'ArcGIS formate automatiquement la chaîne correctement lors de l'appel de cet outil à partir de l'interface utilisateur d'ArcGIS, vous n'avez pas à vous soucier du format. Mais lorsque vous l'invoquez par programmation, veillez à fournir une chaîne correctement formatée. Les numéros sont ordonnés GAUCHE, BAS, DROITE, HAUT. Par exemple, si le raster se trouve dans un système de coordonnées géographiques, il peut être découpé à 10 W, 15 S, 20 E et 25 N avec la chaîne :

Des nombres entiers ou décimaux peuvent être fournis.

Mapper l'expression d'algèbre à exécuter sur le raster.

ATTENTION: Le générateur de modèle de géotraitement ArcGIS peut supprimer de manière aléatoire et silencieuse la valeur de ce paramètre. Il s'agit d'un bogue dans ArcGIS. Avant d'exécuter un modèle que vous avez enregistré, ouvrez cet outil et vérifiez que la valeur du paramètre existe toujours.

L'expression est exécutée une fois que le raster converti est projeté et découpé (si ces options sont spécifiées). Utilisez la chaîne inputRaster sensible à la casse pour représenter le raster sur lequel vous souhaitez maintenant effectuer une algèbre cartographique. Par exemple, pour convertir le raster en raster d'entiers et ajouter 1 à toutes les cellules, utilisez cette expression :

La chaîne inputRaster est sensible à la casse. Avant d'exécuter l'expression d'algèbre de carte, la chaîne est remplacée par le chemin d'accès à un raster temporaire qui représente le raster généré. L'expression finale doit comporter moins de 4 000 caractères, sinon ArcGIS signalera une erreur.

L'outil ArcGIS Single Output Map Algebra est utilisé pour exécuter l'expression d'algèbre cartographique. Vous devez disposer d'une licence pour l'extension ArcGIS Spatial Analyst afin d'effectuer l'algèbre cartographique.

La syntaxe de l'algèbre cartographique peut être très pointilleuse. Voici quelques conseils qui vous aideront à réussir avec cet outil :

Avant d'utiliser cet outil, construisez et testez votre expression d'algèbre cartographique à l'aide de l'outil ArcGIS Single Output Map Algebra. Collez ensuite l'expression dans cet outil et modifiez-la pour utiliser la variable inputRaster plutôt que la valeur de test que vous avez utilisée avec Single Output Map Algebra.

Si vous développez votre expression directement dans cet outil, commencez par une expression très simple. Vérifiez qu'il fonctionne correctement, ajoutez-en un peu et vérifiez à nouveau. Répétez ce processus jusqu'à ce que vous ayez construit l'expression complète.

Séparez toujours les opérateurs mathématiques des chemins raster à l'aide d'espaces. Dans l'exemple ci-dessus, l'opérateur / contient un espace de chaque côté. Suivez ce modèle. Dans certaines circonstances, ArcGIS ne parviendra pas à traiter les expressions d'algèbre raster qui ne séparent pas les chemins raster des opérateurs utilisant des espaces. Le message d'erreur signalé n'indique généralement pas qu'il s'agit du problème, et le retrouver peut être très frustrant.


Remarques

Spécifiez les paramètres d'analyse raster pour chaque tâche du portail d'entreprise dans l'onglet Paramètres de la page Génération de la sortie.

Les paramètres ne sont pas conservés entre les tâches.

Lorsque vous spécifiez l'intersection des entrées comme étendue de traitement, il est possible qu'aucune des entités ou des cellules ne se chevauche et qu'une étendue nulle (largeur et hauteur nulles) puisse en résulter. Dans de tels cas, aucune caractéristique ou cellule ne sera traitée.

Le paramètre Etendue définit les rasters qui seront traités par un outil. Il est utile lorsque vous devez traiter uniquement une partie d'un ensemble de données plus important. Vous pouvez considérer ce paramètre comme un rectangle utilisé pour sélectionner les entités en entrée et les rasters à traiter. L'étendue du jeu de données en sortie sera généralement plus grande que le paramètre Etendue pour tenir compte des entités ou des cellules qui traversent le rectangle d'étendue.

Le paramètre Aligner le raster est utilisé pour aligner ou aligner une étendue pendant l'exécution. Le coin inférieur gauche de l'étendue est accroché à un coin de cellule du raster d'accrochage et le coin supérieur droit est ajusté à l'aide de la taille de cellule en sortie. Par conséquent, lorsque la taille de cellule en sortie est la même que la taille de cellule du raster d'accrochage, les cellules du raster en sortie sont alignées avec les cellules du raster d'accrochage.


Volet Exporter raster

Le volet Exporter le raster vous permet d'exporter l'intégralité du jeu de données raster, de la mosaïque, du service d'imagerie ou la partie de l'affichage.

Fournit des options pour spécifier le jeu de données raster en sortie , le système de coordonnées , les transformations géographiques , la géométrie de découpage , le maintien de l'étendue de découpage , la taille de cellule , la taille de raster , le type de pixel , la valeur de pixel d'échelle , la valeur NoData , les paramètres de rendu , le format de sortie , le type de compression et la qualité de compression .

Vous permet de configurer le raster de capture , la taille de la tuile , le rééchantillonnage , le type de source , les paramètres des pyramides et les paramètres des statistiques.

Vous pouvez également ouvrir ce volet en cliquant sur le bouton Exporter le raster dans l'onglet Données.

Lors du stockage d'un jeu de données raster dans une géodatabase, aucune extension de fichier ne doit être ajoutée au nom du jeu de données raster. Lors du stockage du jeu de données raster dans un format de fichier, vous devez spécifier l'extension appropriée au fichier.

La boîte de dialogue Référence spatiale est contextuelle et répertorie un système de coordonnées xy ou z selon que vous sélectionnez l'option Système de coordonnées XY actuel ou Z actuel.

  1. Sélectionnez la transformation géographique appropriée lorsque vos données sont transformées entre différents systèmes de coordonnées. L'application n'utilisera que les transformations appropriées à la projection, toutes les autres seront ignorées.

Cette option exportera le jeu de données raster à l'aide des spécifications de référencement spatial du jeu de données raster.

L'étendue de l'affichage actuel sera utilisée.

Par exemple, si vous effectuez un zoom avant sur votre zone d'étude particulière, vous pouvez utiliser cette option pour traiter les entités qui se trouvent dans l'étendue d'affichage actuelle.

Vous entrez les coordonnées du rectangle de délimitation minimum : saisissez l'étendue pour Gauche , Droite , Haut et Bas .

Toutes les couches sont répertoriées et vous pouvez en choisir une à utiliser comme étendue.

Comme avec l'option Étendue d'affichage actuelle, l'étendue de la couche est lue et stockée.

Utilisez le bouton Parcourir pour rechercher l'emplacement du dossier de la classe d'entités que vous souhaitez utiliser pour la géométrie de découpage . Une fois qu'une entité en entrée est fournie, une case à cocher apparaît avec l'option Utiliser les entités en entrée pour la géométrie de découpage, avec des options de découpage pour délimiter l'intérieur ou l'extérieur .

Conserve l'alignement des cellules en tant que raster en entrée et ajuste l'étendue en sortie en conséquence.

Ajuste le nombre de colonnes et de lignes et rééchantillonne les pixels pour qu'ils correspondent exactement à l'étendue de découpage spécifiée.

L'option Mettre à l'échelle la valeur de pixel apparaît lorsqu'un type de pixel différent est choisi, qui peut être utilisé pour mettre à l'échelle votre type de pixel d'une profondeur de bits à une autre.

Si le type de pixel est rétrogradé (abaissé), les valeurs raster en dehors de la plage valide pour cette profondeur de pixel seront tronquées et perdues. Par exemple, lorsque la sortie est un type de pixel différent de l'entrée (comme 16 bits à 8 bits), vous pouvez choisir de mettre les valeurs à l'échelle pour s'adapter à la nouvelle plage sinon, les valeurs qui ne rentrent pas dans le nouveau pixel la plage sera supprimée.

En cas de mise à l'échelle, telle que 8 bits à 16 bits, le minimum et le maximum des valeurs 8 bits seront mis à l'échelle au minimum et au maximum dans la plage 16 bits. En cas de réduction, telle que 16 bits à 8 bits, le minimum et le maximum des valeurs 16 bits seront réduits au minimum et au maximum dans la plage 8 bits. Pour en savoir plus sur la capacité de profondeur de bits pour les formats d'exportation pris en charge, voir Liste des capteurs pris en charge.

Les valeurs de pixel resteront les mêmes et ne seront pas mises à l'échelle. Toutes les valeurs qui ne rentrent pas dans la plage de valeurs seront rejetées. C'est la valeur par défaut.

Les valeurs de pixel seront mises à l'échelle au nouveau type de pixel. Lorsque vous mettez à l'échelle la profondeur de vos pixels, votre raster affichera la même chose, mais les valeurs seront mises à l'échelle à la nouvelle profondeur de bits qui a été spécifiée.

Ceci est recommandé si vous exportez vers un jeu de données raster basé sur un fichier et que l'écrêtage graphique est choisi. Lorsqu'un graphique est utilisé pour découper vos données, des pixels NoData existeront très probablement dans la sortie. La spécification de la valeur NoData vous permet de contrôler la profondeur de pixel et la valeur qui stockera NoData, cependant, la valeur des pixels NoData n'est pas stockée pour les données exportées vers une géodatabase ou un format CRF.

Cochez la case Forcer RVB pour exporter le raster en sortie en tant que jeu de données raster RVB à trois canaux avec le moteur de rendu actuel.

De plus, lors de l'exportation vers des formats tels que TIFF, JP2, PNG et MRF qui prennent en charge les bandes alpha, vous pouvez utiliser cette option pour exporter les données en tant que jeu de données raster à quatre canaux avec une bande alpha pour préserver les paramètres de transparence des données d'origine .

La case à cocher Utiliser la palette de couleurs est activée uniquement si votre raster source contient une palette de couleurs ou si le rendu de couche raster actuel est un rendu à valeur unique.

Cochez la case Utiliser le rendu pour exporter le jeu de données raster avec les statistiques et les options actuelles du rendu. Lorsque vous ouvrez le jeu de données raster exporté dans ArcGIS Pro , les règles de rendu par défaut sont appliquées. Pour conserver le même rendu que lorsque vous avez exporté les données, définissez le type d'étirement sur Aucun , car il est déjà étiré.

  • La case à cocher Utiliser le rendu doit être cochée pour utiliser Force RVB ou Utiliser la palette de couleurs .
  • Les options Force RGB et Use Colormap ne peuvent pas être utilisées ensemble.
  • La combinaison d'Utiliser le rendu et de Forcer RVB exporte la sortie en tant que raster RVB avec trois ou quatre bandes (bande alpha le cas échéant) avec l'affichage du rendu de couche raster actuel lors de l'exportation.
  • La combinaison de Use Renderer et Use Colormap exporte le raster avec une palette de couleurs dans les cas suivants :
    • Si le raster source contient une palette de couleurs.
    • Pour les jeux de données qui ne contiennent pas de palette de couleurs, mais peuvent être visualisés dans un moteur de rendu à valeur unique.
    1. Choisissez le type de compression si votre format de sortie le permet.
    2. Choisissez la qualité de compression si votre format de sortie est JP2 ou JPG.

    Lorsque les rasters sont stockés sous forme de blocs de données, les jeux de données raster sont stockés dans un type de données appelé BLOB (binary large object). L'option de taille de tuile vous permet de contrôler le nombre de pixels stockés dans chaque BLOB et, par conséquent, vous permet de contrôler la taille de chaque BLOB. Il est spécifié comme le nombre de pixels dans X (largeur de la tuile) et Y (hauteur de la tuile).

    Le rééchantillonnage est le processus d'interpolation des valeurs de pixels tout en transformant votre jeu de données raster. Ceci est utilisé lorsque l'entrée et la sortie ne s'alignent pas exactement, lorsque la taille des pixels change, lorsque les données sont décalées, ou une combinaison de ceux-ci.

    Effectue une affectation de voisin le plus proche et est la plus rapide des méthodes d'interpolation. Il est principalement utilisé pour les données discrètes, telles qu'une classification de l'utilisation des terres, car il ne modifiera pas les valeurs des pixels. L'erreur spatiale maximale sera d'un demi-pixel.

    Effectue une interpolation bilinéaire et détermine la nouvelle valeur d'un pixel en fonction d'une moyenne de distance pondérée des quatre centres de pixel d'entrée les plus proches. Il est utile pour les données continues et entraînera un certain lissage des données.

    Effectue une convolution cubique et détermine la nouvelle valeur d'un pixel en fonction de l'ajustement d'une courbe lisse passant par les 16 centres de pixels d'entrée les plus proches. Il convient aux données continues, bien qu'il puisse en résulter que le raster en sortie contienne des valeurs en dehors de la plage du raster en entrée. Il est géométriquement moins déformé que le raster obtenu en exécutant l'algorithme de rééchantillonnage du voisin le plus proche. L'inconvénient de l'option Cubic est qu'elle nécessite plus de temps de traitement. I Si le temps de traitement est un problème, utilisez plutôt Bilinéaire.

    Il n'y a pas de type de données spécifié.

    Le raster est un type de données d'altitude.

    Le raster est un type de données thématique qui a des valeurs discrètes, telles que la couverture terrestre.

    La trame a été traitée en couleur et ne doit pas être étirée en contraste.

    Le raster contient des informations scientifiques et sera affiché avec la rampe de couleur bleu à rouge par défaut.

    Le raster est un raster à deux canaux qui contient une composante U et une composante V des données de champ vectoriel.

    Le raster est un raster à deux canaux qui contient l'amplitude et la direction des données de champ vectoriel.

    1. Spécifiez le nombre de niveaux de pyramide . Vous pouvez spécifier le nombre de niveaux à créer, ou vous pouvez laisser la valeur vide pour créer tous les niveaux.
    2. Cochez la case Ignorer le premier pour ignorer le premier niveau de pyramide de votre raster.
    3. Spécifiez la technique de rééchantillonnage de la pyramide utilisée pour construire vos pyramides : Voisin le plus proche , Bilinéaire ou Cubique .
    4. Choisissez le type de compression de la pyramide à utiliser lors de la création des pyramides raster.
    • Par défaut : le système détectera un type de compression approprié. Si les données source sont compressées à l'aide d'une compression par ondelettes, elle construira des pyramides avec le type de compression JPEG sinon, LZ77 sera utilisé.
    • Aucun : aucune compression n'est utilisée lors de la création de pyramides.
    • LZ77 : l'algorithme de compression LZ77 sera utilisé pour construire les pyramides. LZ77 peut être utilisé pour tout type de données.
    • JPEG : l'algorithme de compression JPEG sera utilisé pour créer des pyramides. Seules les données conformes à la spécification de compression JPEG peuvent utiliser ce type de compression. Si JPEG est choisi, vous pouvez alors définir la qualité de compression .
    • JPEG YCbCr : une compression avec perte utilisant les composants de l'espace colorimétrique luma (Y) et chroma (Cb et Cr).
    1. Choisissez si vous souhaitez sauter des pixels entre les échantillons. Les paramètres X Skip Factor et Y Skip Factor représentent respectivement le nombre de pixels horizontaux et verticaux entre les échantillons. La valeur doit être supérieure à zéro et inférieure ou égale au nombre de colonnes ou de lignes du jeu de données raster.
    2. Le paramètre Statistiques ignorer la ou les valeurs vous permet d'ignorer une ou plusieurs valeurs qui ne participeront pas au calcul des statistiques, comme une valeur d'arrière-plan. Les valeurs multiples sont séparées par des points-virgules.

    Une fois l'exportation terminée, le jeu de données raster exporté est ajouté à la carte.


    Syntaxe de script

    ArcGISRasterToLines_GeoEco (inputRaster, outputFeatureClass, backgroundValue, minDangleLength, simplifier, champ, projectedCoordinateSystem, GeographicTransformation, resamplingTechnique, projectedCellSize, registrationPoint, clippingDataset, clippingRectangle, mapAlgebraExpression)

    Le raster sera converti en une classe d'entités linéaires à l'aide de l'outil ArcGIS Raster to Polyline. Pour chaque paire de cellules raster de premier plan adjacentes, l'outil trace une ligne reliant leurs centres. Cet algorithme est approprié pour convertir des entités raster de type ligne, telles que les fronts de température de surface de la mer ou d'autres données de limites, en entités vectorielles.

    Classe d'entités surfaciques en sortie qui contiendra les lignes converties.

    Les répertoires manquants dans ce chemin seront créés s'ils n'existent pas.

    Spécifie la valeur de cellule qui identifiera les cellules d'arrière-plan. Le jeu de données raster est considéré comme un ensemble de cellules de premier plan et de cellules d'arrière-plan. Les caractéristiques linéaires sont formées à partir des cellules de premier plan.

    ZERO - L'arrière-plan est composé de cellules de zéro ou moins ou de NoData. Toutes les cellules avec une valeur supérieure à zéro sont considérées comme une valeur de premier plan.

    NODATA - L'arrière-plan est composé de cellules NoData. Toutes les cellules avec des valeurs valides appartiennent au premier plan.

    Longueur minimale des lignes pendantes qui seront conservées. La valeur par défaut est zéro.

    Si True (valeur par défaut), les lignes de sortie seront lissées selon un algorithme non documenté implémenté par l'outil ArcGIS Raster to Polyline.

    Champ utilisé pour affecter les valeurs des cellules du raster en entrée aux lignes du jeu de données en sortie. Il peut s'agir d'un nombre entier ou d'un champ de chaîne.

    Nouveau système de coordonnées vers lequel projeter le raster.

    Le raster ne peut être projeté dans un nouveau système de coordonnées que si la projection d'origine est définie. Une erreur sera générée si vous spécifiez un nouveau système de coordonnées sans définir le système de coordonnées d'origine.

    L'outil ArcGIS Project Raster est utilisé pour effectuer la projection. La documentation de cet outil recommande de spécifier également une taille de cellule pour le nouveau système de coordonnées.

    J'ai remarqué que pour certains systèmes de coordonnées, l'outil ArcGIS 9.2 Project Raster semble couper le raster projeté dans une mesure arbitraire trop petite. Par exemple, lors de la projection d'une image globale de chlorophylle MODIS Aqua de 4 km en coordonnées géographiques vers Lambert_Azimuthal_Equal_Area avec un méridien central de -60 et une latitude d'origine de -63, l'image résultante est tronquée pour ne montrer qu'un quart de la planète. Ce problème ne se produit pas lorsque Project Raster est appelé de manière interactive à partir de l'interface utilisateur d'ArcGIS, il se produit uniquement lorsque l'outil est appelé par programmation (la méthode ProjectRaster_management du géoprocesseur). Ainsi, vous ne le verrez peut-être pas lorsque vous utilisez Project Raster vous-même, mais cela peut arriver lorsque vous utilisez des outils MGET qui appellent Project Raster dans le cadre de leurs opérations de géotraitement.

    Si vous rencontrez ce problème, vous pouvez le contourner comme ceci :

    Tout d'abord, exécutez cet outil sans spécifier de nouveau système de coordonnées, pour obtenir le raster dans le système de coordonnées d'origine.

    Dans ArcCatalog, utilisez l'outil Projeter un raster pour projeter le raster dans le nouveau système de coordonnées. Vérifiez que tout le raster est présent, qu'il n'a pas été découpé dans une mesure trop petite.

    Dans ArcCatalog, recherchez l'étendue du raster projeté en cliquant dessus avec le bouton droit dans l'arborescence du catalogue, en sélectionnant Propriétés et en faisant défiler jusqu'à Etendue.

    Maintenant, avant d'exécuter l'outil MGET qui projette le raster, définissez le paramètre d'environnement Etendue sur les valeurs que vous avez recherchées. Si vous appelez l'outil MGET de manière interactive à partir d'ArcCatalog ou d'ArcMap, cliquez sur le bouton Environnements dans la boîte de dialogue de l'outil, ouvrez Paramètres généraux, modifiez la liste déroulante Étendue sur "Comme spécifié ci-dessous" et saisissez les valeurs que vous avez recherchées. Si vous l'invoquez à partir d'un modèle de géotraitement, cliquez avec le bouton droit sur l'outil dans le modèle, sélectionnez Créer une variable, À partir de l'environnement, Paramètres généraux, Étendue. Cela placera l'étendue en tant que variable dans votre modèle, attachée à l'outil MGET. Ouvrez la variable Etendue, remplacez-la par "Comme spécifié ci-dessous" et saisissez les valeurs que vous avez recherchées. Si vous appelez l'outil MGET par programme, vous devez définir la propriété Extent du géoprocesseur sur les valeurs que vous avez recherchées. Veuillez consulter la documentation ArcGIS pour plus d'informations à ce sujet et les paramètres d'environnement en général.

    Exécutez l'outil MGET. L'étendue du raster doit maintenant être de la bonne taille.

    Une méthode de transformation utilisée pour convertir entre le système de coordonnées d'origine et le nouveau système de coordonnées.

    Ce paramètre est une nouvelle option introduite par ArcGIS 9.2. Vous devez disposer d'ArcGIS 9.2 pour utiliser ce paramètre.

    Ce paramètre n'est nécessaire que lorsque vous spécifiez que le raster doit être projeté dans un nouveau système de coordonnées et que ce nouveau système utilise un système de référence différent du système de coordonnées d'origine, ou qu'il existe une autre différence entre les deux systèmes de coordonnées qui nécessite une transformation. Pour déterminer si une transformation est nécessaire, je recommande la procédure suivante :

    Tout d'abord, exécutez cet outil sans spécifier de nouveau système de coordonnées, pour obtenir le raster dans le système de coordonnées d'origine.

    Ensuite, utilisez l'outil ArcGIS 9.2 Project Raster sur le raster pour le projeter dans le système de coordonnées souhaité. Si une transformation géographique est nécessaire, cet outil vous en demandera une. Notez le nom exact de la transformation que vous avez utilisée.

    Enfin, si une transformation était nécessaire, saisissez le nom exact dans cet outil, réexécutez-le et vérifiez que le raster a été projeté comme vous le souhaitez.

    Algorithme de rééchantillonnage à utiliser pour projeter le raster d'origine dans un nouveau système de coordonnées. L'outil ArcGIS Project Raster est utilisé pour effectuer la projection et accepte les valeurs suivantes :

    NEAREST - interpolation du voisin le plus proche

    BILINAIRE - interpolation bilinéaire

    Vous devez spécifier l'un de ces algorithmes pour projeter vers un nouveau système de coordonnées. Une erreur sera générée si vous spécifiez un nouveau système de coordonnées sans sélectionner d'algorithme.

    La taille de cellule du système de coordonnées projeté. Bien que ce paramètre soit facultatif, pour obtenir les meilleurs résultats, la documentation ArcGIS vous recommande de toujours le spécifier lors de la projection dans un nouveau système de coordonnées.

    Les coordonnées x et y (dans l'espace de sortie) utilisées pour l'alignement des pixels.

    Ce paramètre est une nouvelle option introduite par ArcGIS 9.2. Vous devez disposer d'ArcGIS 9.2 pour utiliser ce paramètre. Il est ignoré si vous ne spécifiez pas que le raster doit être projeté dans un nouveau système de coordonnées.

    Classe d'entités, raster ou autre jeu de données géographiques existant ayant l'étendue à laquelle le raster doit être découpé.

    AVERTISSEMENT : si vous utilisez cet outil dans un modèle de géotraitement ArcGIS et que vous sélectionnez un jeu de données en cliquant sur l'icône de dossier et en accédant au jeu de données, votre sélection peut mystérieusement disparaître de cette zone de texte après la fermeture de l'outil. Il s'agit d'un bogue dans ArcGIS. Pour contourner ce problème, faites glisser et déposez l'ensemble de données souhaité dans le modèle. Cela créera une couche dans le modèle pour cet ensemble de données. Sélectionnez ensuite ce calque dans cet outil en cliquant sur la liste déroulante au lieu de cliquer sur l'icône du dossier. Le calque sélectionné ne doit pas disparaître lorsque vous fermez l'outil.

    Rectangle auquel le raster doit être découpé.

    Si un système de coordonnées projetées a été spécifié, le découpage est effectué après la projection et les coordonnées du rectangle doivent être spécifiées dans le nouveau système de coordonnées. Si aucun système de coordonnées projetées n'a été spécifié, les coordonnées doivent être spécifiées dans le système de coordonnées d'origine.

    L'outil ArcGIS Clip est utilisé pour effectuer le clip. Le rectangle de découpage doit être transmis à cet outil sous la forme d'une chaîne de quatre nombres séparés par des espaces. L'interface utilisateur d'ArcGIS formate automatiquement la chaîne correctement lors de l'appel de cet outil à partir de l'interface utilisateur d'ArcGIS, vous n'avez pas à vous soucier du format. Mais lorsque vous l'invoquez par programmation, veillez à fournir une chaîne correctement formatée. Les numéros sont ordonnés GAUCHE, BAS, DROITE, HAUT. Par exemple, si le raster se trouve dans un système de coordonnées géographiques, il peut être découpé à 10 W, 15 S, 20 E et 25 N avec la chaîne :

    Des nombres entiers ou décimaux peuvent être fournis.

    Mapper l'expression d'algèbre à exécuter sur le raster.

    ATTENTION: Le générateur de modèle de géotraitement ArcGIS peut supprimer de manière aléatoire et silencieuse la valeur de ce paramètre. Il s'agit d'un bogue dans ArcGIS. Avant d'exécuter un modèle que vous avez enregistré, ouvrez cet outil et vérifiez que la valeur du paramètre existe toujours.

    L'expression est exécutée une fois que le raster converti est projeté et découpé (si ces options sont spécifiées). Utilisez la chaîne inputRaster sensible à la casse pour représenter le raster sur lequel vous souhaitez maintenant effectuer une algèbre cartographique. Par exemple, pour convertir le raster en un raster d'entiers et ajouter 1 à toutes les cellules, utilisez cette expression :

    La chaîne inputRaster est sensible à la casse. Avant d'exécuter l'expression d'algèbre de carte, la chaîne est remplacée par le chemin d'accès à un raster temporaire qui représente le raster généré. L'expression finale doit comporter moins de 4 000 caractères, sinon ArcGIS signalera une erreur.

    L'outil ArcGIS Single Output Map Algebra est utilisé pour exécuter l'expression d'algèbre cartographique. Vous devez disposer d'une licence pour l'extension ArcGIS Spatial Analyst afin d'effectuer l'algèbre cartographique.

    La syntaxe de l'algèbre cartographique peut être très pointilleuse. Voici quelques conseils qui vous aideront à réussir avec cet outil :

    Avant d'utiliser cet outil, construisez et testez votre expression d'algèbre cartographique à l'aide de l'outil ArcGIS Single Output Map Algebra. Collez ensuite l'expression dans cet outil et modifiez-la pour utiliser la variable inputRaster plutôt que la valeur de test que vous avez utilisée avec Single Output Map Algebra.

    Si vous développez votre expression directement dans cet outil, commencez par une expression très simple. Vérifiez qu'il fonctionne correctement, ajoutez-en un peu et vérifiez à nouveau. Répétez ce processus jusqu'à ce que vous ayez construit l'expression complète.

    Séparez toujours les opérateurs mathématiques des chemins raster à l'aide d'espaces. Dans l'exemple ci-dessus, l'opérateur / contient un espace de chaque côté. Suivez ce modèle. In some circumstances, ArcGIS will fail to process raster algebra expressions that do not separate raster paths from operators using spaces. The reported error message usually does not indicate that this is the problem, and tracking it down can be very frustrating.


    Request parameters

    List of input rasters. The input raster can be the Portal Item ID, Image Service URL, cloud raster dataset, or shared multidimensional raster dataset.

    At least one type of input needs to be provided in the JSON object. If multiple inputs are given, the itemIds takes the priority.

    Syntax: List of JSON object describing the input rasters.

    Data identifying positions where you want a sample taken.

    The input can be an image service or a feature service.

    Syntax: JSON object describing the input raster or feature. At least one type of input needs to be provided in the JSON object. If multiple inputs are given, the itemid takes the priority.

    Name of the output table or feature service holding the sampled cell values.

    Syntax: JSON object describing the output table

    • NEAREST —Nearest neighbor assignment. This is the default.
    • BILINEAR —Bilinear interpolation
    • CUBIC —Cubic convolution

    Syntax: A string representing the resamplingType .

    A field containing a different value for every location or feature in the input location raster or point features.

    Syntax: A string representing the field.

    Specify the time, depth or other acquisition data associated with the location features.

    Statistics will be calculated for variables within the dimension range of the following combinations:

    • Dimension + Start field or value
    • Dimension + Start field or value + End field or value
    • Dimension + Start field or value + Relative value or days before + Relative value or days after

    Only non-negative values are supported for the following:

    Syntax: a list of dictionary objects.

    The type of statistic to be calculated.

    • MINIMUM —Finds the minimum value within the specified range.
    • MAXIMUM —Finds the maximum value within the specified range.
    • MEDIAN —Finds the median value within the specified range.
    • MEAN —Calculates the average for the specified range.
    • SUM —Calculates the sum of the variables within the specified range.
    • MAJORITY —Finds the value that occurs most frequently.
    • MINORITY —Finds the value that occurs least frequently.
    • STD —Calculates the standard deviation.
    • PERCENTILE —Calculates a defined percentile within the specified range.

    A string representing the statisticsType .

    The percentile to calculate when the statisticsType parameter is set to PERCENTILE .

    This value can range from 0 to 100. The default is 90.

    Syntax: A double representing the percentileValue .

    The specified distance around the location data features. The buffer distance is specified in the linear unit of the location feature's spatial reference. If the feature uses a geographic reference, the unit will be in degrees.

    Statistics will be calculated within this buffer area.

    Syntax: A value representing the bufferDistance .

    Syntax: A string representing the buffer distance field in the inLocationData .

    • ROW_WISE —Sampled values appear in separate rows in the output table. This is the default.
    • COLUMN_WISE —Sampled values appear in separate columns in the output table. This option is only valid when the input multidimensional raster contains one variable and one dimension, and each slice is a single-band raster.

    Syntax: A string of one of the keywords.

    Boolean value to determine if this tool generates an output feature service containing a feature class with sampled values or only a table with sampled values. The default is false .

      ( extent )—A bounding box that defines the analysis area. ( mask )—Only cells that fall within the analysis mask will be considered in the operation. ( outSR )—The output raster will be projected into the output spatial reference. ( parallelProcessingFactor )—The specified number or percentage of processes will be used for the analysis.
    • True: Samples will be taken for all dimensions (such as time or depth) of a multidimensional dataset.
    • False: Samples will be taken from the current slice of a multidimensional dataset. This is the default.
    • NEAREST—Uses the value of the closest cell to assign a value to the output cell when resampling. This is the default.
    • BILINEAR—Determines the new value of a cell based on a weighted distance average of the four nearest input cell centers.
    • CUBIC—Determines the new value of a cell based on fitting a smooth curve through the 16 nearest input cell centers.

    The response format. The default response format is html.


    Environnements

    Analysis environment settings are additional parameters that affect a tool's results. You can access the tool's analysis environment settings by clicking the gear icon at the top of the tool pane.

    This tool honors the following Analysis Environments :

    • Extent—Specifies the area to be used for analysis.
    • Mask—Specifies a mask layer, where only the cells that fall within the mask area will be used for analysis.
    • Recycle interval of processing workers—Defines how many image sections to process before restarting worker processes.
    • Parallel processing factor—Controls the raster processing CPU or GPU instances.
    • Number of retries on failures—Defines how many retries a worker process will attempt when there is random failure processing a job.

    Convert the units of measurement

    The solar radiation raster uses watt-hours per square meter as its unit of measurement. According to the legend in the Contents pane, some cells have values of over 1 million (expressed with the notation e+06). To reduce the size of these values and make them easier to read, you will convert the raster layer to kilowatt-hours per square meter (kWh/m 2 ).

    1. In the Contents pane, uncheck Building_Footprints to turn it off.
    2. In the Geoprocessing pane, click the Back button (you may need to click it twice). Search for and open Raster Calculator (Spatial Analyst Tools) .

    There are 1,000 watts in a kilowatt, so to convert the units of measurement, you only need to create an expression that divides existing cell values by 1,000.

    If you created your own solar radiation raster layer, the name of the layer in the expression will be Solar_Rad_Whm2_Example instead of Solar_Rad_Whm2 .

    The new raster layer is created and added to the map. It is similar to the original solar radiation layer, but the values are 1,000 times smaller.

    You no longer need the original solar radiation layer, so you will remove it.


    1 réponse 1

    I was able to devise a method to accomplish this task (tho likely not in the most elegant fashion): as described above I used zonal statistics to create a new raster with mean values of the raster (depth) per site (delineated by polygon layer). I then used the same tool to create a raster layer representing the standard deviation per site. Puis dans rater calculator I set all values to NULL which were less than 1 standard deviation from the mean depth per site. SetNull("depth"> "mean_depth" - "sd_depth", "depth")-- and this created a new raster with only those pixel values greater than 1 std from mean depth (ie the deepest habitat per site). note: because depth values were negative, we used > (greater than)


    Convert raster to polygon based on cell value ranges using ArcGIS for Desktop? - Systèmes d'information géographique

    There are many different types of raster analysis available in ArcGIS. Here are just a few of the common analytical functions. More analytical functions on raster surfaces models will be dealt with in 3-D and Surface Modeling.

    Importing data from generic raster files

    Many datasets are available on the World Wide Web. Most of the raster datasets are in a generic format. Those formats that can be imported into ArcGIS are

    • The ASCII raster file format
    • The binary raster file format
    • The USGS Digital Elevation Model (DEM) raster file format
    • The US DMA DTED raster file format

    The most common format you are likely to see is the USGS DEM. There is a page of Washington 10 and 30 m DEMs for USGS 7.5' quad sheet boundaries on a server in Geological Sciences.

    En utilisant ArcToolbox it is possible to import these data sets. The Spatial Analyst Extension must be activated in order to import raster data. Imported raster data will be converted to the ArcInfo raster grid data format.

    Here, a USGS DEM was downloaded from the site in Geological Sciences and unzipped. Importing the grid is very straightforward using ArcGIS 's GUI.

    The same basic process is used to import from the other raster interchange file formats.

    Merging adjacent grids ("mosaicking")

    Merging or mosaicking adjacent grids is used when your study area falls across several grids, and you wish to treat those grids as a single grid. This is commonly used when the data source is the USGS series of DEMs. Because DEMs are created and distributed as tiles, if your study area falls across several tiles, it is often necessary to merge these tiles together.

    In this example, I have downloaded and imported the Elbe DEM as well as the Eatonville DEM. The following images show the grid created from mosaicking the two inputs (before [above], and after [below]).

    Calculating distance surfaces and buffers

    Distance surfaces are grids whose output value is the distance to the closest feature in the input layer. The input layer can be a selected set of any type of feature (point, line, polygon, or grid cell). Distance surfaces are calculated by using the Spatial Analyst> Distance menu choice in the Spatial Analyst toolbar.

    Distance surfaces are similar to buffers in the vector world. The difference between vector buffering and creating distance grids is that the distance surface represents a continu change in distance from the source as you move across the landscape, whereas the buffer analysis changes in user-defined quantized steps.

    Here are the streams of Pack Forest and a distance surface created from them:

    A distance grid is calculated. Every cell in the output dataset is assigned cell value equal to the straight-line distance to the closest stream line feature. Those cells closest to the stream are light yellow in color, and those farthest away are blue (note the southwest corner). Note that this is different from a buffer, which only gives an "inside or outside" encoding of the output data set. Although the image shows "rings" of distance classes, the underlying data are continuous in value.

    Determining proximity

    Determining proximity is similar to calculating a distance surface, but rather than creating a continuous surface whose value is the distance to a feature, the proximity grid contains values in the cells for a corresponding value in the input feature attribute table. Each cell is coded with the closest feature's value from the input layer, rather than for the distance to features.

    Proximity analysis uses as input the selected set of the active layer, and is available from the menu at Spatial Analyst > Distance > Allocation.

    Here, proximity is calculated for some bird nest points. In the output, the value for any given cell in the output Allocation grid layer is the sequential ID number for the closest nest.

    Every cell is encoded for the value of the closest point. This means that, for example, anywhere within grid zone 5 is closer to point 5 than to any other point. The cells on the edge between zone 4 and 5 are equidistant to either point.

    This technique is also known as Thiessen ou alors Voronoi analysis.

    Creating surfaces from point samples

    Frequently point samples are taken to because it is too costly (either in terms of time or money) to sample an entire population. It is possible to generate interpolated surfaces based on point samples. The cells between the sampling points are given a value that represents a smooth transition of value between the sampling points. If you need an estimate of a value somewhere that you do not have a sampling point, you can get a grid value at that spot. Be careful here, because the assumption that values change smoothly across the landscape is not necessarily true! This type of analysis is well-suited to data that definitely do change gradually over a large area, such as precipitation. In any case, if your sampling points are spread too far apart, you may create an interpolated grid that does not capture local variations.

    Here is a surface generated from the Pack Forest CFI plot centers using a Regularized Spline method. Red indicates low standing wood volume and green indicates high standing wood volume for conifer trees in 1994. Plot centers are also displayed here for illustration.

    There are a number of different options for creating surfaces from point samples. If you need to perform surface interpolation from points, you should read the help documents thoroughly.

    If you have data representing a continuous surface, it is possible to create single contour lines for a given grid cell value, or to create a whole group of contour lines at a regular interval. This can be of value if you wish to create a contour map of any continuously changing surface. Although digital vector elevation contours are available for some USGS quad sheets, many areas of the state have not been digitized yet. However, we do have complete statewide coverage for DEMs. These DEMs can be used to create contour lines that can be added to maps.

    Here are contours made from 30 m DEMs:

    Calculating summary attributes for features using a grid layer ("Zonal statistics ")

    Zones in one grid layer can be defined by either polygons or zones of integer grids. For areas within different polygons, or for zones within an integer grid, the input grid values are summarized. The output is a table in which a single record exists for the unique values in the chosen field in the zone-defining layer. Each record in the output table contains the fields Surface, Min, Max, Varier, Signifier, Std, et Somme.

    In this example, the zone-defining layer is Stands. The individual zones are polygons containing the same value for the SITE_INDEX domaine. This means that for every unique occurrence of a site index value in the Stands layer, a new grid zone will be defined (even if the stands are not contiguous). The layer to be summarized is Dém.

    The statistic shown on the graph is the Signifier, that is, the mean elevation within each unique site index zone. In this case, each data marker in the graph signifies the mean elevation for all cells within stands with that site index.

    Based on the graph, the stands with the greatest site index (a proxy measure for productivity) also have the highest mean elevation.

    Cross tabulating areas

    Cross-tabulation allows you to compare the area of one specific value in integer grid layer against one specific value in an another integer grid layer. The input layers and fields are defined in an ArcToolbox tool.

    In this example, the Espèce field in the Stands layer is compared against the Soil.name field in the Sols couche.

    The output table contains a unique value for each Espèce record, and fields representing unique values from the Sols couche.

    Here is the stands grid table showing the values corresponding to the field names in the cross-tabulation table (i.e., Value = 1 corresponds to soil_stand_xtab.VALUE_1):

    The values in the fields are the area (in map units) for the spatial overlap between the classes in the input layers. For example, in Kapowsin soils, there are 2933100 ft^2 of Mixed Redcedar stands.

    If you have two layers representing the same data for a study area at different times, you can use cross-tabulation for change analysis. Tabulation can be used any combination of (integer) grid layers.

    Cross-tabulating areas is a raster analysis technique. When tabulating areas for polygon layers, you need to first convert from polygon to grid. You should select a cell size that will capture the detail of the features in the polygon data. The smaller the cell size, the greater the precision, but the longer processing will take..

    This is a very powerful technique for change analysis. If you have datasets representing two different time slices, you can compare the area of such attributes as land cover or zoning designations.

    "Querying" across multiple grid layers

    While the normal feature attribute table query allows a query only on a single layer, the Raster Calculator allows you to make a complex query based on multiple layers. These types of queries are simple to perform as long as the grid layers representing the properties in question are contained in a single data frame. To do the same query in the vector world requires polygon layers representing the layers (which is in itself a problem, since vector layers are not good at representing continuous phenomena), and the performance of multiple topological overlay operations.

    In this example, I am interested in finding cells closer than 300 ft from a stream, with an elevation > 1500, with greater than 6,000 bd-ft timber volume.

    Those cells displayed in green meet the criteria (coded with a value of 1).

    How would you go about getting the answer to the same query if you only had access to vector data and vector processing?

    Calculating neighborhood statistics

    Neighborhood statistics are the focal functions referred to in Raster Analysis I. The neighborhood is defined as the group of cells for which statistics will be calculated. The neighborhood (a.k.a. noyau ou alors concentrer) can be shaped as a circle, rectangle, ring, or wedge. Statistics available are

    • Minimum
    • Maximum
    • Signifier
    • Médian
    • Somme
    • Varier
    • Écart-type
    • Majorité
    • Minorité
    • Variété

    The processor looks in the neighborhood, identifies cells or point features within that neighborhood, and calculates a single statistic for that neighborhood. That single value is then placed in the output grid in the cell located at the center of the neighborhood. The process is performed for every input cell location in the analysis window.

    It is possible to perform neighborhood statistics on point layers. If a point layer has a numeric field, the process is performed for the entire area within the analysis window, and the statistic is generated for the points located within the kernel at each output cell location.

    A typical use of neighborhood statistics is known as "filtering." A "low pass" filter is nothing more than a 3 by 3 cell focal mean performed for an entire grid. Low pass filters smooth out anomalies and peaks in surfaces.

    A "high pass" filter is also a 3 by 3 focal function, but rather than taking the mean of the 9-cell window, it performs a focal sum of the kernel cells, but first multiplies the cells by these coefficients:

    There are several different coefficients that can be used in a high-pass filter, but they all have the objective of sharpening edges. ArcGIS 's default high-pass filter uses these particular coefficients.

    In this example, an input grid represents several different vegetation zones (stand age).

    The high pass filter makes the zone interiors the same value (0), while the edges get either a high or low value. The edges are most pronounced where the contrast is greatest.

    This analysis is performed using the Filter tool:

    With this grid as the result:

    Edges can be used to define places where animal movement may be hindered, or where species that prefer ecotones may be found.

    Conditional processing

    Conditional processing is a method of creating new grids based on an "if-then" condition. For example, we may be interested in reclassifying cells that have a certain value, but leaving other cells with their original value, this is possible with a reclassification. However, reclassification can be tedious (setting up the output classes), whereas conditional processing can create the new grid based on specific rules rather than simple numerical transformations. The conditions can also include several grids, rather than reclassifying based only on the values within a single grid.

    Going back to the mosaicked Eatonville/Elbe grids, all cells between 500 and 700 m in elevation are multiplied by 100, and anything else is coded with a value of 0.

    The expression means this line-by-line, in English:

    If elevation is greater than 500 and less than 1700, then
    set output value to (elevation * 100), or else
    set output value to 0

    Here is the resultant grid:

    Conditional processing is very useful when you need to select out or analyze a specific group of cells in one way, and another group of cells in another way.

    Converting raster and vector data sources

    It is possible to go back and forth between raster and vector formats. This always is at the expense of the loss, or generalization, of shapes. Any feature layer can be converted to a grid layer.

    Vector to Raster:
    Points are converted to single cells. Lines are converted to groups of cells oriented in a linear arrangement. Polygons are converted to zones. In all cases, only selected features are converted, or all features if no selection is active.

    Raster to Vector:
    Grid layers can only be converted directly to polygon vector layers. Be careful, because a new polygon will be created based on the field that is used for the conversion. If you have an elevation grid layer and you convert this to a polygon feature layer based on the Valeur field, you will get a very large number of very small polygons, and this will take a long time. It is more customary to first reclassify grids to create zones, and then convert these zones to polygon features.

    Here, the Pack Forest dem has been reclassified into 100-ft elevation bands and then saved as a shape file. The value for the new polygon attribute Gridcode matches the original Valeur field from the input grid data source.

    The new polygon layer is displayed in a graduated color classification based on the Gridcode domaine.

    Ici le ruisseaux line feature layer has been converted to a grid layer based on the DNR_TYPE domaine.

    Once a raster dataset has been converted to vector format, all of the vector analysis and overlay tools can be used. Likewise, when a vector datasets is converted to a grid, it can be used in raster analytical techniques.


    Voir la vidéo: ArcGIS - Raster to Polygon conversion (Octobre 2021).