任克彬，趙剛，王博，黃亮

（1.鄭州大學(xué) 水利與環(huán)境學(xué)院，河南 鄭州450000；2.河南省文物建筑保護研究院，河南 鄭州450002）

0 引言

古塔作為中國五千年文明史的載體，具有重要的歷史、藝術(shù)和科學(xué)價值。河南省作為全國現(xiàn)存古塔最多的省份之一，目前有獨立凌空的磚石塔600余座［1］?，F(xiàn)存的古塔在長時間經(jīng)受戰(zhàn)爭、地震、風(fēng)雨侵蝕等人為和自然破壞后，都產(chǎn)生了不同程度的傾斜，特別是磚石結(jié)構(gòu)的古塔，由于其長細(xì)比大、重心高、自重大和強度低，對地基條件很敏感，一旦傾斜，就會造成傾斜一側(cè)地基的應(yīng)力集中，產(chǎn)生不均勻沉降，加速塔體結(jié)構(gòu)的傾斜［2－4］。古塔的傾斜是一個發(fā)展的動態(tài)過程，現(xiàn)有的檢測手段相對單一、成本高、連續(xù)監(jiān)測不足，無法提供連續(xù)系統(tǒng)的監(jiān)測數(shù)據(jù)支撐古塔傾斜的保護加固立項，古塔傾斜監(jiān)測系統(tǒng)的開發(fā)與模型驗證是促進古塔加固保護立項的基礎(chǔ)和依據(jù)。

測量古塔傾斜的方法主要是通過使用全站儀等儀器來測量古塔若干個點的坐標(biāo)，進而推算塔體或各層的傾斜角度［5－6］。由于建造歷史悠久以及受自然和人為因素影響，每層塔體的傾斜變形未必會完全一致，因此，無法確定某個點就是古塔的角點或者中點，故不能保證測量結(jié)果的準(zhǔn)確性。針對此問題，周千探討了確定古塔各層中心位置的通用方法以及古塔的傾斜、彎曲、扭曲等變形情況。對于古塔各層中心的確定，通過最小二乘法擬合出觀測點所在平面，在擬合平面上尋找到各觀測點距離的平方和最小的點作為目標(biāo)層的中心點；對于古塔的傾斜、彎曲和扭曲等變形情況，分別給出了合理的數(shù)學(xué)描述，給出的是一個擬合后的整體傾斜程度，但沒有解決各層的不同傾角問題［7］。李斌等針對實際古塔開展了傾斜測量與加固工作，并監(jiān)控了加固過程中的傾角變化，建立了測量控制網(wǎng)，但成本投入較高［8－9］。總體上說，現(xiàn)有對古塔的傾斜測量多使用全站儀、水準(zhǔn)儀等儀器測量，需要多次轉(zhuǎn)站且棱鏡需要放置在塔體的特定部位，易受周圍環(huán)境和現(xiàn)場條件制約且費時費力。三維激光測量技術(shù)作為高新技術(shù)，以無接觸形式且高精度地獲得研究對象的三維點云數(shù)據(jù)，使三維激光掃描技術(shù)的應(yīng)用領(lǐng)域深入到工程測量、變形監(jiān)測等測繪領(lǐng)域的各個分支，但對現(xiàn)場條件要求高，工作量成倍遞增且費用昂貴，因此多用于小型構(gòu)件的測量，使用三維激光掃描精確的長期監(jiān)測古塔傾斜狀態(tài)存在現(xiàn)實困難。

隨著計算機技術(shù)的高速發(fā)展，數(shù)字圖像技術(shù)的識別方法使得基于數(shù)字圖像的識別技術(shù)變?yōu)榭赡埽?0－12］。李哲等對傾斜問題的數(shù)字圖像識別提出了可行的技術(shù)路線，并取得較好的識別結(jié)果［13－16］。

由于古塔的傾斜危害較大，快速準(zhǔn)確、成本低廉地檢測古塔傾斜角度并進行長期的監(jiān)測，具有十分重要的現(xiàn)實意義，能夠為古塔的保護立項提供數(shù)據(jù)支撐。文章基于數(shù)字圖像技術(shù)開發(fā)了古塔傾斜監(jiān)測系統(tǒng)并進行模型驗證。

1 古塔傾斜測量原理

1.1 傾角測量技術(shù)

現(xiàn)有古塔的構(gòu)型多為正多面體結(jié)構(gòu)，其中每層可提取為等腰棱柱體，每層傾斜角度各異，為取得每層的傾角，可將每層投影提取成等腰梯形，如圖1所示，若能將等腰梯形的兩邊與水平線的夾角α1和α2得到，則每層的傾角可由式（1）、（2）推導(dǎo)為

