De capa 3D de GA a NetCDF (Geostatistical Analyst)

Disponible con una licencia de Geostatistical Analyst.

Resumen

Exporta una o varias capas de estadísticas geográficas 3D creadas con la herramienta Kriging bayesiano empírico 3D a un formato netCDF (archivo *.nc). El archivo de salida se muestra como una capa vóxel en una escena local.

Más información sobre las capas de vóxel

Ilustración

Ilustración de la herramienta Capa 3D de GA a NetCDF
La capa de vóxel muestra los resultados de una interpolación 3D.

Uso

  • Si la herramienta se ejecuta en una escena local con los mismos sistemas de coordenadas horizontales y verticales que las entidades de entrada, se agregará una capa vóxel a la escena que le permitirá explorar interactivamente los resultados. También puede agregar el archivo netCDF de salida como una capa vóxel utilizando la herramienta Crear capa vóxel multidimensional o el cuadro de diálogo Agregar capa vóxel multidimensional.

    Puede convertir el archivo netCDF de salida en un ráster multidimensional utilizando la herramienta Copiar ráster. También puede agregarlo a un mapa como una capa de entidades o ráster con la herramienta Crear capa de entidades NetCDF o con la herramienta Crear capa ráster NetCDF, respectivamente.

  • Todas las capas de estadísticas geográficas de entrada deben ser capas de estadísticas geográficas 3D creadas con Kriging bayesiano empírico 3D.

  • Puede elegir las ubicaciones para predecir y exportar de una de las siguientes formas:

    • Predecir a puntos de cuadrícula en 3D: las predicciones se exportarán a puntos de cuadrícula en 3D. Para usar esta opción, especifique puntos de cuadrícula en 3D para el parámetro Exportar ubicaciones. Los parámetros Espaciado X, Espaciado Y y Espaciado de elevación especifican la distancia entre cada punto de cuadrícula en las dimensiones x, y y z (elevación).
    • Predecir a ubicaciones 3D personalizadas: las predicciones se exportarán a ubicaciones 3D personalizadas definidas por entidades de puntos 3D. Para utilizar esta opción, especifique Puntos 3D personalizados para el parámetro Exportar ubicaciones y proporcione entidades de punto 3D al parámetro Ubicaciones de punto 3D para especificar las ubicaciones que se van a predecir.
      Sugerencia:

      Puede utilizar la herramienta Entidad a 3D por atributo para crear puntos 3D a partir de puntos 2D con un campo de elevación.

  • La extensión de salida predeterminada del archivo netCDF de salida es la combinación de las extensiones de todas las capas de estadísticas geográficas 3D de entrada.

  • Los parámetros Polígonos de área de estudio de entrada, Ráster de recorte de elevación mínima y Ráster de recorte de elevación máxima puede utilizarse para limitar el análisis dentro de un área de estudio específica y entre dos superficies de elevación. Cualquier vóxel fuera de estos límites no tendrá valor y no se mostrará. Por ejemplo, si los puntos se encuentran dentro de una reserva marina, puede crear una capa vóxel que se muestre solo dentro de un polígono de la reserva (área de estudio), por encima del fondo oceánico (ráster de elevación mínima) y por debajo de la termoclina (ráster de elevación máxima).

    Hay varias consideraciones para utilizar las superficies de elevación como rásteres de elevación mínima o máxima. Los servicios de imágenes, las capas de elevación web y las capas de imágenes web tendrán el rendimiento más lento y pueden producirse errores para un gran número de consultas. Los rásteres guardados como archivos locales en disco tendrán el rendimiento más rápido y se recomiendan cuando se crean capas vóxel de alta resolución en grandes extensiones espaciales.

  • Cuando se establece la cuadrícula 3D de puntos que representarán los vóxeles, el primer punto se crea en las coordenadas x, y y z mínimas de la extensión de salida (por defecto, la extensión de las entidades de entrada). Los puntos restantes se crean al iterar las distancias de los parámetros Espaciado X, Espaciado Y y Espaciado de elevación a través de las dimensiones de la extensión de salida. Si alguna de las distancias de espaciado no divide uniformemente la dimensión correspondiente de la extensión de salida, se creará una fila o columna de puntos más allá de la extensión de salida. Por ejemplo, si la extensión de salida para x se especifica como 0 a 10 y el parámetro Espaciado X se especifica como 3, la salida tendrá cinco filas en la extensión x: 0, 3, 6, 9 y 12. Del mismo modo, se creará una fila o columna de puntos adicional si las distancias de espaciado no dividen uniformemente las extensiones y o z.

  • Los rásteres de recorte de elevación mínima y máxima deben estar en un sistema de coordenadas proyectadas. Si los puntos o rásteres tienen un sistema de coordenadas geográficas con coordenadas de latitud y longitud, deben proyectarse a un sistema de coordenadas proyectadas mediante la herramienta Proyecto o Proyectar ráster.

  • Si se proporcionan polígonos de área de estudio de entrada, la extensión del área de estudio se utilizará como la extensión de salida predeterminada y los valores de los parámetros Espaciado X y Espaciado Y se recalcularán según esta extensión. Así se garantiza que la salida rellene toda el área de estudio de forma predeterminada.

