本發(fā)明屬于地理信息技術(shù)應(yīng)用領(lǐng)域，具體涉及一種基于地形地質(zhì)圖，自適應(yīng)判定地表出露巖層產(chǎn)狀的方法。

背景技術(shù)：

巖層的產(chǎn)狀是指巖層在空間的產(chǎn)出狀態(tài)。除水平巖層成水平狀態(tài)產(chǎn)出外，一切傾斜巖層的產(chǎn)狀均以其走向、傾向和傾角表示，稱為巖層產(chǎn)狀三要素。產(chǎn)狀作為地質(zhì)測(cè)量與成圖中常見(jiàn)的地質(zhì)要素之一，在地質(zhì)各領(lǐng)域有著廣泛應(yīng)用，其獲取方法主要包括野外實(shí)測(cè)和間接計(jì)算兩種方法。

野外的巖層產(chǎn)狀測(cè)量方法，通常是使用地質(zhì)羅盤(pán)、坡度儀等儀器，直接在出露的巖層層面上量測(cè)。但是，野外巖石的風(fēng)化、層面和劈理、線理面區(qū)分困難等因素均會(huì)影響產(chǎn)狀的測(cè)量。有時(shí)受限于巖層出露情況，只能借鑒附近區(qū)域的數(shù)據(jù)，再通過(guò)室內(nèi)計(jì)算方法間接獲得產(chǎn)狀數(shù)據(jù)。產(chǎn)狀間接計(jì)算中常用的方法，主要有“三點(diǎn)法”和“相鄰等高線法”，其原理為通過(guò)尋找共面且不共線的多點(diǎn)，確定平面的形態(tài)。在計(jì)算過(guò)程中，一方面，需要人工篩選合適的點(diǎn)和計(jì)算方法，并通過(guò)繪制輔助點(diǎn)和線，對(duì)產(chǎn)狀進(jìn)行計(jì)算；另一方面，需要人為確定產(chǎn)狀符號(hào)位置。相關(guān)計(jì)算工作的人為化處理，導(dǎo)致計(jì)算效率低下，且產(chǎn)狀數(shù)據(jù)的質(zhì)量也因人而異。

為此，本專利針對(duì)地表巖層，擬基于地形地質(zhì)圖，實(shí)現(xiàn)一種地表巖層產(chǎn)狀的自適應(yīng)計(jì)算方法。其重點(diǎn)在于自適應(yīng)完成交點(diǎn)選擇、計(jì)算方法選擇和產(chǎn)狀點(diǎn)位選擇。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明所要解決的技術(shù)問(wèn)題在于，克服現(xiàn)有技術(shù)存在的缺陷，提出了一種地表巖層產(chǎn)狀的自適應(yīng)判定方法，基于地形地質(zhì)圖，自適應(yīng)完成交點(diǎn)選擇、計(jì)算方法選擇和產(chǎn)狀點(diǎn)位選擇。從而實(shí)現(xiàn)地表出露巖層產(chǎn)狀的自適應(yīng)判定。

本專利的主要思路如下：1)提取和預(yù)處理巖層邊界線，計(jì)算巖層界線與等高線的交點(diǎn)；2)遍歷同一巖層界線上的鄰接交點(diǎn)，通過(guò)智能的交點(diǎn)選擇、計(jì)算方法選擇和產(chǎn)狀點(diǎn)位選擇等步驟，實(shí)現(xiàn)產(chǎn)狀的自適應(yīng)判定；3)對(duì)產(chǎn)狀信息進(jìn)行符號(hào)化和標(biāo)注。

基于上述思路，本發(fā)明提出的一種地表出露巖層產(chǎn)狀的自適應(yīng)判定方法，其步驟如下：

(一)交點(diǎn)集合獲取

