[发明专利]一种探地雷达(GPR)数据分解的方法

摘要

本发明公开了一种探地雷达GPR数据分解的方法，首先设计GPR数据分解公式表示为对其进行优化，定义优化公式L(ωk,uk,ρ)=αΣk||∂fL(ωk)||22+||f-Σuk||22+ρ|f-Σuk|,]]>采用ωk计算步骤和uk计算步骤混合迭代优化。本发明的有益效果是提供了一种GPR数据高效分解的方法，使得后期探地雷达采集数据可视化效率提高。

主权项

一种探地雷达GPR数据分解的方法，其特征在于按照以下步骤进行：步骤1：GPR数据分解公式表示为：f(t)=Σk=1Puk(t)=Σk=1PAk(t)ejωkt+η]]>f(t)为GPR数据或观测信号，uk为f(t)的第K个子带，Ak和ωk分别为k子带幅值和频率，为信号的复指数表示，η为噪声；t为变量，t＝1,2,…；j为虚数单位；P为大于1的整数；步骤2：定义优化公式：L(ωk,uk,ρ)=αΣk||∂fL(ωk)||22+||f-Σuk||22+ρ|f-Σuk|]]>L(·)表示经过优化后的函数，ωk为k子带频率，uk为第K个子带；f为探地雷达采集到的阵列数据，范围为0‑1之间；fL表示1D单数据道信息，ρ和α分别为调节参数，表示进行偏微分计算；步骤3：采用ωk计算步骤和uk计算步骤混合迭代优化：ωk计算步骤：解决ωk的最小化问题，ωk的迭代公式为：ωkn+1=argminωkL(ωkn,ukn,ρ);]]>uk步骤：解决uk的最小化问题，迭代公式为：ukn+1=argminukL(ωkn,ukn,ρ).]]>