鐵路大范圍復(fù)雜地質(zhì)區(qū)域環(huán)境三維建模方法研究

呂希奎，白嬌嬌，莊建杰，聶良濤，高崇文

(1.石家莊鐵道大學(xué) 交通運(yùn)輸學(xué)院，河北 石家莊 050043；2.河北雄安軌道快線(xiàn)有限責(zé)任公司，河北 保定 071700；3.中國(guó)鐵路設(shè)計(jì)集團(tuán)有限公司，天津 300308)

山區(qū)普遍構(gòu)造活動(dòng)強(qiáng)烈,工程地質(zhì)條件復(fù)雜,鐵路選線(xiàn)設(shè)計(jì)面臨著巨大難度。如何在地質(zhì)條件復(fù)雜的山區(qū)環(huán)境中科學(xué)合理的規(guī)劃鐵路線(xiàn)位,將鐵路的地質(zhì)災(zāi)害問(wèn)題解決在成災(zāi)之前,已成為國(guó)內(nèi)外學(xué)者關(guān)注的焦點(diǎn)問(wèn)題[1]。全面整合、整理和分析地質(zhì)信息數(shù)據(jù),有效利用地質(zhì)信息進(jìn)行選線(xiàn),解決復(fù)雜地質(zhì)區(qū)域選線(xiàn)難題,研究人員嘗試對(duì)二維地質(zhì)信息進(jìn)行三維可視化建模。但復(fù)雜山區(qū)的地質(zhì)災(zāi)害和地質(zhì)構(gòu)造的復(fù)雜性、多樣性、地質(zhì)數(shù)據(jù)的稀疏性和地質(zhì)數(shù)據(jù)的海量特性,使得鐵路選線(xiàn)大范圍復(fù)雜地質(zhì)環(huán)境的三維建模一直未能真正實(shí)現(xiàn),成為難點(diǎn)問(wèn)題。針對(duì)三維地質(zhì)建模,許多學(xué)者都進(jìn)行了研究。李健等[2]利用鉆孔與地質(zhì)剖面資料進(jìn)行數(shù)據(jù)轉(zhuǎn)換融合,通過(guò)插值算法加密地層矢量點(diǎn),提出一種基于規(guī)則體元分裂的三維地質(zhì)建模方法。文獻(xiàn)[3-5]結(jié)合ArcGIS、SketchUp等三維建模軟件,實(shí)現(xiàn)一種半自動(dòng)的復(fù)雜地質(zhì)體三維建模和可視化方法,并將其與三維地形環(huán)境中的不良地質(zhì)區(qū)域進(jìn)行疊加。陳兵等[6]基于BIM技術(shù)了建立了小范圍的橋址區(qū)三維地質(zhì)環(huán)境模型,用于更精準(zhǔn)地揭示橋址區(qū)的巖溶分布情況。文獻(xiàn)[7-8]基于鉆孔數(shù)據(jù),采用CATIA二次開(kāi)發(fā)技術(shù),實(shí)現(xiàn)了完成了三維地層實(shí)體的建模。Hassen等[9]基于地表地質(zhì)和高程數(shù)據(jù),使用三維隱式方法構(gòu)建了法國(guó)諾曼底Pays d’Auge高原區(qū)域尺度三維地質(zhì)模型。文獻(xiàn)[10-12]應(yīng)用DEM和鉆孔數(shù)據(jù),結(jié)合GIS平臺(tái)和三維建模軟件,分別實(shí)現(xiàn)了杜卡拉地區(qū)、摩洛哥東北部Gareb Bouareg盆地和希臘克里特島Tymbakion盆地的三維地質(zhì)模型構(gòu)建。文獻(xiàn)[13-14]利用既有二維地質(zhì)成果,通過(guò)參數(shù)化,交互式快速建立了地質(zhì)三維模型。Olierook等[15]研究出地表地質(zhì)觀測(cè)與地球物理調(diào)查數(shù)據(jù)相融合可以建立區(qū)域三維地質(zhì)模型方法。