步驟1.1：加載矢量格式的巖層面要素圖層和等高線線圖層，得到所有面要素集合Stra＝{s_i|i＝1,2,3,…,n}(其中：s_i為巖層面要素，n為巖層面要素的數(shù)量)和等高線的線要素集合CourLine＝{cl_i|i＝1,2,3,…,p}(其中：cl_i為等高線的線要素，p為等高線要素的數(shù)量)。

步驟1.2：讀取集合Stra中每個(gè)巖層面要素的邊界，將其預(yù)處理為單線要素，并剔除圖幅邊緣部分、第四紀(jì)巖層邊界部分、斷層界線部分以及不整合接觸界線部分，得到巖層層面邊界線集合SLine＝{sl_i|i＝1,2,3,…,c}，sl_i是巖層層面邊界線要素，c為所有巖層保留的層面出露邊界線的數(shù)量。

步驟1.3：針對(duì)任一邊界線要素sl_i遍歷等高線集合CourLine，求取線要素sl_i與所有等高線的交點(diǎn)集合InterP_i，對(duì)邊界線集合SLine內(nèi)所有元素進(jìn)行交點(diǎn)計(jì)算，得到集合InterP＝InterP₁∪InterP₂∪InterP₃∪…∪InterP_l，l為與等高線相交的邊界線數(shù)量。記錄交點(diǎn)的高程H，源巖層編號(hào)SID和源等高線編號(hào)CID，在同一集合InterP_i內(nèi)，交點(diǎn)的SID屬性均為i。

(二)產(chǎn)狀自適應(yīng)獲取與剔除

針對(duì)交點(diǎn)集合InterP_i＝{pi_j|j＝1,2,3,…,qi}，qi為交點(diǎn)集合InterP_i包含的交點(diǎn)數(shù)量，逐點(diǎn)進(jìn)行條件判斷，獲取產(chǎn)狀a＝{ρ，θ，δ}，并確定產(chǎn)狀點(diǎn)位。其中：ρ為巖層走向，θ是巖層傾向，δ是巖層傾角。

由j＝1開(kāi)始進(jìn)行判斷，具體步驟為：

步驟2.1：讀取交點(diǎn)pi_j、pi_j+1和pi_j+2和pi_j+3(若下標(biāo)超出qi，則使用pi_(j+1)％qi、pi_(j+2)％qi和pi_(j+3)％qi)，判別交點(diǎn)pi_j、pi_j+1、pi_j+2和pi_j+3是否滿足相鄰等高線法的下列選點(diǎn)規(guī)則；如果符合，執(zhí)行步驟2.4，如果不符合執(zhí)行步驟2.2：

a)交點(diǎn)pi_j和pi_j+3高程一致，交點(diǎn)pi_j+1和pi_j+2高程一致；

b)交點(diǎn)pi_j+1和pi_j+2的源等高線編號(hào)屬性CID一致；

c)向量和的最小夾角不大于閾值ε，即判定兩向量方位近似平行；閾值ε的設(shè)定一般不大于6°，兩線夾角在此范圍內(nèi)可視為近似平行。閾值的設(shè)定越小，則平行要求越嚴(yán)格，可根據(jù)需求調(diào)整閾值以滿足不同要求。

步驟2.2：讀取交點(diǎn)pi_j、pi_j+1和pi_j+2(若下標(biāo)超出qi，則使用pi_(j+1)％qi、pi_(j+2)％qi)，判別是否滿足三點(diǎn)法的下列選點(diǎn)規(guī)則；如果符合，執(zhí)行步驟2.6；如果不符合，執(zhí)行步驟2.3:

a)交點(diǎn)pi_j、pi_j+1和pi_j+2中至少兩點(diǎn)的高程不等；

b)交點(diǎn)pi_j、pi_j+1和pi_j+2不共線。

步驟2.3：重新選擇判斷起點(diǎn)。若原線要素sl_i為環(huán)線要素，則取d＝qi為界，若原線要素sl_i非環(huán)線要素，則取d＝qi-3為界。下標(biāo)j增加1之后，若j不超過(guò)d，則返回步驟2.1；若j超過(guò)d，則交點(diǎn)集合InterP_i上的產(chǎn)狀信息獲取完成，執(zhí)行步驟2.8。

