Pure Python library for Ground Penetrating Radar (GPR) data processing: IO, filters and visualization. Tested with Python 3.6.
Either clone the repo and install with setup.py
git clone https://github.com/ahartikainen/siina
cd siina
python setup.py install
or with a pip
python -m pip install git+https://github.com/ahartikainen/siina
Header information is saved as a dictionary: obj.header
Measurement data is saved as a list of ndarrays: obj.data_list
Main channel can be accessed with .data -method
import siina
# create RadarFile object
meas = siina.Radar()
# read in the data
meas.read_file("./example_path/example_file.DZT")
# set the center frequency for GPR (in Hertz) if not done
if meas.header.get('frequency', None) is None:
meas.header['frequency'] = 1e9 # 1 GHz
# print dimensions for the data
print("points in samples={}, samples={}, channels={}".format(meas.nrows, meas.ncols, meas.nchan)
# strip markers (important step with .DZT files)
meas.read_markers()
# center each sample (for each trace do func(trace[500:])
meas.func_dc(start=500)
# apply lowpass filter with cutoff= 6 * frequency
# if cutoff is float -> cutoff = cutoff
# if cutoff is str -> cutoff = float(cutoff) * frequency
meas.func_filter(cutoff='6')
import matplotlib.pyplot as plt
# plot mean function for the first channel
# all channels are found under obj.data_list
plt.plot(meas.data.mean(1))
plt.show()
# plot radargram with plt.imshow
# be careful with the profile size (meas.ncols < 5000)
plt.imshow(meas.data, aspect='auto')
plt.show()