綜上研究成果,目前三維地質(zhì)模型的構(gòu)建主要依賴(lài)大量鉆孔資料基礎(chǔ)上進(jìn)行的三維地質(zhì)建模方法,建模的范圍也以小范圍區(qū)域?yàn)橹鳌ｈF路選線(xiàn)設(shè)計(jì)的分階段勘察特性,使得鐵路預(yù)可研和可研的前期規(guī)劃階段所獲取的地質(zhì)資料往往十分匱乏,更是缺少目前三維地質(zhì)建模模型時(shí)所需的海量勘探鉆孔數(shù)據(jù)[16]這使得前期選線(xiàn)規(guī)劃階段的三維地質(zhì)建模面臨著難題。在前期規(guī)劃設(shè)計(jì)階段,衛(wèi)星遙感是獲取地質(zhì)信息主要的應(yīng)用手段[17],因此,設(shè)計(jì)人員在預(yù)可研階段也主要通過(guò)二維柵格地質(zhì)遙感解譯圖進(jìn)行粗略的線(xiàn)路規(guī)劃方案設(shè)計(jì),但在進(jìn)行地質(zhì)條件分析及方案比選時(shí)缺少直觀性,導(dǎo)致預(yù)可研階段設(shè)計(jì)的方案往往與最終設(shè)計(jì)方案存在較大的偏差。本文針對(duì)鐵路選線(xiàn)前期規(guī)劃設(shè)計(jì)階段三維地質(zhì)建模難題,基于空間數(shù)據(jù)三維可視化建模技術(shù),研究鐵路大范圍復(fù)雜地質(zhì)區(qū)域環(huán)境整體三維可視化方法,解決以下關(guān)鍵問(wèn)題:①選線(xiàn)整體地質(zhì)區(qū)域環(huán)境海量三維地質(zhì)建模數(shù)據(jù)處理和矢量化地質(zhì)對(duì)象生成算法;②解決復(fù)雜山區(qū)大范圍所有地質(zhì)信息一次性整體三維建模和可視化關(guān)鍵問(wèn)題;③實(shí)現(xiàn)在海量三維地質(zhì)區(qū)域環(huán)境中的選線(xiàn)應(yīng)用方法。實(shí)現(xiàn)充分展示整個(gè)鐵路區(qū)域的三維地質(zhì)環(huán)境,實(shí)現(xiàn)三維地質(zhì)區(qū)域環(huán)境和三維地理信息的統(tǒng)一集成和應(yīng)用,為預(yù)可研和可研階段復(fù)雜山區(qū)科學(xué)合理的路規(guī)鐵劃線(xiàn)位提供解決方案。

1 基于遙感地質(zhì)解譯圖的地質(zhì)對(duì)象空間數(shù)據(jù)生成方法

1.1 遙感解譯地質(zhì)對(duì)象的分類(lèi)處理

遙感地質(zhì)解譯圖中存儲(chǔ)著多種地質(zhì)體、不良地質(zhì)和地質(zhì)現(xiàn)象等大量的地質(zhì)信息,如圖1所示。

圖1 地質(zhì)遙感解譯圖示例

這些地質(zhì)信息以柵格圖像格式在地質(zhì)解譯圖上進(jìn)行表達(dá)。解譯的地質(zhì)類(lèi)型可在地質(zhì)解譯圖圖例中獲取,見(jiàn)圖2。

圖2 地質(zhì)解譯圖圖例

為更好地為三維數(shù)字化選線(xiàn)設(shè)計(jì)服務(wù),需要將柵格化的地質(zhì)信息提取并轉(zhuǎn)化為矢量化信息才能被計(jì)算機(jī)識(shí)別。對(duì)于各類(lèi)地質(zhì)對(duì)象的邊界信息,利用遙感解譯軟件AOI功能得到各解譯地質(zhì)對(duì)象的二維邊界坐標(biāo)數(shù)據(jù)(x,y),為實(shí)現(xiàn)三維建模則需要進(jìn)行三維坐標(biāo)(x,y,z)的計(jì)算處理。

依據(jù)各類(lèi)地質(zhì)對(duì)象在遙感解譯圖上的標(biāo)識(shí)形式,將圖1中的地質(zhì)類(lèi)型劃分為點(diǎn)狀、線(xiàn)狀和面狀三類(lèi),詳細(xì)劃分見(jiàn)表1。