步驟2.4：相鄰等高線法計(jì)算產(chǎn)狀?；诮稽c(diǎn)pi_j(x_j,y_j,h₁)、pi_j+1(x_j+1,y_j+1,h₂)、pi_j+2(x_j+2,y_j+2,h₂)和pi_j+3(x_j+3,y_j+3,h₁)使用相鄰等高線法計(jì)算產(chǎn)狀：

a)對(duì)于交點(diǎn)pi_j、pi_j+1、pi_j+2和pi_j+3，計(jì)算平面向量和分

別取和方向上的單位向量和取和向量和的方位角度ρ為巖層走向；

b)根據(jù)公式(1)計(jì)算巖層傾向θ；

c)計(jì)算交點(diǎn)pi_j+1和交點(diǎn)pi_j+2的中點(diǎn)，記為輔助點(diǎn)ap₁(x₁,y₁,h₂)，過(guò)點(diǎn)ap₁，沿θ方向的

直線與向量v₂交于點(diǎn)ap₂(x₂,y₂,h₁)，根據(jù)公式(2)計(jì)算傾角δ；

步驟2.5：相鄰等高線法產(chǎn)狀點(diǎn)位自動(dòng)確定。記直線l為過(guò)輔助點(diǎn)ap₁且方向?yàn)棣鹊闹本€，l與交點(diǎn)pi_j+1與pi_j+2之間的sl_i邊界線部分相交于點(diǎn)ap₃，則以ap₃坐標(biāo)(x₃,y₃)作為產(chǎn)狀標(biāo)注點(diǎn)位坐標(biāo)，存儲(chǔ)產(chǎn)狀信息a。返回步驟2.2。

步驟2.6：三點(diǎn)法計(jì)算產(chǎn)狀?；诮稽c(diǎn)pi_j(x₁,y₁,h₁)，pi_j+1(x₂,y₂,h₂)，pi_j+2(x₃,y₃,h₃)使用三點(diǎn)法計(jì)算產(chǎn)狀：

a)將交點(diǎn)pi_j、pi_j+1和pi_j+2的坐標(biāo)代入公式(3)，得到巖層層面方程(公式(4)所示，

其中D為任意常數(shù))的系數(shù)A、B和C；

Ax+By+Cz+D＝0 (4)

b)根據(jù)系數(shù)A、B，確定巖層傾向線方程(公式(5)所示)，其中M為任意常數(shù)；

Bx+Ay+M＝0 (5)

c)由公式(6)、(7)分別計(jì)算、修正傾向θ；

d)根據(jù)傾向θ，由公式(8)計(jì)算巖層走向；

ρ＝θ±90 (8)

e)根據(jù)系數(shù)A、B和C，由公式(9)計(jì)算巖層傾角δ；

步驟2.7：三點(diǎn)法產(chǎn)狀標(biāo)注點(diǎn)位自動(dòng)確定。以交點(diǎn)pi_j+1的坐標(biāo)作為產(chǎn)狀點(diǎn)位坐標(biāo)，存儲(chǔ)產(chǎn)狀信息a。返回步驟2.3。

步驟2.8：遍歷巖層邊界線sl_i對(duì)應(yīng)的交點(diǎn)集合InterP_i上的交點(diǎn)，得到Att_i＝{ai_j|j＝1,2,3,…,ci}，其中ci為邊界線sl_i上的產(chǎn)狀點(diǎn)數(shù)量。對(duì)邊界線集合SLine內(nèi)所有巖層邊界線sl_i進(jìn)行產(chǎn)狀計(jì)算，得到產(chǎn)狀集合Att＝Att₁∪Att₂∪Att₃∪…Att_s，s為有效獲取了產(chǎn)狀的交點(diǎn)集合數(shù)量。

