Статистика точек LAS как растра (LAS Point Statistics As Raster) (Управление данными)

Уровень лицензии:BasicStandardAdvanced

Резюме

Создает растр, растра, значения ячеек которого отражают о статистическую информацию о размерах из файлов LAS, на которые ссылается набор данных LAS.

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

Синтаксис

LasPointStatsAsRaster_management (in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})
ПараметрОбъяснениеТип данных
in_las_dataset

Входной набор данных LAS.

LAS Dataset Layer
out_raster

Местоположение и имя выходного растра. При хранении набора растровых данных в базе геоданных или в папке, например в Esri Grid, расширение файла к имени набора растровых данных добавлять не нужно. При хранении набора растровых данных в папке для определения формата растра может быть задано расширение файла:

  • .bil – Esri BIL
  • .bip – Esri BIP
  • .bsq – Esri BSQ
  • .dat – ENVI DAT
  • .img – ERDAS IMAGINE
  • .png – PNG
  • .tif – TIFF

Если растр хранится в базе геоданных или в виде файла TIFF, тип и качество сжатия этого растра можно задать с помощью настроек среды геообработки.

Raster Dataset
method
(дополнительно)

Тип собираемых статистических данных по точкам LAS в каждой из ячеек выходного растра.

  • PULSE_COUNT (Число импульсов)Количество последних возвращаемых точек.
  • POINT_COUNTКоличество точек во всех возвращаемых данных.
  • PREDOMINANT_LAST_RETURNНаиболее частое последнее возвращаемое значение.
  • PREDOMINANT_CLASS (Преобладающий класс)Наиболее часто встречаемый код классификации.
  • INTENSITY_RANGE (Диапазон интенсивности)Диапазон значений интенсивности.
  • Z_RANGE (Z диапазон)Диапазон значений рельефа.
String
sampling_type
(дополнительно)

Указывает метод, используемый для интерпретации Значение выборки (Sampling Value) для определения разрешения выходного растра.

  • OBSERVATIONSОпределяет число ячеек, разделяющих самую длинную сторону экстента набора данных LAS.
  • CELLSIZEОпределяет размер ячейки выходного набора растровых данных. Это значение используется по умолчанию.
String
sampling_value
(дополнительно)

Задает значение, используемое вместе с параметром Значение выборки (Sampling Value) для определения разрешения выходного растра.

Double

Пример кода

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

В следующем примере показано использование этого инструмента в окне Python:

import arcpy
from arcpy import env

env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img", 
                             "INTENSITY_RANGE", "CELLSIZE", 15)
LasPointStatsAsRaster — пример 2 (автономный скрипт)

В следующем примере показано использование этого инструмента в автономном скрипте Python:

'''**********************************************************************
Name: LAS Point Statistics As Raster
Description: Identifies the most frequently occurring return value for 
             each pulse in a given set of LAS files. 
             Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy
import exceptions, sys, traceback

# Set Local Variables
lasD = arcpy.GetParameterAsText(0)
inLas = arcpy.GetParameterAsText(1) #input las files
sr = arcpy.GetParameter(2) #spatial reference of las dataset
statsRaster = arcpy.GetParameterAsText(3)

try:
    arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
    # Execute LasPointStatsAsRaster
    arcpy.management.LasPointStatsAsRaster(lasD, statsRaster, 
                                           "PREDOMINANT_RETURNS_PER_PULSE", 
                                           "CELLSIZE", 15)
    arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
    print arcpy.GetMessages()
except:
    # Get the traceback object
    tb = sys.exc_info()[2]
    tbinfo = traceback.format_tb(tb)[0]
    # Concatenate error information into message string
    pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\
          .format(tbinfo, str(sys.exc_info()[1]))
    msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2))
    # Return python error messages for script tool or Python Window
    arcpy.AddError(pymsg)
    arcpy.AddError(msgs)

Параметры среды

Связанные темы

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

ArcGIS for Desktop Basic: Требует 3D Analyst или Spatial Analyst
ArcGIS for Desktop Standard: Да
ArcGIS for Desktop Advanced: Да
9/10/2013