表1 地質(zhì)類(lèi)型分類(lèi)

1.2 點(diǎn)狀空間地質(zhì)對(duì)象生成

點(diǎn)狀范圍地質(zhì)對(duì)象可用三維點(diǎn)數(shù)據(jù)(x,y,z)格式表示,利用解譯的地質(zhì)對(duì)象數(shù)據(jù)二維邊界坐標(biāo)數(shù)據(jù),結(jié)合數(shù)字高程模型,以邊界點(diǎn)Pj(xj,yj)與內(nèi)插點(diǎn)P(x,y)之間距離dj平方的倒數(shù)作為權(quán)重ωj,基于Pj的高程Zj,采用加權(quán)平均法內(nèi)插計(jì)算邊界坐標(biāo)點(diǎn)高程Zp,從而完成三維坐標(biāo)點(diǎn)數(shù)據(jù)的生成,計(jì)算式為

(1)

利用此方法可確定古火山口、溫泉、熱泉等點(diǎn)狀地質(zhì)對(duì)象三維坐標(biāo)數(shù)據(jù),即可定位其空間位置。

1.3 線(xiàn)狀空間地質(zhì)對(duì)象生成

線(xiàn)狀地質(zhì)對(duì)象實(shí)體主要是指在幾何形態(tài)上呈線(xiàn)狀的地質(zhì)空間實(shí)體[18],其空間邊界由一系列有序的控制點(diǎn)連接而成,如斷層、斷裂等。拉格朗日插值曲線(xiàn)原理是從一系列控制點(diǎn)生成一系列點(diǎn)構(gòu)成的插值曲線(xiàn),因此,為更真實(shí)地在三維空間中表達(dá)斷層等具有不同走向或形狀的線(xiàn)性地質(zhì)對(duì)象,對(duì)其控制點(diǎn)集進(jìn)行拉格朗日插值處理,定義如下:

線(xiàn)狀地質(zhì)對(duì)象三維控制點(diǎn)坐標(biāo)集Qi(xi,yi,zi),i=0,1,…,n,取其二維平面坐標(biāo)集構(gòu)造n個(gè)插值基函數(shù),利用插值基函數(shù)和已知節(jié)點(diǎn)坐標(biāo)求得插值多項(xiàng)式,再使用插值多項(xiàng)式獲得插值點(diǎn)坐標(biāo)。

假定二維平面坐標(biāo)集中的點(diǎn)滿(mǎn)足條件

(2)

為滿(mǎn)足條件式(2),構(gòu)建n次插值基函數(shù)為

i=0,1,2,…,n

(3)

則滿(mǎn)足條件式(2)的插值多項(xiàng)式為

(4)

通過(guò)求解線(xiàn)性方程組(4),解得Ln(xp),以此求得插值點(diǎn)的二維平面坐標(biāo),插值點(diǎn)Pi(i=0,1,2,…,n)的高程Zi通過(guò)式(1)求解獲得。最終連接空間三維型值點(diǎn)Pi(xi,yi,zi),即可實(shí)現(xiàn)線(xiàn)狀地質(zhì)對(duì)象在三維空間上的繪制,見(jiàn)圖3。其中,紅色點(diǎn)為地質(zhì)對(duì)象原有型值點(diǎn),藍(lán)色點(diǎn)為拉格朗日插值點(diǎn)。繪制的線(xiàn)狀地質(zhì)對(duì)象,見(jiàn)圖4。

圖4 線(xiàn)狀地質(zhì)對(duì)象繪制效果

1.4 面狀空間地質(zhì)對(duì)象生成

滑坡、泥石流、巖土層界線(xiàn)等特殊巖土體都屬于面狀地質(zhì)實(shí)體,這些地質(zhì)對(duì)象在空間分布上呈封閉多邊形。面狀地質(zhì)實(shí)體邊界由空間型值點(diǎn)序列組成,采用分段三次Hermite插值法對(duì)面狀地質(zhì)對(duì)象邊界曲線(xiàn)進(jìn)行插值,可同時(shí)滿(mǎn)足邊界曲線(xiàn)的連續(xù)性和光滑性要求,能夠很好地描述面狀地質(zhì)對(duì)象,構(gòu)造過(guò)程如下所述。

