滕 飛，周 楊

(信息工程大學(xué) 地理空間信息學(xué)院，鄭州，450052)

科學(xué)技術(shù)發(fā)展進(jìn)入計(jì)算機(jī)應(yīng)用時(shí)代后，全數(shù)字測(cè)繪得到了飛速發(fā)展，測(cè)繪生產(chǎn)作業(yè)效率達(dá)到了量級(jí)的提高，地圖數(shù)據(jù)更新迅速，4D數(shù)字測(cè)繪產(chǎn)品覆蓋世界各個(gè)角落。地圖形式由紙制地圖到數(shù)字地圖，再到三維透視景觀圖，人們對(duì)地圖的直觀性、可讀性要求越來(lái)越高，對(duì)地圖信息量的需求越來(lái)越大。但是三維透視景觀圖制作工藝要求高，地物建模工作量巨大，生產(chǎn)周期相對(duì)較長(zhǎng)，且對(duì)感興趣地表細(xì)部不具備精確可量測(cè)性，因此對(duì)于應(yīng)急測(cè)繪保障來(lái)說(shuō)，不能在第一時(shí)間完成任務(wù)并交給上級(jí)決策層進(jìn)行應(yīng)急指揮決策，就會(huì)對(duì)其他各項(xiàng)工程任務(wù)的可靠性、準(zhǔn)確性產(chǎn)生一定的影響，甚至帶來(lái)不可挽救的嚴(yán)重?fù)p失。因此，我們利用正射影像平面精度高、像幅可拼接性和與立體輔助影像組成的立體像對(duì)的高精度量測(cè)功能，對(duì)整個(gè)測(cè)區(qū)進(jìn)行三維環(huán)境工程快速構(gòu)建，從而實(shí)現(xiàn)具有可量測(cè)性三維可視化功能的應(yīng)急測(cè)繪工程設(shè)計(jì)。

一、立體正射影像對(duì)的基本原理

數(shù)字微分糾正技術(shù)是數(shù)字?jǐn)z影測(cè)量領(lǐng)域里的一項(xiàng)成熟技術(shù)，主要用于數(shù)字正射影像的生成。數(shù)字正射影像由于其既具有精確的平面位置，又保持著豐富的影像信息，因此其直觀可讀性比較強(qiáng)，但是它是二維的，不含有三維信息，不具備三維高程量測(cè)的能力。雖然將等高線與正射影像套合起來(lái)，可以部分的彌補(bǔ)這個(gè)缺點(diǎn)，但卻不可能取代通過(guò)立體觀察所獲得的立體感覺(jué)。為此，可制作一個(gè)具有人工視差的輔助影像即立體輔助片[1]，將正射影像和立體輔助影像共同構(gòu)成立體正射影像對(duì)，從而實(shí)現(xiàn)高精度的三維立體觀測(cè)環(huán)境。圖1描述了立體正射影像對(duì)的基本原理。

圖1 立體正射影像對(duì)的基本原理

如圖所示，數(shù)字高程模型DEM規(guī)則格網(wǎng)高程Z及其平面坐標(biāo)XY，左邊的圖形是制作正射影像，直接將DEM格網(wǎng)上的高程值通過(guò)共線條件方程變換到像片上去，從而獲取制作正射影像。構(gòu)造虛擬的立體模型，我們按照數(shù)字高程模型所描述的地形情況引入人工視差，最簡(jiǎn)單的方法是用投影角為α的平行光法，此人工視差的大小反映地形的起伏情況，從而形成一個(gè)立體輔助片，并與正射影像構(gòu)成立體模型。該立體模型由于是在左右像片的正射影像基礎(chǔ)上構(gòu)建的，所以正射像片和立體輔助片的同名像點(diǎn)坐標(biāo)具有只有左右視差而沒(méi)有上下視差的性質(zhì)，并且同時(shí)又分別來(lái)自立體像對(duì)的左片和右片，且比例尺相同，這就滿足了立體觀測(cè)的條件，即可對(duì)地表進(jìn)行高差和坐標(biāo)量測(cè)。

立體正射影像對(duì)通過(guò)人工引入視差的方法可以進(jìn)行立體觀察，看到與實(shí)地相似的立體幾何模型，但是它的更大的意義在于可量測(cè)性。數(shù)字高程模型DEM具有一定的高程信息，但只是對(duì)整個(gè)測(cè)區(qū)大的地勢(shì)起伏信息的描述，關(guān)于地表細(xì)部高度的信息無(wú)法全部進(jìn)行詳細(xì)表達(dá)，但這些地表細(xì)部高度都包含在立體像對(duì)模型中[3]，通過(guò)立體觀測(cè)就可以達(dá)到對(duì)高度信息量測(cè)的目的。下面將討論幾種引入視差的方法。

二、對(duì)數(shù)投影法人工視差引入的方法

