[发明专利]一种基于快速独立分量分析的地震属性融合方法

权利要求书

1.一种基于快速独立分量分析的地震属性融合方法，其特征在于包括：

（1）对待融合属性数据分块：选取大小为N×N的窗W对各大小均为M1×M2待融合属性切片平滑分块，把分出来的每一个小块存储成N²×1的向量，放到各属性块的集合{Ii}里，得到一个块的集合C，C={{I1},{I2},…,{In}}；

（2）选块：根据集合C的维数和各小块的均值，从集合C里面随机的选取P块属性块作为ICA算法的输入来求取ICA变换的变换核T{·}；

（3）求取ICA域的变换核：由步骤2得到一个由P块属性块组成的矩阵mixedX,其维数为N²×P，把mixedX作为快速独立分量分析算法的输入，求出mixedX的协方差矩阵CovX，CovX的维数为N²×N²，经过特征分解得到白化矩阵WhiteMattix和去白化矩阵deWhiteMattix，去白化矩阵与白化矩阵是互逆矩阵，通过求出的白化矩阵，求出白化信号Whitesig，把求得的白化信号作为快速独立分量核心算法中的x，求出ICA域的变换核T{·}即分离矩阵W和ICA域的变换核的逆T^-1{·}即混合矩阵A；

（4）空间域映射到ICA域：求出各属性对应块的独立分量，即把空域的块映射到ICA域；

（5）ICA域融合：由步骤4把各属性的所有块均映射到ICA域里面，用如下的融合规则得到第i块在ICA域的融合结果：

Fi(l)=Σk=1n|SSki(l)|SSki(l)Σk=1n|SSki(l)|,l=1,2,...,N2---(2)]]>

依照公式(2)求出所有块在ICA域的融合结果；

（6）ICA域映射到空域：第i块在空域融合的结果,

MFi=A×Fi+1nΣs=1nmsi---(3)]]>

依照公式(3)得到所有块在空域的融合结果，然后把(M1-N+1)*(M2-N+1)*n块属性块的空间域融合结果组合成M1*M2大小的融合属性，对得到的融合结果采用求平均值法分块去重叠处理，得到最终的融合结果。