[發(fā)明專利]一種地層界線圖層的自動(dòng)化生成方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201610023762.2
• 申請(qǐng)日: 2016-01-14
• 公開（公告）號(hào): CN105701848B
• 公開（公告）日: 2018-07-20
• 發(fā)明（設(shè)計(jì)）人: 姚蒙蒙;李安波;陳楹;李夢(mèng)圓;吳春瑩
• 申請(qǐng)（專利權(quán)）人: 南京師范大學(xué)
• 主分類號(hào): G06T11/00
• 分類號(hào): G06T11/00
• 代理公司: 南京知識(shí)律師事務(wù)所 32207
• 代理人: 李媛媛

本發(fā)明公開了一種地層界線圖層的自動(dòng)化生成方法。該方法包括如下步驟：(1)基于沉積、火山地層界線圖層和變質(zhì)地層單位圖層，搜尋出圖層中所有鄰接地層，構(gòu)建地層鄰接關(guān)系圖；(2)基于地層鄰接關(guān)系圖提取相鄰地層的公共邊界，生成地層界線；(3)基于斷層圖層，識(shí)別斷層接觸關(guān)系；(4)基于產(chǎn)狀圖層，識(shí)別基本接觸關(guān)系，完成地層界線圖層的自動(dòng)化生成。與人工制作的方法相較，本發(fā)明方法生成的地層界線圖層不僅高效、快捷，且通過降低人工參與程度而盡可能避免人為誤差的引入，滿足了地質(zhì)圖空間數(shù)據(jù)庫地層界線圖層制作的需要。

技術(shù)領(lǐng)域

本發(fā)明屬于地理信息技術(shù)與應(yīng)用領(lǐng)域，具體涉及一種基于沉積和火山地層單位圖層、變質(zhì)地層單位圖層、產(chǎn)狀圖層和斷層圖層，自動(dòng)化生成地層界線圖層的方法。

背景技術(shù)

在地質(zhì)圖空間數(shù)據(jù)庫建設(shè)過程中，地層界線圖層的制作是一項(xiàng)必不可少的工作。在中國地質(zhì)調(diào)查局所制定的行業(yè)標(biāo)準(zhǔn)“數(shù)字化地質(zhì)圖圖層及屬性文件格式[DZ/T 0197-1997]”中，明確地對(duì)地層界線圖層的屬性數(shù)據(jù)做出了規(guī)定。

目前，地層界線圖層的制作主要通過人工制作來實(shí)現(xiàn)。制作方法主要包括兩方面工作：一是掃描并矢量化已有的地質(zhì)圖件得到地層界線，并對(duì)界線圖元進(jìn)行編輯，完成地層界線的提取；二是依據(jù)不同地層接觸關(guān)系在地質(zhì)圖上的表現(xiàn)，人工識(shí)別接觸關(guān)系并進(jìn)行屬性信息的錄入。該制作方法不僅存在工作量大、人工識(shí)別效率低，而且易引入人為誤差，如接觸關(guān)系識(shí)別錯(cuò)誤、屬性信息錯(cuò)誤、與相關(guān)地層單位圖層界線不一致等缺點(diǎn)。

在地質(zhì)圖空間數(shù)據(jù)庫建設(shè)過程中，著重制作的沉積和火山地層單位圖層、變質(zhì)地層單位圖層中，已蘊(yùn)含了相關(guān)地層界線位置和接觸關(guān)系信息。其中，地層界線的位置在沉積和火山地層單位圖層和變質(zhì)地層單位圖層中表現(xiàn)為具有相鄰?fù)負(fù)潢P(guān)系的兩地層要素的公共邊界。接觸關(guān)系可以基于地層先后順序、產(chǎn)狀數(shù)據(jù)、斷層數(shù)據(jù)等進(jìn)行自動(dòng)化識(shí)別。

發(fā)明內(nèi)容

為此，本發(fā)明主要針對(duì)人工制作地層界線圖層時(shí)出現(xiàn)的工作量大、識(shí)別效率低，以及易引入人為誤差等問題，提出了一種基于沉積和火山地層單位圖層、變質(zhì)地層單位圖層、產(chǎn)狀圖層和斷層圖層，進(jìn)行自動(dòng)化生成地層界線圖層的方法。

本發(fā)明的技術(shù)方案是：

一種地層界線圖層的自動(dòng)化生成方法，包括如下步驟：(1)基于沉積、火山地層界線圖層和變質(zhì)地層單位圖層，搜尋出圖層中所有鄰接地層，構(gòu)建地層鄰接關(guān)系圖(StratumAdjacency Relation Graph，簡稱SARG)；(2)基于地層鄰接關(guān)系圖提取相鄰地層的公共邊界，生成地層界線；(3)基于斷層圖層，識(shí)別斷層接觸關(guān)系；(4)基于產(chǎn)狀圖層，識(shí)別基本接觸關(guān)系，完成地層界線圖層的自動(dòng)化生成。

所述步驟(1)的具體步驟如下：

步驟11：加載沉積和火山地層單位圖層、變質(zhì)地層單位圖層數(shù)據(jù)，并將相應(yīng)的地層面要素讀入集合S＝{s_i|i＝1，2，…，n}，n為地層面要素的個(gè)數(shù)；

步驟12：讀出每個(gè)地層面要素s_i的頂點(diǎn)集v_i＝{P_j|j＝1，2，…，m}，其中P_j表示構(gòu)成地層面要素s_i的頂點(diǎn)，m為頂點(diǎn)的個(gè)數(shù)；并對(duì)每個(gè)地層面要素s_i的頂點(diǎn)集v_i進(jìn)行如下的預(yù)處理：

d)將頂點(diǎn)集v_i中的點(diǎn)P_j按順時(shí)針排序；

e)判斷每個(gè)頂點(diǎn)集v_i中是否存在P₁＝P_m，即首節(jié)點(diǎn)等于尾節(jié)點(diǎn)，若存在，則跳轉(zhuǎn)下一步

c)，對(duì)頂點(diǎn)集v_i進(jìn)行處理；否則，跳轉(zhuǎn)步驟13；