Parámetros

EtiquetaExplicaciónTipo de datos
Capas de estadísticas geográficas 3D de entrada

Las capas de estadísticas geográficas 3D que se exportarán a un archivo netCDF de salida.

Geostatistical Layer
Archivo netCDF de salida

El archivo netCDF de salida que contiene los valores exportados de las capas de estadísticas geográficas de entrada. Los resultados de cada capa de estadísticas geográficas se guardan como variables diferentes en el archivo netCDF.

File
Exportar ubicaciones
(Opcional)

Especifica las ubicaciones que se exportarán desde el valor del parámetro Capas de estadísticas geográficas 3D de entrada. Si elige la opción puntos de cuadrícula en 3D, debe proporcionar valores para los parámetros Espaciado X, Espaciado Y y Espaciado de elevación que representen la distancia entre cada punto de cuadrícula de todas las dimensiones. Si elige la opción Puntos 3D personalizados, debe proporcionar entidades de puntos 3D en el parámetro Ubicaciones de puntos 3D que representen las ubicaciones que se van a exportar.

  • Puntos de cuadrícula en 3DLas ubicaciones de predicción son puntos de cuadrícula en 3D. Esta es la opción predeterminada.
  • Puntos 3D personalizadosLas ubicaciones de predicción se definen mediante entidades de puntos 3D personalizadas.
String
Espaciado X
(Opcional)

El espaciado entre cada punto de cuadrícula de la dimensión x. El valor predeterminado crea 40 puntos a lo largo de la extensión x de salida.

Linear Unit
Espaciado Y
(Opcional)

El espaciado entre cada punto de cuadrícula de la dimensión y. El valor predeterminado crea 40 puntos a lo largo de la extensión y de salida.

Linear Unit
Espaciado de elevación
(Opcional)

El espaciado entre cada punto de cuadrícula de la dimensión (z) de elevación. El valor predeterminado crea 40 puntos a lo largo de la extensión z de salida.

Linear Unit
Ubicaciones de puntos 3D
(Opcional)

Las entidades de puntos 3D que representan las ubicaciones que se van a exportar. Las entidades de puntos deben tener sus elevaciones almacenadas en el atributo de geometría Shape.Z.

Feature Layer
Variables de salida
(Opcional)

Especifica los tipos de salida para los valores de las Capas de estadísticas geográficas 3D de entrada. Puede especificar uno o más tipos de salida para cada una de las capas o puede aplicar un tipo de salida a todas las capas de estadísticas geográficas de entrada. De forma predeterminada, se exportarán las predicciones de todas las capas.

Para exportar otros tipos de salida, especifique la capa a exportar (o elija Todas para especificar todas las capas) en la primera entrada de la tabla de valores. Especifique el tipo de salida en la segunda entrada de la tabla de valores. Si elige Probabilidad o Cuantil como tipo de salida, especifique el valor de umbral (para probabilidad) o el valor de cuantil (para cuantil) en la tercera entrada de la tabla de valores. Si elige Predicción o Error estándar de predicción como tipo de salida, puede dejar la tercera entrada de la tabla de valores vacía.

Más información sobre los tipos de salida de estadísticas geográficas

Value Table
Polígonos de área de estudio de entrada
(Opcional)

Las entidades poligonales que representan el área de estudio. Solo los puntos que se encuentran dentro del área de estudio se guardan en el archivo netCDF de salida. Cuando se visualizan como una capa de vóxel, solo se mostrarán en la escena los vóxeles dentro del área de estudio. Se determina si los puntos están dentro o fuera del área de estudio usando solo sus coordenadas x e y.

Feature Layer
Ráster de recorte de elevación mínima
(Opcional)

El ráster de elevación que se utilizará para recortar la parte inferior de la capa vóxel. Solo se asignarán predicciones a los vóxeles que estén por encima de este ráster de elevación. Por ejemplo, si se utiliza un ráster de elevación del terreno, la capa vóxel solo se mostrará por encima de la superficie. También puede utilizarse para las superficies de roca firme o el fondo de un depósito de lutita.

El ráster debe estar en un sistema de coordenadas proyectadas y los valores de elevación deben estar en la misma unidad que la unidad vertical del ráster.

Raster Layer
Ráster de recorte de elevación máxima
(Opcional)

