He empezado a aprender cómo manipular LAS datos en python y quería ver cómo los demás manejar archivos LAS. Me gustaría leer los puntos (estoy usando una colección de la matriz), y filtrar las clases 1 y 2 (sin clasificar y tierra) a una matriz independiente. Tengo el siguiente código, pero no parecen llegar a los puntos de filtrado.
# Import modules
from liblas import file
import numpy as np
if __name__=="__main__":
'''Read LAS file and create an array to hold X, Y, Z values'''
# Get file
las_file = r"E:\Testing\ground_filtered.las"
# Read file
f = file.File(las_file, mode='r')
# Get number of points from header
num_points = int(f.__len__())
# Create empty numpy array
PointsXYZIC = np.empty(shape=(num_points, 5))
# Load all LAS points into numpy array
counter = 0
for p in f:
newrow = [p.x, p.y, p.z, p.intensity, p.classification]
PointsXYZIC[counter] = newrow
counter += 1
He visto arcpy.da.featureClassToNumpyArray, pero yo no quería importar arcpy ni tiene que convertir a shapefile. Cualquier sugerencia sobre el filtrado o la manera más fácil leer LAS datos en una colección de la matriz son bienvenidos.