import laspy

# Cargar el archivo .laz
file_path = "/home/data13/LiDAR_3a_cobertura/santander_laz/santander_copc/santander.copc.laz"
las = laspy.read(file_path)

# Obtener los valores mínimos y máximos de las coordenadas (bounding box)
min_x, max_x = las.x.min(), las.x.max()
min_y, max_y = las.y.min(), las.y.max()
min_z, max_z = las.z.min(), las.z.max()

# Imprimir el bounding box
print(f"Bounding box: \n"
      f"X: [{min_x}, {max_x}] \n"
      f"Y: [{min_y}, {max_y}] \n"
      f"Z: [{min_z}, {max_z}]")
x_med = (min_x+max_x)/2
y_med = (min_y+max_y)/2
z_med = (min_z+max_z)/2
print(x_med)
print(y_med)
print(z_med)