李張翼，張群喜，徐 棟，任秀芬

（1.西安融達(dá)多媒體技術(shù)有限公司，陜西 西安 710077；2.陜西歷史博物館，陜西 西安 710061；3.云南省文物考古研究所，云南 昆明 650000；4.銀川西夏陵區(qū)管理處 文物管理科，寧夏 銀川 750021）

0 引 言

我國(guó)歷史悠久，擁有豐富的文化遺產(chǎn)，但目前我國(guó)正在經(jīng)歷巨大的文化生態(tài)變化。隨著經(jīng)濟(jì)全球化趨勢(shì)和現(xiàn)代化進(jìn)程的加快，很多古建筑、古遺址的整體風(fēng)貌已經(jīng)遭到破壞，甚至一些重要的物質(zhì)文化遺產(chǎn)因?yàn)樗亮魇Ш惋L(fēng)化裂變，已經(jīng)消亡或失傳。留存延續(xù)的危機(jī)降臨在了這些不可再生的物質(zhì)文化遺產(chǎn)上。為了保護(hù)這些文化遺產(chǎn)，我們需要采取比搶救性修復(fù)更有效果的預(yù)防性保護(hù)措施。近15 年來(lái)我國(guó)開(kāi)展了館藏文物保存環(huán)境監(jiān)測(cè)預(yù)警和預(yù)防性保護(hù)實(shí)踐，成效顯著，文物保護(hù)工作也已逐步從搶救性保護(hù)向預(yù)防性保護(hù)轉(zhuǎn)變，預(yù)防性保護(hù)的理念是通過(guò)控制引起文物劣化的主要影響因子（比如文物保存環(huán)境）來(lái)避免建筑遺產(chǎn)的損毀或破壞[1]。不同于以往的應(yīng)急性保護(hù)工程，文物預(yù)防性保護(hù)的重點(diǎn)在于通過(guò)信息的收集、精密勘察、風(fēng)險(xiǎn)和價(jià)值評(píng)估等方法來(lái)確定文化遺產(chǎn)面臨的風(fēng)險(xiǎn)影響因子，基于系統(tǒng)監(jiān)測(cè)和定期檢測(cè)的方法分析物質(zhì)文化遺產(chǎn)的劣化規(guī)律，在上述研究基礎(chǔ)上及時(shí)降低或者消除面臨的風(fēng)險(xiǎn)，主要采取災(zāi)害預(yù)防、日常維護(hù)、科學(xué)管理等措施，從而實(shí)現(xiàn)遺產(chǎn)的全面保護(hù)[2]。

目前變化檢測(cè)方法在遙感影像領(lǐng)域得到了越來(lái)越多的關(guān)注和發(fā)展，并在自然災(zāi)害的災(zāi)情監(jiān)測(cè)、土地利用/覆蓋動(dòng)態(tài)研究、環(huán)境變化的監(jiān)視、森林覆蓋和水域檢測(cè)以及農(nóng)作物生長(zhǎng)狀況的監(jiān)測(cè)和評(píng)估等[3-4]方面發(fā)揮了重要的作用。同時(shí)關(guān)于圖像變化檢測(cè)的方法有很多，如最大后驗(yàn)概率估計(jì)法、模糊吉布斯馬爾科夫隨機(jī)場(chǎng)模型法[5]、非下采樣Contourlet 變換和脈沖耦合神經(jīng)網(wǎng)絡(luò)結(jié)合法[6]、期望值最大化法、局部擬合模型法[7]、基于分形維數(shù)和支持向量機(jī)法[8]、基于馬爾科夫隨機(jī)場(chǎng)和水平集方法的交互式分割法[9]和主成分分析法等[10-11]。

