王式太 張博宇 殷 敏

1 桂林理工大學(xué)測(cè)繪地理信息學(xué)院，桂林市雁山街319號(hào)，541006 2 廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室，桂林市雁山街319號(hào)，541006

1 動(dòng)態(tài)閾值陰影匹配算法

陰影匹配算法的原理是以依據(jù)信噪比判定標(biāo)準(zhǔn)獲得的衛(wèi)星可見性為參考基準(zhǔn)，對(duì)接收機(jī)在城市峽谷中的可能位置(搜索范圍)上依據(jù)已有的3D建筑模型和衛(wèi)星方位角、仰角等信息獲得的衛(wèi)星空間可見性進(jìn)行評(píng)分，所有可能位置中評(píng)分(可見性匹配度)最高者作為接收機(jī)定位位置。傳統(tǒng)算法定位步驟為：

1)移動(dòng)終端初始化定位。通過GNSS接收機(jī)獲取單點(diǎn)定位結(jié)果，得到精度較低的P0點(diǎn)位置信息，將P0點(diǎn)設(shè)為可能位置的搜索原點(diǎn)。

2)設(shè)定搜索區(qū)域。以P0點(diǎn)為中心設(shè)定搜索區(qū)域，通常情況下可設(shè)置為圓形或者矩形，搜索半徑可根據(jù)不同終端的定位精度設(shè)為10～50 m。

3)預(yù)測(cè)搜索區(qū)域內(nèi)每個(gè)候選位置的衛(wèi)星可見性。在每個(gè)候選位置上，通過衛(wèi)星坐標(biāo)和候選位置坐標(biāo)計(jì)算該時(shí)刻衛(wèi)星方位角和高度角，同時(shí)結(jié)合周圍建筑物三維模型判定該候選位置上的衛(wèi)星可見性。

4)依據(jù)信噪比判斷衛(wèi)星可見性。一般認(rèn)為，接收機(jī)直接接收到衛(wèi)星信號(hào)時(shí)其SNR值較高，當(dāng)接收到經(jīng)過建筑物衍射和反射的衛(wèi)星信號(hào)時(shí)，SNR值會(huì)變低。設(shè)定固定的信噪比閾值作為判斷標(biāo)準(zhǔn)，當(dāng)信噪比低于閾值時(shí)則判定衛(wèi)星不可見，否則為可見。

5)候選位置評(píng)分。對(duì)比步驟3)和4)評(píng)估衛(wèi)星可見度預(yù)測(cè)值與觀測(cè)值之間的匹配度，匹配較好的候選位置將被賦予更高的權(quán)重。利用評(píng)分方案，根據(jù)預(yù)測(cè)和觀測(cè)到的可見度，采用打分模板(圖1)對(duì)每顆衛(wèi)星進(jìn)行打分。

圖1 打分模板示意圖

每個(gè)候選位置的總得分為：

(1)

式中，fpos(j)為侯選位置j的最終得分，fsat(i,j)為衛(wèi)星i在侯選位置j的得分，n為大于截止高度角的衛(wèi)星顆數(shù)。

對(duì)于我國來說，制定促進(jìn)農(nóng)戶適應(yīng)新型農(nóng)業(yè)食品體系的政策是非常必要的。政策旨在解決農(nóng)戶參與超市供應(yīng)鏈時(shí)可能遇到的營銷壁壘，以及促進(jìn)集體行動(dòng)，從而盡可能地降低農(nóng)戶參與過程中的交易成本。

6)定位解算。由打分方案可知，在衛(wèi)星預(yù)測(cè)和觀測(cè)準(zhǔn)確的情況下，高分值所在的侯選位置更接近真實(shí)位置。高分侯選位置通常不止一個(gè)，則取坐標(biāo)均值作為陰影匹配的定位結(jié)果。