面狀地質(zhì)對(duì)象邊界函數(shù)構(gòu)造方法。設(shè)有P1,P2,…,Pn共n個(gè)空間型值點(diǎn),Pi(i=0,1,…,n)的型值點(diǎn)坐標(biāo)為(xi,yi),相鄰兩點(diǎn)之間通過(guò)三次Hermite曲線(xiàn)段連接,插值多項(xiàng)式H(x)可計(jì)算為

(5)

插值多項(xiàng)式求解。

Hi(x),x∈[xi-1,xi]滿(mǎn)足如下條件:

Hi(xi-1)=f(xi-1)=yi-1Hi(xi)=f(xi)=yi

Hi′(xi-1)=f′(xi-1)=y′i-1Hi′(xi)=f′(xi)=y′i

通過(guò)φi-1(x)、φi(x)、ψi-1(x)、ψi(x)四個(gè)三次多項(xiàng)式作為基函數(shù)構(gòu)造Hi(x),得

Hi(x)=yi-1φi-1(x)+yiφi(x)+y′i-1ψi-1(x)+y′iψi(x)

(6)

則由插值條件可以分別給出基函數(shù)滿(mǎn)足的條件為

令hi=xi-xi-1,解此方程組得

由此得到

(7)

這樣便求出滿(mǎn)足面狀地質(zhì)對(duì)象邊界的分段三次Hermite插值多項(xiàng)式。將面狀地質(zhì)對(duì)象的邊界輪廓等分為m個(gè)區(qū)間,對(duì)每個(gè)區(qū)間分別構(gòu)造三次Hermite插值多項(xiàng)式,利用方程(7)計(jì)算各插值點(diǎn)的坐標(biāo)并添加至原型值點(diǎn)集,然后按式(1)的加權(quán)平均法計(jì)算各型值點(diǎn)的高程值,最終在三維空間將該區(qū)間的點(diǎn)依次連接,繪制出面狀地質(zhì)對(duì)象,見(jiàn)圖5。

圖5 面狀地質(zhì)對(duì)象繪制效果

2 矢量化地質(zhì)對(duì)象三維模型創(chuàng)建方法

2.1 地質(zhì)對(duì)象矢量化模型創(chuàng)建

以插值處理后的地質(zhì)對(duì)象空間點(diǎn)集為數(shù)據(jù)源,基于CADC++接口,開(kāi)發(fā)地質(zhì)遙感解譯圖矢量化地質(zhì)對(duì)象三維模型創(chuàng)建模塊,實(shí)現(xiàn)點(diǎn)、線(xiàn)和面狀矢量化地質(zhì)對(duì)象模型的自動(dòng)創(chuàng)建。生成的矢量化地質(zhì)對(duì)象模型見(jiàn)圖6。

圖6 生成的矢量模式地質(zhì)對(duì)象模型

2.2 地質(zhì)對(duì)象矢量化顏色填充模式模型創(chuàng)建

為更充分展示各地質(zhì)對(duì)象模型,在矢量化基礎(chǔ)上實(shí)現(xiàn)顏色填充模式建模。

(8)

式中:m為地質(zhì)對(duì)象i的邊界數(shù)據(jù)點(diǎn)數(shù)量。

建立色系,其中線(xiàn)狀質(zhì)對(duì)象實(shí)測(cè)斷裂和推薦斷裂為紅色,解譯斷裂為藍(lán)色,標(biāo)注文字為黑色。

對(duì)面狀地質(zhì)對(duì)象建立RGB填充色系,顏色由淺到深,見(jiàn)表2。

表2 面狀地質(zhì)對(duì)象填充顏色色系

統(tǒng)計(jì)各類(lèi)地質(zhì)對(duì)象所占范圍面積總和,按照面積由大到小順序排列,進(jìn)行填充顏色確定,所占比例越高,則其配色越淺,以達(dá)到更為清晰的地質(zhì)展現(xiàn)效果。

