Sensibilidad de semivariograma (Geostatistical Analyst)

Disponible con una licencia de Geostatistical Analyst.

Resumen

Esta herramienta realiza un análisis de sensibilidad sobre los valores predichos y los errores estándar asociados cambiando los parámetros del semivariograma del modelo (los rangos nugget, meseta parcial y rangos principal/secundario) dentro de un porcentaje de los valores originales.

La herramienta toma una fuente de modelo geoestadístico para poblar estos valores iniciales de los rangos nugget, meseta parcial y mayor/menor. La salida de la herramienta es una tabla que indica qué valores de parámetro se han utilizado y cuáles eran los valores de error estándar y predichos resultantes. Si hay grandes fluctuaciones en la salida con pequeños cambios en los valores de parámetro del modelo, no se puede confiar mucho en el resultado. Por otro lado, si los cambios en el resultado son pequeños, puede confiar en las predicciones del modelo y tomar decisiones basadas en su salida.

Uso

  • La fuente del modelo de estadísticas geográficas es una capa de estadísticas geográficas o un modelo de estadísticas geográficas (XML).

  • En la tabla de salida, PtID corresponde al valor OID de la clase de entidad Ubicaciones de observación de puntos de entrada y el campo RndParam identifica cuál de los parámetros del semivariograma se cambió aleatoriamente.

  • Esta herramienta utiliza un generador de números aleatorios en sus operaciones. El Valor de inicialización utilizado se puede controlar en el entorno Generador de números aleatorios.

    • Si se utiliza un valor de inicialización de 0 (el valor predeterminado), cada vez que se ejecute la herramienta se utilizará un conjunto diferente de números aleatorios y se generará una salida diferente.
    • Si el valor de inicialización de número aleatorio se establece en un número fijo mayor que 0, la herramienta generará la misma salida cada vez que se ejecute, hasta que se cambie el valor de inicialización.

    Nota:

    Solo se admite el tipo de generador de números aleatorios Mersenne Twister; si se elige Algoritmo 599 capturado por ACM o Rand C estándar, se utilizará Mersenne Twister en su lugar.

  • En la secuencia de comandos Python, la clase GeostatisticalDatasets de ArcPy será útil para rellenar el parámetro Dataset de entrada (s ).

  • Para los formatos de datos que admiten valores nulos, como clases de entidad de geodatabase de archivos, se utilizará un valor nulo para indicar que no se pudo hacer una predicción para esa ubicación o que el valor se debe ignorar cuando se usa como entrada. Para los formatos de datos que no admiten valores nulos, como los shapefiles, se utiliza el valor -1.7976931348623158e+308 (es el valor negativo de la constante de C++ DBL_MAX definida) para indicar que no se pudo realizar una predicción para esa ubicación.

Parámetros

EtiquetaExplicaciónTipo de datos
Fuente del modelo de estadísticas geográficas de entrada

La fuente del modelo de estadísticas geográficas que se analizará.

File; Geostatistical Layer
Datasets de entrada

El nombre de los datasets de entrada y los nombres de campo que se utilizan en la creación de la capa de salida.

Geostatistical Value Table
Ubicaciones de observación de puntos de entrada

Ubicaciones de puntos en los que se realiza el análisis de sensibilidad.

Feature Layer
Alcance de Nugget (% del valor del modelo)
(Opcional)

Porcentaje restado y sumado al parámetro Nugget para crear un rango para una próxima selección aleatoria de parámetros Nugget.

Double
Número de cálculos para Nugget
(Opcional)

Número de valores aleatorios de Nugget muestreados aleatoriamente del intervalo Nugget.

Long
Alcance de Meseta parcial (% del valor del modelo)
(Opcional)

Porcentaje restado y sumado al parámetro Meseta parcial para crear un rango para una selección aleatoria de Meseta parcial.

Double
Número de cálculos para Meseta parcial
(Opcional)

Número de valores de Meseta parcial muestreadas aleatoriamente del intervalo Meseta parcial.

Long
Alcance de Rango mayor (% del valor del modelo)
(Opcional)

Porcentaje restado y sumado al parámetro Rango principal para crear un rango para una selección aleatoria de Rango principal.

Double
Número de cálculos para Rango mayor
(Opcional)

Número de valores de Rango principal muestreados aleatoriamente del intervalo Rango principal.

Long
Alcance de Rango menor (% del valor del modelo)
(Opcional)

Porcentaje restado y sumado al parámetro Rango secundario para crear un rango para una selección aleatoria de Rango secundario.

Double
Número de cálculos para Rango menor
(Opcional)

