中巴地球資源一號衛(wèi)星紅外多光譜掃描irmss零級圖像質(zhì)量評價(jià)

中巴地球資源一號衛(wèi)星紅外多光譜掃描irmss零級圖像質(zhì)量評價(jià)

0衛(wèi)星屬性定量評價(jià)的必要性1999年10月14日，中國和巴基斯坦聯(lián)合開發(fā)出了中國的地球資源衛(wèi)星。中巴地球資源1號衛(wèi)星（cns-1）發(fā)射成功，并進(jìn)入了預(yù)定義的路線。從那時(shí)起，大量衛(wèi)星數(shù)據(jù)不斷地傳輸?shù)降孛?，這表明中國的空間遙感應(yīng)用進(jìn)入了一個(gè)新階段。此后，中國結(jié)束了資源衛(wèi)星的歷史，但迫切需要客觀地評估衛(wèi)星數(shù)據(jù)的質(zhì)量和應(yīng)用潛力，并為應(yīng)用提供指導(dǎo)。該操作僅接收到cbs-1中紅外多光譜采樣（iriss）中的四種破碎波圖像數(shù)據(jù)，并對應(yīng)用潛力進(jìn)行了全面的定量評估。1紅外信號特性研究使用的數(shù)據(jù)是由中巴地球資源一號衛(wèi)星(CBERS-1)IRMSS傳感器于2000年3月23日獲取的一景零級圖象,軌道編號為1,像幅編號為56.圖象共有4個(gè)波段,即1個(gè)全色波段(B6:0.50～0.90μm),2個(gè)短波紅外波段(B7:1.55～1.75μm,B8:2.08～2.35μm),1個(gè)熱紅外波段(B9:10.40～12.50μm);圖象覆蓋區(qū)地處冀中保定—廊坊地區(qū),經(jīng)緯度范圍為115°13′44″E～116°51′32″E,38°32′25″N～39°48′19″N;星下點(diǎn)分辨率:6、7、8波段為77.8m,9波段為156m;圖象大小:6、7、8波段為1712×1536pixels,9波段為864×768pixels,地面幅寬為119km.2正、噪聲濾除數(shù)據(jù)預(yù)處理工作分錯(cuò)位校正、輻射校正、噪聲濾除(包括去條紋、條帶噪聲,隨機(jī)噪聲等),另外,對圖象還進(jìn)行了直方圖修正、數(shù)據(jù)融合等增強(qiáng)處理.2.1奇偶掃描帶間的位置錯(cuò)動IRMSS系統(tǒng)在通過掃描鏡前向和后向振蕩掃描來獲取地物輻射能量實(shí)現(xiàn)成像的過程中,由于前、后掃描速率不同及衛(wèi)星姿態(tài)的不一樣等原因,可能導(dǎo)致圖象掃描帶間或掃描線間的錯(cuò)位等畸變現(xiàn)象,必須予以校正.(1)掃描線配準(zhǔn)在6、7、8波段原始圖象上奇偶掃描帶間(8行為一掃描帶)有約172個(gè)像元的位置錯(cuò)動,9波段原始圖象的奇偶掃描帶間(4行為一掃描帶)有84個(gè)像元的位置錯(cuò)動,因此處理時(shí),以奇數(shù)掃描帶為基準(zhǔn),將偶數(shù)掃描帶向左平移相應(yīng)的錯(cuò)動像元數(shù),即可完成掃描帶的配準(zhǔn).在經(jīng)掃描帶配準(zhǔn)后的圖象上,奇偶掃描線間還存在2個(gè)像元的位置錯(cuò)動,處理方法與掃描帶配準(zhǔn)類似,即以奇數(shù)掃描線為基準(zhǔn),對偶數(shù)掃描線向左平移2個(gè)像元.(2)掃描非線性配準(zhǔn)經(jīng)過掃描線配準(zhǔn)后的圖象上還存在奇偶掃描帶間的細(xì)微錯(cuò)動,處理方法是以奇數(shù)掃描帶為基準(zhǔn),在相鄰的偶數(shù)掃描帶上選取多個(gè)控制點(diǎn),建立奇偶掃描帶間的位置關(guān)系,把偶數(shù)掃描帶校正到與奇數(shù)掃描帶配準(zhǔn)的位置上.配準(zhǔn)公式為:Y=1E-05X2+0.9924X-0.2818(X,Y分別為奇、偶掃描帶的列坐標(biāo))2.2波前定標(biāo)值l由于沒有設(shè)定合適的輻射定標(biāo)場,這里僅對圖象進(jìn)行了相對定標(biāo).通過相對定標(biāo)可以消除圖象的邊緣輻射現(xiàn)象,并且也有利于波段間的輻射匹配.IRMSS6、7、8波段以定標(biāo)燈為內(nèi)定標(biāo)源,而9波段則以高低溫黑體為內(nèi)定標(biāo)源,輻射定標(biāo)過程須區(qū)別對待,具體作法如下:(1)前3個(gè)波段:定標(biāo)燈的輸出值分為4個(gè)級別,即最大可能輸出(這里設(shè)定為255)的100%,66%,33%和0%.在定標(biāo)圖象上,對每個(gè)級別選取比較均勻穩(wěn)定的區(qū)域,提取出10～20個(gè)連續(xù)亮條,每個(gè)亮條為8行,其分別對應(yīng)8個(gè)探測器,對每個(gè)行取平均得到8個(gè)值,4個(gè)級別共32個(gè)值,由每個(gè)探測器的4個(gè)定標(biāo)值可擬合一條定標(biāo)曲線,共獲得8條曲線,來對相應(yīng)圖象逐探測器行進(jìn)行定標(biāo).表1給出了7波段的定標(biāo)公式,可見定標(biāo)曲線均為顯著線性關(guān)系,相關(guān)系數(shù)R2均大于0.99.(2)第9波段:基于普朗克函數(shù)在一定溫度范圍內(nèi)近似為線性關(guān)系的原理,根據(jù)星上配備的高低溫黑體的兩組溫度和DN(Digitalnumber)值,可以逐波段、逐像元地將DN值轉(zhuǎn)換為像元亮度溫度值.轉(zhuǎn)換公式為Tij=1DNBB2(i)?DNBB1(i)[(TBB2?TBB1)×DNij+DNBB2(i)×TBB1?DNBB1(i)×TBB2]Τij=1DΝBB2(i)-DΝBB1(i)[(ΤBB2-ΤBB1)×DΝij+DΝBB2(i)×ΤBB1-DΝBB1(i)×ΤBB2]其中,Tij為某波段圖象上i行、j列像元亮度溫度值(K);DNBB2(i),DNBB1(i)為圖象上第i行所對應(yīng)的高溫黑體溫度TBB2和低溫黑體溫度TBB1的平均灰度值.2.3壞行間斷性噪聲的校正在IRMSS圖象上,7波段間斷性噪聲最明顯,而條帶現(xiàn)象在各個(gè)波段圖象上表現(xiàn)得都很明顯,7波段和9波段的隨機(jī)噪聲較嚴(yán)重,7、8、9波段圖象上存在壞行現(xiàn)象,為此采用不同的方法濾除噪聲:周期性噪聲用二維傅立葉變換法提取噪聲分量;壞行用相鄰兩行或幾行的平均值代替來修補(bǔ);部分壞行(間斷性噪聲)模板卷積和域值法;去條帶單個(gè)探測器分圖象的累積直方圖與總的累積直方圖相匹配的非線性校正方法;隨機(jī)噪聲(斑點(diǎn))中值濾波法.2.4圖像幾何的校正和配置(1)基于envi的幾何校正由于IRMSS是掃描成像,在水平方向上會產(chǎn)生橫向拉伸,需要進(jìn)行正切校正.但是,橫向的拉伸不足2個(gè)像元,所以只用控制點(diǎn)進(jìn)行幾何校正.根據(jù)圖象區(qū)域范圍及特性,選取1∶100000地形圖9張,對照地形圖,選取地面控制點(diǎn)共18個(gè),以高斯克呂格為投影標(biāo)準(zhǔn),在標(biāo)準(zhǔn)遙感圖象處理軟件ENVI上對6波段圖象進(jìn)行幾何精糾正,均方根誤差(RMS)為0.3851.(2)橫向錯(cuò)動的相對穩(wěn)定性以圖象質(zhì)量較好的第6波段圖象為基準(zhǔn),選取控制點(diǎn)對7、8、9波段圖象進(jìn)行波段間的配準(zhǔn).結(jié)果發(fā)現(xiàn),波段間的橫向錯(cuò)動比較大,縱向錯(cuò)動相對較小.7波段相對6波段向右偏移了18個(gè)像元,向下偏移了8個(gè)像元;8波段相對于6波段向右偏移了43個(gè)像元,向下偏移了11個(gè)像元.經(jīng)過以上預(yù)處理后的各波段圖象基本上消除了條紋、條帶、壞行以及隨機(jī)噪聲等的影響(見圖版Ⅰ圖1和圖2).2.5圖像信息的擴(kuò)展和澄清(1)灰度變換為了增強(qiáng)圖象的顯示效果以便于識別和分析,需對圖象的灰度信息進(jìn)行變換處理,這里采用的是分段線性拉伸或直方圖均衡化的方法.(2)處理后的構(gòu)圖采用直方圖匹配的方法來提高彩色合成圖象的分辨力.處理后的圖象見圖版Ⅰ圖1和圖2,從處理后的圖象上,可以清晰地識別出村莊、水體、道路、林地等多種地物,可識別能力大大提高.2.6表象融合的方法利用高空間分辨率第6波段圖象的細(xì)節(jié)信息和低空間分辨率第9波段圖象的光譜信息進(jìn)行融合處理,使第9波段圖象的質(zhì)量得到明顯改善,其具體處理是基于小波變換進(jìn)行的,處理流程如下:(1)以6波段圖象為基準(zhǔn),對9波段圖象進(jìn)行配準(zhǔn),并且重采樣至相同大小.(2)對兩圖象進(jìn)行直方圖匹配.(3)對兩圖象進(jìn)行n次小波分解,通常n取2或3.結(jié)果為I=∑i=1nIpi+IrΙ=∑i=1nΙpi+Ιr其中,I代表原圖象,Ipi表示圖象經(jīng)小波變換分解得到的高頻細(xì)節(jié),Ir表示圖象經(jīng)n次小波變換分解得到的低頻輪廓部分.(4)以6波段圖象中的細(xì)節(jié)部分∑i=1nIpi∑i=1nΙpi代替9波段圖象的細(xì)節(jié)部分.(5)對9波段圖象進(jìn)行小波逆變換得到融合后的結(jié)果圖象.融合前后的結(jié)果見圖版Ⅰ圖3,通過對比不難發(fā)現(xiàn),融合后圖象清晰度明顯增強(qiáng),可分辨力提高,這表明小波變換是有效的,能夠較好地提高圖象的質(zhì)量.3圖像質(zhì)量分析3.1數(shù)字遙感影像綜合評價(jià)指標(biāo)圖象質(zhì)量評價(jià)最常用的一個(gè)指標(biāo)是圖象的分解力,它表示某一遙感成像系統(tǒng)對微小細(xì)節(jié)反差表達(dá)的能力.但它只是表示對微小地物的極限分辨能力,不能從分解力的值去了解遙感器對較大地物的表達(dá)質(zhì)量.清晰度表示影像邊界的清晰程度,數(shù)字遙感影像的清晰程度,可用相鄰像元之間在x,y兩個(gè)方向的灰度變化速率的加權(quán)平均值來表達(dá),清晰度的數(shù)值本身并無絕對的意義,只是用于相互比較的一種指標(biāo),并且由于地物種類的不同,清晰度有很大的差異.圖象的反差能在一定程度上反映遙感影像的目視效果,通常都要求影像的反差適中.以上這些評價(jià)標(biāo)準(zhǔn)都局限于某一方面,因此要全面反映圖象的質(zhì)量,必須將這些指標(biāo)統(tǒng)一起來,進(jìn)行綜合評價(jià).具體指標(biāo)為:地面分解力、清晰度、信噪比、反差、輻射精度.(1)實(shí)際地面調(diào)制度采用調(diào)制傳遞函數(shù)(MTF)計(jì)算地面分解力.首先定義以下3個(gè)變量:反差:u=(ro?rb)rb(ro?rbu=(ro-rb)rb(ro?rb:目標(biāo)與背景的輻射度)調(diào)制度:M=(u-1)/(u+1)調(diào)制函數(shù):T(f)=Mi(f)/Mb(f),其中Mb(f)和Mi(f)分別為在頻率f的條件下,目標(biāo)物的地面真實(shí)調(diào)制度和對應(yīng)圖象上像元的調(diào)制度.地面分解力計(jì)算的具體步驟如下:①找到圖象中有明顯變化的邊緣作樣區(qū),沿垂直該邊緣掃取N條數(shù)據(jù)平均后,獲得刀刃曲線G(k),其中,k表示曲線上離散點(diǎn)的序號;②采用“五點(diǎn)三次法”對G作平滑處理;③三次樣條內(nèi)插;④求導(dǎo),進(jìn)行平滑處理,得到曲線A(k),其中,k=1,…,n-1;⑤進(jìn)行傅立葉變換x(L)=∑k=0N?1A(k)e?i2πkLN(k=0,1,2,?,n?1,n為離散數(shù)據(jù)個(gè)數(shù))；x(L)=∑k=0Ν-1A(k)e-i2πkLΝ(k=0,1,2,?,n-1,n為離散數(shù)據(jù)個(gè)數(shù))；⑥求調(diào)制傳遞函數(shù):T(L)=x(L)max(|x(L)|)Τ(L)=x(L)max(|x(L)|);⑦當(dāng)T(L)=0.05時(shí),反求L;⑧地面分解力為Sp×W/L,其中,Sp為像元尺寸,Sp=30m,W為掃描數(shù)據(jù)寬度,單位為pixel.(2)方向總復(fù)雜度v采用邊緣銳度值的“后驗(yàn)?zāi)Ｐ突惴ā眮碛?jì)算清晰度(EAV)EAV=∑ab(df/dx)2|f(b)?f(a)|)EAV=∑ab(df/dx)2|f(b)-f(a)|)其中,df/dx為垂直于邊緣的灰度變化率,f(b)-f(a)為該方向總對比度.(3)求象方差極大值lmsdm利用局域平均值與方差比來計(jì)算SNR,具體步驟如下:①求圖象平均值M;②圖象分塊(4×4,8×8…),求局域均值與方差LM,LSD;③求出局域方差極大值LSDm;④SNR=M/LSDm(4)比較反差利用調(diào)制傳遞函數(shù)(MTF)曲線所圍總面積來衡量:面積越大,圖象越清晰.(5)輻射精度評價(jià)輻射精度是反映圖象信息量豐富程度的指標(biāo).同一地區(qū)的不同圖象,若灰度分布范圍越大,方差越大,則表明圖象信息越豐富.采用如下8個(gè)變量來綜合評價(jià)輻射精度:指標(biāo)均值m=∑(i×p(i))方差d=∑((i-m)2p(i))偏斜S=∑((i-m)3p(i))/d3陡度K=∑((i-m)4p(i))/d4其中,i為圖象像元可能取值,p(i)為圖象像元值為i的概率.邊緣輻射畸變R由行方向均值矢量的均值和方差來衡量.增益調(diào)整畸變P由列方向均值矢量的均值和方差來衡量.3.2合成構(gòu)圖的差異根據(jù)上述評價(jià)指標(biāo)和步驟,得到中巴地球資源一號衛(wèi)星IRMSS數(shù)據(jù)的評價(jià)結(jié)果見表2,為便于比較,表中同時(shí)給出TM數(shù)據(jù)的評價(jià)結(jié)果.需要說明的是,在所統(tǒng)計(jì)的指標(biāo)中,地面分解力、圖象清晰度、反差、信噪比具有一定的普遍性和可比性.而其他指標(biāo)隨不同地域和季節(jié)的不同,具有較大的可變性.由表2分析得出,在圖象的信噪比方面,TM具有絕對優(yōu)勢;在圖象偏斜度方面,TM好于IRMSS;雖然IRMSS7、8波段的地面分解力稍微高于TM,但是與此對應(yīng)的空間分辨率(地面分解力×像元空間分辨率)還是遜于TM;IRMSS8、9波段的圖象清晰度稍好于TM,但7波段比TM差;圖象反差方面,IRMSS6、7波段高于TM,但8波段低于TM;輻射精度方面,IRMSS8波段方差好于TM,但7、9波段明顯低于TM.另外,根據(jù)上述評價(jià)指標(biāo)和步驟,同時(shí)評價(jià)了處理前后圖象的質(zhì)量(表3).其中,處理前圖象指的是只做過錯(cuò)位糾正的圖象,而處理后圖象指的是經(jīng)綜合處理后的圖象.由表3可以看到,處理后圖象的地面分解力明顯提高;清晰度也有一定程度改善;反差增強(qiáng);信噪比大幅度提高;圖象方差增強(qiáng),圖象反映的信息更加豐富,從而圖象質(zhì)量大大改善.第6波段圖象相對于其他3個(gè)波段的信噪比為最高,反差最大,地面分解力最優(yōu),這也是它從直觀上看,圖象質(zhì)量遠(yuǎn)遠(yuǎn)好于其他波段的原因.4ss:6、7、8段IRMSS數(shù)據(jù)影像雖然在空間分辨率上低于TM影像,但其信息量大,不同波段合成的假彩色影像圖對地物、地貌景觀反映十分清楚.單波段影像清晰,色調(diào)一致,反差適中.IRMSS6、7、8波段合成圖象色彩鮮艷,層次分明,對河流、溝渠、水庫、坑塘反映較好,影像清晰;水體的深淺及濕地的賦水程度也有所反映;特別對干河谷反映更好,在視覺上比TM圖象清晰;對地貌、構(gòu)造,特別是對沖洪積扇反映效果好于TM;對礦產(chǎn)、地?zé)岬慕庾g尚有一定的難度,但由于IRMSS空間分辨率較低,譜帶寬,對土地資源解譯只能達(dá)到一級分類水平,對林地解譯也不能細(xì)分.總的來說,中巴地球資源一號衛(wèi)星IRMSS數(shù)據(jù),在多個(gè)領(lǐng)域的應(yīng)用前景還是很廣闊的,可為國土資源調(diào)查、能源、環(huán)保、城市規(guī)劃等領(lǐng)域提供服務(wù).5視角方面的改善通過對中巴地球資源一號衛(wèi)星IRMSS圖象零級數(shù)