在實(shí)際運(yùn)算過程中，傳統(tǒng)方法通過衛(wèi)星信噪比經(jīng)驗(yàn)閾值來判定衛(wèi)星的可見性不夠準(zhǔn)確，因?yàn)橥瑸榭梢娀蛲瑸檎趽醯牟煌l(wèi)星，其SNR值由于信號(hào)穿過不同大氣/云層會(huì)存在差異，采用相同閾值可能會(huì)存在誤判。在高分值區(qū)域內(nèi)，如果閾值設(shè)置正確，理論上當(dāng)搜索位置與實(shí)際觀測(cè)位置一致時(shí)，預(yù)測(cè)可見性與觀測(cè)可見性將完全一致，但實(shí)際上每個(gè)歷元的候選位置中通常不存在滿分位置。據(jù)此，本文提出在每個(gè)歷元中找到可使最優(yōu)區(qū)域產(chǎn)生滿分的一組動(dòng)態(tài)閾值，閾值組隨時(shí)間繼承，直至因非滿分再進(jìn)行局部閾值改變。動(dòng)態(tài)閾值算法計(jì)算流程如圖2所示，具體步驟為：

圖2 動(dòng)態(tài)閾值計(jì)算流程

1)設(shè)定初始閾值組為k(k1=k2=…=kn)(n為觀測(cè)到的衛(wèi)星個(gè)數(shù))，采用此閾值組進(jìn)行該歷元衛(wèi)星觀測(cè)可見性判斷，與預(yù)測(cè)可見性對(duì)比獲得最高分區(qū)域。

2)在該歷元的最高分區(qū)域內(nèi)統(tǒng)計(jì)每個(gè)候選位置上預(yù)測(cè)可見性與觀測(cè)可見性不符的衛(wèi)星，設(shè)分差為d，全部不符衛(wèi)星數(shù)為m(m≥d)，統(tǒng)計(jì)m顆衛(wèi)星在最高分區(qū)域內(nèi)的不符次數(shù)。

3)選擇不符衛(wèi)星中出現(xiàn)次數(shù)由高到低的前d顆衛(wèi)星。

4)依據(jù)衛(wèi)星預(yù)測(cè)可見性對(duì)d顆衛(wèi)星依次改進(jìn)SNR閾值Si：

①依據(jù)幾何關(guān)系判斷為遮擋，SNR閾值判斷為不遮擋：

Si=SNRi+1

(2)

②依據(jù)幾何關(guān)系判斷為不遮擋，SNR閾值判斷為遮擋：

Si=SNRi-1

(3)

5)重新對(duì)搜索區(qū)域內(nèi)的候選位置進(jìn)行打分，根據(jù)最高分區(qū)域計(jì)算改正后的點(diǎn)位，將閾值組傳入下一歷元。

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

本實(shí)驗(yàn)采用某品牌智能手機(jī)作為接收機(jī)，同時(shí)使用RTK對(duì)實(shí)驗(yàn)區(qū)周圍建筑物進(jìn)行測(cè)量，包括建筑物各角點(diǎn)測(cè)量以及樓層高度測(cè)量，以此建立三維模型。實(shí)驗(yàn)區(qū)為桂林理工大學(xué)雁山校區(qū)7A、7B、6B、6C四棟樓之間的過道，峽谷總長(zhǎng)為130 m，寬度為16 m，進(jìn)行靜態(tài)模式測(cè)量。

實(shí)驗(yàn)中同時(shí)采用RTK測(cè)定坐標(biāo)。由于RTK獲得的坐標(biāo)精度遠(yuǎn)高于手機(jī)定位的坐標(biāo)精度，因此將其作為分析實(shí)驗(yàn)精度的參考值。數(shù)據(jù)采集時(shí)間為2020-05-10 14：13～14：36，采樣間隔為1 s。以3 s時(shí)刻采集的坐標(biāo)數(shù)據(jù)為例，以采集的坐標(biāo)為中心建立20 m×20 m矩陣搜索區(qū)域，候選區(qū)域邊長(zhǎng)為0.5 m，共1 600個(gè)候選區(qū)域。

