Estadísticas de dataset LAS (Administración de datos)
Resumen
Calcula o actualiza estadísticas para un dataset LAS y genera un informe estadístico opcional.
Uso
Las estadísticas calculan el espaciado de punto promedio para cada archivo LAS al usar un método de colocación en bin que evalúa en aumento las áreas pequeñas del archivo para determinar un cálculo.
Las estadísticas habilitan las opciones de filtrado para una capa de dataset LAS para visualizar automáticamente los códigos de clase y devolver valores que se encuentran en los archivos LAS a los que hace referencia. Las opciones de filtrado se puede especificar mediante el cuadro de diálogo Propiedades de capa en ArcMap y ArcScene.
-
La siguiente información acerca de los archivos LAS se almacenará en el dataset LAS y se puede revisar por medio del cuadro de diálogo Propiedades del dataset de LAS en ArcCatalog:
- Rango de valores de elevación e intensidad
- Número de puntos de datos por código de clase y valor devuelto
- Existencia de valores RGB
- Versión LAS e ID del formato de punto
- ID del proyecto e ID de la fuente de archivo
- Nombre del software usado para generar el archivo LIDAR
El archivo de informe estadístico opcional proporciona una vista general de las propiedades de todos los archivos LAS en el dataset LAS o cada archivo LAS en un formato de texto que se pueden importar en cualquier número de aplicaciones.
Sintaxis
| Parámetro | Explicación | Tipo de datos |
in_las_dataset |
El dataset de LAS de entrada. | LAS Dataset Layer |
calculation_type | Especifica si las estadísticas se calcularán para todos los archivos LIDAR o solamente para aquellos que no tienen estadísticas:
| Boolean |
out_file (Opcional) |
El archivo de texto de salida que incluirá el resumen de las estadísticas de dataset LAS. | Text File |
summary_level (Opcional) |
Especifica el tipo de resumen incluido en el informe.
| String |
delimiter (Opcional) |
El delimitador de campo que se utiliza en el archivo de texto.
| String |
decimal_separator (Opcional) | El carácter decimal usado en el archivo de texto para diferenciar al entero de un número de su parte fraccionaria.
| String |
Ejemplo de código
El siguiente ejemplo muestra cómo usar de esta herramienta en la ventana Python:
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.LASDatasetStatistics_3d("test.lasd", "NO_FORCE", "LAS_FILE",
"DECIMAL_POINT", "SPACE", "LAS_summary.txt")
El siguiente ejemplo muestra cómo usar esta herramienta en una secuencia de comandos independiente de Python:
'''*********************************************************************
Name: Modify Files in LAS Dataset& Calculate Stats for LASD
Description: Adds files & surface constraints to a LAS dataset, then
calculates statistics and generates report.
*********************************************************************'''
# Import system modules
import arcpy
from arcpy import env
import exceptions, sys, traceback
try:
# Script variables
env.workspace = 'C:/data'
lasd = 'sample.lasd'
oldLas = ['2006', '2007/file2.las']
newLas = ['2007_updates_1', '2007_updates_2']
oldSurfaceConstraints = ['boundary.shp', 'streams.shp']
newSurfaceConstraints = [['sample.gdb/boundary', '<None>',
'Soft_Clip']
['sample.gdb/streams', 'Shape.Z',
'Hard_Line']]
arcpy.management.RemoveFilesFromLasDataset(lasd, oldLas,
oldSurfaceConstraints)
arcpy.management.AddFilesToLasDataset(lasd, newLas, 'RECURSION',
newSurfaceConstraints)
arcpy.management.LasDatasetStatistics(lasd, "UPDATED_FILES",
"lasd_stats.txt",
"LAS_FILE", "DECIMAL_POINT",
"SPACE", "LAS_summary.txt")
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)