3D-Слой GA в NetCDF (Geostatistical Analyst)

Доступно с лицензией Geostatistical Analyst.

Сводка

Экспортирует один или несколько 3D геостатистических слоев, созданных с помощью инструмента Эмпирический байесовский кригинг 3D, в формат netCDF (файл *.nc). Основная цель этого инструмента – подготовить 3D геостатистические слои для визуализации в виде слоя вокселов в локальной сцене.

Более подробно о слоев вокселов

Иллюстрация

3D геостатистический слой экспортируется в слой вокселов
Слой вокселов отображает результаты 3D-интерполяции.

Использование

  • Параметр Выходной файл netCDF не будет автоматически добавлять слои на карту или сцену. Чтобы отобразить выходные данные в виде слоя вокселов, они должны быть добавлены в локальную сцену с той же горизонтальной и вертикальной системой координат, что и файл netCDF. Он также может быть добавлен на карту или сцену в качестве векторного или растрового слоя с помощью инструмента Создать векторный слой NetCDF или Создать растровый слой NetCDF соответственно.

  • Чтобы отобразить Выходной файл netCDF как слой вокселов, укажите для параметра Экспорт местоположений значение 3D-точки на сетке. Для визуализации слоя вокселов требуются местоположения на сетке.

  • Все Входные 3D геостатистические слои должны быть 3D геостатистическими слоями, созданными при помощи Эмпирического байесовского кригинга 3D.

  • Вы можете выбрать местоположения для прогнозирования и экспорта одним из следующих способов:

    • Прогнозирование в 3D-точек на сетке – прогнозированные значения будут экспортированы в 3D-точки на сетке. Чтобы использовать эту опцию, укажите 3D-точки на сетке для параметра Экспорт местоположений. Параметры Интервал X, Интервал Y и Интервал высоты задают расстояние между каждой точкой на сетке в x, y и z (высота) измерениях.
    • Прогнозирование для пользовательских 3D-местоположений – прогнозированные значения будут экспортированы в пользовательские 3D-местоположения, определенные 3D-точечными объектами. Чтобы использовать эту опцию, укажите Пользовательские 3D-точки для параметра Экспорт местоположений, и добавьте 3D-точечные объекты в параметр Местоположения 3D-точек, чтобы указать местоположения для экспорта. Местоположения 3D-точек должны иметь данные высоты, которые хранятся в атрибуте геометрии Shape.Z. Если высоты хранятся в другом поле атрибутов, не в Shape.Z, вы можете преобразовать их в 3D-точки с помощью инструмента Объекты в 3D по атрибуту.

  • Выходной экстент по умолчанию для параметра Выходной файл netCDF объединяет экстенты всех Входных геостатистических 3D-слоев.

  • Если выбрана опция 3D-точки на сетке для параметра Экспорт местоположений, первая точка создается с минимальной x-, минимальной y- и минимальной z-координатой для параметра среды Выходной экстент. Оставшиеся точки создаются с использованием итерации расстояний параметров Интервал X, Интервал Y и Интервал высоты в размерах выходного экстента. Если какое-либо из интервалов расстояний не делит равномерно соответствующий размер выходного экстента, то за пределами выходного экстента будет создана одна строка или столбец точек. Например, если выходной экстент для x указан как значения от 0 к 10, а параметр Интервал X задан 3, то в выходных данных будет 4 строки в экстенте x: 0, 3, 6, 9 и 12. Аналогично, будет создана дополнительная строка или столбец точек, если расстояния интервалов между ними не будут равномерно разделять y- или z-экстенты.

  • Параметр Полигоны входной изучаемой области используется для ограничения результатов интерполяции конкретной изучаемой областью вместо показа полного куба. В выходном файле netCDF только для местоположений, расположенных в пределах изучаемой области, будут созданы точки.

  • Если заданы Входные полигоны изучаемой области, в качестве выходного экстента по умолчанию будет использоваться экстент изучаемой области, а значения параметров Интервал X, Интервал Y и Интервал высоты будут пересчитаны на основе этого экстента. Это позволит быть уверенными в том, что выходные данные по умолчанию заполнят всю изучаемую область.

