| Étiquette | Explication | Type de données |
Zone d’étude en entrée | Zone d’étude en entrée dans laquelle les localisations d’échantillonnage sont créées. La zone d’étude doit correspondre à des polygones ou à un raster entier (catégoriel). Pour les rasters, les cellules dont les valeurs sont nulles ne sont pas incluses dans la zone d’étude. | Feature Layer; Raster Layer |
Entités en sortie | Entités en sortie représentant les localisations d’échantillonnage. Pour un échantillonnage aléatoire simple et l’échantillonnage stratifié, les entités en sortie sont des points. Pour un échantillonnage par agrégat, la sortie correspond à des polygones. Pour un échantillonnage systématique, la sortie peut correspondre à des points ou des polygones. | Feature Class |
Méthode d'échantillonnage (Facultatif) | Spécifie la méthode d’échantillonnage utilisée pour créer les localisations d’échantillonnage.
| String |
Champ d’ID de strate (Facultatif) | Pour un échantillonnage stratifié par champ d’ID de strate, champ d’ID de strate définissant la strate. | Field |
Strata Sample Count Allocation Method (Méthode d’allocation du nombre d’échantillonnage des strates) (Facultatif) | Pour un échantillonnage stratifié, spécifie la méthode utilisée pour déterminer le nombre de localisations d’échantillonnage à créer dans chaque strate.
| String |
Bin Shape (Forme du groupe) (Facultatif) | Pour un échantillonnage systématique et par agrégat, spécifie la forme de chaque polygone de la tessellation maillée.
| String |
Bin Size (Taille du groupe) [nombre ou surface] (Facultatif) | Pour un échantillonnage systématique et par agrégat, taille de chaque polygone de la tessellation. Cette valeur peut être spécifiée sous forme de nombre (nombre total de polygones tesselés créés dans la zone d’étude) ou de surface (surface de chaque polygone tesselé). Pour une entrée sous forme de nombre, la valeur par défaut est de 100. Pour une entrée sous forme de surface, une valeur doit être spécifiée. Si un nombre est spécifié, l’outil tente de créer le nombre spécifié de localisations d’échantillonnage. Si le nombre exact ne peut pas être créé, un avertissement est renvoyé. | Areal Unit; Long |
Résolution H3 (Facultatif) | Pour un échantillonnage systématique ou par agrégat avec des groupes d’hexagones H3, spécifie la résolution H3 des hexagones. À chaque augmentation de la valeur de la résolution, la surface des polygones est divisée par sept.
| Long |
Nombre d’échantillons (Facultatif) | Nombre de localisations d’échantillonnage à créer. Ce paramètre s’applique toujours aux échantillonnages aléatoires simples et par agrégat. Pour un échantillonnage stratifié, ce paramètre s’applique si le nombre d’échantillons est proportionnel à la surface de la strate ou à un champ de population. Pour un échantillonnage aléatoire simple et un échantillonnage stratifié, la valeur par défaut est de 100. Pour un échantillonnage par agrégat, la valeur par défaut est de 10. | Long |
Nombre d’échantillons par strate (Facultatif) | Pour un échantillonnage stratifié avec un nombre d’échantillons égal dans chaque strate, nombre de localisations d’échantillonnage créées dans chaque strate. Le nombre total d’échantillons correspond à cette valeur multipliée par le nombre de strates. La valeur par défaut est 100. | Long |
Population Field (Facultatif) | Champ de population d’un échantillonnage stratifié si le nombre d’échantillons est égal ou proportionnel à un champ de population. | Field |
Type de géométrie en sortie (Facultatif) | Pour un échantillonnage systématique, indique si les localisations d’échantillonnage correspondent à des polygones tesselés ou des centroïdes (points) de polygones tesselés.
| String |
Distance minimale entre les points d’échantillonnage (Facultatif) | Pour un échantillonnage aléatoire simple et un échantillonnage stratifié, plus petite distance autorisée entre les localisations d’échantillonnage. Pour un échantillonnage aléatoire simple, tous les points sont séparés d’au moins cette distance. Pour un échantillonnage stratifié, les points d’une même strate sont séparés d’au moins cette distance, mais les points dans les strates voisines peuvent être plus proches que cette distance. Pour les grandes distances, il est possible de créer moins de localisations d’échantillonnage que prévu pour que les localisations restent suffisamment distantes les unes des autres. Dans ce cas, un message d’avertissement est renvoyé. | Linear Unit |
Relation spatiale (Facultatif) | Spécifie les polygones d’une tessellation en arrière-plan à inclure comme localisations d’échantillonnage. Ce paramètre s’applique à l’échantillonnage par agrégat et à l’échantillonnage systématique si le type de géométrie en sortie est le polygone.
| String |
Synthèse
Créez des localisations d’échantillonnage dans une zone d’étude continue à l’aide de conceptions d’échantillonnage aléatoire simple, stratifé, systématique (maillé) ou par agrégat.
L’échantillonnage est le processus de sélection d’individus au sein d’une population afin de les étudier et de tirer des conclusions sur l’ensemble de cette population. L’échantillonnage spatial continu traite la population comme une zone continue à partir de laquelle toute localisation ou zone peut être échantillonnée. Par exemple, cet outil permet de créer des localisations d’échantillonnage pour les arbres d’une forêt dense ou de collecter des mesures d’humidité du sol dans un champ de culture. Cet outil ne convient pas à l’échantillonnage des populations discrètes, telles que les foyers, les animaux ou les villes.
Illustration
Utilisation
La zone d’étude en entrée doit être une classe d’entités surfaciques ou un raster entier (catégoriel). Vous pouvez également dessiner la zone d’étude sur une carte à l’aide d’une entrée d’entité interactive. Pour les rasters, les cellules dont les valeurs sont nulles ne sont pas prises en compte dans cette zone d’étude.
Des localisations d’échantillonnage peuvent être créées pour les principales conceptions d’échantillonnage suivantes :
- Échantillonnage aléatoire simple : crée des points d’échantillonnage de manière aléatoire dans la zone d’étude. Chaque localisation de la zone d’étude a autant de chance d’être sélectionnée comme localisation d’échantillonnage. La zone d’étude est traitée comme une zone unique et toutes les limites entre les polygones ou les catégories de raster sont ignorées (par exemple, une classe d’entités surfaciques de tous les comtés d’un État définit la même zone d’étude qu’un seul polygone de tout l’État). L’échantillonnage aléatoire simple est utile si vous souhaitez étudier toute la zone d’étude, mais qu’aucune localisation n’est plus importante qu’une autre pour l’échantillonnage. Pour l’échantillonnage aléatoire simple, sélectionnez l’option Simple random (Aléatoire simple) du paramètre Sampling Method (Méthode d’échantillonnage).
- Exemple d’application : si la zone d’étude est une forêt dense dans laquelle il est supposé que toute localisation contient un arbre, l’échantillonnage aléatoire simple peut être utilisé pour échantillonner aléatoirement les arbres dans cette forêt.
- Échantillonnage aléatoire stratifié : crée des points d’échantillonnage en divisant la zone d’étude en strates distinctes (comme des classes de sols ou des types d’utilisation du sol) et en effectuant un échantillonnage aléatoire simple séparément dans chaque strate. L’échantillonnage aléatoire stratifié est utile si vous souhaitez vous assurer que toutes les strates sont représentées dans l’échantillon. Pour effectuer un échantillonnage aléatoire stratifié, sélectionnez l’une des trois options de stratification du paramètre Sampling Method (Méthode d’échantillonnage) (pour des informations sur chaque type de stratification, consultez le prochain conseil d’utilisation).
- Exemple d’application : si un parc national est divisé en classes d’altitudes, l’échantillonnage aléatoire stratifié peut être utilisé pour collecter des échantillons de sol séparément pour chaque classe d’altitude. De la sorte, le nombre d’échantillons de sol est suffisant pour toutes les altitudes du parc.
- Échantillonnage systématique : crée des localisations d’échantillonnage dans un modèle maillé non aléatoire à l’intérieur de la zone d’étude. La grille est créée par une tessellation de polygones de forme régulière (tels que des hexagones, carrés ou triangles). Les localisations d’échantillonnage peuvent être renvoyées comme polygones tesselés ou points (centroïdes des polygones tesselés). L’échantillonnage systématique est utile pour s’assurer qu’aucune section de la zone d’étude n’est échantillonnée plus que les autres, ce qui est souvent souhaitable si l’objectif consiste à créer une carte des échantillons et non de tirer des conclusions sur l’ensemble de la zone d’étude. Pour un échantillonnage systématique, sélectionnez l’option Systematic (Systématique) du paramètre Sampling Method (Méthode d’échantillonnage).
- Exemple d’application : pour étudier le plancher océanique d’une zone maritime, vous pouvez créer une grille hexagonale des localisations d’échantillonnage pour échantillonner les espèces de plantes maritimes.
- Échantillonnage par agrégat : crée des polygones d’échantillonnage en créant un échantillon systématique et en sélectionnant aléatoirement certains des polygones de la tessellation. Les polygones résultants sont appelés agrégats et, généralement, ces agrégats sont étudiés de manière exhaustive, par un échantillonnage aussi poussé que possible dans chaque agrégat. L’échantillonnage par agrégat est utile si vous vous intéressez particulièrement au mode d’interaction entre les échantillons à de courtes distances et que vous acceptez volontier que de grandes sections de la zone d’étude ne contiennent pas d’échantillon. Pour un échantillonnage par agrégat, sélectionnez l’option Cluster (Agrégat) du paramètre Sampling Method (Méthode d’échantillonnage).
- Exemple d’application : lors de l’échantillonnage de colonies d’insectes, l’échantillonnage par agrégat peut être utilisé pour créer de petites zones d’une parcelle. Toutes les colonies d’insectes à l’intérieur de ces agrégats seront alors échantillonnées.
- Échantillonnage aléatoire simple : crée des points d’échantillonnage de manière aléatoire dans la zone d’étude. Chaque localisation de la zone d’étude a autant de chance d’être sélectionnée comme localisation d’échantillonnage. La zone d’étude est traitée comme une zone unique et toutes les limites entre les polygones ou les catégories de raster sont ignorées (par exemple, une classe d’entités surfaciques de tous les comtés d’un État définit la même zone d’étude qu’un seul polygone de tout l’État). L’échantillonnage aléatoire simple est utile si vous souhaitez étudier toute la zone d’étude, mais qu’aucune localisation n’est plus importante qu’une autre pour l’échantillonnage. Pour l’échantillonnage aléatoire simple, sélectionnez l’option Simple random (Aléatoire simple) du paramètre Sampling Method (Méthode d’échantillonnage).
Pour un échantillonnage stratifié, vous pouvez définir les strates de trois manières. Chacune est disponible comme option du paramètre Sampling Method (Méthode d’échantillonnage) :
- Stratify by individual polygon (Stratifier par polygone) : chaque enregistrement d’une classe d’entités surfaciques correspond à une strate différente. Par exemple, si la zone d’étude est un champ avec des sous-parcelles stockées comme polygones distincts, des points d’échantillonnage sont créés séparément pour chaque sous-parcelle. La zone d’étude en entrée doit être des polygones.
- Stratify by contiguous raster region (Stratifier par région de raster contigu) : chaque région d’un raster entier (catégoriel) correspond à une strate. Une région de raster est un bloc contigu de cellules de même valeur (du champ Value) connectées par des contours de cellule partagés. Si deux régions possèdent la même valeur, mais sont déconnectées l’une de l’autre, elles se trouvent dans des strates différentes. La zone d’étude en entrée doit être un raster.
- Stratify by strata ID field (Stratifier par champ d’ID de strate) : tous les polygones ou cellules raster possédant la même valeur d’ID de strate correspondent à une strate. Il n’est pas nécessaire que les polygones ou les cellules raster soient contigus pour être dans la même strate. Spécifiez le champ contenant les valeurs des ID de strate dans le paramètre Strata ID Field (Champ d’ID de strate). Le champ doit être un entier ou du texte.
Vous pouvez spécifier le nombre d’échantillons à créer dans chaque strate à l’aide de l’une des options suivantes du paramètre Strata Sample Count Allocation Method (Méthode d’allocation du nombre d’échantillonnage des strates) :
- Equal count in each stratum (Nombre égal dans chaque strate) : un nombre égal d’échantillons sont créés dans chaque strate. Spécifiez la valeur du paramètre Number of Sample Per Stratum (Nombre d’échantillons par strate).
- Count proportional to stratum area (Nombre proportionnel à la surface de la strate) : le nombre d’échantillons de la strate est proportionnel à la taille de cette dernière. Spécifiez le nombre global d’échantillons dans le paramètre Number of Samples (Nombre d’échantillons) ; le nombre total est distribué dans chaque strate proportionnellement à sa surface.
- Count equal to population field (Nombre égal au champ de population) : le nombre d’échantillons dans chaque strate est égal aux valeurs d’un champ de population. Spécifiez le champ dans le paramètre Population Field (Champ de population). Le champ ne doit pas contenir de valeur négative et doit être de type entier.
- Count proportional to population field (Nombre proportionnel au champ de population) : le nombre d’échantillons dans chaque strate est proportionnel aux valeurs d’un champ de population. Spécifiez ce champ dans le paramètre Population Field (Champ de population) et le nombre global d’échantillons dans le paramètre Number of Samples (Nombre d’échantillons).
- Stratify by individual polygon (Stratifier par polygone) : chaque enregistrement d’une classe d’entités surfaciques correspond à une strate différente. Par exemple, si la zone d’étude est un champ avec des sous-parcelles stockées comme polygones distincts, des points d’échantillonnage sont créés séparément pour chaque sous-parcelle. La zone d’étude en entrée doit être des polygones.
L’outil permet alors de créer des conceptions d’échantillonnage plus avancées, qui ne sont pas disponibles comme options explicites du paramètre Sampling Method (Méthode d’échantillonnage).
- Échantillonnage par agrégat en deux étapes : crée des agrégats de points dans l’ensemble de la zone d’étude en créant un échantillon par agrégat, puis en créant des points (aléatoires simples, stratifiés ou systématiques) dans chaque agrégat. Cette conception d’échantillonnage est utile si un échantillon par agrégat est requis, mais qu’il n’est pas possible d’étudier de manière exhaustive chaque polygone d’agrégat. Elle s’avère également utile si vous vous intéressez principalement au mode d’interaction entre les échantillons à de courtes distances. Pour un échantillonnage par agrégat en deux étapes, utilisez l’outil pour créer un échantillon par agrégat, puis utilisez les polygones d’agrégat comme zone d’étude en entrée dans une conception d’échantillonnage aléatoire simple, stratifié ou systématique.
- Échantillonnage mixte (composite) : crée des localisations d’échantillonnage séparément à partir de différentes conceptions d’échantillonnage, puis les fusionne dans un même jeu de données. Par exemple, la combinaison d’un échantillon aléatoire simple et d’un échantillon d’agrégat en deux étapes génère des localisations d’échantillonnage dans l’ensemble de la zone d’étude (aléatoire simple), mais inclut également de petites parcelles avec davantage de points (agrégat en deux étapes). Cela est utile car l’échantillonnage aléatoire simple seul peut ne pas permettre de déterminer le mode d’interaction entre les échantillons à de courtes distances, tandis que l’échantillonnage par agrégat en deux étapes laisse des surfaces importantes de la zone d’étude sans localisation d’échantillonnage. En combinant les deux, vous pouvez vous assurer que l’intégralité de la zone d’étude est représentée, mais étudier quand même l’interaction entre les échantillons à de courtes distances.
Un avertissement est renvoyé si le nombre spécifié de localisations d’échantillonnage ne peuvent pas être créés. Cela peut survenir dans les cas suivants :
- La valeur du paramètre Minimum Distance Between Sample Points (Distance minimale entre les points d’échantillonnage) est suffisamment grande pour que le nombre spécifié de localisations d’échantillonnage ne puissent pas être créées dans la zone d’étude (ou la strate) sans que des points ne soient plus proches les uns des autres que la distance minimale admise. Dans ce cas, le nombre de localisations créées est inférieur à celui spécifié.
- Si la valeur du paramètre Bin Size (Taille du groupe) est spécifiée sous forme de nombre, il n’est pas toujours possible de créer le nombre spécifié de localisations d’échantillonnage dans la zone d’étude. L’outil essaye diverses valeurs de surface et utilise la surface qui crée un nombre d’échantillons le plus proche de la valeur spécifiée. La surface (dans l’unité du système de coordonnées en sortie) et le nombre résultant de localisations d’échantillonnage sont renvoyés comme messages de géotraitement.
Si les paramètres spécifiés ne créent pas de localisation d’échantillonnage (par exemple, si une étendue en sortie n’intersectant pas la zone d’étude est utilisée), une erreur est renvoyée.
Pour un échantillonnage systématique et par agrégat et toute forme de groupe à l’exception des hexaones H3, le centroïde du premier polygone de la tessellation est créé dans le coin inférieur gauche de l’étendue en sortie. Les hexagones H3 se trouvent à des localisations fixes. Pour toutes les formes de groupe, vous pouvez utiliser le paramètre Spatial Relationship (Relation spatiale) pour renvoyer les polygones qui intersectent la zone d’étude, ceux qui se trouvent complètement à l’intérieur ou ceux dont les centroïdes se trouvent dans la zone d’étude.
Si vous stratifiez par champ d’ID de strate et utilisez un champ de population (égal ou proportionnel), la population de chaque strate correspond à la somme des valeurs des champs de population de chacun des polygones ou catégories de raster de la strate.
Si vous stratifiez par région de raster contigu, vous ne pouvez pas utiliser de champ de population. En effet, chaque champ de population représente la population totale d’une catégorie de raster même si cette catégorie est composée de plusieurs régions disjointes. Pour utiliser des champs de population lors d’une stratification par région de raster contigu, utilisez l’outil Raster To Polygon (Raster vers polygones) pour convertir le raster en polygones et affecter des valeurs de population à chaque polygone (par exemple, en allouant la population de chaque catégorie proportionnellement au nombre de cellules dans chacune de ses régions).
Pour un échantillonnage stratifié avec un nombre d’échantillons de strate proportionnel à la surface ou à un champ de population, la méthode du plus grand reste est utilisée afin que le nombre global d’échantillons ne soit pas modifié par l’arrondi.
Paramètres
arcpy.management.CreateSpatialSamplingLocations(in_study_area, out_features, {sampling_method}, {strata_id_field}, {strata_count_method}, {bin_shape}, {bin_size}, {h3_resolution}, {num_samples}, {num_samples_per_strata}, {population_field}, {geometry_type}, {min_distance}, {spatial_relationship})| Nom | Explication | Type de données |
in_study_area | Zone d’étude en entrée dans laquelle les localisations d’échantillonnage sont créées. La zone d’étude doit correspondre à des polygones ou à un raster entier (catégoriel). Pour les rasters, les cellules dont les valeurs sont nulles ne sont pas incluses dans la zone d’étude. | Feature Layer; Raster Layer |
out_features | Entités en sortie représentant les localisations d’échantillonnage. Pour un échantillonnage aléatoire simple et l’échantillonnage stratifié, les entités en sortie sont des points. Pour un échantillonnage par agrégat, la sortie correspond à des polygones. Pour un échantillonnage systématique, la sortie peut correspondre à des points ou des polygones. | Feature Class |
sampling_method (Facultatif) | Spécifie la méthode d’échantillonnage utilisée pour créer les localisations d’échantillonnage.
| String |
strata_id_field (Facultatif) | Pour un échantillonnage stratifié par champ d’ID de strate, champ d’ID de strate définissant la strate. | Field |
strata_count_method (Facultatif) | Pour un échantillonnage stratifié, spécifie la méthode utilisée pour déterminer le nombre de localisations d’échantillonnage à créer dans chaque strate.
| String |
bin_shape (Facultatif) | Pour un échantillonnage systématique et par agrégat, spécifie la forme de chaque polygone de la tessellation maillée.
| String |
bin_size (Facultatif) | Pour un échantillonnage systématique et par agrégat, taille de chaque polygone de la tessellation. Cette valeur peut être spécifiée sous forme de nombre (nombre total de polygones tesselés créés dans la zone d’étude) ou de surface (surface de chaque polygone tesselé). Pour une entrée sous forme de nombre, la valeur par défaut est de 100. Pour une entrée sous forme de surface, une valeur doit être spécifiée. Si un nombre est spécifié, l’outil tente de créer le nombre spécifié de localisations d’échantillonnage. Si le nombre exact ne peut pas être créé, un avertissement est renvoyé. | Areal Unit; Long |
h3_resolution (Facultatif) | Pour un échantillonnage systématique ou par agrégat avec des groupes d’hexagones H3, spécifie la résolution H3 des hexagones. À chaque augmentation de la valeur de la résolution, la surface des polygones est divisée par sept.
| Long |
num_samples (Facultatif) | Nombre de localisations d’échantillonnage à créer. Ce paramètre s’applique toujours aux échantillonnages aléatoires simples et par agrégat. Pour un échantillonnage stratifié, ce paramètre s’applique si le nombre d’échantillons est proportionnel à la surface de la strate ou à un champ de population. Pour un échantillonnage aléatoire simple et un échantillonnage stratifié, la valeur par défaut est de 100. Pour un échantillonnage par agrégat, la valeur par défaut est de 10. | Long |
num_samples_per_strata (Facultatif) | Pour un échantillonnage stratifié avec un nombre d’échantillons égal dans chaque strate, nombre de localisations d’échantillonnage créées dans chaque strate. Le nombre total d’échantillons correspond à cette valeur multipliée par le nombre de strates. La valeur par défaut est 100. | Long |
population_field (Facultatif) | Champ de population d’un échantillonnage stratifié si le nombre d’échantillons est égal ou proportionnel à un champ de population. | Field |
geometry_type (Facultatif) | Pour un échantillonnage systématique, indique si les localisations d’échantillonnage correspondent à des polygones tesselés ou des centroïdes (points) de polygones tesselés.
| String |
min_distance (Facultatif) | Pour un échantillonnage aléatoire simple et un échantillonnage stratifié, plus petite distance autorisée entre les localisations d’échantillonnage. Pour un échantillonnage aléatoire simple, tous les points sont séparés d’au moins cette distance. Pour un échantillonnage stratifié, les points d’une même strate sont séparés d’au moins cette distance, mais les points dans les strates voisines peuvent être plus proches que cette distance. Pour les grandes distances, il est possible de créer moins de localisations d’échantillonnage que prévu pour que les localisations restent suffisamment distantes les unes des autres. Dans ce cas, un message d’avertissement est renvoyé. | Linear Unit |
spatial_relationship (Facultatif) | Spécifie les polygones d’une tessellation en arrière-plan à inclure comme localisations d’échantillonnage. Ce paramètre s’applique à l’échantillonnage par agrégat et à l’échantillonnage systématique si le type de géométrie en sortie est le polygone.
| String |
Exemple de code
Le script Python suivant illustre l’utilisation de la fonction CreateSpatialSamplingLocations.
# Create 50 sampling locations in the dissolved California counties.
import arcpy
arcpy.management.CreateSpatialSamplingLocations(
in_study_area="CA_counties",
out_features="outputSamplingLocations"
sampling_method="RANDOM",
strata_id_field=None,
strata_count_method="EQUAL",
bin_shape="HEXAGON",
bin_size=None,
h3_resolution=7,
num_samples=50,
num_samples_per_strata=100,
population_field=None,
geometry_type="POINT",
min_distance="15 NauticalMilesInt",
spatial_relationship = "HAVE_THEIR_CENTER_IN"
)Le script Python suivant illustre l’utilisation de la fonction CreateSpatialSamplingLocations.
# Simple random sampling
# Create 50 sample points in a polygon study area.
# Import system modules.
import arcpy
# Allow overwriting output.
arcpy.env.overwriteOutput = True
# Define study area and output features.
inputStudyArea = "C:/samplingdata/inputs.gdb/study_area_polygons"
outputFeatures = "C:/samplingdata/outputs.gdb/out_samples_SRS"
# Define the sampling method and number of samples.
samplingMethod = "RANDOM"
numSamples=50
# Define the minimum distance between any two points.
minDistance= "15 NauticalMilesInt"
# Run tool.
try:
arcpy.management.CreateSpatialSamplingLocations(inputStudyArea, outputFeatures,
samplingMethod, "", "", "", "", "", numSamples, "", "", "",
minDistance)
except arcpy.ExecuteError:
# If an error occurred when running the tool, print the error message.
print(arcpy.GetMessages())Le script Python suivant illustre l’utilisation de la fonction CreateSpatialSamplingLocations.
# Stratify by individual polygons
# Create 100 sample points in each polygon.
# Import system modules.
import arcpy
# Allow overwriting output.
arcpy.env.overwriteOutput = True
# Define the study area and output features.
inputStudyArea = "C:/samplingdata/inputs.gdb/study_area_polygons"
outputFeatures = "C:/samplingdata/outputs.gdb/out_samples_SBIP"
# Define the sampling method.
samplingMethod = "STRAT_POLY"
# Create 100 samples in each polygon.
strataCountMethod = "EQUAL"
numSamplesPerStrata=100
# Define the minimum distance between any two points in the same polygon.
minDistance= "15 Meters"
# Run tool.
try:
arcpy.management.CreateSpatialSamplingLocations(inputStudyArea, outputFeatures,
samplingMethod, "", strataCountMethod, "", "", "", "",
numSamplesPerStrata, "", "", minDistance)
except arcpy.ExecuteError:
# If an error occurred when running the tool, print the error message.
print(arcpy.GetMessages())Le script Python suivant illustre l’utilisation de la fonction CreateSpatialSamplingLocations.
# Stratify by contiguous raster region
# Create 100 points in a raster study area with number of samples in
# each region proportional to the area of the region.
# Import system modules.
import arcpy
# Allow overwriting output.
arcpy.env.overwriteOutput = True
# Define the study area and output features.
inputStudyArea = "C:/samplingdata/raster_study_area.tif"
outputFeatures = "C:/samplingdata/outputs.gdb/out_samples_SBCRR"
# Define the sampling method.
samplingMethod = "STRAT_RAST"
# Create 100 points and allocate proportionally to the area of the regions.
strataCountMethod = "PROP_AREA"
numSamples=100
# Run tool.
try:
arcpy.management.CreateSpatialSamplingLocations(inputStudyArea, outputFeatures,
samplingMethod, "", strataCountMethod, "", "", "", numSamples)
except arcpy.ExecuteError:
# If an error occurred when running the tool, print the error message.
print(arcpy.GetMessages())Le script Python suivant illustre l’utilisation de la fonction CreateSpatialSamplingLocations.
# Stratify by strata ID field
# Create sample points in each land use category of a raster.
# Use a population field to define the number of samples in each category.
# Import system modules.
import arcpy
# Allow overwriting output.
arcpy.env.overwriteOutput = True
# Define the study area and output features.
inputStudyArea = "C:/samplingdata/land_use_raster.tif"
outputFeatures = "C:/samplingdata/outputs.gdb/out_samples_SBSIDF"
# Define the sampling method.
samplingMethod = "STRAT_ID"
# All raster cells with the same value are in the same stratum.
strataIDField = "LandUse"
# Define the number of samples using a population field.
strataCountMethod = "FIELD"
populationField="Population"
# Run tool.
try:
arcpy.management.CreateSpatialSamplingLocations(inputStudyArea, outputFeatures,
samplingMethod, strataIDField, strataCountMethod, "", "", "",
"", "", populationField)
except arcpy.ExecuteError:
# If an error occurred when running the tool, print the error message.
print(arcpy.GetMessages())Le script Python suivant illustre l’utilisation de la fonction CreateSpatialSamplingLocations.
# Systematic sampling
# Create sample points in a hexagonal tessellation in a polygon study area.
# Import system modules.
import arcpy
# Allow overwriting output.
arcpy.env.overwriteOutput = True
# Define the study area and output features.
inputStudyArea = "C:/samplingdata/inputs.gdb/study_area_polygons"
outputFeatures = "C:/samplingdata/outputs.gdb/out_samples_SYS"
# Define the sampling method.
samplingMethod = "SYSTEMATIC"
# Create points in a hexagonal tessellation.
binShape = "HEXAGON"
binSize = "10000 SquareFeet"
outputGeometryType = "POINT"
# Run tool.
try:
arcpy.management.CreateSpatialSamplingLocations(inputStudyArea, outputFeatures,
samplingMethod, "", "", binShape, binSize, "", "", "", "",
outputGeometryType)
except arcpy.ExecuteError:
# If an error occurred when running the tool, print the error message.
print(arcpy.GetMessages())Le script Python suivant illustre l’utilisation de la fonction CreateSpatialSamplingLocations.
# Cluster sampling
# Create 100 cluster polygons that are diamond shaped.
# Import system modules.
import arcpy
# Allow overwriting output.
arcpy.env.overwriteOutput = True
# Define the study area and output features.
inputStudyArea = "C:/samplingdata/inputs.gdb/study_area_polygons"
outputFeatures = "C:/samplingdata/outputs.gdb/out_samples_CLUST"
# Define the sampling method.
samplingMethod = "CLUSTER"
# Create a diamond tessellation and randomly choose 100 polygons.
binShape = "DIAMOND"
binSize = "1000000 SquareFeet"
numSamples=100
spatialRelationship = "INTERSECT"
# Run tool.
try:
arcpy.management.CreateSpatialSamplingLocations(inputStudyArea, outputFeatures,
samplingMethod, "", "", binShape, binSize, "", numSamples, "",
"", "", "", spatialRelationship)
except arcpy.ExecuteError:
# If an error occurred when running the tool, print the error message.
print(arcpy.GetMessages())Environnements
Rubriques connexes
Vous avez un commentaire à formuler concernant cette rubrique ?