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

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

Краткая информация

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

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

Иллюстрация

Результаты 3D-интерполяции визуализированы в слое вокселов.

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

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

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

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

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

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

    Подсказка:

    Вы можете использовать инструмент Объекты в 3D по атрибуту для создания 3D-точек из 2D-точек на основании поля высот.

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

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

  Есть много информации об использовании поверхностей высот в качестве минимального или максимального растра высот. Сервисы изображений, веб-слои высот и веб-слои изображений будут медленнее обрабатываться и могут вызывать ошибки при большом количестве запросов. Растры, сохраненные в виде файлов на локальных дисках, обрабатываются быстрее и рекомендуются при создании слоев вокселов высокого разрешения на больших пространственных экстентах.

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

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

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

Параметры

Диалог

Подпись	Описание	Тип данных
Входные геостатистические 3D-слои	Геостатистические 3D-слои, которые будут экспортированы в Выходной файл netCDF.	Geostatistical Layer
Выходной файл netCDF	Выходной файл netCDF, содержащий экспортированные значения из входных геостатистических слоев. Результаты каждого геостатистического слоя сохраняются как разные переменные в файле netCDF.	File
Экспорт местоположений (Дополнительный)	Задает местоположения для экспорта данных, указанных в параметре Входные геостатистические 3D-слои. Если вы выбираете опцию 3D-точки на сетке, вам требуется предоставить значения для параметров Интервал X, Интервал Y и Интервал высот, которые представляют расстояние от каждой точки на сетке во всех измерениях. Если вы выбираете опцию Пользовательские 3D-точки, вы должны указать точечные 3D-объекты в параметре Местоположения 3D-точек, представляющие местоположения для экспорта. 3D-точки на сетке—Прогнозируемыми значениями являются 3D-точки на сетке. Это значение по умолчанию Пользовательские 3D-точки—Прогнозируемые значения указываются пользовательскими объектами 3D-точек.	String
Интервал X (Дополнительный)	Расстояние между каждой точкой сетки в x-измерении. Значение по умолчанию создает 40 точек вдоль выходного x-экстента.	Linear Unit
Интервал Y (Дополнительный)	Расстояние между каждой точкой сетки в y-измерении. Значение по умолчанию создает 40 точек вдоль выходного y-экстента.	Linear Unit
Интервал высоты (Дополнительный)	Расстояние между каждой точкой сетки в z-измерении. Значение по умолчанию создает 40 точек вдоль выходного z-экстента.	Linear Unit
Местоположения 3D-точек (Дополнительный)	Точечные 3D-объекты, представляющие местоположения для экспорта. Точечные объекты должны иметь данные высоты, которые хранятся в атрибуте геометрии Shape.Z.	Feature Layer
Выходные переменные (Дополнительный)	Задает выходные типы для Входных геостатистических 3D-слоев. Вы можете указать один или несколько выходных типов для каждого из слоев или применить выходной тип ко всем входным геостатистическим слоям. По умолчанию будут экспортироваться прогнозируемые значения для всех слоев. Чтобы экспортировать в другие типы выходных данных, укажите слой для экспорта (или выберите Все, чтобы взять все слои) в первой записи таблица значений. Укажите выходной тип во второй записи таблицы значений. Если вы выберете Вероятность или Квантиль в качестве выходного типа, укажите предельное значение (для вероятности) или значение квантиля (для квантиля) в третьей записи таблицы значений. Если вы выберете Прогнозирование или Стандартная ошибка прогнозирования в качестве выходного типа, вы можете не заполнять третью запись таблицы значений. Более подробно о выходных типах геостатистики	Value Table
Входные области изучаемой области (Дополнительный)	Полигональные объекты области изучения. В выходном файле netCDF сохраняются только точки, находящиеся в пределах изучаемой области. При визуализации в виде слоя вокселов в сцене будут отображаться только вокселы в пределах изучаемой области. Нахождение точек внутри или вне изучаемой области определяется только их координатами x и y.	Feature Layer
Минимальная высота вырезания растра (Дополнительный)	Растр высот будет использован для вырезания нижнего края слоя воксела. Прогнозированные значения получают только те вокселы, высоты которых выше растра высот. Например, если вы используете растр высот земной поверхности, слой вокселов будет отображаться только выше поверхности. Такой растр также можно использовать для обозначения уровня поверхности коренных пород или дна сланцевого месторождения. Растр должен быть в системе координат проекции, и значения высот должны быть в тех же единицах измерения, что и единицы растра по вертикали.	Raster Layer
Максимальная высота вырезания растра (Дополнительный)	Растр высот будет использован для вырезания верхнего края слоя воксела. Прогнозированные значения получают только те вокселы, высоты которых ниже растра высот. Например, если вы используете растр высот земной поверхности, слой вокселов будет отображаться только ниже поверхности. Его также можно использовать для отсечения вокселов в верхней части ограниченного воздушного пространства. Растр должен быть в системе координат проекции, и значения высот должны быть в тех же единицах измерения, что и единицы растра по вертикали.	Raster Layer

Производные выходные данные

Подпись	Описание	Тип данных
Выходной слой вокселов	Слой вокселов прогнозируемых значений.	Voxel Layer

Python

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

Имя	Описание	Тип данных
in_3d_geostat_layers [in_3d_geostat_layers,...]	Геостатистические 3D-слои, которые будут экспортированы в Выходной файл netCDF.	Geostatistical Layer
out_netcdf_file	Выходной файл netCDF, содержащий экспортированные значения из входных геостатистических слоев. Результаты каждого геостатистического слоя сохраняются как разные переменные в файле netCDF.	File
export_locations (Дополнительный)	Задает местоположения для экспорта данных, указанных в параметре in_3d_geostat_layers. Если вы выбираете опцию 3D_GRIDDED_POINTS, вам требуется предоставить значения для параметров x_spacing, u_spacing и elevation_spacing, которые представляют расстояние от каждой точки на сетке во всех измерениях. Если вы выбираете опцию CUSTOM_3D_POINTS, вы должны указать точечные 3D-объекты в параметре in_points_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
min_elev_raster (Дополнительный)	Растр высот будет использован для вырезания нижнего края слоя воксела. Прогнозированные значения получают только те вокселы, высоты которых выше растра высот. Например, если вы используете растр высот земной поверхности, слой вокселов будет отображаться только выше поверхности. Такой растр также можно использовать для обозначения уровня поверхности коренных пород или дна сланцевого месторождения. Растр должен быть в системе координат проекции, и значения высот должны быть в тех же единицах измерения, что и единицы растра по вертикали.	Raster Layer
max_elev_raster (Дополнительный)	Растр высот будет использован для вырезания верхнего края слоя воксела. Прогнозированные значения получают только те вокселы, высоты которых ниже растра высот. Например, если вы используете растр высот земной поверхности, слой вокселов будет отображаться только ниже поверхности. Его также можно использовать для отсечения вокселов в верхней части ограниченного воздушного пространства. Растр должен быть в системе координат проекции, и значения высот должны быть в тех же единицах измерения, что и единицы растра по вертикали.	Raster Layer

Производные выходные данные

Имя	Описание	Тип данных
out_voxel_layer	Слой вокселов прогнозируемых значений.	Voxel 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