मैंने सीखना शुरू कर दिया है कि अजगर में एलएएस डेटा में हेरफेर कैसे करना है और यह देखना चाहता है कि अन्य लोग एलएएस फ़ाइलों को कैसे संभालते हैं। मैं अंक पढ़ना चाहता हूं (मैं एक संख्यात्मक सरणी का उपयोग कर रहा हूं), और एक अलग सरणी के लिए कक्षा 1 और 2 (अवर्गीकृत और जमीन) को फ़िल्टर करता हूं। मेरे पास निम्नलिखित कोड है, लेकिन अंक फ़िल्टर किए गए प्रतीत नहीं हो सकते।
# 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
मैंने arcpy.da.featureClassToNumpyArray देखा है, लेकिन मैं आर्कपी आयात नहीं करना चाहता था और न ही इसे आकार में बदलना है।
मैं कैसे LAS डेटा को एक सुन्न सरणी में फ़िल्टर / पढ़ सकता हूं?