[发明专利]一种基于迭代收缩的快速稀疏Radon变换方法

权利要求书

1.一种基于迭代收缩的快速稀疏Radon变换方法，其特征在于，包括如下步骤：

S1：设置变量初始值；

S2：构造变换算子L并计算变换算子L的广义逆(L^TL)^-1L^T；

S3：利用步骤S2得到的变换算子L的广义逆(L^TL)^-1L^T对待处理地震道集d进行处理；

S4：判断地震数据体中所有道集是否全部处理完毕，如果否，返回步骤S3，如果全部处理完毕，则结束。

2.如权利要求1所述的基于迭代收缩的快速稀疏Radon变换方法，其特征在于，在所述步骤S1中，所述变量包括阈值系数α，迭代步长β，最大迭代次数N和RT变换参数。

3.如权利要求1所述的基于迭代收缩的快速稀疏Radon变换方法，其特征在于，所述步骤S3包括如下步骤：

S31：输入一个待处理的地震道集d，所述地震道集d为n维，所述n为大于1的正整数；

S32：设迭代数t＝0，利用步骤S2得到的变换算子L的广义逆(L^TL)^-1L^T求取n维模型m(t)=(L^TL)^-1L^Td；

S33：令t=t+1，得到更新后模型m(t+1)，

m(t+1)=T_α{m(t)+β(LTL)^-1L^T[d-Lm(t)]}

其中，T_α:Rⁿ→Rⁿ是收缩算子，定义为：

Tα{m}ij=(|mij|-αm~ij)+sgn(mij)]]>

其中，α是阈值,且0<α＜1,m={m_ij},将一个n维均值滤波器应用于|m|，得到滤波结果且

(x)+=xx≥00x<0;]]>

S34：判断迭代次数t是否达到最大迭代次数N，如没达到，返回步骤S33；如果达到，输出当前地震道集的处理结果。

4.如权利要求1所述的基于迭代收缩的快速稀疏Radon变换方法，其特征在于，在频域内实施或者在时域内实施。