ArcMap

Suivi de particules

Disponible avec une licence Spatial Analyst.

Résumé

Calcule le chemin d'une particule via un champ de vitesse et retourne un fichier ASCII de données de suivi de particules et, le cas échéant, une classe d'entités d'informations de suivi.

Pour en savoir plus sur le suivi de particules

Utilisation

  • Les rasters de grandeur et de direction en entrée doivent provenir de la même exécution de l'outil Flux de Darcy.

  • Le fichier de chemin généré par cet outil est un fichier texte ASCII contenant des informations relatives à la position, à la direction et à l'amplitude de la vitesse locale, ainsi qu'à la longueur et au temps cumulés du déplacement le long du chemin. Ce fichier est utilisé comme entrée par Concentration. Le format de ce fichier est le suivant :

    time         x            y            length       flow dir     flow mag
0.000000000  0.000000000  482.8400000  0.000000000  90.00000000  0.04418909563
113.1648712  4.999804443  482.7957786  5.000000000  91.01366126  0.04418332249
226.2741353  9.998043277  482.6630814  10.00000000  92.02765240  0.04420504404
339.3574334  14.99315255  482.4419855  15.00000000  93.04094157  0.04421519432
452.3447720  19.98356700  482.1325285  20.00000000  94.05521317  0.04425274599
565.2657591  24.96772671  481.7348453  25.00000000  95.06807622  0.04427874865
678.0514031  29.94406931  481.2490323  30.00000000  96.08254679  0.04433188322
790.7309576  34.91104149  480.6752838  35.00000000  97.09488082  0.04437362239

  • Aucun système particulier d'unités n'est spécifié par Suivi de particules. Il est important que toutes les données soient dans un ensemble cohérent d'unités, utilisant les mêmes unités de temps (secondes, jours, années) et de longueur (pieds, mètres).

  • L'emplacement source doit être dans la limite des rasters en entrée et ne peut pas être dans une zone de valeur NoData.

  • Le fichier de suivi se termine si la piste atteint l'extérieur de la zone d'étude et n'a pas rencontré le temps de suivi maximal spécifié.

  • Si la particule faisant l'objet d'un suivi a atteint la limite de la zone d'étude au moment indiqué et que le point de prédiction est à l'extérieur de la zone d'étude, le fichier de suivi se termine.

  • Si la particule faisant l'objet d'un suivi migre vers une dépression au moment indiqué, le fichier de suivi se termine. Une dépression peut être créée par un puits de décharge ou une autre cuvette.

  • Les deux sorties de cet outil sont :

    1. un fichier ASCII de suivi de particules utilisant le nom spécifié en tant que fichier de suivi de particules en sortie ;
    2. une classe d'entités polylignes facultative.

  • Pour plus d’informations sur les environnements de géotraitement qui s’appliquent à cet outil, reportez-vous à la rubrique Environnements d’analyse et Spatial Analyst.

Syntaxe

ParticleTrack(in_direction_raster, in_magnitude_raster, source_point, out_track_file, {step_length}, {tracking_time}, {out_track_polyline_features})
ParamètreExplicationType de données
in_direction_raster

Raster en entrée dans lequel chaque valeur de cellule représente la direction du vecteur de vitesse d'infiltration (vitesse linéaire moyenne) au centre de la cellule.

Les directions sont exprimées dans les coordonnées du compas, en degrés dans le sens horaire à partir du Nord. Elles peuvent être créées par l'outil Flux de Darcy.

Les valeurs de direction doivent être des valeurs à virgules flottantes.

Raster Layer
in_magnitude_raster

Raster en entrée dans lequel chaque valeur de cellule représente la magnitude du vecteur de vitesse d'infiltration (vitesse linéaire moyenne) au centre de la cellule.

Les unités correspondent à une longueur/temps. Elles peuvent être créées par l'outil Flux de Darcy.

Raster Layer
source_point

Un objet de classe Python Point présente l’emplacement du point source, en unités cartographiques, à partir duquel commencer le suivi des particules.

L'objet a la forme suivante :

  • point(x,y)
Point
out_track_file

Fichier texte ASCII en sortie qui contient les données de suivi de particules.

File
step_length
(Facultatif)

Longueur de pas à utiliser pour le calcul du suivi des particules.

La valeur par défaut correspond à la moitié de la taille des cellules. Les unités correspondent à une longueur.

Double
tracking_time
(Facultatif)

Temps écoulé maximal pour le suivi des particules.

L'algorithme suit la piste tant que ce temps n'est pas écoulé ou que la particule ne sort pas du raster ou n'entre pas dans une dépression.

La valeur par défaut correspond à l'infini. Les unités correspondent à un temps.

Double
out_track_polyline_features
(Facultatif)

Classe d'entités linéaires en sortie facultative qui contient le suivi de particules.

Cette classe d'entités contient une série d'arcs avec des attributs pour la position, la direction et l'amplitude de la vitesse locale, ainsi que la longueur cumulée et la durée de transport le long du chemin.

Feature Class

Exemple de code

Exemple 1 d'utilisation de l'outil ParticleTrack (fenêtre Python)

Cet exemple exécute l'outil sur les entrées requises et génère un fichier ASCII de données de suivi de particules et une classe d'entités de fichiers de formes du suivi de particules.

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
ParticleTrack("gwdir", "gwmag", arcpy.Point(-200,-200), 
              "C:/sapyexamples/output/trackfile.txt",10, 100000, 
              "C:/sapyexamples/output/trackpolyline.shp")
Exemple 2 d'utilisation de l'outil ParticleTrack (script autonome)

Cet exemple exécute l'outil sur les entrées requises et génère un fichier ASCII de données de suivi de particules et une classe d'entités de fichiers de formes du suivi de particules.

# Name: ParticleTrack_Ex_02.py
# Description: Calculates the path of a particle through a velocity field.
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "C:/sapyexamples/data"

# Set local variables
inDirectionRaster = "gwdir"
inMagnitudeRaster = "gwmag"
sourcePoint = arcpy.Point(-200, -200)
outTrackFile = "C:/sapyexamples/output/trackfile.txt"
stepLength = 10
trackingTime = 10000000
outTrackPolylineFeatures = "C:/sapyexamples/output/trackpolyline.shp"

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Execute ParticleTrack
ParticleTrack(inDirectionRaster, inMagnitudeRaster, sourcePoint, outTrackFile,
              stepLength, trackingTime, outTrackPolylineFeatures)

Environnements

Informations de licence

  • Basic: Requiert Spatial Analyst
  • Standard: Requiert Spatial Analyst
  • Advanced: Requiert Spatial Analyst

Rubriques connexes