Exemplo n.º 1
0
    align_file_name = 'Align_coeff.bin'  # Имя файла с текущими коэффициентами выравнивания АЧХ
    # current_catalog = r'2021/Results'           # Текущий каталог (за определенный период, здесь - год)

    current_data_dir = '2022'
    primary_data_dir = 'Primary_data'  # Каталог исходных данных (за определенный период, здесь - год)
    converted_data_dir = 'Converted_data'  # Каталог для записи результатов конвертации данных и заголовков
    data_treatment_dir = 'Data_treatment'  # Каталог для записи результатов обработки, рисунков

    current_primary_dir = '2022_01_27calibr'
    current_converted_dir = current_primary_dir + '_conv'
    current_converted_path = Path(converted_data_dir, current_converted_dir)
    current_treatment_dir = current_primary_dir + '_treat'
    current_treatment_path = Path(data_treatment_dir, current_treatment_dir)

    current_primary_file = '2022-01-27_13'
    converted_data_file_path, head_path = path_to_data(current_data_dir,
                                                       current_converted_path)
    data_treatment_file_path, head_path = path_to_data(current_data_dir,
                                                       current_treatment_path)

    folder_align_path = Path(head_path, 'Alignment')
    date = current_primary_file[0:10]

    # !!!! ******************************************* !!!!
    # ****** Блок исходных параметров для обработки *******

    freq_res = 8  # Установка разрешения по частоте в МГц
    kt = 1  # Установка разрешения по времени в единицах минимального разрешения 8.1925e-3 сек
    delta_t = 8.3886e-3
    delta_f = 7.8125
    att_val = [i * 0.5 for i in range(64)]
    att_dict = {s: 10**(s / 10) for s in att_val}
Exemplo n.º 2
0
            if num >= 10:
                add_pass0 = add_pass0[:_l - 2] + num_str
            else:
                add_pass0 = add_pass0[:_l - 2] + '0' + num_str
            add_pass1 = add_pass0 + '.' + format
    _path_to_picture = Path(_file_path, add_pass1)
    return _path_to_picture


if __name__ == '__main__':

    # ********************************** Путь к файлу данных ****************************
    current_catalog = r'2022\Results'  # Текущий каталог (за определенный период, здесь - год)
    current_data_dir = r'2022_01_20test'  # Папка с текущими даннымиH:\Fast_Acquisition\2021\Results\2021_09_22test\2021-09-22_01_14bit_pm20
    current_data_file = r'2022-01-20_02test'  # Имя файла с исходными текущими данными без расширенияH:\Fast_Acquisition\2022\Results\2022_01_20test
    file_path_data, head_path = path_to_data(current_catalog, current_data_dir)

    #           ************** Загрузка матрицы спектров *************
    spectrum = np.load(Path(file_path_data,
                            current_data_file + '_spectrum.npy'),
                       allow_pickle=True)
    #               **********************************************
    spectrum0 = spectrum[0]
    spectrum1 = spectrum[1]
    spectrum2 = spectrum[2]
    spectrum3 = spectrum[3]
    delta_t = 8.1925e-3
    spectrum_current = np.transpose(spectrum1[:, 20:22])
    spectrum_signal_av = low_freq_noise_spectrum(spectrum_current, 16384)
    legend = ['freq20', 'freq21']
    plot_low_freq_spec(spectrum_signal_av, delta_t,
if __name__ == '__main__':

    start = datetime.now()

    current_data_dir = '2022'
    primary_data_dir = 'Primary_data'           # Каталог исходных данных (за определенный период, здесь - год)
    converted_data_dir = 'Converted_data'       # Каталог для записи результатов конвертации данных и заголовков
    data_treatment_dir = 'Data_treatment'       # Каталог для записи результатов обработки, рисунков

    current_primary_dir = '2022_01_27test'
    current_primary_path = Path(primary_data_dir, current_primary_dir)
    current_converted_dir = current_primary_dir + '_conv'
    current_converted_path = Path(converted_data_dir, current_converted_dir)

    current_primary_file = '2022-01-27_02'
    primary_data_file_path, head_path = path_to_data(current_data_dir, current_primary_path)
    converted_data_file_path, head_path = path_to_data(current_data_dir, current_converted_path)

    align_file_name = 'Align_coeff.bin'         # Имя файла с текущими коэффициентами выравнивания АЧХ
    folder_align_path = Path(head_path, 'Alignment')

    date = current_primary_file[0:10]

    # !!!! ******************************************* !!!!
    # ****** Блок исходных параметров для обработки *******
    band_size_init = 'whole'
    num_of_polar = 2
    shift = 0
    # band_size = 'whole'   Параметр 'whole' означает работу в диапазоне 1-3 ГГц, 'half' - диапазон 1-2 или 2-3 ГГц
    # polar = 'both'        Принимает значения поляризаций: 'both', 'left', 'right'