由于變化檢測(cè)方法在文物保護(hù)監(jiān)測(cè)領(lǐng)域有很少的應(yīng)用，面對(duì)當(dāng)前文物預(yù)防性保護(hù)的緊急性和必要性，本文針對(duì)土遺址風(fēng)化裂變和水土流失問(wèn)題，考慮到遺址保存環(huán)境光照變化情況，設(shè)計(jì)了一種基于光照變化的三色差值主成分分析方法，該方法通過(guò)架設(shè)在土遺址周邊的固定攝像機(jī)拍攝的實(shí)時(shí)圖像與初始圖像，能夠識(shí)別出其風(fēng)化和水土流失部分。該方法計(jì)算成本低；結(jié)合通信和云平臺(tái)，可嵌入到預(yù)防性保護(hù)在線監(jiān)測(cè)系統(tǒng)。

1 方法流程

本文通過(guò)光照對(duì)RGB 三色影像進(jìn)行分析，為減少光照對(duì)變化檢測(cè)的干擾，在建立差值圖像過(guò)程中加入三色權(quán)重系數(shù)，結(jié)合圖像預(yù)處理、主成分分析和閾值分割等方法，設(shè)計(jì)了一種土遺址輪廓變化方法，流程如圖1 所示。

1.1 圖像預(yù)處理

雖然攝像機(jī)是架設(shè)在固定位置進(jìn)行定期拍攝的，但由于室外天氣變化和風(fēng)力影響，設(shè)備在拍攝期間會(huì)略有擺動(dòng)，因此需對(duì)目標(biāo)圖像進(jìn)行幾何校正以避免出現(xiàn)誤差。本文利用參考圖像和目標(biāo)圖像的特征點(diǎn)進(jìn)行匹配，并用RANSAC 隨機(jī)抽樣方法進(jìn)行配對(duì)篩選，計(jì)算出體現(xiàn)兩幅圖映射關(guān)系的變換矩陣，并用該矩陣對(duì)目標(biāo)圖像進(jìn)行幾何校正。

在實(shí)際監(jiān)測(cè)過(guò)程中，監(jiān)測(cè)主體為原位保存的文化遺址，為避免遺址周邊植物生長(zhǎng)、人類活動(dòng)變化等影響，須人工設(shè)定監(jiān)測(cè)區(qū)域。針對(duì)本文劃定的多邊形輪廓，利用射線法建立影像掩膜，如圖2 所示。

圖2 建立掩膜圖像

1.2 主成分分析法

主成分分析（Principal Component Analysis, PCA）是一種多維正交線性變換，是建立在統(tǒng)計(jì)特征基礎(chǔ)上的一種離散的K-L 變換[12]，被廣泛應(yīng)用于遙感圖像處理中的變化檢測(cè)任務(wù)。

主成分分析算法如下：設(shè)有向量集X={Xi,i=1, 2,...,N}∈Rn，E(X)為X的數(shù)學(xué)期望，C是X的協(xié)方差矩陣，U是C特征向量按照特征值由大到小的順序排列而成的變換矩陣，則有

這樣就構(gòu)成了主成分分析算法，其中：

差異主成分變換法和主成分差異法是當(dāng)前被廣泛應(yīng)用于變化檢測(cè)的主成分分析方法。

差異主成分變換法是先對(duì)需要對(duì)比的兩幅圖像作相差并取絕對(duì)值處理，從而得到一個(gè)差值矩陣，這個(gè)差值矩陣集中了不同時(shí)段拍攝的原兩幅圖像中絕大部分的變化信息，過(guò)濾掉了影像中相同的背景部分，再對(duì)差值矩陣進(jìn)行主成分變換。由其原理可知，變換輸出的第一分量集中了圖像的主要信息，因此差值影像進(jìn)行主成分變換后的第一分量應(yīng)該集中了原兩時(shí)相影像的主要差異信息。這個(gè)分量可以認(rèn)為是變化信息而被單獨(dú)提取出來(lái)。