步驟2.9：對(duì)具明顯錯(cuò)誤的產(chǎn)狀值進(jìn)行剔除。以產(chǎn)狀集合Att_i為單位，基于構(gòu)造特征結(jié)合地質(zhì)知識(shí)對(duì)錯(cuò)誤的產(chǎn)狀值進(jìn)行剔除，得到篩選后的產(chǎn)狀集合Att_i＝{ai_j|j＝1,2,3,…,di}，其中di為邊界線sl_i上的糾錯(cuò)后保留的產(chǎn)狀點(diǎn)數(shù)量。錯(cuò)誤產(chǎn)狀值剔除步驟如下：

a)褶皺具有“兩翼產(chǎn)狀較穩(wěn)定，中間部分產(chǎn)狀過(guò)渡”的構(gòu)造特征。結(jié)合“V字型法則”，整體判斷兩翼產(chǎn)狀傾向。具體表現(xiàn)為：同斜褶皺兩翼產(chǎn)狀同向，直立向斜兩翼產(chǎn)狀指向外部，直立背斜兩翼產(chǎn)狀指向內(nèi)部。剔除褶皺上不滿足上述褶皺構(gòu)造特征的產(chǎn)狀值；

b)若巖層厚度均勻，則兩側(cè)產(chǎn)狀基本為同向。對(duì)于巖層形狀簡(jiǎn)單，兩側(cè)界線基本平行的巖層，首先可根據(jù)周圍構(gòu)造或結(jié)合“V字型法則”大致確定一側(cè)界線上的產(chǎn)狀；繼而推測(cè)另一側(cè)界線上的產(chǎn)狀。剔除簡(jiǎn)單巖層上不滿足其產(chǎn)狀特征的產(chǎn)狀值；

c)形態(tài)較平直的界線部分，若界線上產(chǎn)狀傾向不穩(wěn)定，統(tǒng)計(jì)界線上產(chǎn)狀傾向在不同方向區(qū)間上的頻率，保留高頻并連續(xù)的產(chǎn)狀區(qū)間數(shù)據(jù)。區(qū)間大小根據(jù)實(shí)際情況而定，在產(chǎn)狀個(gè)數(shù)較少的情況下，可取較大的區(qū)間。取一平直界線上的的16個(gè)產(chǎn)狀數(shù)據(jù)為例(如下表所示)，以10°為間隔對(duì)產(chǎn)狀傾向做雷達(dá)圖(如圖2所示)。其產(chǎn)狀數(shù)據(jù)傾向連續(xù)且較多分布于90°到130°之間，則保留90°到130°之間的產(chǎn)狀數(shù)據(jù)，剔除其他低頻、離散的數(shù)據(jù)值

步驟3、產(chǎn)狀數(shù)據(jù)符號(hào)化

步驟3.1：對(duì)于產(chǎn)狀點(diǎn)ai_j，過(guò)其坐標(biāo)(axi_j,ayi_j)點(diǎn)，沿其傾向θ繪制單位長(zhǎng)度的線，垂直于傾向θ繪制單位長(zhǎng)度的線，完成產(chǎn)狀符號(hào)的繪制。

步驟3.2：將傾向θ由極坐標(biāo)系轉(zhuǎn)為大地坐標(biāo)系，對(duì)產(chǎn)狀要素傾角δ進(jìn)行標(biāo)注。

本發(fā)明方法針對(duì)地表出露巖層，基于地形地質(zhì)圖，自適應(yīng)完成交點(diǎn)選擇、計(jì)算方法選擇和產(chǎn)狀點(diǎn)位選擇的產(chǎn)狀判定，實(shí)現(xiàn)了地表出露巖層產(chǎn)狀的自適應(yīng)判定，實(shí)用性強(qiáng)，適用性廣。

附圖說(shuō)明

圖1本發(fā)明方法流程圖。

圖2示例產(chǎn)狀傾向雷達(dá)圖

圖3實(shí)施例實(shí)驗(yàn)數(shù)據(jù)。