Синтаксис

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})
ParameterОбъяснениеТип данных
in_3d_geostat_layers
[in_3d_geostat_layers,...]

3D геостатистические слои, которые будут экспортированы в Выходной файл netCDF. Если предоставлено более одного слоя, то в результате получится многомерный файл netCDF.

Geostatistical Layer
out_netcdf_file

Выходной файл netCDF, содержащий экспортированные значения из Входных геостатистических 3D-слоев.

File
export_locations
(Дополнительный)

Задает местоположения для экспорта из Входных геостатистических 3D-слоев. Вы можете экспортировать в 3D-точки на сетке или указать пользовательские точечные 3D-объекты для представления местоположений для экспорта. Если вы выбираете 3D-точки на сетке, вам требуется предоставить значения для параметров Интервал X, Интервал Y и Интервал высоты, которые представляют расстояние от каждой точки на сетке во всех измерениях. Если вы выбираете Пользовательские 3D-точки, вы должны указать точечные 3D-объекты в параметре Местоположения 3D-точек, представляющие местоположения для экспорта.

  • 3D_GRIDDED_POINTSПрогнозируемыми значениями являются 3D-точки на сетке. Это значение по умолчанию.
  • CUSTOM_3D_POINTSПрогнозируемые значения указываются пользовательскими объектами 3D-точек.
String
x_spacing
(Дополнительный)

Расстояние между каждой точкой сетки в x-измерении. Значение по умолчанию создает 40 точек вдоль выходного x-экстента.

Linear Unit
y_spacing
(Дополнительный)

Расстояние между каждой точкой сетки в y-измерении. Значение по умолчанию создает 40 точек вдоль выходного y-экстента.

Linear Unit
elevation_spacing
(Дополнительный)

Расстояние между каждой точкой сетки в z-измерении. Значение по умолчанию создает 40 точек вдоль выходного z-экстента.

Linear Unit
in_points_3d
(Дополнительный)

Точечные 3D-объекты, представляющие местоположения для экспорта. Точечные объекты должны иметь данные высоты, которые хранятся в атрибуте геометрии Shape.Z.

Feature Layer
output_variables
[[layer_name, output_type, quantile_probability_value],...]
(Дополнительный)

Задает выходные типы для Входных геостатистических 3D-слоев. Вы можете указать один или несколько выходных типов для каждого из слоев или применить выходной тип ко всем входным геостатистическим слоям. По умолчанию будут экспортироваться прогнозируемые значения для всех слоев.

Чтобы экспортировать в другие типы выходных данных, укажите слой для экспорта (или выберите Все, чтобы взять все слои) в первой записи таблица значений. Укажите выходной тип во второй записи таблицы значений. Если вы выберете Вероятность или Квантиль в качестве выходного типа, укажите предельное значение (для вероятности) или значение квантиля (для квантиля) в третьей записи таблицы значений. Если вы выберете Прогнозирование или Стандартная ошибка прогнозирования в качестве выходного типа, вы можете не заполнять третью запись таблицы значений.

Более подробно о выходных типах геостатистики

Value Table
in_study_area
(Дополнительный)

Полигональные объекты области изучения. В выходном файле netCDF сохраняются только точки, находящиеся в пределах изучаемой области. При визуализации в виде слоя вокселов в сцене будут отображаться только вокселы в пределах изучаемой области. Нахождение точек внутри или вне изучаемой области определяется только их координатами x и y.

Feature Layer

Пример кода

GALayer3DToNetCDF пример 1 (окно Python)

Дважды интерполирует 3D-точки и конвертирует выходные данные в многомерный файл netCDF.

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 #")
GALayer3DToNetCDF, пример 2 (автономный скрипт)

Дважды интерполирует 3D-точки и конвертирует выходные данные в многомерный файл netCDF. Экспортирует пользовательские 3D-точки и 3D-точки сетки.

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

Информация о лицензиях

  • Basic: Требуется Geostatistical Analyst
  • Standard: Требуется Geostatistical Analyst
  • Advanced: Требуется Geostatistical Analyst

Связанные разделы