聯(lián)立式（1）和式（2）可求得式（3）為

所以，傾角的測量問題可轉(zhuǎn)化為圖1中兩邊界線的傾角提取問題。

圖1 單層古塔傾角計算示意圖

1.2 圖像識別技術(shù)

邊界線傾角問題可以通過數(shù)字圖像技術(shù)加以解決。（1）將圖像中需要求傾角的部分進行剪裁；（2）將其轉(zhuǎn)化為灰度模式，變?yōu)榛叶葓D像后，只剩下一個維度，相當(dāng)于將三維向量投影為一維向量；設(shè)置閾值將灰度圖像轉(zhuǎn)變?yōu)槎祱D像?；叶葓D像的像素值為0～255，而二值圖像的像素值只有0和255。二值化過程為掃描整個圖片，確定亮度最高和最低點，分別作為1和0，從1和0中間依次按一個比較小的臺階，依次作為二值化的分界，大于最高亮度的算1，小于的算0，計算所有點實際亮度與二值化以后的亮度的方差，依次遍歷計算所有的分界點，計算所有的方差和，取方差和最小的分界點作為最終的分界點，深度高于最終分界點值的則是1，低于此值則是0，將整個圖像呈現(xiàn)出明顯的黑白效果。（3）識別二值圖像中的直線，并去除其他非直線段的干擾，對圖像進行邊緣檢測。但此步驟面臨2個問題，（1）通常來說古塔的照片圖像文件是比較大的，在處理過程中會耗時較長，（2）復(fù)雜的處理算法可能造成內(nèi)存的溢出。采用分塊處理的方法，先把一副圖像分割成為一批固定大小的子圖片，然后在處理完后進行整合，從圖像的對比來看，此方法能夠有效地去除圖像的干擾。在2個具有不同灰度值的相鄰區(qū)域之間總存在邊緣，邊界分割就是基于亮度值的不連續(xù)變化來進行的。采用Hough變換函數(shù)來檢測圖像中的直線［17］，求各直線與坐標(biāo)軸的夾角，并用顏色予以標(biāo)記，并把相應(yīng)的角度值標(biāo)注于直線旁邊。

基于上述邊界線傾角算法，采用Matlab數(shù)值軟件編制了計算程序，并通過等腰直角三角板（斜邊45°）進行圖像識別加以驗證，如圖2所示。

圖2 邊界線傾角算法驗證圖

對比顯示，從圖像中提取出的斜邊角度為45°，說明程序的計算結(jié)果是可靠的。

2 古塔傾斜測量系統(tǒng)開發(fā)

2.1 硬件開發(fā)

所取照片需滿足拍攝裝置在目標(biāo)物正面水平位置，因此，開發(fā)了成套數(shù)字取像裝置［18］，如圖3所示。

圖3 古塔傾斜測量取像系統(tǒng)圖

取像裝置由相機、支架、側(cè)盤等組成。相機固定連接在設(shè)于支架上端的托盤上，相機外部設(shè)有盒子，盒子內(nèi)設(shè)有與相機固定連接的橫軸，盒子上固定設(shè)置俯仰角測盤刻度面，橫軸上固定設(shè)有俯仰角測盤指針；盒子底部固定連接在與托盤轉(zhuǎn)動連接的縱軸；托盤上設(shè)置水平測盤刻度面，縱軸上設(shè)置水平測盤指針。通過轉(zhuǎn)動縱軸調(diào)節(jié)相機的水平轉(zhuǎn)角，通過轉(zhuǎn)動橫軸調(diào)節(jié)相機的俯仰角，對準(zhǔn)目標(biāo)物體的一個側(cè)面，即可拍攝滿足正拍要求的圖像。

2.2 軟件開發(fā)

古塔傾斜測量系統(tǒng)TTMS（Tower Tilt Measurement System），目前開發(fā)版本為 V1.1，Window為平臺軟件，.exe為可執(zhí)行性文件，由Matlab語言編寫開發(fā)，針對我國現(xiàn)有古塔的基本形式，對塔身平面不同形狀進行分項處理，如圖4所示。

圖4 TTMS 1.1界面示意圖

界面左側(cè)用來填寫待測古塔的基本信息。界面右則是根據(jù)需要測量古塔的平面形式進行選擇。如待測古塔平面是四邊形，點擊四邊形按鈕即可進入測量界面，如圖5所示。

圖5 TTMS分部界面圖