最終得到的矢量化地質(zhì)填充模式建模效果可直觀地展現(xiàn)各類(lèi)地質(zhì)對(duì)象的分布情況,見(jiàn)圖7。圖7(a)為原始地質(zhì)遙感解譯柵格圖,7(b)為對(duì)應(yīng)的矢量化填充模式地質(zhì)建模效果。

圖7 填充模式的地質(zhì)建模效果

3 地質(zhì)信息三維空間定位算法

3.1 基于七參數(shù)的坐標(biāo)系轉(zhuǎn)換方法

為實(shí)現(xiàn)大范圍矢量化地質(zhì)模型在3DGIS中集成和展示,需要解決二者的坐標(biāo)系匹配問(wèn)題,本文采用基于七參數(shù)的坐標(biāo)系轉(zhuǎn)換方法(以下簡(jiǎn)稱(chēng)七參數(shù)法)解決,其基于橢球間三維轉(zhuǎn)換,是由兩個(gè)坐標(biāo)系之間的平移值Δx0、y2平移值Δy0、z2平移值Δz0、x軸旋轉(zhuǎn)角度ωx、y軸旋轉(zhuǎn)角度ωy、z軸旋轉(zhuǎn)角度ωz,縮放尺度比m組成,數(shù)學(xué)模型為

(9)

七參數(shù)法坐標(biāo)轉(zhuǎn)換見(jiàn)圖8。由七參數(shù)法數(shù)學(xué)模型可知,要求解七參數(shù)最少需要知道原坐標(biāo)系與目標(biāo)坐標(biāo)系中位置相對(duì)應(yīng)的三組控制點(diǎn)。可將目標(biāo)坐標(biāo)系中的三維地形坐標(biāo)賦予地質(zhì)矢量模型中對(duì)應(yīng)的控制點(diǎn),計(jì)算出坐標(biāo)轉(zhuǎn)換七參數(shù),進(jìn)而計(jì)算得到地質(zhì)矢量模型中各點(diǎn)的三維坐標(biāo),實(shí)現(xiàn)地質(zhì)矢量模型在3DGIS環(huán)境的集成。

圖8 七參數(shù)法示意

3.2 控制點(diǎn)選取策略研究

七參數(shù)法中的控制點(diǎn)選取精確與否是影響匹配精度的主要因素,其誤差主要體現(xiàn)在空間偏移距離與俯瞰圖像匹配度兩方面。本文擬定控制點(diǎn)的選取范圍參數(shù)(控制點(diǎn)三點(diǎn)所圍成的三角形區(qū)域面積與正射影像圖(DOM)影像整體面積的比值r)與圖像匹配精度參數(shù)(兩個(gè)標(biāo)記區(qū)域的空間距離誤差L、原三維地理環(huán)境的標(biāo)記區(qū)域圖像(基準(zhǔn)圖)與地質(zhì)矢量模型上相對(duì)應(yīng)的標(biāo)記區(qū)域圖像(匹配圖)各灰度狀態(tài)比值的均值A(chǔ)vg),研究通過(guò)對(duì)不同范圍內(nèi)的控制點(diǎn)組的結(jié)果進(jìn)行對(duì)比,得到兩類(lèi)參數(shù)之間的相關(guān)關(guān)系,從而確定能夠得到高匹配度結(jié)果的控制點(diǎn)選取策略。

3.2.1 基于灰度直方圖的圖像匹配度分析方法

本文以冰磧湖為例,將部分冰磧湖區(qū)域設(shè)定為標(biāo)記區(qū)域,其中紅色邊界線(xiàn)為冰磧湖邊界輪廓,黃色邊界為矢量化模型邊界。通過(guò)調(diào)整地質(zhì)矢量模型的透明度,將地形覆蓋前后的標(biāo)記區(qū)域范圍進(jìn)行對(duì)比,利用圖像處理與分離技術(shù),得到基準(zhǔn)圖與匹配圖,圖像分離過(guò)程見(jiàn)圖9。

圖9 標(biāo)記區(qū)域圖像分離