主成分差異法的不同之處在于：對(duì)兩幅不同時(shí)段拍攝的圖像做差值之前，先對(duì)影像進(jìn)行主成分變換，然后取差值的絕對(duì)值為處理結(jié)果。雖然在對(duì)不同時(shí)段拍攝的圖像分別進(jìn)行主成分變換時(shí)，靠前的維度集中了圖像的主成分信息，但后面維度的差值有時(shí)也能體現(xiàn)出原始圖像的變化信息。本文經(jīng)過(guò)測(cè)試后采用第二種方法，即主成分差異法對(duì)不同時(shí)段拍攝圖像的RGB 三色通道分別進(jìn)行主成分分析后作差值，在獲得三個(gè)通道差值的絕對(duì)值后，進(jìn)行加權(quán)疊加以及后續(xù)處理。

1.3 建立差值圖像

將圖像進(jìn)行主成分分析后，圖像中不同的三色通道對(duì)應(yīng)的信息會(huì)有所不同，因而可以從不同色道的顯示情況中研究建立差值圖像。通過(guò)分析發(fā)現(xiàn)，戶外土遺址在光照變化的情況下，圖像中RGB 三色的平均值變化波動(dòng)不同，本文以西夏王陵土遺址不同天氣情況圖像為例，如圖3 所示。

圖3 不同天氣情況下西夏王陵土遺址圖像

假定以多云圖像圖3（a）為參考圖像，以圖3（b）、圖3（c）、圖3（d）為不同天氣情況下的目標(biāo)圖像，對(duì)圖像進(jìn)行重點(diǎn)區(qū)域劃定后，分別獲取重點(diǎn)區(qū)域RGB 三色通道的平均值，并將其與參考圖像圖3（a）的值作比值，獲取三色通道因光照變化而產(chǎn)生的不同變化，數(shù)據(jù)見(jiàn)表1 所列。

表1 RGB 三色通道均值比值

從表1 中可以看出，在陰雨天氣沒(méi)有光照時(shí)，R 通道和G 通道的均值比值有明顯變化，而B(niǎo) 通道則有較小浮動(dòng)；同時(shí)在出現(xiàn)風(fēng)化等裂變的情況中，由于背景圖像比重增多，且多為藍(lán)、白等冷色調(diào)，因此B 通道會(huì)有明顯變化。基于上述情況，本文設(shè)計(jì)出一種基于光照變化的三色差值圖像建立方法，根據(jù)比值減少R 通道與G 通道所屬權(quán)重，降低光照不同所帶來(lái)的影響。首先計(jì)算出兩幅圖像R 通道和G 通道的均值，再得出各通道的比值，如式（3）所示：

式中：α為求出的比值；分別代表經(jīng)過(guò)預(yù)處理、重點(diǎn)區(qū)域劃分的參考圖像和目標(biāo)圖像的單通道均值；R、G 分別代表圖像的R 通道和G 通道。在求出αR、αG后，取其與1 的差的絕對(duì)值，通過(guò)該值可以判斷光照影響強(qiáng)度，該值越大，則權(quán)重應(yīng)越小，具體計(jì)算公式為：

式中，C為常量。當(dāng)|1-α|大于C時(shí)，說(shuō)明該圖像與參考圖像相比光照變化較大，對(duì)應(yīng)須降低相關(guān)影響，須根據(jù)現(xiàn)場(chǎng)監(jiān)測(cè)主體、背景的不同對(duì)C值進(jìn)行細(xì)微調(diào)整。在本文的實(shí)驗(yàn)中，C值取0.15。

在計(jì)算出差值權(quán)重后，根據(jù)式（4）建立差值圖像：

式中：D為所建立的差值圖像矩陣；DB、DG和DR分別為通過(guò)主成分差異法形成的RGB 三通道差值結(jié)果矩陣。

1.4 閾值分割與形態(tài)學(xué)處理

