def write_xyz_file(file_name):
points = read_raw_las_data(file_name)
pc = create_vtkpc_from_array(points)
writer = vtk.vtkSimplePointsWriter()
writer.SetFileName(file_name[:-4] + ".xyz")
writer.SetInputData(pc.vtkPolyData)
print("Writing to %s" % file_name[:-4] + ".xyz")
writer.Write()
def main():
input1 = "Data/MantecaFiltered/MantecaRailSlice.las"
with File(input1, mode='r') as f:
input_header = f.header
to_plot = [] # list of VTK Point Clouds to plot
points = subsample_frac_from_las_data(input1, .01)
pc = create_vtkpc_from_array(points)
to_plot.append(pc)
create_point_cloud_plot_qt(to_plot,
input_header=input_header,
axes_on=True)
print("Finding histogram took %d seconds" % end)
# del point_set
print(max(H))
return max(H)
with File(input1, mode='r') as f:
input_header = f.header
to_plot = []
# points = subsample_frac_from_las_data(input1, .0625)
# pc = create_vtkpc_from_array(points)
# to_plot.append(pc)
points2 = subsample_frac_from_las_data(input2, .0625)
pc2 = create_vtkpc_from_array(points2)
to_plot.append(pc2)
theta = cp.Variable(1)
objective = cp.Maximize(max_bin(points2, theta))
constraints = [0 <= theta, theta <= np.pi / 2.]
prob = cp.Problem(objective, constraints)
result = prob.solve()
# print(max_bin(points2, theta))
print(theta.value)
pc3 = create_vtkpc_from_array(points2)
to_plot.append(pc3)
# pdb.set_trace()
create_point_cloud_plot_qt(to_plot,
# pc2 = create_vtkpc_from_array(points2)
# to_plot.append(pc2)
# points = read_raw_las_data(input1)
# warehouse = subsample_frac_from_las_data(input2, .1)
#
# # points = subsample_frac_from_las_data(input1, .01)
# pc = create_vtkpc_from_array(points+warehouse)
# to_plot.append(pc)
points = naive_slice_from_las(input1,
AxisAlignedBox3D([7, -13, -.5], [9, -11, 3]))
# points = subsample_frac_from_las_data(input1, .01)
points = ann_guided_filter(points, num_neighbors=50, filter_eps=.07)
# points = subsample_frac_from_las_data(input1, .1)
pc = create_vtkpc_from_array(points)
to_plot.append(pc)
# points2 = threshold_filter(points, threshold=.01)
# pc2 = create_vtkpc_from_array(points2)
# to_plot.append(pc2)
# points2 = read_raw_las_data(input2)
# pc2 = create_vtkpc_from_array(points2)
# to_plot.append(pc2)
# points2 = threshold_filter(points, threshold=.015)
# pc2 = create_vtkpc_from_array(points2)
# to_plot.append(pc2)
# points3 = rounding_filter(points2, decimal_places=3)