El ráster de elevación que se utilizará para recortar la parte superior de la capa vóxel. Solo se asignarán predicciones a los vóxeles que estén por debajo de este ráster de elevación. Por ejemplo, si se utiliza un ráster de elevación del terreno, la capa vóxel solo se mostrará por debajo de la superficie. También puede utilizarse para recortar vóxeles en la parte superior de un espacio aéreo restringido.

El ráster debe estar en un sistema de coordenadas proyectadas y los valores de elevación deben estar en la misma unidad que la unidad vertical del ráster.

Raster Layer

Salida derivada

EtiquetaExplicaciónTipo de datos
Capa de vóxel de salida

Una capa vóxel de los valores predichos.

Voxel Layer

arcpy.ga.GALayer3DToNetCDF(in_3d_geostat_layers, out_netcdf_file, {export_locations}, {x_spacing}, {y_spacing}, {elevation_spacing}, {in_points_3d}, {output_variables}, {in_study_area}, {min_elev_raster}, {max_elev_raster})
NombreExplicaciónTipo de datos
in_3d_geostat_layers
[in_3d_geostat_layers,...]

Las capas de estadísticas geográficas 3D que se exportarán a un archivo netCDF de salida.

Geostatistical Layer
out_netcdf_file

El archivo netCDF de salida que contiene los valores exportados de las capas de estadísticas geográficas de entrada. Los resultados de cada capa de estadísticas geográficas se guardan como variables diferentes en el archivo netCDF.

File
export_locations
(Opcional)

Especifica las ubicaciones que se van a exportar desde el valor del parámetro in_3d_geostat_layers. Si elige esta opción 3D_GRIDDED_POINTS, deberá proporcionar valores para los parámetros x_spacing, u_spacing y elevation_spacing que representan la distancia entre cada punto de la cuadrícula en todas las dimensiones. Si la opción CUSTOM_3D_POINTS, debe proporcionar entidades de puntos 3D en el parámetro in_points_3d que representen las ubicaciones que se van a exportar.

  • 3D_GRIDDED_POINTSLas ubicaciones de predicción son puntos de cuadrícula en 3D. Esta es la opción predeterminada.
  • CUSTOM_3D_POINTSLas ubicaciones de predicción se definen mediante entidades de puntos 3D personalizadas.
String
x_spacing
(Opcional)

El espaciado entre cada punto de cuadrícula de la dimensión x. El valor predeterminado crea 40 puntos a lo largo de la extensión x de salida.

Linear Unit
y_spacing
(Opcional)

El espaciado entre cada punto de cuadrícula de la dimensión y. El valor predeterminado crea 40 puntos a lo largo de la extensión y de salida.

Linear Unit
elevation_spacing
(Opcional)

El espaciado entre cada punto de cuadrícula de la dimensión (z) de elevación. El valor predeterminado crea 40 puntos a lo largo de la extensión z de salida.

Linear Unit
in_points_3d
(Opcional)

Las entidades de puntos 3D que representan las ubicaciones que se van a exportar. Las entidades de puntos deben tener sus elevaciones almacenadas en el atributo de geometría Shape.Z.

Feature Layer
output_variables
[[layer_name, output_type, quantile_probability_value],...]
(Opcional)

Especifica los tipos de salida para los valores de las Capas de estadísticas geográficas 3D de entrada. Puede especificar uno o más tipos de salida para cada una de las capas o puede aplicar un tipo de salida a todas las capas de estadísticas geográficas de entrada. De forma predeterminada, se exportarán las predicciones de todas las capas.

Para exportar otros tipos de salida, especifique la capa a exportar (o elija Todas para especificar todas las capas) en la primera entrada de la tabla de valores. Especifique el tipo de salida en la segunda entrada de la tabla de valores. Si elige Probabilidad o Cuantil como tipo de salida, especifique el valor de umbral (para probabilidad) o el valor de cuantil (para cuantil) en la tercera entrada de la tabla de valores. Si elige Predicción o Error estándar de predicción como tipo de salida, puede dejar la tercera entrada de la tabla de valores vacía.

Más información sobre los tipos de salida de estadísticas geográficas

Value Table
in_study_area
(Opcional)

Las entidades poligonales que representan el área de estudio. Solo los puntos que se encuentran dentro del área de estudio se guardan en el archivo netCDF de salida. Cuando se visualizan como una capa de vóxel, solo se mostrarán en la escena los vóxeles dentro del área de estudio. Se determina si los puntos están dentro o fuera del área de estudio usando solo sus coordenadas x e y.

Feature Layer
min_elev_raster
(Opcional)

El ráster de elevación que se utilizará para recortar la parte inferior de la capa vóxel. Solo se asignarán predicciones a los vóxeles que estén por encima de este ráster de elevación. Por ejemplo, si se utiliza un ráster de elevación del terreno, la capa vóxel solo se mostrará por encima de la superficie. También puede utilizarse para las superficies de roca firme o el fondo de un depósito de lutita.

