[發(fā)明專利]一種地球化學數(shù)據(jù)異常信息的提取方法及系統(tǒng)

在審

• 申請?zhí)枺?/td>: 201710321886.3
• 申請日: 2017-05-09
• 公開（公告）號: CN107291798A
• 公開（公告）日: 2017-10-24
• 發(fā)明（設計）人: 曹夢雪;路來君
• 申請（專利權）人: 吉林大學
• 主分類號: G06F17/30
• 分類號: G06F17/30
• 代理公司: 北京輕創(chuàng)知識產(chǎn)權代理有限公司11212
• 代理人: 楊立,朱毅

1.一種地球化學數(shù)據(jù)異常信息的提取方法，其特征在于，包括：

S1、從地球化學數(shù)據(jù)集合中抽取多個樣品數(shù)據(jù)，每個樣品數(shù)據(jù)均具有相同的變量組，所述變量組包含多個變量，將所述變量組設置為原始矩陣；

S2、將所述原始矩陣構造成反映所述地球化學數(shù)據(jù)集合中數(shù)據(jù)的實際空間分布信息的典型變量組合，將所述多個樣品數(shù)據(jù)的坐標構造成反映所述地球化學數(shù)據(jù)集合中數(shù)據(jù)的模擬空間分布信息的典型坐標組合；

S3、當所述典型變量組合與所述典型坐標組合滿足第一預設條件時，計算得到所述典型變量組合的第一系數(shù)與所述典型坐標組合的第二系數(shù)；

S4、基于所述第一系數(shù)和所述第二系數(shù)，當所述典型變量組合與所述典型坐標組合滿足第二預設條件時，計算得到第一常數(shù)向量與第二常數(shù)向量之間的關系式；

S5、基于所述第一系數(shù)、所述第二系數(shù)和所述第一常數(shù)向量與所述第二常數(shù)向量之間的關系式，當所述典型變量組合與所述典型坐標組合滿足第三預設條件時，得到所述地球化學數(shù)據(jù)集合的組合異常信息。

2.根據(jù)權利要求1所述的一種地球化學數(shù)據(jù)異常信息的提取方法，其特征在于，步驟S1中，從地球化學數(shù)據(jù)集合中抽取n個樣品數(shù)據(jù)，所述變量組包括m個變量，所述原始矩陣為X＝(β₁,β₂,…,β_m)，m和n均為正整數(shù)。

3.根據(jù)權利要求2所述的一種地球化學數(shù)據(jù)異常信息的提取方法，其特征在于，步驟S2中，將所述原始矩陣構造成典型變量組合u＝a₁β₁+a₂β₂+…+a_mβ_m，設所述第一系數(shù)a＝(a₁,a₂,…,a_m)，則u＝Xa。

4.根據(jù)權利要求3所述的一種地球化學數(shù)據(jù)異常信息的提取方法，其特征在于，步驟S2中，設第i個樣品數(shù)據(jù)的坐標(x_i,y_i)，i＝1,2,…,n，j＝1,2,…,m，設n個樣品數(shù)據(jù)的橫坐標x＝(x₁,x₂,…,x_n)和縱坐標y＝(y₁,y₂,…,y_n)，用n個樣品數(shù)據(jù)的坐標的傅里葉有限項級數(shù)表達另一典型變量v：

v=Σk=0MΣl=0Nλkl(αklcoskπxLcoslπyH+bklsinkπxLcoslπyH+cklcoskπxLsinlπyH+dklsinkπxLsinlπyH)]]>

將上式展開，得到將該公式中的各項表達成z₁,z₂,…,z_p，n個樣品數(shù)據(jù)的坐標調(diào)和多項式構成矩陣Z＝(z₁,z₂,…,z_p)，將典型變量v轉(zhuǎn)化為典型坐標組合設所述第二系數(shù)b＝(b₁,b₂,…,b_p)，則υ＝Zb；

其中，k為x方向上的調(diào)和階數(shù)，l為y方向上的調(diào)和階數(shù)，α_kl為第k階、l階調(diào)和的余弦-余弦項系數(shù)，為b_kl第k階、l階調(diào)和的正弦-余弦項系數(shù)，c_kl為第k階、l階調(diào)和的余弦-正弦項系數(shù)，d_kl為第k階、l階調(diào)和的正弦-正弦項系數(shù)，M為x方向上的最高調(diào)和階數(shù)，N為y方向上的最高調(diào)和階數(shù)，L為x方向上取樣長度之半，H為y方向上取樣長度之半，p為正整數(shù)。