采用S．Collins教授提出的對(duì)數(shù)函數(shù)的方法[2]來(lái)作為視差引入函數(shù)，可有效的解決人工視差與天然視差不一致的問(wèn)題。對(duì)數(shù)投影法引入視差的原理如圖2所示。以地面上的A點(diǎn)為例，它相對(duì)于投影面的高差為Z，該點(diǎn)的正射投影為P1，對(duì)數(shù)投影為P0，正射投影得到正射像片，對(duì)數(shù)投影得到立體輔助片。立體觀測(cè)所得到的左右視差為P。

圖2 對(duì)數(shù)投影人工視差引入法

對(duì)數(shù)投影法引入視差的函數(shù)為

式中，B為攝影基線;H為航高;Z地面點(diǎn)高程;P為高程Z對(duì)應(yīng)的視差。

則得高差為

三、立體正射影像對(duì)的高程精度評(píng)估分析

立體正射影像對(duì)的量測(cè)精度主要取決于生成的正射影像的平面精度和引入人工視差所生成的立體正射影像對(duì)的高程量測(cè)精度。根據(jù)張劍清等人的研究，立體正射影像對(duì)的高程量測(cè)精度比用來(lái)制作正射影像和立體輔助影像的數(shù)字高程模型DEM的高程精度還要高3倍左右[2]。從理論上講，如果正射影像和立體匹配片的像素大小與原始影像的像素大小相同，利用上述方法進(jìn)行坐標(biāo)量測(cè)可以獲得與用原始影像量測(cè)同樣的精度[4]。在正射立體像對(duì)中，不僅能夠看到由數(shù)字地面模型所描述的大的地勢(shì)起伏情況，而且也能夠真實(shí)的在立體環(huán)境中看到數(shù)字地面模型未采集到的許多地物碎部，如樹(shù)高、建筑物高和微型地貌等。但是原始立體像對(duì)上產(chǎn)生的視差和通過(guò)視差函數(shù)引入的人工視差其原理不同，因而對(duì)于地物碎部的高程量測(cè)應(yīng)該予以仔細(xì)分析。

如下圖所示，假設(shè)地表某地物A，其所在地形表面在小范圍內(nèi)是保持水平的，且該地物在數(shù)字高程模型DEM中未被采集高程信息。設(shè)地物A頂在左右航片上的位移在X方向的大小為d X1和d X2，位移d X1在正射影像上的位移量為d P1，位移d X2在輔助片上的位移量為d P2。由圖3的幾何關(guān)系可以得出天然視差為

即

圖3 斜平行投影法碎部高程測(cè)量[5]

即

因此，對(duì)數(shù)投影人工視差引入函數(shù)能夠很好地使人工視差和天然視差保持一致，理論嚴(yán)密。高程、高差量測(cè)精度高，能夠很好地恢復(fù)原始立體像對(duì)立體觀測(cè)的視覺(jué)效果和量測(cè)功能。因此，在制作立體輔助影像時(shí)采用對(duì)數(shù)投影的人工視差引入方法。

四、試驗(yàn)與結(jié)果

立體正射影像對(duì)制作流程如圖4所示。

具體實(shí)施步驟如下：

1)利用原始遙感影像、定向參數(shù)對(duì)和一定的地面控制點(diǎn)，對(duì)影像模型進(jìn)行空三定位;

2)利用已定位的影像模型進(jìn)行自動(dòng)匹配，生成數(shù)字高程模型DEM;

3)采用數(shù)字微分糾正的方法，利用生成的DEM對(duì)測(cè)區(qū)所有影像進(jìn)行正射糾正，生成正射影像;

4)以每個(gè)立體像對(duì)的左片為左正射影像，右片為右正射影像，并且以每個(gè)立體模型的重疊區(qū)域?yàn)橛行ц偳秴^(qū)域，在此區(qū)域內(nèi)選取共同鑲嵌線分別對(duì)左正射影像和右正射影像進(jìn)行鑲嵌;

5)利用對(duì)數(shù)投影的方法，對(duì)右正射影像進(jìn)行人工視差引入，生成立體輔助影像;

6)將左正射影像和立體輔助影像組成立體正射影像對(duì)，進(jìn)行真三維立體環(huán)境量測(cè)。

圖4 立體正射影像對(duì)制作流程

其中，在“引入人工視差”步驟之前的流程在實(shí)際生產(chǎn)作業(yè)中已經(jīng)非常成熟了，這里就制作立體輔助影像進(jìn)行著重研究，并對(duì)在立體正射影像對(duì)立體環(huán)境下和原始立體像對(duì)立體環(huán)境下高程與高差量測(cè)精度進(jìn)行誤差分析。

圖5是利用地面分辨率為1m的IKONOS衛(wèi)星影像數(shù)據(jù)及其RPC定向參數(shù)，使用自動(dòng)匹配得到的數(shù)字高程模型DEM的格網(wǎng)間距為20 m，采用數(shù)字微分糾正生成左(右)正射影像，并使用對(duì)數(shù)投影的投影函數(shù)對(duì)右正射影像進(jìn)行人工視差的引入，生成右立體輔助影像。