將圖像分離完畢后,利用計(jì)算機(jī)獲取基準(zhǔn)圖與匹配圖的灰度直方圖,將基準(zhǔn)圖H1與匹配圖H2的直方圖相似度數(shù)值d作為評(píng)價(jià)匹配精度的指標(biāo)之一。采用OpenCV計(jì)算直方圖相似度的數(shù)值方法如下:

采用相關(guān)對(duì)比標(biāo)準(zhǔn)(CV_COMP_CORREL)計(jì)算d(H1,H2),d(H1,H2)越接近1,表示兩個(gè)圖像匹配度越高,計(jì)算式為

計(jì)算流程如下:

1)裝載基準(zhǔn)圖像和匹配圖像進(jìn)行對(duì)比。

2)轉(zhuǎn)換圖像為HSV圖像格式。

3)計(jì)算所有圖像的H-S直方圖并歸一化。

4)對(duì)比基準(zhǔn)圖像與匹配圖像的直方圖。

5)計(jì)算直方圖相似度數(shù)值。

若完全匹配,則兩邊界完全重合,直方圖相似度數(shù)值為1,圖10 中則具有一定的偏差。標(biāo)記區(qū)域空間距離誤差L越小,則重合度越高。

圖10 兩個(gè)多邊形區(qū)域

3.2.2 標(biāo)記區(qū)域空間距離誤差計(jì)算方法

Hausdorff距離是一種基于輪廓的匹配方法,用來(lái)衡量?jī)蓚€(gè)點(diǎn)集間的距離[19]。設(shè)A={a1,a2,…,an},ai(xi,yi)∈R2,i=1,2,…,n;B={b1,b2,…,bm},bj(xj,yj)∈R2,j=1,2,…,m為兩個(gè)多邊形區(qū)域有限集,輪廓邊緣由連續(xù)的點(diǎn)表示,見(jiàn)圖10。

A,B之間的雙向Hausdorff距離定義為

H(A,B)=max{h(A,B),h(B,A)}

(11)

式中:‖·‖為點(diǎn)集A和B的歐氏距離。

Hausdorff距離H(A,B)取h(A,B)和h(B,A)的最大值,H(A,B)越小,A,B點(diǎn)集之間匹配程度越高[20]。因此,Hausdorff距離匹配方法合適計(jì)算標(biāo)記區(qū)域與疊加區(qū)域的距離L值,將L取值為雙向Hausdorff距離。

取D(a,b)為a與b之間的歐幾里得距離,采用窮舉算法計(jì)算集合A和集合B兩個(gè)點(diǎn)集距離,算法流程見(jiàn)圖11。

3.2.3 控制點(diǎn)的選取策略

通過(guò)五組控制點(diǎn)組得到的控制點(diǎn)選取范圍參數(shù)與匹配度參數(shù)見(jiàn)表3。由表3可以看出,控制點(diǎn)三點(diǎn)所圍成的三角形區(qū)域面積與DOM影像整體面積的比值r,影響著圖像直方圖相似度系數(shù)d與標(biāo)記區(qū)域空間距離誤差L:r越大,L越小,d 越大,圖像匹配度越高。

表3 控制點(diǎn)選取范圍參數(shù)與匹配度參數(shù)

為進(jìn)一步研究?jī)深?lèi)參數(shù)的關(guān)系,采用Matlab對(duì)表3中的控制點(diǎn)數(shù)據(jù)擬合,得到了在95%的置信區(qū)間下,控制點(diǎn)選取范圍參數(shù)與圖像匹配參數(shù)的函數(shù)關(guān)系表達(dá)式與關(guān)系圖像,見(jiàn)圖12和圖13。根據(jù)圖像匹配度要求,可反算出控制點(diǎn)的選取范圍。

圖12 控制點(diǎn)區(qū)域占比與圖像灰度比均值擬合

圖13 控制點(diǎn)區(qū)域占比與空間距離誤差擬合

r與d的函數(shù)關(guān)系為

d=0.105 2×r1.561+73.09

r與L的函數(shù)關(guān)系為

L=6 922×r-0.014 2-6 593

