Multi-frequency CRTomo data import

Import CRTomo-style multiy-frequency measurement data into a sEIT container

imports

import os
import subprocess
import reda

show the directory/file structure of the input data

# structure of the input data:
input_dir = 'data_crtomo_format'
for filename in sorted(os.listdir(input_dir)):
    print(filename)

# the frequencies.dat filename contains the frequencies, corresponding to the
# volt* files
print('content of frequencies.dat')
subprocess.call('cat {}/frequencies.dat'.format(input_dir), shell=True)

frequencies.dat
volt_01_0.462963Hz.crt
volt_02_1.0Hz.crt
volt_03_2.12766Hz.crt
volt_04_4.54545Hz.crt
volt_05_10.0Hz.crt
content of frequencies.dat

0

create the container

seit = reda.sEIT()
seit.import_crtomo('data_crtomo_format/')
seit.frequencies

Summary:
                 a            b  ...    frequency         rpha
count  4125.000000  4125.000000  ...  4125.000000  4125.000000
mean     19.608000    22.526788  ...     3.715244    -0.337727
std      10.681428    10.646378  ...     3.512372     1.339764
min       1.000000     2.000000  ...     0.462963    -7.241622
25%      10.000000    14.000000  ...     1.000000    -0.250827
50%      21.000000    26.000000  ...     2.127660    -0.037169
75%      29.000000    31.000000  ...     4.545454     0.208738
max      37.000000    38.000000  ...    10.000000     2.983013

[8 rows x 7 columns]

[0.46296296, 1.0, 2.1276596, 4.5454545, 10.0]