為了能夠確定每個(gè)變化的像素點(diǎn)，我們需從差值圖像中確定變化閾值，基于閾值標(biāo)識(shí)出二值圖像。這一步需要不斷調(diào)整影像閾值范圍才能夠較好體現(xiàn)出正確的變化信息。本文選擇大津閾值分割方法，提取差值圖像D的變化主體。此外，差值圖像往往會(huì)出現(xiàn)很多零散雜亂的圖斑，為了更好地規(guī)整變化信息，體現(xiàn)出土遺址風(fēng)化、水土流失的規(guī)整性，本文使用形態(tài)學(xué)、篩選最小連通區(qū)域等方法對(duì)差值影像進(jìn)行處理，消除掉像元值較為零散的像素。

2 實(shí)驗(yàn)分析

2.1 研究目標(biāo)概況

本文以西夏王陵6 號(hào)陵為實(shí)驗(yàn)對(duì)象，選取不同時(shí)期和不同光照情況下的圖像進(jìn)行檢測(cè)。目前西夏王陵文化遺址保存狀況良好，沒(méi)有明顯的裂變情況。為了更好地驗(yàn)證本方法，本文對(duì)目標(biāo)圖像進(jìn)行了人工處理，模擬風(fēng)化裂變等情況進(jìn)行實(shí)驗(yàn)，如圖4 所示。

圖4 通過(guò)人工處理的目標(biāo)圖像

實(shí)驗(yàn)中本文以文獻(xiàn)[10-11]的實(shí)驗(yàn)方法作為對(duì)比，其中文獻(xiàn)[10]提出了一種主成分分析法和變換向量法相結(jié)合的變化檢測(cè)方法，文獻(xiàn)[11]提出了一種基于多波段差值和主成分分析相結(jié)合的變化檢測(cè)方法。

2.2 實(shí)驗(yàn)結(jié)果

在選取固定的參考圖像后，選擇圖4 中處理后的兩圖進(jìn)行變化檢測(cè)比對(duì)，并與文獻(xiàn)[10-11]的結(jié)果進(jìn)行對(duì)比。

圖5（a）為參考圖像，圖5（b）、圖5（f）的模擬圖像為人工修改后的目標(biāo)圖像。圖5（b）由人工修改了大塊的主體信息，并存在不同的光照強(qiáng)度；圖5（f）的主體頂端至左側(cè)有少部分人工模擬變化痕跡，并且光照強(qiáng)度有很大的變化。對(duì)比實(shí)驗(yàn)結(jié)果可以看出，本文方法在圖像大塊主體信息隨著光照變化的情況下仍能夠較好地獲取整體變化情況；并在光照情況相差較大時(shí)，也能夠減少光照對(duì)檢測(cè)結(jié)果的影響，同時(shí)識(shí)別出變化部分。相比文獻(xiàn)[10]和文獻(xiàn)[11]，本文設(shè)計(jì)的檢測(cè)方法有較好的效果。

圖5 實(shí)驗(yàn)結(jié)果

3 結(jié) 語(yǔ)

本文以原位保存的文化遺址為監(jiān)測(cè)目標(biāo)，設(shè)計(jì)了一種基于光照變化的三色差值主成分分析方法。結(jié)合圖像幾何校正、建立掩膜影像等預(yù)處理手段和大津閾值分割等核心算法，提出了一種面向文化遺址風(fēng)化裂變的變化檢測(cè)方法。該方法能夠有效減少光照對(duì)于變化檢測(cè)結(jié)果的影響，并可以通過(guò)多相機(jī)多角度的拍攝和互聯(lián)網(wǎng)可視化功能，形成較為完整的輪廓監(jiān)測(cè)系統(tǒng)，契合文物預(yù)防性保護(hù)中的定期檢測(cè)和系統(tǒng)監(jiān)測(cè)的技術(shù)要求。

本文所提出的文化遺址輪廓變化檢測(cè)方法的不足之處在于：圖像最終以二維平面的結(jié)果呈現(xiàn)，如果只以面積變化比為量化標(biāo)準(zhǔn)，無(wú)法真正體現(xiàn)遺址文物裂化狀態(tài)。因此在滿足無(wú)損監(jiān)測(cè)、易于施工的要求下，后期需針對(duì)三維配準(zhǔn)方向開(kāi)展研究。