從圖12和圖13中可以得出,面積比值r越大,圖像直方圖相似度值d越大,標(biāo)記區(qū)域空間距離誤差L越小,二者分別呈指數(shù)增長(zhǎng)和指數(shù)遞減關(guān)系。

3.3 地質(zhì)信息三維呈現(xiàn)方法

本文選用WGS84坐標(biāo)系(代碼為L(zhǎng)L84),將DEM與DOM影像均以光柵文件在InfraworksBIM平臺(tái)中重建三維場(chǎng)景,生成項(xiàng)目所在地真實(shí)的三維地形環(huán)境,局部地形效果見(jiàn)圖14。

圖14 局部地形效果

采用最優(yōu)控制點(diǎn)的選取策略進(jìn)行控制點(diǎn)選取和空間定位,采用地形覆蓋形式,將整體矢量化地質(zhì)模型在Infraworks中集成,實(shí)現(xiàn)大范圍復(fù)雜地質(zhì)信息一次性整體三維建模和可視化,圖15為示例控制點(diǎn)組,圖16為實(shí)現(xiàn)方法流程。

圖15 控制點(diǎn)組選取

圖16 整體三維建模和可視化方法實(shí)現(xiàn)流程

整體一次性的大范圍地質(zhì)區(qū)域環(huán)境三維建模方法,能夠快速建立選線(xiàn)區(qū)域完整的復(fù)雜地質(zhì)環(huán)境,能夠更加清晰直觀地對(duì)地質(zhì)情況進(jìn)行展示。在建立的三維地質(zhì)區(qū)域環(huán)境中,可導(dǎo)入已有設(shè)計(jì)方案。為更直觀地展示設(shè)計(jì)方案三維效果,將Revit參數(shù)化生成的線(xiàn)路BIM模型沿線(xiàn)路中心線(xiàn)進(jìn)行空間定位,得到在大范圍三維地質(zhì)區(qū)域環(huán)境中設(shè)計(jì)方案BIM模型集成展示。同時(shí)可更清晰直觀地了解和查詢(xún)線(xiàn)路方案周邊區(qū)域的地質(zhì)信息。

4 結(jié)論

1)以衛(wèi)星遙感地質(zhì)解譯影像為主要數(shù)據(jù)源,實(shí)現(xiàn)了基于遙感地質(zhì)解譯圖的矢量化信息提取與處理方法。采用拉格朗日插值和三次Hermite曲線(xiàn)插值算法,實(shí)現(xiàn)了數(shù)據(jù)驅(qū)動(dòng)的點(diǎn)狀、線(xiàn)狀和面狀矢量化地質(zhì)對(duì)象三維模型創(chuàng)建方法。

2)以插值處理后的地質(zhì)對(duì)象空間點(diǎn)集為數(shù)據(jù)源,對(duì)各類(lèi)地質(zhì)對(duì)象所占范圍面積總和比進(jìn)行配色選取,實(shí)現(xiàn)了顏色填充模式的矢量化地質(zhì)環(huán)境建模效果。

3)以面積比值、雙向Hausdorff距離的空間距離誤差和灰度狀態(tài)比值的均值,建立了基于灰度直方圖的圖像匹配度分析方法和最優(yōu)控制點(diǎn)選取策略,實(shí)現(xiàn)了地質(zhì)信息三維空間定位算法。

4)基于七參數(shù)法以地形覆蓋形式實(shí)現(xiàn)了完整矢量化地質(zhì)模型在InfraworksBIM平臺(tái)的空間定位、集成和三維展示,實(shí)現(xiàn)了大范圍復(fù)雜地質(zhì)區(qū)域環(huán)境的三維可視化建模。

5)為進(jìn)一步深入應(yīng)用,將設(shè)計(jì)方案和鐵路線(xiàn)路BIM設(shè)計(jì)結(jié)果在InfraworksBIM平臺(tái)集成,實(shí)現(xiàn)了在三維環(huán)境中無(wú)損展示三維線(xiàn)路方案BIM模型設(shè)計(jì)效果。研究成果為預(yù)可研階段的大范圍三維地質(zhì)建模和選線(xiàn)應(yīng)用提供了可行的解決方案。