print("num_returns shape",num_returns.shape)
print("First number of returns:",num_returns[0])
unique_elements, counts_elements = np.unique(num_returns, return_counts=True)
print("Frequency of unique values of the said array:")
print(np.asarray((unique_elements, counts_elements)))
src_ids = inFile.get_pt_src_id()
print("First number of returns:",num_returns[0])
unique_elements, counts_elements = np.unique(src_ids, return_counts=True)
print("Frequency of", "source ids", "of the said array:")
print(np.asarray((unique_elements, counts_elements)))
#get_scan_dir_flag()
flight_lines = inFile.get_edge_flight_line()
raw_classes = inFile.get_raw_classification()
classes = inFile.get_classification()
#scanner_channel = inFile.get_scanner_channel()
synthetics = inFile.get_synthetic()
key_points = inFile.get_key_point()
withehelds_points = inFile.get_withheld()
try:
overlap = inFile.get_overlap()
pointformat = inFile.point_format
for spec in inFile.point_format:
print(spec.name)
print("==========")
headerformat = inFile.header.header_format
for spec in headerformat:
print(spec.name)
print(inFile.header.offset)
print(inFile.header.scale)
print("==========")
print(len(inFile.points))
print(inFile.points[0])
print(inFile.X[0])
print(inFile.Y[0])
print(inFile.Z[0])
print(inFile.intensity[0])
print(inFile.flag_byte[0])
print(inFile.raw_classification[0])
print(inFile.user_data[0])
print(inFile.get_edge_flight_line()[0])
print(inFile.pt_src_id[0])
print(inFile.gps_time[0])
print("==========")
print(inFile.header.min)
longitude = inFile.X[0] * inFile.header.scale[0] + inFile.header.offset[0]
latitude = inFile.Y[0] * inFile.header.scale[1] + inFile.header.offset[1]
waterLevel = inFile.Z[0] * inFile.header.scale[2] + inFile.header.offset[2]
print(twd97.fromwgs84(latitude, longitude), waterLevel)