圖4巖層邊界線和等高線交點(diǎn)圖。

圖5示例交點(diǎn)圖。

圖6示例三點(diǎn)法自適應(yīng)選點(diǎn)結(jié)果圖。

圖7示例相鄰等高線法自適應(yīng)選點(diǎn)結(jié)果圖。

圖8實(shí)施例全圖三點(diǎn)法自適應(yīng)選點(diǎn)結(jié)果圖。

圖9實(shí)施例全圖相鄰等高線法自適應(yīng)選點(diǎn)結(jié)果圖。

圖10實(shí)施例全圖產(chǎn)狀符號(hào)化和標(biāo)注結(jié)果圖。

圖11示例實(shí)測(cè)產(chǎn)狀數(shù)據(jù)與計(jì)算產(chǎn)狀數(shù)據(jù)對(duì)比圖。

具體實(shí)施方式

下面結(jié)合實(shí)施例和附圖，對(duì)本發(fā)明方法作進(jìn)一步詳細(xì)說(shuō)明。

實(shí)施例：選用某地1:50000地形地質(zhì)圖(如圖2所示)數(shù)據(jù)，下面結(jié)合附圖，基于shp格式的地質(zhì)體面圖層和等高線圖層，描述本方法方法實(shí)施過(guò)程并進(jìn)一步說(shuō)明本發(fā)明的效果。

(一)交點(diǎn)集合獲取

步驟1.1：加載巖層面要素圖層和等高線圖層，得到所有面要素集合Stra＝{s_i|i＝1,2,3,…,29}，等高線要素集合CourLine＝{cl_i|i＝1,2,3,…,21}。

步驟1.2：讀取集合Stra中巖層面要素的邊界線分別進(jìn)行預(yù)處理，去除第四紀(jì)地層邊界線及地層與第四紀(jì)接觸的邊界線部分，并去除圖幅邊界的部分，得到邊界線集合SLine＝{sl_i|i＝1,2,3,…,47}。

步驟1.3：對(duì)邊界線集合SLine和等高線集合CourLine求取交點(diǎn)，邊界線要素sli上的交點(diǎn)存入一個(gè)集合InterP_i，得到交點(diǎn)集合InterP₁＝{(400875.27,3271721.62),(400957.22,3271801.49),(401056.83,3271892.77)}，InterP₂＝{(401235.55,3272047.50),(401334.05,3272128.06),(401434.51,3272198.33),(401551.99,3272129.90),(401497.99,3272039.65),(401398.55,3272039.69),(401376.85,3271795.00),(401606.74,3271790.62)}，…，InterP₄₄＝{(400144.21,3274334.45),(400102.31,3274263.53),(400051.63,3274186.33),(400050.79,3273975.12),(400131.35,3274013.82),(400196.16,3274044.96),(400278.58,3274088.92),(400354.36,3274135.09)}，交點(diǎn)記錄等高線高程H，源巖層編號(hào)SID和源等高線編號(hào)CID，在同一集合InterP_i內(nèi)，交點(diǎn)的SID屬性一致為i。所有交點(diǎn)集合InterP＝InterP₁∪InterP₂∪InterP₃∪…∪InterP₄₄，包括交點(diǎn)398個(gè)(如圖3所示)。

(二)產(chǎn)狀自適應(yīng)計(jì)算與剔除

針對(duì)交點(diǎn)集合自適應(yīng)選點(diǎn)、計(jì)算及點(diǎn)位確定，在此設(shè)定向量平行判斷的角度閾值ε為6°。以巖層面要素s₇的邊界線為例，其邊界包括巖層邊界線sl₁₁、sl₁₂、sl₁₃和sl₁₄，與等高線分別存在16個(gè)、3個(gè)、0個(gè)和9個(gè)交點(diǎn)，則邊界線sl₁₁、sl₁₂和sl₁₄分別對(duì)應(yīng)交點(diǎn)集合InterP₉＝{p9_j|j＝1,2,3,…,16}、InterP₁₀＝{p10_j|j＝1,2,3}和InterP₁₁＝{p11_j|j＝1,2,3,…,9}(如圖4所示)。以交點(diǎn)集合InterP₉為例，交點(diǎn)坐標(biāo)值如下表所示。由交點(diǎn)p9₁開(kāi)始進(jìn)行判斷和計(jì)算；