El ráster debe estar en un sistema de coordenadas proyectadas y los valores de elevación deben estar en la misma unidad que la unidad vertical del ráster.

Raster Layer
max_elev_raster
(Opcional)

El ráster de elevación que se utilizará para recortar la parte superior de la capa vóxel. Solo se asignarán predicciones a los vóxeles que estén por debajo de este ráster de elevación. Por ejemplo, si se utiliza un ráster de elevación del terreno, la capa vóxel solo se mostrará por debajo de la superficie. También puede utilizarse para recortar vóxeles en la parte superior de un espacio aéreo restringido.

El ráster debe estar en un sistema de coordenadas proyectadas y los valores de elevación deben estar en la misma unidad que la unidad vertical del ráster.

Raster Layer

Salida derivada

NombreExplicaciónTipo de datos
out_voxel_layer

Una capa vóxel de los valores predichos.

Voxel Layer

Muestra de código

Ejemplo 1 de GALayer3DToNetCDF (ventana de Python)

Interpola puntos 3D dos veces y convierte las salidas a un archivo netCDF multivariante.

import arcpy
arcpy.ga.EmpiricalBayesianKriging3D("my3DPoints1", "Shape.Z",
                                    "myValueField1", "my3DGALayer1")
arcpy.ga.EmpiricalBayesianKriging3D("my3DPoints2", "Shape.Z",
                                    "myValueField2", "my3DGALayer2")
arcpy.ga.GALayer3DToNetCDF("my3DGALayer1;my3DGALayer2", "outputNCDF.nc",
                           "3D_GRIDDED_POINTS", "50 Meters", "50 Meters", "5 Meters",
                           "", "<ALL> PREDICTION #")
Ejemplo 2 de GALayer3DToNetCDF (script independiente)

Interpola puntos 3D dos veces y convierte las salidas a un archivo netCDF multivariante. Exporte a puntos 3D personalizados y de cuadrícula.

# Name: GALayer3DToNetCDF_Example_02.py
# Description: Interpolates 3D points and exports to a netCDF file.
# Requirements: Geostatistical Analyst Extension
# Author: Esri

# Import system modules
import arcpy

# Allow overwriting output
arcpy.env.overwriteOutput = True

# Set up variables
in3DPoints1 = "C:/gapydata/inputs.gdb/my3DPoints1"
in3DPoints2 = "C:/gapydata/inputs.gdb/my3DPoints2"
elevationField1 = "Shape.Z"
elevationField2 = "Shape.Z"
valueField1 = "myValueField1"
valueField2 = "myValueField2"
outGALayer1 = "myGALayer1"
outGALayer2 = "myGALayer2"


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

# Execute Empirical Bayesian Kriging 3D twice
arcpy.ga.EmpiricalBayesianKriging3D(in3DPoints1, elevationField1, valueField1, outGALayer1)
arcpy.ga.EmpiricalBayesianKriging3D(in3DPoints2, elevationField2, valueField2, outGALayer2)


# Export predictions for first model and probability that second model exceeds 10
# Export to gridded 3D points

# Set up variables
in_3d_ga_layers = outGALayer1+";"+outGALayer2
out_ncdf = "C:/gapydata/outputs/outputNetCDF1.nc"
export_locations = "3D_GRIDDED_POINTS"
x_spacing = "50 Meters"
y_spacing = "50 Meters"
elev_spacing = "5 Meters"
custom_points = ""
out_vars = "myGALayer1 PREDICTION #;myGALayer2 PROBABILITY 10"

# Run tool.
arcpy.ga.GALayer3DToNetCDF(in_3d_ga_layers, out_ncdf, export_locations,
                           x_spacing, y_spacing, elev_spacing, custom_points, out_vars)


# Export standard errors for first model and 75th quantile for second model
# Export to custom 3D points

# Set up variables
in_3d_ga_layers = outGALayer1+";"+outGALayer2
out_ncdf = "C:/gapydata/outputs/outputNetCDF2.nc"
export_locations = "CUSTOM_3D_POINTS"
x_spacing = ""
y_spacing = ""
elev_spacing = ""
custom_points = "C:/gapydata/inputs.gdb/myCustom3DPoints"
out_vars = "myGALayer1 PREDICTION_STANDARD_ERROR #;myGALayer2 QUANTILE 0.75"

# Run tool.
arcpy.ga.GALayer3DToNetCDF(in_3d_ga_layers, out_ncdf, export_locations,
                           x_spacing, y_spacing, elev_spacing, custom_points, out_vars)

Información de licenciamiento

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

Temas relacionados