def _get_file_info(file_paths):
"""
Getting file info from file paths
:param file_paths: list of file paths
:return:
"""
if len(file_paths) == 0:
return None
files_data = dict()
for path in file_paths:
file_name = os.path.basename(path)
try:
name, extension = file_name.split('.')
except ValueError:
continue
if extension in ('00', 'xx', 'bin'):
bin_data = BinaryFile()
bin_data.path = path
file_info = dict()
file_info['path'] = path
file_info['file_type'] = bin_data.device_type
file_info['frequency'] = bin_data.signal_frequency
file_info['datetime_start'] = bin_data.datetime_start
file_info['datetime_stop'] = bin_data.datetime_stop
files_data[file_name] = file_info
if len(files_data.keys()) == 0:
return None
else:
return files_data
def run(self):
params = self.parameters
for index_i, item_i in enumerate(params):
export_name, file_stat = item_i
export_folder = os.path.dirname(file_stat[0][0])
if export_name[-3:] != '.00':
export_name += '.00'
output_file = open(os.path.join(export_folder, export_name), 'wb')
for index_j, item_j in enumerate(file_stat):
file_path, dt_start, dt_stop = item_j
bin_data = BinaryFile()
bin_data.path = file_path
bin_data.use_avg_values = False
if index_j == 0:
output_file.write(bin_data.main_header.get_binary_format())
for k in range(3):
output_file.write(
bin_data.get_channel_header(k).get_binary_format())
# channel_headers_hole = np.array([0] * 54, dtype=np.int32)
# channel_headers_hole.tofile(output_file)
signal = bin_data.signals
signal.astype(np.int32).tofile(output_file)
self.finished_part_stitching.emit(
(index_j + 1) / len(file_stat) * 100)
self.finished_file.emit((index_i + 1) / len(params) * 100)
def select_files(self):
self._ui.teFilesList.clear()
file_dialog = QFileDialog()
file_dialog.setFileMode(QFileDialog.ExistingFiles)
names = file_dialog.getOpenFileNames()[0]
folder_path = None
min_dt = None
max_dt = None
max_frequency = None
for path in names:
if folder_path is None:
folder_path = os.path.dirname(path)
self._ui.leInputFolder.setText(folder_path)
file_name = os.path.basename(path)
t = file_name.split('.')
if len(t) != 2:
continue
name, extension = t
if extension not in ('00', 'xx'):
continue
bin_data = BinaryFile()
bin_data.path = path
bin_data.record_type = 'XYZ'
is_correct, _ = bin_data.check_correct()
if not is_correct:
continue
if min_dt is None:
min_dt = bin_data.datetime_start
max_dt = bin_data.datetime_stop
max_frequency = bin_data.signal_frequency
else:
min_dt = max(min_dt, bin_data.datetime_start)
max_dt = min(max_dt, bin_data.datetime_stop)
max_frequency = max(max_frequency, bin_data.signal_frequency)
self._ui.teFilesList.append(file_name)
if min_dt is None:
return None
if min_dt >= max_dt:
self._ui.teFilesList.clear()
msg = QMessageBox()
msg.setIcon(QMessageBox.Information)
msg.setText("Отсутствует перекрытие по времени")
msg.exec_()
else:
self._ui.dtMinDt.setMinimumDateTime(min_dt)
self._ui.dtMinDt.setMaximumDateTime(max_dt)
self._ui.dtMaxDt.setMinimumDateTime(min_dt)
self._ui.dtMaxDt.setMaximumDateTime(max_dt)
self._ui.dtMinDt.setDateTime(min_dt)
self._ui.dtMaxDt.setDateTime(min_dt + datetime.timedelta(hours=1))
self._ui.sbResampleFrequency.setMinimum(0)
self._ui.sbResampleFrequency.setMaximum(max_frequency)
self._ui.sbResampleFrequency.setValue(max_frequency)
def collect_parameters(self):
result = list()
grid_data = self.ui.tGrid
for i in range(grid_data.rowCount()):
if grid_data.item(i, 0) is None:
continue
export_name = grid_data.item(i, 0).text()
folder = grid_data.item(i, 1).text()
file_statistics = list()
freq = None
for index, file in enumerate(os.listdir(folder)):
name_parse = file.split('.')
extension = name_parse[-1]
if extension != '00part':
continue
path = os.path.join(folder, file)
bin_data = BinaryFile()
bin_data.path = path
if freq is not None:
if bin_data.signal_frequency != freq:
continue
else:
freq = bin_data.signal_frequency
file_statistics.append(
(path, bin_data.datetime_start, bin_data.datetime_stop))
file_statistics = sorted(file_statistics, key=lambda x: x[1])
is_right = True
for j in range(len(file_statistics) - 1):
current_t_stop = file_statistics[j][2]
next_t_start = file_statistics[j + 1][1]
dt = (next_t_start - current_t_stop).total_seconds()
discrete_amount = int(dt * freq)
if discrete_amount > 1:
is_right = False
break
if is_right:
result.append((export_name, file_statistics))
color = QtGui.QColor(0, 207, 0)
else:
color = QtGui.QColor(248, 207, 0)
for j in range(grid_data.columnCount()):
grid_data.item(i, j).setBackground(color)
return result
def load_file(self):
path = self.ui.leFilePath.text()
if len(path) == 0:
return
bin_data = BinaryFile()
bin_data.path = path
dt_start = bin_data.datetime_start
dt_stop = bin_data.datetime_stop
freq = bin_data.signal_frequency
self.__dt_file_start = dt_start
self.__dt_file_stop = dt_stop
self.ui.leFilePath.setText(path)
self.ui.sbResampleFrequency.setValue(freq)
self.ui.dtStartTime.setDateTime(dt_start)
self.ui.dtStopTime.setDateTime(dt_stop)
def load_origin_signal(self):
if self.__parameters is None:
self.__origin_signal = None
self.__origin_spectrogram = None
if self.__parameters is not None and self.__origin_signal is None:
params = self.__parameters
bin_data = BinaryFile()
bin_data.path = params['file_path']
bin_data.use_avg_values = False
bin_data.resample_frequency = params['resample_frequency']
bin_data.read_date_time_start = params['dt_start']
bin_data.read_date_time_stop = params['dt_stop']
if bin_data.signals is None:
self.__origin_signal = None
self.__origin_spectrogram = None
else:
index_channel = 'XYZ'.index(params['component'])
signal = bin_data.ordered_signal_by_components[:,
index_channel]
self.__origin_signal = signal
self.__origin_spectrogram = specgram(
signal_data=signal,
frequency_of_signal=params['resample_frequency'])
import os
import numpy as np
from SeisCore.BinaryFile.BinaryFile import BinaryFile
files = os.listdir('/media/michael/Transcend/1')
export_folder = '/media/michael/Data/TEMP'
for file in files:
if file in os.listdir(export_folder):
continue
path = os.path.join('/media/michael/Transcend/1', file)
bin_data = BinaryFile()
bin_data.path = path
signals = bin_data.signals
main_header = bin_data.main_header
bin_main_h = main_header.get_binary_format()
ch_header = [
bin_data.get_channel_header(channel_index=0),
bin_data.get_channel_header(channel_index=1),
bin_data.get_channel_header(channel_index=2)
]
bin_channels = [x.get_binary_format() for x in ch_header]
export_path = os.path.join(export_folder, file)
with open(export_path, 'wb') as handle:
handle.write(bin_main_h)
handle.write(bin_channels[0])
from SeisCore.BinaryFile.BinaryFile import BinaryFile
bin_data = BinaryFile()
bin_data.path=r'/media/michael/Data/Projects/GRP/DemkinskoeDeposit' \
r'/Demkinskoye_data/4772-4773/Mail/DM04_063_122_2019-07-30_07-14-11/DM04_063_122_2019-07-30_07-14-11.xx'
# bin_data.signal_frequency=1000
bin_data.record_type = 'ZXY'
# bin_data.read_date_time_start=datetime.datetime(year=2019, month=7, day=28,
# hour=13, minute=0, second=0)
# bin_data.read_date_time_stop=bin_data.read_date_time_start+datetime\
# .timedelta(hours=1)
# signal=bin_data.signals
print(bin_data.datetime_start, bin_data.datetime_stop)
print(bin_data.signals.shape)
overlap_size = 4093
median_param = 7
marmett_param = 7
time_step_minutes = 60
components = ['X', 'Y', 'Z']
f_min = 1
f_max = 10
# ----------------------- Don't touch ---------------------------------------
for root, folders, files in os.walk(root_folder):
for file in files:
name, extension = file.split('.')
if extension not in ['00', 'xx', 'bin']:
continue
bin_data = BinaryFile()
bin_data.path = os.path.join(root, file)
bin_data.use_avg_values = False
i = 0
while True:
left_time_lim = bin_data.datetime_start + timedelta(
minutes=i * time_step_minutes)
if left_time_lim > bin_data.datetime_stop:
break
right_time_lim = bin_data.datetime_start + timedelta(
minutes=(i + 1) * time_step_minutes)
if right_time_lim > bin_data.datetime_stop:
right_time_lim = bin_data.datetime_stop
bin_data.read_date_time_start = left_time_lim
bin_data.read_date_time_stop = right_time_lim
from SeisCore.BinaryFile.BinaryFile import BinaryFile
from SeisCore.Functions.Spectrogram import create_spectrogram
from datetime import timedelta
bin_file = r'/media/michael/77E9-B4C6/YM_002_K14_2019-12-08_04-30-40.xx'
bin_data = BinaryFile()
bin_data.path = bin_file
bin_data.resample_frequency = 250
bin_data.read_date_time_start = bin_data.datetime_start + timedelta(minutes=30)
bin_data.read_date_time_stop = bin_data.datetime_start + timedelta(minutes=60)
signal_z = bin_data.signals[:, 0]
create_spectrogram(signal_z, 250, '/media/michael/77E9-B4C6', 'test2', 1, 30,
0)
def run(self):
finished_parts_count = 0
params = self.parameters
for index, key in enumerate(self.files_info):
cur_file_data = self.files_info[key]
dt_start = cur_file_data['datetime_start']
dt_stop = cur_file_data['datetime_stop']
file_path = cur_file_data['path']
file_name = os.path.basename(file_path).split('.')[0]
cross_dt_start = max(params['dt_start'], dt_start)
cross_dt_stop = min(params['dt_stop'], dt_stop)
bin_data = BinaryFile()
bin_data.path = file_path
bin_data.use_avg_values = True
bin_data.resample_frequency = params['resample_frequency']
part_count = int(round(
(cross_dt_stop - cross_dt_start).total_seconds() // (
params['interval_size'] * 60)))
for i in range(part_count):
read_dt_start = cross_dt_start + timedelta(
minutes=params['interval_size'] * i)
read_dt_stop = read_dt_start + timedelta(
minutes=params['interval_size'])
read_dt_start_label = datetime.strftime(read_dt_start,
'%Y-%m-%d_%H-%M-%S')
bin_data.read_date_time_start = read_dt_start
bin_data.read_date_time_stop = read_dt_stop
signals = bin_data.signals
record_type=bin_data.record_type
for component in params['components']:
component_signal=signals[:, record_type.index(component)]
output_folder = os.path.join(params['output_folder'],
file_name,
f'{component}Component')
if not os.path.exists(output_folder):
os.makedirs(output_folder)
output_name = f'SP_{component}Component_{file_name}_{read_dt_start_label}'
create_spectrogram(signal_data=component_signal,
frequency=params['resample_frequency'],
output_folder=output_folder,
output_name=output_name,
min_frequency=params['visual_min_frequency'],
max_frequency=params['visual_max_frequency'],
time_start=0)
output_name = f'SIG_{component}Component_{file_name}_{read_dt_start_label}'
plot_signal(time_start_sec=0,
frequency=params['resample_frequency'],
signal=component_signal, label=output_name,
output_folder=output_folder,
output_name=output_name)
finished_parts_count += 1
percent_value = finished_parts_count / self.parts_count * 100
self.finished_percent.emit(percent_value)
from SeisCore.BinaryFile.BinaryFile import BinaryFile
from datetime import datetime, timedelta
from matplotlib import pyplot as plt
bin_data = BinaryFile()
bin_data.path = '/media/michael/Data/Projects/HydroFracturing/NorthDeposit' \
'/Perforation/Binary/12_SigmaN001_2020-04-11_12-15-15.bin'
signal = bin_data.signals[1000:2000]
plt.plot(signal)
plt.show()
print(signal.shape[0])
signal = signal[
(signal[:, 0] != 0) * (signal[:, 1] != 0) * (signal[:, 2] != 0)]
print(signal.shape[0])
print("")
exit(1)
bin_data.read_date_time_start = datetime(2019, 7, 28, 8, 0, 0) + timedelta(
seconds=-1 * 3)
bin_data.read_date_time_stop = datetime(2019, 7, 28, 8, 0, 2)
print(bin_data.discrete_amount)
print(bin_data.check_correct())
exit(0)
bin_data.read_date_time_start = datetime(2019, 7, 28, 8, 0, 0) + timedelta(
seconds=-1 * 3)
bin_data.read_date_time_stop = datetime(2019, 7, 28, 8, 0, 2)
signal = bin_data.signals[:, 0]
def run(self):
params = self.parameters
dt = (params['dt_stop'] - params['dt_start']).total_seconds()
discrete_amount = int(round(dt * params['resample_frequency']))
components_amount = len(params['components'])
files_amount = len(self.files_info.keys())
frequencies_list = rfftfreq(params['window_size'],
1 / params['resample_frequency'])
frequency_count = frequencies_list.shape[0]
joined_signal_array = np.empty(
shape=(components_amount, discrete_amount, files_amount),
dtype=np.int)
average_spectrum_data = np.empty(
shape=(components_amount, 2, frequency_count, files_amount),
dtype=np.float)
bin_files_list = list(self.files_info.keys())
for index in range(len(bin_files_list)):
bin_files_list[index] = bin_files_list[index].split('.')[0]
finished_parts_count = 0
for index_a, key in enumerate(self.files_info):
cur_file_data = self.files_info[key]
file_path = cur_file_data['path']
bin_data = BinaryFile()
bin_data.path = file_path
bin_data.use_avg_values = True
bin_data.resample_frequency = params['resample_frequency']
bin_data.read_date_time_start = params['dt_start']
bin_data.read_date_time_stop = params['dt_stop']
signals = bin_data.signals
if signals.shape[0]>discrete_amount:
signals=signals[:discrete_amount]
if params['detrend_frequency']!=0:
for i in range(signals.shape[1]):
signals[:,i]=detrend(signal=signals[:,i],
frequency=bin_data.resample_frequency,
edge_frequency=params['detrend_frequency'])
record_type=bin_data.record_type
for index_b, component in enumerate(params['components']):
compoment_signal = signals[:, record_type.index(component)]
joined_signal_array[index_b, :, index_a] = compoment_signal
# calculating average unsmoothed spectrum
av_spec_simple_component = average_spectrum(
signal=compoment_signal,
frequency=params['resample_frequency'],
window=params['window_size'],
offset=params['overlap_size'],
med_filter=None,
marmett_filter=None)
average_spectrum_data[index_b, 0, :, index_a] = \
av_spec_simple_component[:, 1]
# calculating average smoothed spectrum
av_spec_smooth_component = average_spectrum(
signal=compoment_signal,
frequency=params['resample_frequency'],
window=params['window_size'],
offset=params['overlap_size'],
med_filter=params['median_filter'],
marmett_filter=params['marmett_filter'])
average_spectrum_data[index_b, 1, :, index_a] = \
av_spec_smooth_component[:, 1]
finished_parts_count += 1
percent_value = finished_parts_count / self.parts_count * 100
self.finished_percent.emit(percent_value)
result_correlate_matrix = np.empty(
shape=(components_amount, files_amount, files_amount),
dtype=np.float)
for component in range(components_amount):
result_correlate_matrix[component, :, :] = cross_correlation(
frequency=frequencies_list,
f_min_analysis=params['correlation_min_frequency'],
f_max_analysis=params['correlation_max_frequency'],
amplitudes=average_spectrum_data[component, 1, :, :])
# result data exporting
read_dt_start_label = datetime.strftime(params['dt_start'],
'%Y-%m-%d_%H-%M-%S')
for index_a, bin_file_name in enumerate(self.files_info):
bin_file_name = bin_file_name.split('.')[0]
file_processing_result_folder = \
os.path.join(params['output_folder'], bin_file_name)
if not os.path.exists(file_processing_result_folder):
os.mkdir(file_processing_result_folder)
for index_b, component_label in enumerate(params['components']):
if params['is_selection_signal_file']:
dat_file_name = f'{bin_file_name}_ClearSignal_' \
f'{component_label}_Component.dat'
export_path = os.path.join(file_processing_result_folder,
dat_file_name)
np.savetxt(fname=export_path,
X=joined_signal_array[index_b, :, index_a],
fmt='%i')
if params['is_selection_signal_graph']:
png_file_name = \
f'{bin_file_name}_ClearSignal_{component_label}_{read_dt_start_label}'
plot_signal(time_start_sec=0,
frequency=params['resample_frequency'],
signal=joined_signal_array[index_b, :, index_a],
label=png_file_name,
output_folder=file_processing_result_folder,
output_name=png_file_name)
if params['is_spectrum_graph']:
png_file_name = '{}_AverageSpectrum_{}_Component_' \
'Graph'.format(bin_file_name,
component_label)
plot_average_spectrum(
frequency=frequencies_list,
origin_amplitudes=average_spectrum_data[
index_b, 0, :, index_a],
smooth_amplitudes=average_spectrum_data[
index_b, 1, :, index_a],
f_min=params['visual_min_frequency'],
f_max=params['visual_max_frequency'],
output_folder=file_processing_result_folder,
output_name=png_file_name)
if params['is_unsmoothed_spectrum_file']:
file_name = \
f'{bin_file_name}_{component_label}_Component.sc'
export_path = os.path.join(
file_processing_result_folder, file_name)
temp_array = np.empty(
shape=(frequencies_list.shape[0], 2),
dtype=np.float)
temp_array[:, 0] = frequencies_list
temp_array[:, 1] = average_spectrum_data[index_b, 0, :,
index_a]
np.savetxt(fname=export_path, X=temp_array, fmt='%f',
delimiter='\t')
if params['is_smoothed_spectrum_file']:
file_name = \
f'{bin_file_name}_{component_label}_Component.ssc'
export_path = os.path.join(
file_processing_result_folder, file_name)
temp_array = np.empty(
shape=(frequencies_list.shape[0], 2),
dtype=np.float)
temp_array[:, 0] = frequencies_list
temp_array[:, 1] = \
average_spectrum_data[index_b, 1, :, index_a]
np.savetxt(fname=export_path, X=temp_array, fmt='%f',
delimiter='\t')
if params['is_correlation_graph']:
png_file_name = \
f'{bin_file_name}_Separate_Correlation_{component_label}_Component_Graph'
plot_single_correlation(
devices=bin_files_list,
correlation_data=result_correlate_matrix[
index_b, index_a, :],
output_folder=file_processing_result_folder,
output_name=png_file_name)
finished_parts_count += 1
percent_value = finished_parts_count / self.parts_count * 100
self.finished_percent.emit(percent_value)
colors = list()
for i in range(len(self.files_info.keys())):
colors.append(generate_random_color())
for index_a, component_label in enumerate(params['components']):
if params['is_general_smooth_spectrum_graph']:
file_name = f'SmoothSpectrums_{component_label}_Component'
plot_all_smooth_spectrums(
spectrums_name_list=bin_files_list,
frequency=frequencies_list,
spectrum_data=average_spectrum_data[index_a, 1, :, :],
f_min_visualize=params['visual_min_frequency'],
f_max_visualize=params['visual_max_frequency'],
colors=colors, output_folder=params['output_folder'],
output_name=file_name)
if params['is_correlation_matrix_file']:
file_name = f'CorrelationMatrix_{component_label}_Component.dat'
header = 'NULL\t' + '\t'.join(bin_files_list) + '\n'
write_lines = list()
write_lines.append(header)
corr_matrix = result_correlate_matrix[index_a, :, :]
for i in range(corr_matrix.shape[0]):
t = list()
t.append(bin_files_list[i])
for j in range(corr_matrix.shape[1]):
t.append(str(corr_matrix[i, j]))
s = '\t'.join(t) + '\n'
write_lines.append(s)
export_path = os.path.join(params['output_folder'], file_name)
f = open(export_path, 'w')
for line in write_lines:
f.write(line)
f.close()
if params['is_general_correlation_graph']:
file_name = f'Correlations_{component_label}_Component'
plot_correlation(devices=bin_files_list,
colors=colors,
correlation_matrix=result_correlate_matrix[
index_a, :, :],
output_folder=params['output_folder'],
output_name=file_name)
from SeisCore.BinaryFile.BinaryFile import BinaryFile
import struct
import numpy as np
file = r'/media/michael/GSP1RMCSTFREO_RU_DVD/YM_013_K07_2019-11-25_00-00-00.xx'
bin_data=BinaryFile()
bin_data.path=file
bin_data.use_avg_values=False
signal=bin_data.signals
signal=signal[:5000000][:,0]
a=np.average(signal)
signal=signal-a
print(signal[0], signal[-1])
np.savetxt(r'/media/michael/MAHCYP32GB/190907_55207_Sigma№002/signal.dat',
signal, '%i', '\t', header='Z\tX\tY', comments='')
input()
#
# bin_data = open(file, 'rb')
# ch_count = _binary_read(bin_data, 12, 'I', 1)
# version = _binary_read(bin_data,0,'I',1)
# resolution = _binary_read(bin_data,0,'I',1)
# frequency = _binary_read(bin_data, 0, 'I', 1)
# # name=_binary_read(bin_data, 0, 's',12)
# latitude=_binary_read(bin_data, 12, 's',8)
# longitude=_binary_read(bin_data,0, 's', 9)
# start_date=_binary_read(bin_data,3,'I',1)
# start_time=_binary_read(bin_data,0,'I',1)
# print(frequency, latitude, longitude, start_date, start_time)
read_limits = [(datetime.datetime(2019, 11, 29,
13), datetime.datetime(2019, 11, 29, 13,
30))]
output_folder = r'/home/michael/Temp/Signals'
# read_dt_start=datetime.datetime(2019,12,5,14)
# read_dt_stop=datetime.datetime(2019,12,5,14,10)
bandpass_filter_limits = (0, 0.4)
bandpass_step = 0.1
spectrograms_freq_limits = (0, 250)
for file in files:
file_path = os.path.join(root, file)
bin_data = BinaryFile()
bin_data.path = file_path
bin_data.use_avg_values = False
for dt_lims in read_limits:
bin_data.read_date_time_start = dt_lims[0]
bin_data.read_date_time_stop = dt_lims[1]
signal = bin_data.signals[:, 0]
file_name = os.path.basename(file_path).split('.')[0]
for band_f_min in np.arange(bandpass_filter_limits[0],
bandpass_filter_limits[1], bandpass_step):
# band_f_max=band_f_min+bandpass_step
band_f_max = 500
exp_folder_name = f'{band_f_min}-{band_f_max}_Hz'
from SeisCore.BinaryFile.BinaryFile import BinaryFile
import numpy as np
origin_file = r'/home/michael/Загрузки/K42/YM_027_K42_2019-11-30_09-49.00'
test_file = r'/home/michael/Загрузки/K42/parts/file.00'
origin_bin = BinaryFile()
origin_bin.path = origin_file
origin_bin.use_avg_values = False
print(origin_bin.discrete_amount, origin_bin.signals.shape[0])
test_bin = BinaryFile()
test_bin.path = test_file
test_bin.use_avg_values = False
for i in range(3):
diff = np.abs(origin_bin.signals[:, i] - test_bin.signals[:, i])
print(np.max(diff))
import os
import numpy as np
from SeisCore.BinaryFile.BinaryFile import BinaryFile
root = r'/home/michael/Загрузки/K42/parts'
output_file_name = 'file.00'
file_statistics = list()
for file in os.listdir(root):
path = os.path.join(root, file)
bin_data = BinaryFile()
bin_data.path = os.path.join(root, file)
file_statistics.append(
(path, bin_data.datetime_start, bin_data.datetime_stop))
print(bin_data.discrete_amount)
file_statistics = sorted(file_statistics, key=lambda x: x[1])
is_right = True
for i in range(len(file_statistics) - 1):
t_stop = file_statistics[i][2]
t_start = file_statistics[i + 1][1]
dt = (t_start - t_stop).total_seconds()
if dt != 0.001:
is_right = False
break
if not is_right:
print('Different parts. The long time delay between signal parts')
exit()
join_signal_data = np.zeros(shape=(0, 3), dtype=np.int32)
import os
import numpy as np
from SeisCore.MSICore.CalcFunctions.Wavelett import detrend
from SeisCore.BinaryFile.BinaryFile import BinaryFile
folder = '/media/michael/Садыков_РР/mihael'
for root, folders, files in os.walk(folder):
for file in files:
name, extension = file.split('.')
if extension not in ('00', 'xx'):
continue
print(name)
bin_data = BinaryFile()
bin_data.path = os.path.join(root, file)
bin_data.record_type = 'ZXY'
signals = bin_data.ordered_signal_by_components
main_header = bin_data.main_header
channel_header_a = bin_data.get_channel_header(channel_index=0)
channel_header_b = bin_data.get_channel_header(channel_index=1)
channel_header_c = bin_data.get_channel_header(channel_index=2)
for i in range(3):
detrend_signal = detrend(signal=signals[:, i],
frequency=bin_data.signal_frequency,
edge_frequency=2)
signals[:, i] = detrend_signal[:signals[:, i].shape[0]]
output_file_path = os.path.join(
def calculation(self):
folder_path = self._ui.leInputFolder.text()
export_folder_path = self._ui.leOutputFolder.text()
file_list = self._ui.teFilesList.toPlainText()
dt_start = self._ui.dtMinDt.dateTime()
dt_stop = self._ui.dtMaxDt.dateTime()
if folder_path is None or file_list is None or export_folder_path is None:
return None
files = file_list.split('\n')
extract_data = None
for index, file in enumerate(files):
full_path = os.path.join(folder_path, file)
bin_data = BinaryFile()
bin_data.path = full_path
bin_data.record_type = 'XYZ'
bin_data.use_avg_values = False
bin_data.read_date_time_start = dt_start.toPyDateTime()
bin_data.read_date_time_stop = dt_stop.toPyDateTime()
bin_data.resample_frequency = self._ui.sbResampleFrequency.value()
if bin_data.signals is None:
continue
if extract_data is None:
discrete_count = bin_data.signals.shape[0]
extract_data = np.empty(shape=(discrete_count, 3 * len(files)),
dtype=np.int)
extract_data[:, 3 * index:3 * (index + 1)] = bin_data.signals
if extract_data is None:
return None
sum_trace = np.empty(shape=(extract_data.shape[0], 3),
dtype=np.int)
channel_plots = [self.channel_a_plot, self.channel_b_plot, self.channel_c_plot]
spectrum_plots = [self.channel_a_spectrum, self.channel_b_spectrum,
self.channel_c_spectrum]
window_size = self._ui.sbWindowSize.value()
overlap_size = self._ui.sbOverlapSize.value()
median_filter = self._ui.sbMedianFilter.value()
marmett_filter = self._ui.sbMarmettFilter.value()
resample_frequency = self._ui.sbResampleFrequency.value()
if marmett_filter == 0:
marmett_filter = None
if median_filter == 0:
median_filter = None
spectrums = None
for i in range(3):
channel_signals = extract_data[:, i:3 * len(files):3]
sum_trace[:, i] = np.sum(channel_signals, axis=1)
channel_plot = channel_plots[i]
spectrum_plot = spectrum_plots[i]
channel_plot.clear()
channel_plot.plot(sum_trace[:, i], pen=(196, 255, 18))
spectrum_plot.clear()
spectrum_data = average_spectrum(signal=sum_trace[:, i],
frequency=resample_frequency,
window=window_size,
overlap=overlap_size,
med_filter=median_filter,
marmett_filter=marmett_filter)
if spectrums is None:
spectrums = np.empty(shape=(spectrum_data.shape[0], 4),
dtype=np.float)
spectrums[:, :2] = spectrum_data
else:
spectrums[:, 1 + i] = spectrum_data[:, 1]
spectrum_plot.plot(spectrum_data, pen=(255, 0, 0))
tmp_file_a = os.path.join(export_folder_path, 'sum_trace.npy')
tmp_file_b = os.path.join(export_folder_path, 'spectrums.npy')
np.save(tmp_file_a, sum_trace)
np.save(tmp_file_b, spectrums)
output_file_sum_trace = os.path.join(export_folder_path, 'SumTraces.dat')
output_file_spectrums = os.path.join(export_folder_path, 'Spectrums.dat')
np.savetxt(output_file_sum_trace, sum_trace, fmt='%i', delimiter='\t',
header='Channel_1\tChannel_2\tChannel_3', comments='')
np.savetxt(output_file_spectrums, spectrums, fmt='%f', delimiter='\t',
header='Frequency\tChannel_1\tChannel_2\tChannel_3',
comments='')
