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)