一種對合成孔徑雷達(dá)影像中的道路進(jìn)行半自動(dòng)提取的方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于遙感圖像處理技術(shù)領(lǐng)域，設(shè)及一種對圖像中道路進(jìn)行半自動(dòng)提取的方 法，具體設(shè)及一種改進(jìn)剖面匹配和擴(kuò)展卡爾曼濾波對合成孔徑雷達(dá)（SyntheticAperture Radar,SAR)影像進(jìn)行道路半自動(dòng)提取的方法。
【背景技術(shù)】
[0002] 合成孔徑雷達(dá)（SyntheticApertureRadar,SAR)技術(shù)是一種重要的對地觀測手 段，廣泛應(yīng)用于軍事、民用等眾多領(lǐng)域。利用SAR影像進(jìn)行線目標(biāo)提取是遙感影像目標(biāo)解譯 的研究熱點(diǎn)，在眾多線目標(biāo)（道路、橋梁、海岸線、河流等）提取中道路提取占有重要地位， 從SAR影像中提取道路是典型線狀目標(biāo)提取的研究熱點(diǎn)，廣泛應(yīng)用于道路制圖、城市規(guī)劃、 構(gòu)建城市道路網(wǎng)和GIS數(shù)據(jù)庫更新等諸多領(lǐng)域。SAR影像道路提取根據(jù)是否需要人機(jī)交互， 可W分為自動(dòng)和半自動(dòng)的提取方法，但是由于道路的復(fù)雜性和多樣性使得SAR影像道路在 進(jìn)行全自動(dòng)提取后依然會(huì)進(jìn)行大量的人工后處理，因此道路全自動(dòng)提取近期還難W實(shí)現(xiàn)。 而道路半自動(dòng)提取的方法能夠結(jié)合操作員可靠的識別能力和計(jì)算機(jī)高效的計(jì)算能力，是對 當(dāng)前SAR影像道路提取的折中處理，具有重要的應(yīng)用潛力。
[000引 目前，針對SAR影像的道路半自動(dòng)提取方法有很多，主要的是對半自動(dòng)道路提取 方法的初始化，預(yù)測，觀測，更新等步驟進(jìn)行改進(jìn)。但是復(fù)雜的初始化使得算法的復(fù)雜度增 加，預(yù)測階段獲取的預(yù)測值不準(zhǔn)確會(huì)影響后續(xù)觀測值和擴(kuò)展卡爾曼濾波最優(yōu)估計(jì)值獲取的 準(zhǔn)確度和效率，觀測階段獲取觀測值得不精確會(huì)使得道路提取偏離道路中屯、位置產(chǎn)生錯(cuò)誤 的道路點(diǎn)。最優(yōu)觀測值得更新則需要用到預(yù)測值和觀測值。預(yù)測值和觀測值的好壞直接影 響了最優(yōu)觀測值的精度和人工干預(yù)次數(shù)的多少，如何在人工干預(yù)盡可能少的同時(shí)能夠?qū)Φ?路進(jìn)行精準(zhǔn)提取是當(dāng)前需要解決的難點(diǎn)問題。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明的目的在于從根本上克服道路半自動(dòng)提取存在的兩種缺陷，在分析SAR影 像道路特性的基礎(chǔ)上，構(gòu)建了SAR影像道路提取模型，提出利用改進(jìn)剖面匹配和擴(kuò)展卡爾 曼濾波對SAR影像進(jìn)行道路半自動(dòng)提取方法：在獲取道路觀測值上，對在道路面上的像素 設(shè)置權(quán)重為2,而對非道路面上的像素設(shè)置權(quán)重為1 ;在獲取道路最優(yōu)估計(jì)值上，利用擴(kuò)展 的卡爾曼濾波的方法；從而實(shí)現(xiàn)在人工干預(yù)盡可能少的同時(shí)對道路進(jìn)行精準(zhǔn)提取。 陽0化]本發(fā)明所采用的技術(shù)方案是：一種對合成孔徑雷達(dá)影像中的道路進(jìn)行半自動(dòng)提取 的方法，其特征在于，包括W下步驟：
[0006] 步驟1 :打開一幅SAR影像，并預(yù)處理；
[0007] 步驟2 :人工初始化道路參數(shù)；通過在SAR影像道路兩側(cè)垂直于道路方向上人工隨 機(jī)選取兩個(gè)點(diǎn)，道路中屯、點(diǎn)坐標(biāo)（r。，C。)，道路方向4。，道路寬度W，參考剖面長度L;
[0008] 步驟3 :構(gòu)建SAR影像道路提取模型，獲取道路預(yù)測狀態(tài)并計(jì)算目標(biāo)剖面；通過道 路提取模型得到下一時(shí)刻道路的預(yù)測狀態(tài)；設(shè)定捜索夾角f和步長化，W系統(tǒng)狀態(tài)預(yù)測值 為中屯、利用參考剖面長度，獲得目標(biāo)剖面的長度LT，同時(shí)獲取擴(kuò)展卡爾曼濾波的先驗(yàn)概率 密度；
[0009] 步驟4 :對參考剖面和目標(biāo)剖面進(jìn)行優(yōu)化處理；每隔一個(gè)像素生成平行于參考剖 面且長度與參考剖面相同的若干剖面，對生成剖面中相同位置像素進(jìn)行多視處理，得到經(jīng) 過優(yōu)化處理后的參考剖面；對目標(biāo)剖面采用與參考剖面優(yōu)化處理相同的方法進(jìn)行優(yōu)化；
[0010] 步驟5 :獲取道路觀測值及判斷；將參考剖面中的像素與目標(biāo)剖面上對應(yīng)的像素 進(jìn)行灰度差平方處理，并利用改進(jìn)剖面匹配算法，將參考剖面中道路面上的像素灰度差平 方的權(quán)值設(shè)定為2,不在道路上的像素灰度差平方的權(quán)值設(shè)定為1，使得改進(jìn)剖面匹配算法 有效的抑制了道路周邊弱散射地物對道路提取的影響；判斷能否準(zhǔn)確的獲取道路觀測值， 若能夠準(zhǔn)確的獲取，則執(zhí)行所述的步驟6 ;若不能夠獲取，則需要加入人工干預(yù)，判斷是否 到達(dá)圖像邊界，若到達(dá)圖像邊界，則道路提取結(jié)束；若沒有到達(dá)圖像邊界，貝陽日入人工干預(yù)， 并回轉(zhuǎn)執(zhí)行所述的步驟2;
[0011] 步驟6 :道路最優(yōu)估計(jì)值的獲取及判斷；利用步驟3獲取當(dāng)前時(shí)刻道路先驗(yàn)概率密 度和所述步驟5獲取的道路觀測值進(jìn)行處理，計(jì)算當(dāng)前時(shí)刻的道路最優(yōu)估計(jì)值；判斷能否 獲得準(zhǔn)確的最優(yōu)估計(jì)值；若能夠獲得則將當(dāng)前時(shí)刻的最優(yōu)觀測值作為下一時(shí)刻道路的初始 值，并回轉(zhuǎn)執(zhí)行所述的步驟3,直至完成整條道路的提??；若不能夠準(zhǔn)確獲取最優(yōu)估計(jì)值， 則需要加入人工干預(yù)，判斷是否到達(dá)圖像邊界，若到達(dá)圖像邊界，則道路提取結(jié)束；若沒有 到達(dá)圖像邊界，貝陽日入人工干預(yù)，并回轉(zhuǎn)執(zhí)行所述的步驟2 ;
[0012] 作為優(yōu)選，步驟1中所述的對SAR影像進(jìn)行預(yù)處理，其具體實(shí)現(xiàn)過程為：對待提取 SAR影像進(jìn)行精致Lee濾波處理來濾除斑點(diǎn)噪聲對道路提取的影響。
[0013] 作為優(yōu)選，步驟3中所述的構(gòu)建SAR影像道路提取模型為：
[0014] Xk=f(Xki)+Wk
[0015] 其中，f(?)表示為狀態(tài)轉(zhuǎn)移策略的關(guān)系式，Xk表示道路提取的第k個(gè)時(shí)刻道路系 統(tǒng)狀態(tài)，Wk代表系統(tǒng)噪聲。
[0016] 作為優(yōu)選，步驟3中所述的目標(biāo)剖面的長度LT的計(jì)算公式為：
[0017]
[001引作為優(yōu)選，步驟3中所述的擴(kuò)展卡爾曼濾波的先驗(yàn)概率密度，即預(yù)測狀態(tài)，通過道 路的預(yù)測狀態(tài)，計(jì)算擴(kuò)展卡爾曼濾波的第k時(shí)刻的先驗(yàn)概率密度和先驗(yàn)概率的協(xié)方差 矩陣Pkiki:
[0022]
[0023] 其中表示第k-1時(shí)刻道路點(diǎn)的最優(yōu)估計(jì)值，Qk表示和道路寬度有關(guān)的系統(tǒng) 噪聲，第A表示的線性化的過程通過雅可比行列式獲得，PkIlk1表示后驗(yàn)概率的協(xié)方差矩陣。
[0024] 作為優(yōu)選，步驟5中所述的改進(jìn)剖面匹配算法，其具體實(shí)現(xiàn)包括W下子步驟：
[00巧]步驟5. 1 :構(gòu)建道路觀測模型： 陽026]Zk=HkXk|ki;
[0027] 其中Zk表示道路的觀測值，Hk表示道路的觀測矩陣；
[00測步驟5. 2 :剖面匹配的加權(quán)策略計(jì)算：
[0029]
[0030] 其中W表示的是道路的寬度，L表示的是參考剖面的長度。
[0031] 步驟5. 3 :改進(jìn)剖面匹配算法的計(jì)算公式為：
[0032]
[003引其中g(shù)i(Xr，yr)表示參考剖面的像素灰度值，g'w(XT，yT)表示目標(biāo)剖面的像素灰 度值，j代表捜索夾角和步長所確定的長度，j=LT-L+1 ;
[0034] 作為優(yōu)選，步驟6中所述的獲取道路最優(yōu)估計(jì)值，其具體實(shí)現(xiàn)包括W下子步驟：
[0035] 步驟6. 1 :計(jì)算擴(kuò)展卡爾曼濾波增益：
[0036]
[0037] 其中Kk表示擴(kuò)展卡爾曼濾波增益，Hk表示觀測矩陣，Rk表示系統(tǒng)觀測噪聲；
[0038] 步驟6. 2 :計(jì)算擴(kuò)展卡爾曼濾波的后驗(yàn)概率密度，即道路的最優(yōu)估計(jì)值，獲取道路 的后驗(yàn)概率密度為巧后驗(yàn)概率協(xié)方差矩陣Pkik;
[0039]
[0040] Pk|k= (I-KkHk)Pkiki;
[0041] 其中Zk表示改進(jìn)剖面匹配獲取的觀測值，I表示單位矩陣；通過更新階段的計(jì)算 公式得出通過擴(kuò)展卡爾曼濾波提取的道路點(diǎn)即是一個(gè)最優(yōu)估計(jì)值。
[0042] 作為優(yōu)選，步驟6中所述的人工干預(yù)，即人工的在道路中屯、選擇一個(gè)道路點(diǎn)，來初 始化道路參數(shù)。
[0043] 本發(fā)明的創(chuàng)新之處在于：
[0044] (1)提出SAR影像道路提取模型，能夠準(zhǔn)確的獲取道路的預(yù)測值，并為利用擴(kuò)展卡 爾曼濾波算法提供理論條件； W45] 似提出了改進(jìn)剖面匹配的算法獲取觀測值，利用對SAR影像中道路部分像素設(shè) 置權(quán)值為2,非道路部分像素設(shè)置權(quán)值為1的方法，有效的抑制了道路周邊地物對道路提取 的影響，提高算法的魯棒性，能夠獲得較為準(zhǔn)確的觀測值；
[0046] (3)利用擴(kuò)展卡爾曼濾波的方法獲取道路最優(yōu)估計(jì)值，擴(kuò)展的卡爾曼濾波方法能 夠在對非線性系統(tǒng)進(jìn)行道路提取，并能夠在道路面被障礙物（周圍高大地物遮擋或出現(xiàn)隔 斷）干擾時(shí)，能夠通過已觀測獲得的數(shù)據(jù)對當(dāng)前時(shí)刻進(jìn)行預(yù)測，從而能夠跨過障礙物，繼續(xù) 道路提取，從而很大程度上提高了道路提取的效率，同時(shí)擴(kuò)展卡爾曼濾波有較為嚴(yán)格的檢 查機(jī)制，能夠在及時(shí)發(fā)現(xiàn)錯(cuò)誤的觀測值并停止道路提取等待操作員進(jìn)行人機(jī)交互，從而大 大的減少道路提取出錯(cuò)的可能，保證了道路提取算法的準(zhǔn)確性。
【附圖說明】
[0047] 圖1 :是本發(fā)明實(shí)施例的流程圖； W48] 圖2 :是本發(fā)明實(shí)施例中用到的道路模型圖； W例圖3 :是本發(fā)明實(shí)施例中用到的改進(jìn)剖面匹配算法。
【具體實(shí)施方式】
[0050] 為了便于本領(lǐng)域普通技術(shù)人員理解和實(shí)施本發(fā)明，下面結(jié)合附圖及實(shí)施例對本發(fā) 明作進(jìn)一步的詳細(xì)描述，應(yīng)當(dāng)理解，此處所描述的實(shí)施示例僅用于說明和解釋本發(fā)明，并不 用于限定本發(fā)明。
[0051] 請見圖1，本發(fā)明提供的一種對合成孔徑雷達(dá)影像中的道路進(jìn)行半自動(dòng)提取的方 法，包括W下步驟：
[0052] 步驟1 :打開一幅SAR影像，并對SAR影像進(jìn)行濾波，如：精致Lee濾波，來濾除斑 點(diǎn)噪聲對道路提取的影響；
[0053] 步驟2 :人工初始化道路參數(shù)；通過在SAR影像道路兩側(cè)垂直于道路方向上人工隨 機(jī)選取兩個(gè)點(diǎn)，后續(xù)獲取其它參數(shù)均由計(jì)算機(jī)自動(dòng)完成；道路中屯、點(diǎn)坐標(biāo)（r。，C。），道路方 向4。，道路寬度W，參考剖面長度L;
[0054] 步驟3 :構(gòu)建SAR影像道路提取模型，獲取道路預(yù)測狀態(tài)并計(jì)算目標(biāo)剖面；通過道 路提取模型得到下一時(shí)刻道路的預(yù)測狀態(tài)；設(shè)定捜索夾角口;和步長化，W系統(tǒng)狀態(tài)預(yù)測值 為中屯、利用參考剖面長度，獲得目標(biāo)剖面的長度LT，同時(shí)獲取擴(kuò)展卡爾曼濾波的先驗(yàn)概率 密度； 陽化5] 請見圖2,是本發(fā)明實(shí)施例中用到的道路模型圖；其中實(shí)屯、圓表示道路最優(yōu)估計(jì) 值，空屯、圓表示道路獲取的觀測值，X型表示預(yù)測值。（Cki，rki)，（Ck，rk)分別為第k-1時(shí) 刻和第k時(shí)刻道路最優(yōu)估計(jì)值XkIlk1和Xkik的坐標(biāo)，4k1和4k分別為第k-1時(shí)刻和第k 時(shí)刻道路的方向。
[0056] 本實(shí)施例的構(gòu)建SAR影像道路提取模型為：
[0057]
[0058] 其中，f( ?)表示為狀態(tài)轉(zhuǎn)移策略的關(guān)系式，Xk表示道路提取的第k個(gè)時(shí)刻道路系 統(tǒng)狀態(tài)，Wk代表系統(tǒng)噪聲。
[0059] 道路當(dāng)前時(shí)刻狀態(tài)可W通過一個(gè)四維向量表示，
[0060]
[OOW] 其中，（rk，Ck)表示第k個(gè)時(shí)刻道路中屯、點(diǎn)坐標(biāo)、(K表示第k個(gè)時(shí)刻道路方向、^ 表示第k個(gè)時(shí)刻道路方向變化角。
[0062] 目標(biāo)剖面的長度LT的計(jì)算公式為：