步驟2.1：由交點(diǎn)p9₁開(kāi)始，先對(duì)交點(diǎn)p9₁、p9₂、p9₃和p9₄進(jìn)行相鄰等高線法的選點(diǎn)規(guī)則判斷，并不符合相鄰等高線法的選點(diǎn)規(guī)則；則對(duì)交點(diǎn)p9₁、p9₂和p9₃進(jìn)行三點(diǎn)法的選點(diǎn)規(guī)則判斷，并不符合三點(diǎn)法的選點(diǎn)規(guī)則，則繼續(xù)判斷交點(diǎn)p9₂。

步驟2.2：由交點(diǎn)p9₂開(kāi)始，對(duì)交點(diǎn)p9₂、p9₃、p9₄和p9₅進(jìn)行相鄰等高線法的選點(diǎn)規(guī)則判斷，并不符合相鄰等高線法的選點(diǎn)規(guī)則；則對(duì)交點(diǎn)p9₂、p9₃和p9₄進(jìn)行三點(diǎn)法的選點(diǎn)規(guī)則判斷，三點(diǎn)滿足如下條件：

1)三點(diǎn)的高程分別為250、250和300，高程不全相等；

2)三點(diǎn)不共線；

因此，三點(diǎn)滿足三點(diǎn)法選點(diǎn)規(guī)則，進(jìn)行三點(diǎn)法計(jì)算。

步驟2.3：將交點(diǎn)p9₂、p9₆和p9₇代入公式計(jì)算，得到所在平面方程的系數(shù)：

計(jì)算得到產(chǎn)狀結(jié)果：

步驟2.4：產(chǎn)狀符號(hào)位置取交點(diǎn)p9₃的坐標(biāo)，繼續(xù)判斷交點(diǎn)p9₃。

步驟2.5：由交點(diǎn)p9₃開(kāi)始，對(duì)交點(diǎn)p9₃、p9₄、p9₅和p9₆進(jìn)行相鄰等高線法的選點(diǎn)規(guī)則判斷，并不符合相鄰等高線法的選點(diǎn)規(guī)則；則對(duì)交點(diǎn)p9₃、p9₄和p9₅進(jìn)行三點(diǎn)法的選點(diǎn)規(guī)則判斷，三點(diǎn)滿足三點(diǎn)法選點(diǎn)規(guī)則，進(jìn)行三點(diǎn)法計(jì)算，得到產(chǎn)狀結(jié)果：

產(chǎn)狀符號(hào)位置取交點(diǎn)p9₄的坐標(biāo)，繼續(xù)判斷交點(diǎn)p9₄。

由交點(diǎn)p9₄開(kāi)始，對(duì)交點(diǎn)p9₄、p9₅、p9₆和p9₇進(jìn)行相鄰等高線法的選點(diǎn)規(guī)則判斷，并不符合；則對(duì)交點(diǎn)p9₄、p9₅和p9₆進(jìn)行三點(diǎn)法的選點(diǎn)規(guī)則判斷，三點(diǎn)滿足三點(diǎn)法選點(diǎn)規(guī)則，進(jìn)行三點(diǎn)法計(jì)算，得到產(chǎn)狀結(jié)果：

產(chǎn)狀符號(hào)位置取交點(diǎn)p9₅的坐標(biāo)，繼續(xù)判斷交點(diǎn)p9₅。

步驟2.6：由交點(diǎn)p9₅開(kāi)始，對(duì)交點(diǎn)p9₅、p9₆、p9₇和p9₈進(jìn)行相鄰等高線法的選點(diǎn)規(guī)則判斷，四點(diǎn)滿足：

