Beispiel #1
0
# calibration.experiment.plot(
#     x_type='energy',
#     source_to_detector_m=source_to_detector_m,
#     offset_us=offset_us,
#     logx=True,
#     fmt='-'
# )

calibrate_result = calibration.calibrate(
    source_to_detector_m=source_to_detector_m,
    offset_us=offset_us,
    vary='none',
    each_step=each_step)
calibration.index_peak(thres_exp=0.05,
                       min_dist_exp=10,
                       min_dist_map=5,
                       thres_map=0.05,
                       rel_tol=0.05)
# calibration.analyze_peak()
calibration.plot(
    y_type='attenuation',
    # y_type='transmission',
    x_type='energy',
    # t_unit='ms',
    # before=True,
    # interp=True,
    mixed=True,
    # peak_exp='all',
    table=False,
    peak_exp='indexed',
    peak_height=True,
                          energy_min=energy_min,
                          energy_max=energy_max,
                          energy_step=energy_step,
                          folder=folder,
                          baseline=baseline)

calibration.experiment.norm_to(norm_to_file, norm_factor=norm_factor)
calibration.experiment.slice(start=image_start, end=image_end)

calibrate_result = calibration.calibrate(
    source_to_detector_m=source_to_detector_m,
    offset_us=offset_us,
    vary='all',
    each_step=each_step)
calibration.index_peak(thres_exp=0.05,
                       min_dist_exp=2,
                       min_dist_map=5,
                       thres_map=0.05)
# calibration.analyze_peak()
calibration.plot(
    y_type='attenuation',
    # y_type='transmission',
    x_type='energy',
    # t_unit='ms',
    # before=True,
    # interp=True,
    # mixed=True,
    # peak_exp='all',
    # table=False,
    peak_exp='indexed',
    peak_height=False,
    index_level='iso',
Beispiel #3
0
                          baseline=baseline,
                          baseline_deg=baseline_deg,
                          x_type=x_type,
                          y_type=y_type)

calibration.experiment.norm_to(file=norm_to_file, norm_factor=norm_factor)
calibration.experiment.slice(start=image_start, end=image_end)

calibrate_result = calibration.calibrate(
    source_to_detector_m=source_to_detector_m,
    offset_us=offset_us,
    vary='all',
    each_step=each_step)
calibration.index_peak(thres_exp=0.15,
                       min_dist_exp=20,
                       min_dist_map=15,
                       thres_map=0.12,
                       rel_tol=0.01)
calibration.analyze_peak(report=False,
                         fit_model='Lorentzian')  # ['Gaussian', 'Lorentzian']

calibration.plot(
    y_type=y_type,
    x_type=x_type,
    # t_unit='ns',
    # before=True,
    # interp=True,
    mixed=True,
    table=True,
    peak_exp='all',
    peak_height=True,