圖5 立體正射影像對(duì)

在數(shù)字?jǐn)z影測(cè)量軟件平臺(tái)上，對(duì)原始立體像對(duì)上15個(gè)建筑物比高進(jìn)行量測(cè)，同時(shí)在生成的立體正射影像對(duì)上對(duì)這15個(gè)建筑物比高點(diǎn)再進(jìn)行量測(cè)，量測(cè)結(jié)果及誤差見(jiàn)表1。

表1 地表建筑物比高量測(cè)結(jié)果 m

從上面的試驗(yàn)結(jié)果的精度可以看出，在一定的精度要求范圍內(nèi)立體正射影像對(duì)完全可以用來(lái)量測(cè)地物碎部(如建筑物比高、樹(shù)高等)的高度，使立體三維環(huán)境具有精確可量測(cè)性?；诹Ⅲw正射影像對(duì)的三維環(huán)境構(gòu)建技術(shù)特別是針對(duì)大區(qū)域三維環(huán)境構(gòu)建速度快、效率高，兼顧三維環(huán)境觀測(cè)直觀性和立體量測(cè)精度的可靠性，在應(yīng)急測(cè)繪任務(wù)中能夠發(fā)揮巨大的作用。

五、總結(jié)與展望

基于立體正射影像對(duì)的三維可視化環(huán)境構(gòu)建技術(shù)，相對(duì)于三維透視景觀可視化技術(shù)具有以下幾個(gè)特點(diǎn)：

1)立體正射影像對(duì)無(wú)須進(jìn)行地表地物三維建模，三維環(huán)境構(gòu)建迅速，能夠確保應(yīng)急測(cè)繪保障工作順利開(kāi)展。通過(guò)對(duì)左(右)正射影像和右(左)立體輔助影像進(jìn)行真三維立體觀測(cè)，因此具有地物高差精確可量測(cè)性。

2)由于正射影像的可拼接性，立體正射影像對(duì)可對(duì)整個(gè)測(cè)區(qū)進(jìn)行整體拼接，形成覆蓋整個(gè)測(cè)區(qū)的大區(qū)域的一個(gè)立體像對(duì)，既可以進(jìn)行地表細(xì)部精確量測(cè)，也可以針對(duì)全測(cè)區(qū)進(jìn)行整體規(guī)劃和工程設(shè)計(jì)。

3)立體正射影像對(duì)可將地表碎部信息的高度或深度等一些非GIS基本地物地貌采集數(shù)據(jù)在立體像對(duì)中顯示出來(lái)，使影像地理信息量最大化得到利用。

4)立體正射影像對(duì)的立體量測(cè)設(shè)備簡(jiǎn)單，現(xiàn)行主流的計(jì)算機(jī)配置都可以配備相應(yīng)的立體顯示和觀測(cè)設(shè)備，并且整個(gè)測(cè)量過(guò)程可由非攝影測(cè)量專業(yè)人員操作完成，其適用性比較廣泛。

基于立體正射影像對(duì)的三維可視化環(huán)境所具有的可量測(cè)性可以使攝影測(cè)量突破測(cè)繪學(xué)科的局限，廣泛的深入到其他領(lǐng)域。例如，林業(yè)部門(mén)對(duì)各類樹(shù)木的高感興趣，規(guī)劃部門(mén)對(duì)建筑物的高度感興趣，地質(zhì)部對(duì)陡坎高度或小溝深度感興趣，而這些信息在遙感攝影時(shí)已經(jīng)包含在立體像對(duì)中，這樣就可以自行在立體模型上進(jìn)行量測(cè)，從而達(dá)到各個(gè)行業(yè)對(duì)不同地物地貌高度或深度量測(cè)的需求。

[1]李德仁，王密，潘俊，等．無(wú)縫立體正射影像數(shù)據(jù)庫(kù)的概念、原理及其實(shí)現(xiàn)[J]．武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2007(11)：951-954．

[2]COLLINSSH．Stereoscopic Orthophoto Maps[J]．Canadian Surveyor，1968(22)：167-176．

[3]張劍清、潘勵(lì)、王樹(shù)根．?dāng)z影測(cè)量學(xué)[M]．武漢：武漢大學(xué)出版社．2003．

[4]張祖勛，張劍清．三維可視化工程設(shè)計(jì)的研究[J]武漢大學(xué)學(xué)報(bào)，2002(4)：337-341．

[5]王密．大型無(wú)縫影像數(shù)據(jù)庫(kù)系統(tǒng)(GeoImageDB)的研制與可量測(cè)虛擬現(xiàn)實(shí)(MVR)的可行性研究[D]．武漢：武漢大學(xué)，2001．