1.一種基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像系統(tǒng),其特征在于:包括,點云文件輸入模塊、點云數(shù)據(jù)管理模塊、點云數(shù)據(jù)投影處理模塊、激光雷達(dá)影像生成模塊;
所述點云文件輸入模塊將地面激光雷達(dá)獲取的帶有與隧道骨架線坐標(biāo)系一致的原始三維點云數(shù)據(jù)集RD(Nk,Ek,Hk,RLk),通過隧道的骨架線里程投影法得到點k的隧道骨架線里程Mk;然后按照Mk分段的里程長度范圍將原始三維點云數(shù)據(jù)集RD分成多個點云分段數(shù)據(jù)集NDj(Mi,Ni,Ei,Hi,RLi);其中,Nk為隧道壁返回的第k個激光點對應(yīng)的骨架線坐標(biāo)系北坐標(biāo),Ek為隧道壁返回的第k個激光點對應(yīng)的骨架線坐標(biāo)系東坐標(biāo),Hk為隧道壁返回的第k個激光點對應(yīng)的骨架線坐標(biāo)系高程坐標(biāo),RLk為第k個激光點反射回來的激光反射強度;k表示原始三維點云數(shù)據(jù)集中點的編號,j表示分段后的點云數(shù)據(jù)集的編號,i表示每個分段后的點云數(shù)據(jù)集中點的編號,NDj表示第j個分段后的點云數(shù)據(jù)集,其中Mi表示第j個分段后的點云數(shù)據(jù)集中第i個點的隧道骨架線里程,Ni表示第j個分段后的點云數(shù)據(jù)集中隧道壁返回的第i個激光點對應(yīng)的骨架線坐標(biāo)系北坐標(biāo);Ei表示第j個分段后的點云數(shù)據(jù)集中隧道壁返回的第i個激光點對應(yīng)的骨架線坐標(biāo)系東坐標(biāo);Hi表示第j個分段后的點云數(shù)據(jù)集中隧道壁返回的第i個激光點對應(yīng)的骨架線坐標(biāo)系高程坐標(biāo);RLi表示第j個分段后的點云數(shù)據(jù)集中第i個激光點反射回來的激光反射強度;
所述點云數(shù)據(jù)管理模塊按照原始文件名結(jié)合每個分段點云NDj所覆蓋的隧道骨架線里程Mi來管理點云文件輸入模塊處理好的各個點云文件NDj;所述點云數(shù)據(jù)管理模塊用于提供一個用戶界面的入口;
所述點云數(shù)據(jù)投影處理模塊根據(jù)斷面輪廓線投影法將點云文件輸入模塊中獲得的多個點云分段NDj數(shù)據(jù)集中的每個點按照其對應(yīng)的骨架線里程Mi,分別投影到該點對應(yīng)的隧道骨架線里程處的隧道二維斷面坐標(biāo)系OXY中,投影后點i的斷面坐標(biāo)系坐標(biāo)記為(Mi,Pxi,Pyi);同時記錄點i距離對應(yīng)的Mi處斷面坐標(biāo)系原點的法線長度Pzi,點i對應(yīng)的Mi處斷面坐標(biāo)系原點的法線與隧道設(shè)計輪廓線的交點DPi到該斷面坐標(biāo)系原點的法線長度DPzi和點i對應(yīng)的Mi處斷面坐標(biāo)系原點的法線和設(shè)計斷面輪廓線的交點DPi沿設(shè)計斷面輪廓線到達(dá)DP0點的斷面輪廓線最短長度Ppi,根據(jù)ΔPzi=Pzi-DPzi計算實測斷面上點i沿法線方向與設(shè)計斷面輪廓線的差異值ΔPzi,其中,DP0點為設(shè)計斷面輪廓線與斷面坐標(biāo)系Y軸正方向的交點;點云數(shù)據(jù)投影處理模塊將最后獲得的數(shù) 據(jù)結(jié)果組成骨架線斷面坐標(biāo)系數(shù)據(jù)集Pj(Mi,Pxi,Pyi,Pzi,ΔPzi,Ppi,RLi),更新到數(shù)據(jù)庫中存儲;
所述激光雷達(dá)影像生成模塊將點云數(shù)據(jù)投影處理模塊生成的每個點i的骨架線斷面坐標(biāo)系數(shù)據(jù)集Pj(Mi,Pxi,Pyi,Pzi,Ppi,RLi)轉(zhuǎn)換成影像坐標(biāo)系坐標(biāo)數(shù)據(jù)集Ij(IRi,ICi,IGi)從而獲取隧道正射影像,其中IRi為點i對應(yīng)的像素在影像中所在行數(shù),ICi為點i對應(yīng)的像素在影像中所在列數(shù),IGi為點i對應(yīng)的像素灰度;Pj和Ij之間的換算關(guān)系為ICi=(Mi-Ma)/L/RS,IRi=Ppi/RS,IGi=GS×RLi;其中,Ma為起點里程,L為每幅影像覆蓋的隧道長度,RS表示隧道正射影像的分辨率,分辨率用像素的寬度表示,單位為米,GS為激光反射強度值與影像灰度值之間的一個固定比例常數(shù)。
2.根據(jù)權(quán)利要求1所述的基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像系統(tǒng),其特征在于:還包括激光雷達(dá)平整度模塊;所述激光雷達(dá)平整度模塊根據(jù)點云文件輸入模塊通過隧道的骨架線里程投影法分段成新的點云數(shù)據(jù)集NDj(Mi,Ni,Ei,Hi,RLi)計算隧道內(nèi)壁點i的平整度數(shù)值Ui,并將計算出的平整度數(shù)值Ui添加到骨架線斷面坐標(biāo)系數(shù)據(jù)集Pj中,形成新的骨架線斷面坐標(biāo)系數(shù)據(jù)集Pj(Mi,Pxi,Pyi,Pzi,Ppi,RLi,Ui)。
3.根據(jù)權(quán)利要求2所述的基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像系統(tǒng),其特征在于:還包括影像成果報告模塊,所述影像成果報告模塊根據(jù)激光雷達(dá)平整度模塊生成的平整度數(shù)值Ui的數(shù)值結(jié)合不同的范圍區(qū)間,將不同平整度數(shù)值范圍內(nèi)的點采用不同色彩或者不同形狀表示,將點i的平整度數(shù)值Ui以半透明彩色分析圖的方式疊加到激光雷達(dá)隧道正射影像中對應(yīng)位置上;進行隧道里程分段的斷面變形量、平整度,是否有滲水或裂縫分析。
4.根據(jù)權(quán)利要求1所述的基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像系統(tǒng),其特征在于:所述骨架線里程投影法是將點云RD中的每個點的隧道骨架坐標(biāo)系三維坐標(biāo)(Nk,Ek,Hk)投影到隧道骨架線的垂足上,對應(yīng)垂足點距離骨架線起點里程的距離為點k所對應(yīng)的隧道骨架線里程Mk。
5.根據(jù)權(quán)利要求1所述的基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像系統(tǒng),其特征在于:所述Mk分段的里程長度范圍是根據(jù)分段后每個點云單元的數(shù)據(jù)量大小與計算機硬件的最優(yōu)分布處理能力所對應(yīng)的最大數(shù)據(jù)塊粒度大小匹配的原則來選取的。6、一種采用權(quán)利要求1所述的基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像系統(tǒng)的隧道正射影像獲取方法,其特征在于:包括以下步驟:
步驟1:點云文件輸入模塊采集每個原始點云文件數(shù)據(jù)RD(Nk,Ek,Hk,RLk),其中(Nk,Ek,Hk)為點k在隧道骨架線控制坐標(biāo)系下坐標(biāo),k表示點云中每個點的標(biāo)號;
步驟2:讀入骨架線三維設(shè)計數(shù)據(jù)中的各個幾何線型,分段投影到平面中的幾何線型要素分段YSg(Mg,Ng,Eg,Tg,Rg,Lg),Mg為幾何線型要素起點的里程,Ng為在隧道骨架線控制坐標(biāo)系下的北坐標(biāo),Eg為在隧道骨架線控制坐標(biāo)系下的東坐標(biāo),Tg表示要素類型,Rg表示幾何線型要素曲率半徑值,Lg表示幾何線型長度;其中,要素類型為直線、緩和曲線和圓;
步驟3:點云文件輸入模塊根據(jù)骨架線里程投影法計算出原始點云RD中的每個點k投影到隧道骨架線垂足的里程Mk;
步驟4:點云文件輸入模塊根據(jù)每個點的Mk數(shù)值和投影參數(shù)設(shè)置中骨架線里程投影分段的里程長度范圍,將每個讀入的原始三維點云數(shù)據(jù)集RD(Nk,Ek,Hk,RLk)進行分割,形成多個點云分段數(shù)據(jù)集NDj(Mi,Ni,Ei,Hi,RLi),并在點云數(shù)據(jù)管理模塊中顯示拆分后的新點云數(shù)據(jù)集NDj(Mi,Ni,Ei,Hi,RLi);
步驟5:點云數(shù)據(jù)投影處理模塊按照點云數(shù)據(jù)管理模塊界面中用戶選擇的待處理點云數(shù)據(jù)集序列,批量讀入點云數(shù)據(jù)集NDj(Mi,Ni,Ei,Hi,RLi);
步驟6:將步驟5中讀取的每個點云NDj中的每個點i逐個按照骨架線里程投影法投影到相應(yīng)里程對應(yīng)的隧道斷面二維坐標(biāo)系,投影后點i的斷面坐標(biāo)系記為(Mi,Pxi,Pyi),同時記錄點i距離對應(yīng)的Mi處斷面坐標(biāo)系原點的法線長度Pzi,點i對應(yīng)的Mi處斷面坐標(biāo)系原點的法線與隧道設(shè)計輪廓線的交點DPi到該斷面坐標(biāo)系原點的法線長度DPzi和點i對應(yīng)的Mi處斷面坐標(biāo)系原點的法線和設(shè)計斷面輪廓線的交點DPi沿設(shè)計斷面輪廓線到達(dá)DP0點的斷面輪廓線最短長度Ppi,根據(jù)ΔPzi=Pzi-DPzi計算實測斷面上點i沿法線方向與設(shè)計斷面輪廓線的差異值ΔPzi,其中,DP0點為設(shè)計斷面輪廓線與斷面坐標(biāo)系Y軸正方向的交點;將最后獲得的數(shù)據(jù)結(jié)果組成骨架線斷面坐標(biāo)系數(shù)據(jù)集Pj(Mi,Pxi,Pyi,Pzi,ΔPzi,Ppi,RLi),更新到數(shù)據(jù)庫中存儲;
步驟7:點云數(shù)據(jù)投影處理模塊根據(jù)點云數(shù)據(jù)集NDj(Mi,Ni,Ei,Hi,RLi)中點i的三維空間坐標(biāo)(Ni,Ei,Hi)、激光反射強度RLi和骨架線斷面坐標(biāo)系數(shù)據(jù)集Pj(Mi,Pxi,Pyi,Pzi,Ppi,RLi)建立一個數(shù)據(jù)集對應(yīng)關(guān)系,以Mi作為關(guān)聯(lián)字段和檢索關(guān)鍵字段,并更新存儲到點云數(shù)據(jù)庫中;
步驟8:對激光雷達(dá)影像生成模塊進行影像分幅的隧道里程長度和影像分辨率參數(shù) 設(shè)置,按照骨架線里程Mi為關(guān)鍵字讀入需要生成影像的點云j數(shù)據(jù)集Aj(Mi,Pxi,Pyi,Pzi,Ppi,RLi,Ni,Ei,Hi);Aj是對已有的NDj和Pj數(shù)據(jù)集進行聯(lián)合查詢后獲得的新數(shù)據(jù)集;
步驟9:激光雷達(dá)影像生成模塊采用斷面輪廓線投影法將點云數(shù)據(jù)投影處理模塊生成的個分段后的點云數(shù)據(jù)集中點i的骨架線斷面坐標(biāo)系數(shù)據(jù)集Pj(Mi,Pxi,Pyi,Pzi,Ppi,RLi)轉(zhuǎn)換成影像坐標(biāo)系坐標(biāo)數(shù)據(jù)集Ij(IRi,ICi,IGi)從而獲取隧道正射影像,其中IRi為點i對應(yīng)的像素在影像中所在行數(shù),ICi為點i對應(yīng)的像素在影像中所在列數(shù),IGi為點i對應(yīng)的像素灰度;Pj和Ij之間的換算關(guān)系為ICi=(Mi-Ma)/L/RS,IRi=Ppi/RS,IGi=GS×RLi;其中,Ma為起點里程,L為每幅影像覆蓋的隧道長度,RS表示隧道正射影像的分辨率,分辨率用像素的寬度表示,單位為米,GS為激光反射強度值與影像灰度值之間的一個固定比例常數(shù)。
6.根據(jù)權(quán)利要求5所述的基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像方法,其特征在于:以所述步驟9中獲取的隧道正射影像圖為底圖,在隧道正射影像圖上疊加半透明彩色目標(biāo)數(shù)值分析圖;或者疊加半透明不同圖例目標(biāo)數(shù)值分析圖。
7.根據(jù)權(quán)利要求5所述的基于激光雷達(dá)LIDAR點云數(shù)據(jù)獲取隧道正射影像方法,其特征在于:所述目標(biāo)數(shù)值包括斷面變形量和平整度。