Recorrido de una partícula (Spatial Analyst)

Disponible con una licencia de Spatial Analyst.

Resumen

Calcula la ruta de una partícula mediante un campo de velocidad y devuelve un archivo ASCII de los datos del recorrido de la partícula y, opcionalmente, una clase de entidad de información del recorrido.

Más información sobre cómo funciona Recorrido de una partícula

Uso

  • Los rásteres de dirección y magnitud de entrada deben proceder de la misma ejecución de la herramienta Flujo Darcy.

  • El archivo de ruta generado por esta herramienta es un archivo de texto ASCII que contiene información sobre la posición, la magnitud y la dirección de velocidad local y la longitud y el tiempo de viaje acumulativos a lo largo de la ruta. Dispersión tipo puff usa este archivo como entrada. El formato de este archivo es el siguiente:

    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
  • Recorrido de una partícula no especifica ningún sistema de unidades concreto. Es importante que todos los datos estén en un conjunto de unidades constante, con la misma unidad para el tiempo (segundos, días, años) y la longitud (pies, metros).

  • La ubicación del origen debe estar entre los límites de los rásteres de entrada y no puede estar en un área de NoData.

  • El archivo de recorrido finalizará si el recorrido llega fuera del área de estudio y no ha cumplido con el tiempo de recorrido máximo especificado.

  • Si la partícula que está realizando el recorrido ha alcanzado el borde del área de estudio en el tiempo indicado y el punto del indicador está fuera del área de estudio, el archivo de recorrido finalizará.

  • Si la partícula que está realizando el recorrido migra a una depresión en el momento indicado, el archivo de recorrido finalizará. Se puede crear una depresión mediante un pozo de descarga u otro sumidero.

  • Las dos salidas de esta herramienta son:

    1. Un archivo ASCII de recorrido de partícula que utiliza el nombre especificado como archivo de recorrido de partícula de salida
    2. Una clase de entidad de polilínea opcional

  • Consulte Entornos de análisis y Spatial Analyst para obtener detalles adicionales sobre los entornos de geoprocesamiento que se aplican a esta herramienta.

Sintaxis

ParticleTrack(in_direction_raster, in_magnitude_raster, source_point, out_track_file, {step_length}, {tracking_time}, {out_track_polyline_features})
ParámetroExplicaciónTipo de datos
in_direction_raster

Un ráster de entrada donde cada valor de celda representa la dirección del vector de velocidad de filtrado (velocidad lineal promedio) en el centro de la celda.

Las direcciones se expresan en coordenadas de brújula, en grados en el sentido de las agujas del reloj desde el norte. Esto se puede crear con la herramienta Flujo Darcy.

Los valores de dirección deben ser puntos flotantes.

Raster Layer
in_magnitude_raster

Un ráster de entrada donde cada valor de celda representa la magnitud del vector de velocidad de filtrado (velocidad lineal promedio) en el centro de la celda.

Las unidades son de longitud/ tiempo. Esto se puede crear con la herramienta Flujo Darcy.

Raster Layer
source_point

Un objeto de la clase Point de Python indica la ubicación del punto de origen, en unidades de mapa, desde la que debe comenzar el recorrido de la partícula.

El formato del objeto es:

  • point(x,y)
Point
out_track_file

El archivo de texto ASCII de salida que contiene los datos del recorrido de la partícula.

File
step_length
(Opcional)

La longitud de paso que se utilizará para calcular el recorrido de la partícula.

La opción predeterminada es la mitad del tamaño de celda. Las unidades son de longitud.

Double
tracking_time
(Opcional)

Tiempo máximo transcurrido para el recorrido de la partícula.

El algoritmo seguirá el recorrido hasta que se cumpla este tiempo o hasta que la partícula migre fuera del ráster o a una depresión.

El valor predeterminado es el infinito. Las unidades son de tiempo.

Double
out_track_polyline_features
(Opcional)

La clase de entidad de línea de salida opcional que contiene el recorrido de la partícula.

Esta clase de entidad contiene una serie de arcos con atributos para la posición, la magnitud y la dirección de velocidad local y la longitud y el tiempo de viaje acumulativos a lo largo de la ruta.

Feature Class

Muestra de código

Ejemplo 1 de ParticleTrack (ventana de Python)

En este ejemplo se ejecuta la herramienta en las entradas necesarias y genera un archivo ASCII de datos de rastreo de partículas y una clase de entidad de shapefile del rastro de partículas.

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")
Ejemplo 2 de ParticleTrack (secuencia de comandos independiente)

En este ejemplo se ejecuta la herramienta en las entradas necesarias y genera un archivo ASCII de datos de rastreo de partículas y una clase de entidad de shapefile del rastro de partículas.

# 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"

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

Información de licenciamiento

  • Basic: Requiere Spatial Analyst
  • Standard: Requiere Spatial Analyst
  • Advanced: Requiere Spatial Analyst

Temas relacionados