進入界面之后，鑒于不同的古塔具有不同的層數(shù)，故需要定義古塔的層數(shù)，然后點擊“系統(tǒng)初始化”按鈕，對古塔數(shù)據(jù)進行初始化處理，此時古塔各面的各層數(shù)據(jù)都被初始為零，可以對各個面進行測量。如點擊進入界面，如圖6所示。

測量過程如下：（1）確定要測量古塔的層數(shù)和位置，點擊“選擇照片并處理按鈕”，選擇目標(biāo)照片，程序讀取之后，對照片進行縮放，在合適大小的窗口，對相應(yīng)的邊界進行框選；（2）雙擊鼠標(biāo)確認(rèn)，程序進行計算，并返回結(jié)果；（3）返回上一層，依次操作可以得到其他面的數(shù)據(jù)。在分部界面點擊顯示結(jié)果，可以顯示對應(yīng)的計算結(jié)果以及相應(yīng)的傾斜模擬圖，當(dāng)所有面的數(shù)據(jù)計算完畢之后，點擊結(jié)果匯總，顯示所有計算結(jié)果，整個古塔的傾斜數(shù)據(jù)就全部計算完畢，點擊“保存結(jié)果”，形成計算結(jié)果文件。對一個古塔進行長時間的監(jiān)測，就可以得到一系列數(shù)據(jù)，從而得到古塔傾斜的發(fā)展歷程。

圖6 單面圖像的處理與計算圖

3 古塔傾斜測量系統(tǒng)模型驗證

圖7 驗證模型示意圖

為驗證古塔監(jiān)測系統(tǒng)的可行性與正確性，設(shè)計了一個模型進行試驗。制作一個四邊形塔模型，將模型放在水平面上，模擬古塔未產(chǎn)生傾斜的狀況，模型4個面的照片，如圖7所示。確定測量的拍攝點，如圖8所示，面ABCD為古塔的待測面，在古塔待測面底邊兩端點C、D各用一定長度r1為半徑做弧線，圓弧交點E1即在中垂線上。重復(fù)上述過程，取另一長度r2半徑做弧線，圓弧交點E2也在中垂線上，取E1—E2線上任一合適點作為確定的拍攝點。

圖8 拍攝點的確定圖

確定拍攝點之后，采用開發(fā)的硬件系統(tǒng)進行圖像采樣，再使用開發(fā)的軟件系統(tǒng)進行識別，結(jié)果如圖9所示。

圖9 驗證識別結(jié)果圖

由圖9可知，計算結(jié)果總體上符合較好，和實際模型的對比，見表1。

表1 計算結(jié)果與實際值的對比／%

從結(jié)果對比來看，最大誤差絕對值為0.24，出現(xiàn)在第3層的面1，并且由表1可以看出，主要的誤差均出現(xiàn)在面1的下面3層，其主要原因是模型的制作誤差。由于模型材料是泡沫板，在制作過程中的線條切割精度不足導(dǎo)致；總體上看，非理想情況下的誤差控制在正負(fù)0.05°以內(nèi)，測量精度超過99.5%，驗證了古塔監(jiān)測系統(tǒng)技術(shù)的可行性。

4 工程技術(shù)經(jīng)濟比較分析

河南省駐馬店市西平縣城東關(guān)的寶嚴(yán)寺塔，又被稱為“東關(guān)塔”，為六角七級樓閣式磚塔，高28.8 m。為進行技術(shù)經(jīng)濟比較分析，同時采用所述測量系統(tǒng)和全站儀進行測量。

經(jīng)過古塔監(jiān)測系統(tǒng)一系列如圖10所示的處理過程，可以得到最后如圖11所示各層的傾角。

由圖11可以看出，對塔體各層傾斜角度根據(jù)其各層高度進行加權(quán)求和求得其平均值為1.1256°。

全站儀網(wǎng)點測量綜合擬合結(jié)果為（1°16′），折合角度為1.266°，可以看出，兩者測量結(jié)果較好符合。但與全站儀只能給出一個整體結(jié)果不同，古塔監(jiān)測系統(tǒng)可以給出古塔在不同位置的不同傾角。由于材料非均勻性塔身各層的傾斜角度并不相同，塔身下部傾斜角度比較小，傾斜最大的是第4和5層。因此，古塔監(jiān)測系統(tǒng)計算結(jié)果全面且具有較高的準(zhǔn)確度，對于古塔保護加固更有實際價值。

圖10 寶嚴(yán)寺塔圖像識別處理過程圖

圖11 寶嚴(yán)寺塔傾斜計算結(jié)果圖／°

技術(shù)經(jīng)濟分析：根據(jù)寶嚴(yán)寺塔的現(xiàn)狀和特點，對此塔監(jiān)測周期為3年，每季度1次。將古塔監(jiān)測系統(tǒng)的方法與傳統(tǒng)全站儀測量方法綜合經(jīng)濟評價對比，見表2。