圖3為整個(gè)觀測(cè)過程中部分衛(wèi)星信噪比值變化，從圖中可以看出，衛(wèi)星信噪比值在每個(gè)歷元均會(huì)發(fā)生波動(dòng)，不同衛(wèi)星信噪比所處區(qū)間具有較大差異。為避免采用單一閾值可能造成的衛(wèi)星可見性誤判，引入動(dòng)態(tài)閾值來減弱該影響(圖4)，對(duì)每顆衛(wèi)星在不同歷元均采用動(dòng)態(tài)閾值后，閾值與相應(yīng)衛(wèi)星的信噪比區(qū)間能較好對(duì)應(yīng)。對(duì)于單顆衛(wèi)星，在SNR值波動(dòng)較小的情況下閾值變化也非常穩(wěn)定。

圖3 衛(wèi)星信噪比值變化

圖4 衛(wèi)星動(dòng)態(tài)閾值

根據(jù)文獻(xiàn)[7]，基于高分值加權(quán)的改進(jìn)陰影匹配定位算法研究中推薦使用38作為單閾值，本文所采集數(shù)據(jù)的信噪比值偏低，經(jīng)驗(yàn)證后采用32作為單閾值。以3 s時(shí)刻為例，手機(jī)原始觀測(cè)接收到16顆衛(wèi)星的信號(hào)，每顆衛(wèi)星的信噪比值如圖5(a)所示，灰色標(biāo)記的衛(wèi)星為5°截止高度角過濾的衛(wèi)星，不參與匹配過程，紅色橫線為單閾值(32閾值)，紅色部分為該時(shí)刻從上一時(shí)刻繼承及本時(shí)刻改進(jìn)后的動(dòng)態(tài)閾值。從圖中可以看出，采用動(dòng)態(tài)閾值后在衛(wèi)星可見性判斷方面與單一閾值相比具有明顯差異，紅色字體衛(wèi)星為采用動(dòng)態(tài)閾值后可見性發(fā)生變化的衛(wèi)星。圖5(b)為每個(gè)歷元采用單閾值和動(dòng)態(tài)閾值判斷的衛(wèi)星可見性與SNR判斷一致的衛(wèi)星個(gè)數(shù)(得分)。

圖5 依據(jù)32閾值、動(dòng)態(tài)閾值判斷第3 s衛(wèi)星可見性及與SNR判斷一致衛(wèi)星個(gè)數(shù)隨時(shí)間的變化

圖6為分別依據(jù)單一閾值法和動(dòng)態(tài)閾值法進(jìn)行陰影匹配計(jì)算后的得分，單一閾值法搜索區(qū)域中最高得分的位置較多，采用鄰近點(diǎn)KNN-SM算法確定的位置精度存在偏差。通過動(dòng)態(tài)閾值處理后原高分區(qū)域內(nèi)出現(xiàn)分化，存在相對(duì)單一閾值的更高分區(qū)域，面積相對(duì)變小，計(jì)算的點(diǎn)位位置相對(duì)更加精確。圖中紅色“*”代表真實(shí)位置，紫色“○”代表動(dòng)態(tài)閾值陰影匹配計(jì)算的位置，黑色“+”代表32閾值陰影匹配計(jì)算的位置。從圖中可以看出，動(dòng)態(tài)閾值較32單閾值在匹配結(jié)果精度上具有明顯提升。通過對(duì)比可知，第3 s時(shí)刻動(dòng)態(tài)閾值算法的改正值相較于32閾值在點(diǎn)位誤差上降低1.08 m，相比于原始觀測(cè)值降低8.37 m。如圖7所示，在第96 s時(shí)刻動(dòng)態(tài)閾值改進(jìn)后的點(diǎn)位誤差比32閾值SM算法降低3.46 m，定位精度具有較大提升。

圖6 32閾值和動(dòng)態(tài)閾值下第3 s打分圖

圖7 第96 s動(dòng)態(tài)閾值打分圖

