GPR 背景去除:擦掉公共水平条带
整张 GPR 剖面上常有一层「每道都长一样」的横向条带,像照片蒙了灰雾。擦掉它,局部目标才显出来。
从哪来
有些成分在全测线每一道都几乎相同:
- 直达波:发射天线直接窜到接收天线,每道开头都有;
- 地表或天线面的反射、耦合响应;
- 大范围的均匀层状反射。
它们在每个位置出现的时间和强度都差不多,画出来就是水平横条,盖住局部异常。
平均道就是背景
对每个时间样点,把所有道的值求平均,得到一条「平均道」。水平背景在每道都一样,在平均里相干保留;局部目标位置相位不同,互相抵消。
从每道减去这条平均道,背景没了,局部目标活下来。
统一公式:
$$D_{out}[t,x]=D[t,x]-\alpha\cdot\hat B[t,x],\qquad \alpha\in[0,1]$$
$\alpha=1$ 全减;$\alpha<1$ 部分减,适合既看点目标又想看连续界面的情况。
背景估计的几种方法
- 全局 mean/median:所有道一起求均值或中位数。
- 滑窗 mean/median:只取当前道附近窗口内的道求平均,适合背景沿测线缓慢变化。
- SVD:取数据矩阵前几个奇异图,背景通常低秩。
- f-k 滤波:在二维频谱里把低空间波数成分当背景去掉。
mean 还是 median?
一个 3×4 的迷你 B-scan:第 0 行是直达波(每道都是 50),第 1 行第 3 道有一根钢筋(30):
道1 道2 道3 道4 |
用 mean:行 1 平均被钢筋拉到 9,减完后其他道出现 −7/−8 的凹坑。
用 median:行 1 中位数是 2.5,减完后钢筋 +27.5,其他道几乎不受影响。
有孤立强目标时,优先用 median。
最大的陷阱
如果你要找的目标本身就是连续水平界面,比如层间脱空、连续裂隙,它在每道都出现,也会被当成背景减掉。
- 找钢筋、管线这种点目标:背景去除是好工具;
- 找脱空、连续裂隙这种面目标:关掉背景去除,或者把窗调得很大。
判读要点:顶部直达波亮纹变淡、局部双曲线更清楚,就对了。同时要警惕横向连续反射有没有莫名消失。
本博客所有文章除特别声明外,均采用 CC BY-NC-SA 4.0 许可协议。转载请注明来自 CYK's Blog!
评论
