GPRPy - Ground Penetrating Radar Processing

Quick Reference

import gprpy.gprpy as gp
import matplotlib.pyplot as plt

# Load and display
data = gp.gprpyProfile()
data.importdata('profile.DZT')
data.showProfile()
plt.show()

# Access data
print(f"Traces: {data.data.shape[1]}")
print(f"Samples: {data.data.shape[0]}")
print(f"Time range: {data.twtt.max():.1f} ns")

Supported Formats

Essential Operations

Basic Processing

data = gp.gprpyProfile()
data.importdata('profile.DZT')

data.dewow(window=10)           # Remove low-frequency drift
data.remMeanTrace(ntraces=50)   # Remove background ringing
data.tpowGain(power=1.5)        # Time-power gain
data.agcGain(window=25)         # Automatic gain control

data.showProfile()

Apply Filters

data.bandpassFilter(minfreq=100, maxfreq=800)  # MHz
data.lowpassFilter(maxfreq=500)
data.highpassFilter(minfreq=50)

Time-to-Depth Conversion

velocity = 0.1  # m/ns (typical for dry sand)
data.setVelocity(velocity)
data.showProfile(yrng=[0, 5])  # Top 5 meters

Topographic Correction

data.topoCorrect(topofile='topography.txt', velocity=0.1)
# File format: x_position, elevation

Export Results

data.exportFig('processed.png', dpi=300)
data.exportSEGY('processed.sgy')
data.exportASCII('processed.txt')

Velocity Analysis (CMP)

cmp = gp.gprpyCMP()
cmp.importdata('cmp_survey.DZT')
cmp.showCMP()

cmp.semblance(vmin=0.05, vmax=0.15, vstep=0.01)
cmp.showSemblance()

Processing Parameters

Material Velocities

When to Use vs Alternatives

Use gprpy when you need open-source GPR processing in Python, batch processing of survey lines, or integration with other Python geoscience tools.

Consider alternatives when you need forward modelling of GPR responses (use GPRMax), advanced migration or commercial-grade processing (use REFLEXW), or highly custom signal processing algorithms (use scipy directly).

Common Workflows

Process raw GPR profile for interpretation

• Import raw data with gp.gprpyProfile() and importdata()
• Apply dewow filter to remove low-frequency drift
• Remove mean trace to eliminate background ringing
• Apply time-power gain or AGC for depth equalization
• Apply bandpass filter to remove noise
• Determine velocity from CMP analysis or material tables
• Convert time axis to depth with setVelocity()
• Apply topographic correction if survey has elevation changes
• Export processed profile as image and/or SEG-Y

References

• [Processing Steps](references/processing_steps.md) - Complete processing workflow guide
• [Material Velocities](references/processing_steps.md#velocity-selection) - Velocity selection by material type

Scripts

• [scripts/processgpr.py](scripts/processgpr.py) - Batch process GPR files with standard workflow