對(duì)實(shí)驗(yàn)中采用不同算法獲得的點(diǎn)位誤差進(jìn)行分析發(fā)現(xiàn)，動(dòng)態(tài)閾值法相對(duì)于單閾值法(32閾值)和原始觀測(cè)值在過街方向上均有一定程度的改進(jìn)(圖8)，點(diǎn)位誤差也有一定程度改善(圖9)。然而，也存在部分位置的改進(jìn)效果不理想甚至變差的情況，為該時(shí)段衛(wèi)星數(shù)較少所致。衛(wèi)星數(shù)越少，動(dòng)態(tài)閾值所能改正的衛(wèi)星閾值越少，同樣單閾值時(shí)能使用的衛(wèi)星也就越少，導(dǎo)致誤差增大，改正效果不理想。

圖8 過街方向平均絕對(duì)誤差

圖9 點(diǎn)位誤差

從整體上看，動(dòng)態(tài)閾值法(DTSM)相對(duì)于32閾值陰影匹配(SM)算法在計(jì)算精度上有所改進(jìn)，動(dòng)態(tài)閾值相比于固定閾值，點(diǎn)位平均絕對(duì)誤差(mean absolute error，MAE)降低約15.1%，均方根誤差(root mean square error，RMSE)降低約14.9%。在過街方向上，相對(duì)于固定閾值，動(dòng)態(tài)閾值MAE降低約26.7%，RMSE降低約19.7%；相對(duì)于原始值，MAE降低約28.9%，RMSE降低約16.4%。

對(duì)原始觀測(cè)值進(jìn)行卡爾曼濾波(KF)后再進(jìn)行陰影匹配，采用動(dòng)態(tài)閾值和32閾值計(jì)算的坐標(biāo)點(diǎn)如圖10所示，卡爾曼濾波后采用動(dòng)態(tài)閾值陰影匹配(KF+DTSM)的RMSE為3.09 m，相比于卡爾曼濾波的4.09 m降低24.4%。圖11為不同方法的點(diǎn)位誤差，動(dòng)態(tài)閾值相比于其他兩種方法更加穩(wěn)定，卡爾曼濾波的MAE為4.77 m，KF+32閾值SM的MAE為4.58 m，KF+DTSM的MAE為3.55 m。KF+DTSM相比于KF在誤差上降低1.22 m，精度提升25.6%；相比于KF+32閾值SM降低1.03 m，精度提升22.4%。實(shí)驗(yàn)中不同算法的點(diǎn)位誤差統(tǒng)計(jì)見表1。

圖10 點(diǎn)位分布

圖11 加入卡爾曼濾波后的點(diǎn)位誤差對(duì)比

表1 不同算法點(diǎn)位誤差統(tǒng)計(jì)

3 結(jié) 語

陰影匹配算法對(duì)于城市峽谷中智能移動(dòng)設(shè)備GNSS導(dǎo)航定位的精度具有一定提升，但采用單一閾值進(jìn)行衛(wèi)星可見性判斷會(huì)增加部分衛(wèi)星誤判的可能。針對(duì)該問題，本文提出采用動(dòng)態(tài)閾值法進(jìn)行衛(wèi)星可見性判斷，在每個(gè)歷元中對(duì)基于衛(wèi)星幾何判斷可見性與信噪比判斷可見性進(jìn)行一致性分析，采用動(dòng)態(tài)閾值使一致性最大，從而改進(jìn)陰影匹配算法的定位精度。實(shí)驗(yàn)分析表明，相對(duì)于單一閾值法，動(dòng)態(tài)閾值陰影匹配算法對(duì)于手機(jī)GNSS定位結(jié)果的改正效果更加明顯，點(diǎn)位精度更加穩(wěn)定。在實(shí)驗(yàn)中也發(fā)現(xiàn)，初始定位精度對(duì)算法的有效性具有較大影響，利用卡爾曼濾波處理后再進(jìn)行動(dòng)態(tài)閾值陰影匹配，初始定位精度會(huì)有明顯改進(jìn)，原因可能為初始定位偏離真實(shí)值較大，搜索區(qū)域無法覆蓋到真實(shí)值，進(jìn)而導(dǎo)致算法改進(jìn)不理想。沿街方向由于直視信號(hào)較好，其定位更準(zhǔn)確，陰影匹配的改進(jìn)效果不夠理想。