1)交點(diǎn)p9₅和交點(diǎn)p9₈高程一致為350m；交點(diǎn)p9₆和交點(diǎn)p9₇高程一致為400m；

2)交點(diǎn)p9₆和交點(diǎn)p9₇的源等高線編號(hào)屬性CID一致；

3)交點(diǎn)p9₅指向交點(diǎn)p9₈的方位角為214.75°，交點(diǎn)p9₆指向交點(diǎn)p9₇的方位角為217.06°，兩方位角相差2.31°，小于閾值ε。

因此，四點(diǎn)滿足相鄰等高線法，利用相鄰等高線法進(jìn)行產(chǎn)狀計(jì)算。計(jì)算得到產(chǎn)狀結(jié)果：

步驟2.7：計(jì)算得到產(chǎn)狀符號(hào)位置坐標(biāo)(400450.92,3272412.92)，繼續(xù)判斷交點(diǎn)p9₅、p9₆和p9₇是否滿足三點(diǎn)法選點(diǎn)規(guī)則，直至遍歷完成。

步驟2.8：交點(diǎn)集InterP₉產(chǎn)狀計(jì)算完成，得到邊界線sl₁₁的產(chǎn)狀集合Att₆＝{a6_j|j＝1,2,3,…,15}。其中，使用三點(diǎn)法計(jì)算得到的有13個(gè)，其三點(diǎn)組合形式如圖5所示；使用相鄰等高線法計(jì)算得到的有2個(gè)，其四點(diǎn)組合形式如圖6所示。

步驟2.9：對(duì)所有InterP_i計(jì)算產(chǎn)狀，除去因不符合條件而沒(méi)有計(jì)算產(chǎn)狀值的交點(diǎn)集，最后得到Att＝Att₁∪Att₂∪…∪Att₂₆。考慮位置重復(fù)且產(chǎn)狀值一致的產(chǎn)狀點(diǎn)為同一個(gè)，則共包含產(chǎn)狀點(diǎn)151個(gè)。其中，使用三點(diǎn)法計(jì)算得到的有132個(gè)，其三點(diǎn)組合形式如圖7所示；使用相鄰等高線法計(jì)算得到的有19個(gè)，其四點(diǎn)組合形式如圖8所示。通過(guò)對(duì)錯(cuò)誤的產(chǎn)狀計(jì)算值進(jìn)行剔除，得到114個(gè)產(chǎn)狀計(jì)算值(如圖9所示)。

(三)產(chǎn)狀數(shù)據(jù)符號(hào)化

加載實(shí)測(cè)產(chǎn)狀數(shù)據(jù)與計(jì)算產(chǎn)狀數(shù)據(jù)進(jìn)行對(duì)比(如圖10所示)，實(shí)測(cè)數(shù)據(jù)包含32個(gè)實(shí)測(cè)產(chǎn)狀值。取實(shí)測(cè)值300m以內(nèi)的計(jì)算值，或與實(shí)測(cè)值在同一巖層界線上的計(jì)算值與實(shí)測(cè)值作比較，對(duì)比結(jié)果如下表：

基于上表計(jì)算產(chǎn)狀值與實(shí)測(cè)產(chǎn)狀值的誤差值，將傾向、傾角劃分為“正確”、“尚可”和“錯(cuò)誤”三個(gè)等級(jí)，統(tǒng)計(jì)其數(shù)量。

通過(guò)32個(gè)產(chǎn)狀實(shí)測(cè)值與其附近的計(jì)算值的對(duì)比，存在傾向或傾角計(jì)算結(jié)果誤差較大的產(chǎn)狀點(diǎn)有5個(gè)，傾向傾角完全計(jì)算錯(cuò)誤的點(diǎn)位0個(gè)，其余數(shù)據(jù)誤差在可接受范圍內(nèi)，總體正確率為84.4％。