print('Constructing the dataframe...')
#Construct the consolidated processed dataframe
filt_sig_ldr = pd.DataFrame(columns=[
'timestamp', 'datetime', 'rlng_accl', 'quant_sig', 'trnry_sig', 'binry_sig'
])
filt_sig_ldr['timestamp'] = ldr_timestamp
filt_sig_ldr['datetime'] = ldr_datetime
filt_sig_ldr['rlng_accl'] = ldr_cordev
filt_sig_ldr['quant_sig'] = ldr_quant
filt_sig_ldr['trnry_sig'] = ldr_trnry
filt_sig_ldr['binry_sig'] = ldr_binry
#Detect PCBs and obtain PCB processing durations
pcb_data_ldr = routine_get_pcb_dur.get_pcb_dur(filt_sig_ldr.binry_sig)
#Eliminate PCB loading events
filt_sig_ldr['cor_binry_sig'] = routine_elim_ldng_pcb.elim_ldng_pcb(
ldr_binry, pcb_data_ldr, 1000)
#Recalculate PCB processing durations
pcb_data_ldr = routine_get_pcb_dur.get_pcb_dur(filt_sig_ldr.cor_binry_sig)
#Obtain local modes for PCB active duration, and also display vital statistics about machine
THRESHOLD = 200000
HISTBIN_WINDOW = 50
ldr_mode_df = routine_get_histogram.get_mode_df(pcb_data_ldr, THRESHOLD,
HISTBIN_WINDOW)
print('Mode calculated at :', datetime.datetime.now())
pp1_vibr_file = pp1_vibr_file.reset_index()
pp1_vibr_file = pp1_vibr_file.iloc[:, 1:len(pp1_vibr_file.columns)]
pp1_vibr_file = pp1_vibr_file.assign(
net_accl=np.sqrt(pp1_vibr_file.data_ax**2 + pp1_vibr_file.data_ay**2 +
pp1_vibr_file.data_az**2))
print('Size of pp1_pxstr_file =', pp1_pxstr_file.shape)
print('Size of pp1_vibr_file =', pp1_vibr_file.shape)
#-------------------------------------------------------------------------
#Raw-acceleration to binary sequene conversion
filt_sig_pp1 = routine_rawaccl_to_binryseq.rawaccl_to_binryseq(
pp1_vibr_file, 1, 'Y', 'Y', 0.6, 1.1)
#Detect PCBs and obtain PCB processing durations
pcb_data_pp1 = routine_get_pcb_dur.get_pcb_dur(filt_sig_pp1.binry_sig)
#Eliminate spuriois PCB detections
filt_sig_pp1['cor_binry_sig'] = routine_elim_spur_pcb.elim_spur_pcb(
filt_sig_pp1.binry_sig, pcb_data_pp1, 500)
#Recalculate PCB processing durations
pcb_data_pp1 = routine_get_pcb_dur.get_pcb_dur(filt_sig_pp1.cor_binry_sig)
#Obtain PCB level dataframe with weightages
pcb_level_pp1 = routine_get_pcb_level_data.get_pcb_level_data(
filt_sig_pp1, pcb_data_pp1, 50)
#-----------------------------------------------------------------------------------------------------------------------
#Merging process
pcb_merge_pp1 = pcb_level_pp1