表2 綜合技術(shù)經(jīng)濟評價／萬元

從表2可以看出，采用古塔監(jiān)測系統(tǒng)的方法監(jiān)測費用較傳統(tǒng)方法降低85%以上，隨著監(jiān)測周期和頻率的提高，成本低的優(yōu)勢將更加明顯。

5 結(jié)論

通過上述研究可知：

（1）古塔傾斜傳統(tǒng)測量方法成本高，無法系統(tǒng)建立古塔傾斜長期監(jiān)測體系，嚴(yán)重影響古塔保護加固立項?；跀?shù)字圖像技術(shù)，開發(fā)古塔傾斜檢測／監(jiān)測系統(tǒng)的軟硬件設(shè)施，形成了古塔長期監(jiān)測信息系統(tǒng)；在此基礎(chǔ)上進行了模型試驗，對比結(jié)果表明，開發(fā)的系統(tǒng)對于傾斜的檢測準(zhǔn)確度達(dá)到99.5%以上。

（2）針對河南省駐馬店市西平縣城東關(guān)的寶嚴(yán)寺塔，利用古塔傾斜檢測系統(tǒng)進行了現(xiàn)場測試和傳統(tǒng)全站儀測試，結(jié)果對比表明，兩者差值在0.14°，再一次驗證了系統(tǒng)的可靠性；技術(shù)經(jīng)濟評價的結(jié)果表明：基于圖像識別的古塔傾斜檢測系統(tǒng)能夠減少成本85%以上。

［1］楊煥成.楊煥成古建筑文集［M］.北京：文物出版社，2009.

［2］張煒，徐磊.陜西省古塔現(xiàn)狀調(diào)查及研究［J］.文博，2012（2）：64－70.

［3］白占微.福建古塔文化研究［D］.福州：福建師范大學(xué)，2011.

［4］馬海霞.河南古塔景觀價值的開發(fā)利用研究［D］.福州：福建農(nóng)林大學(xué)，2014.

［5］胡志曉.古塔傾斜觀測和數(shù)據(jù)分析［J］.江蘇建筑，2011（6）：34－35，44.

［6］李金鎖.古塔傾斜變形的分析和預(yù)測［J］.電子制作，2014（8）：262－264.

［7］周千.基于數(shù)學(xué)模型的古塔變形問題研究［J］.計算機與數(shù)字工程，2014，42（8）：1346－1348，1525.

［8］李斌，馮柯.開封杞縣大云寺古塔糾偏加固方案［J］.河南大學(xué)學(xué)報（自然科學(xué)版），2013，43（4）：469－472.

［9］劉殿華，袁建力，樊華，等.虎丘塔加固工程的監(jiān)控和觀測技術(shù)［J］.土木工程學(xué)報，2004，37（7）：51－58.

［10］袁建力，劉殿華，李勝才，等.虎丘塔的傾斜控制和加固技術(shù)［J］.土木工程學(xué)報，2004，37（5）：44－49，91.

［11］曹志剛，馮仲科，龍春玲.數(shù)字化測繪技術(shù)在源影寺古磚塔測繪中的應(yīng)用［J］.測繪科學(xué)，2006，31（5）：141－143.

［12］周曉敏，孟曉林，張雪萍，等.傾斜攝影測量的城市真三維模型構(gòu)建方法［J］.測繪科學(xué)，2016，41（9）：159－163.

［13］李哲，周成傳奇，閆宇，等.低空攝影測量及三維數(shù)據(jù)處理在磚石塔測繪中的實踐應(yīng)用——以天津獨樂寺塔為例［J］.中國文化遺產(chǎn)，2016（3）：61－67.

［14］張曉娟，樊東燕.傾斜車牌圖像邊緣模糊特征識別方法研究［J］.計算機仿真，2017，34（1）：372－375.

［15］王倩文.一種快速準(zhǔn)確的傾斜圖像校正方法［J］.機電工程技術(shù)，2017，2（S2）：179－183.

［16］馮雷，耿英楠.基于hough變換的傾斜文檔圖像快速校正［J］.內(nèi)蒙古師范大學(xué)學(xué)報（自然科學(xué)漢文版），2017，46（2）：219－222.

［17］劉春閣.基于自適應(yīng)中值濾波和改進Hough變換的直線提?。跩］.測繪科學(xué)，2011，36（6）：267－269.

［18］馬清文，任克彬，李威，等.一種用于傾角測量的取像裝置：中國，ZL201420100022.0［P］.2014－07－23.