Número de valores para Rango menor muestreados aleatoriamente del intervalo Rango secundario.

Si se ha configurado Anisotropy en la fuente del modelo geoestadístico de entrada, se requiere un valor.

Long
Tabla de salida

Tabla que almacena los resultados de sensibilidad.

Table

arcpy.ga.GASemivariogramSensitivity(in_ga_model_source, in_datasets, in_locations, {nugget_span_percents}, {nugget_calc_times}, {partialsill_span_percents}, {partialsill_calc_times}, {range_span_percents}, {range_calc_times}, {minrange_span_percents}, {minrange_calc_times}, out_table)
NombreExplicaciónTipo de datos
in_ga_model_source

La fuente del modelo de estadísticas geográficas que se analizará.

File; Geostatistical Layer
in_datasets

Un objeto GeostatisticalDatasets.

También puede ser una cadena de elementos delimitada por punto y coma. Cada elemento consta de los siguientes componentes:

  • La ruta de catálogo y el nombre de un dataset o el nombre de una capa en la tabla de contenido actual, seguido de un espacio.
  • Una secuencia de nombres de campo, cada nombre de campo separado por un espacio. En el caso de un ráster, se utilizarán los valores de celda.
Geostatistical Value Table
in_locations

Ubicaciones de puntos en los que se realiza el análisis de sensibilidad.

Feature Layer
nugget_span_percents
(Opcional)

Porcentaje restado y sumado al parámetro Nugget para crear un rango para una próxima selección aleatoria de parámetros Nugget.

Double
nugget_calc_times
(Opcional)

Número de valores aleatorios de Nugget muestreados aleatoriamente del intervalo Nugget.

Long
partialsill_span_percents
(Opcional)

Porcentaje restado y sumado al parámetro Meseta parcial para crear un rango para una selección aleatoria de Meseta parcial.

Double
partialsill_calc_times
(Opcional)

Número de valores de Meseta parcial muestreadas aleatoriamente del intervalo Meseta parcial.

Long
range_span_percents
(Opcional)

Porcentaje restado y sumado al parámetro Rango principal para crear un rango para una selección aleatoria de Rango principal.

Double
range_calc_times
(Opcional)

Número de valores de Rango principal muestreados aleatoriamente del intervalo Rango principal.

Long
minrange_span_percents
(Opcional)

Porcentaje restado y sumado al parámetro Rango secundario para crear un rango para una selección aleatoria de Rango secundario.

Double
minrange_calc_times
(Opcional)

Número de valores para Rango menor muestreados aleatoriamente del intervalo Rango secundario.

Si se ha configurado Anisotropy en la fuente del modelo geoestadístico de entrada, se requiere un valor.

Long
out_table

Tabla que almacena los resultados de sensibilidad.

Table

Muestra de código

Ejemplo SemivariogramSensitivity 1 (ventana de Python)

Realizar un análisis de sensibilidad variando los valores de Nugget, Meseta parcial y Rango.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GASemivariogramSensitivity_ga(
    "C:/gapyexamples/data/kriging.lyr", "C:/gapyexamples/data/ca_ozone_pts.shp OZONE",
    "C:/gapyexamples/data/obs_pts.shp", "", "", "", "", "", "", "", "", 
    "C:/gapyexamples/output/outtabSS")
Ejemplo de SemivariogramSensitivity 2 (script independiente)

Realizar un análisis de sensibilidad variando los valores de Nugget, Meseta parcial y Rango.

# Name: SemivariogramSensitivity_Example_02.py
# Description: The semivariogram parameters Nugget, Partial Sill and Range can
#   be varied to perform a sensitivity analysis.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

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

# Set local variables
inLayer = "C:/gapyexamples/data/kriging.lyr"
inData = "C:/gapyexamples/data/ca_ozone_pts.shp OZONE"
inObs = "C:/gapyexamples/data/obs_pts.shp"
nugPercents = ""
nugCalc = ""
sillPercents = ""
sillCalc = ""
rangePercents = ""
rangeClac = ""
minrangePercent = ""
midrangeCalc = ""
outTable = "C:/gapyexamples/output/outtabSS"

# Execute SemivariogramSensitivity
arcpy.GASemivariogramSensitivity_ga(inLayer, inData, inObs, nugPercents, nugCalc,
                                    sillPercents, sillCalc, rangePercents, rangeClac,
                                    minrangePercent, midrangeCalc, outTable)

Información de licenciamiento

  • Basic: Requiere Geostatistical Analyst
  • Standard: Requiere Geostatistical Analyst
  • Advanced: Requiere Geostatistical Analyst

Temas relacionados