杜偉娜， 徐愛(ài)功， 宋耀鑫， 孫華生

(1.遼寧工程技術(shù)大學(xué),阜新 123000； 2.公安部交通管理科學(xué)研究所,無(wú)錫 214151)

?

新型SAR傳感器一級(jí)地距產(chǎn)品絕對(duì)輻射定標(biāo)方法

杜偉娜1， 徐愛(ài)功1， 宋耀鑫2， 孫華生1

(1.遼寧工程技術(shù)大學(xué),阜新 123000； 2.公安部交通管理科學(xué)研究所,無(wú)錫 214151)

針對(duì)目前新型SAR傳感器數(shù)據(jù)產(chǎn)品預(yù)處理軟件缺乏的現(xiàn)狀，全面細(xì)致地論述了ENVISAT ASAR,Radarsat2,Cosmoskymed,TerraSAR-X和Sentinel1等5種新型SAR傳感器的一級(jí)地距產(chǎn)品(level-1 detected ground range product,L-1 DGRP)數(shù)據(jù)的絕對(duì)輻射定標(biāo)方法以及軟件化過(guò)程中所涉及的后向散射系數(shù)等參數(shù)獲取方法； 并以Sentinel1傳感器的L-1 DGRP數(shù)據(jù)為例，采用C++語(yǔ)言編程實(shí)現(xiàn)了絕對(duì)輻射定標(biāo)過(guò)程； 最后將本文方法軟件處理結(jié)果和歐空局S1 ToolBox軟件處理結(jié)果進(jìn)行對(duì)比，二者得到的后向散射系數(shù)值基本相同，證實(shí)了本文所介紹方法的正確性。

新型SAR； 絕對(duì)輻射定標(biāo)； 一級(jí)地距產(chǎn)品(L-1 DGRP)； 后向散射系數(shù)； Sentinel1

0 引言

隨著國(guó)外高分辨率合成孔徑雷達(dá)(synthetic aperture Radar,SAR)數(shù)據(jù)對(duì)我國(guó)各領(lǐng)域的開放以及“龍計(jì)劃”項(xiàng)目[1]的深入開展，ENVISAT ASAR,Radarsat1/2,Cosmoskymed,TerraSAR-X,ALOS PALSAR和Sentinel1等高分辨率星載SAR影像數(shù)據(jù)受到廣大關(guān)注。然而，當(dāng)前能夠支持這些星載SAR數(shù)據(jù)預(yù)處理的商業(yè)軟件還比較少； 因此，研究這些SAR傳感器各級(jí)產(chǎn)品數(shù)據(jù)的輻射定標(biāo)算法并進(jìn)行軟件實(shí)現(xiàn)，對(duì)于促進(jìn)新型高分辨率合成孔徑成像雷達(dá)數(shù)據(jù)的定量應(yīng)用具有重要意義。與SAR數(shù)據(jù)的其他各級(jí)產(chǎn)品相比，一級(jí)地距產(chǎn)品(level-1 detected ground range product,L-1 DGRP)數(shù)據(jù)在陸地資源環(huán)境遙感中的應(yīng)用更為廣泛，也是衛(wèi)星地面接收站和數(shù)據(jù)生產(chǎn)商主要提供的數(shù)據(jù)產(chǎn)品。L-1 DGRP數(shù)據(jù)記錄的是雷達(dá)微波后向散射信號(hào)的振幅，而非地球生物物理參數(shù)定量化研究所需的地物后向散射系數(shù)[2]。在地表參數(shù)反演中，輻射定標(biāo)處理是不可或缺的，該過(guò)程是將傳感器接收的地物后向散射強(qiáng)度信息轉(zhuǎn)化為地物后向散射系數(shù)的唯一途徑。因此，本文根據(jù)SAR系統(tǒng)輻射定標(biāo)的基本原理，介紹了ENVISAT ASAR,Radarsat2,Cosmoskymed,TerraSAR-X和Sentinel1等5種新型星載SAR傳感器的L-1 DGRP數(shù)據(jù)的絕對(duì)輻射定標(biāo)方法； 并以Sentinel1數(shù)據(jù)產(chǎn)品的輻射定標(biāo)為例，采用C++語(yǔ)言編程實(shí)現(xiàn)了其絕對(duì)輻射定標(biāo)過(guò)程。

1 SAR輻射定標(biāo)基本原理

星載合成孔徑雷達(dá)在探測(cè)地物目標(biāo)時(shí)，雷達(dá)系統(tǒng)主動(dòng)向目標(biāo)發(fā)射無(wú)線電脈沖，其發(fā)射能量部分被目標(biāo)吸收，其余大部分能量經(jīng)目標(biāo)散射后被雷達(dá)天線接收，形成回波信號(hào)。這一能量傳輸過(guò)程可用雷達(dá)方程定量表示[3-4]為

(1)

式中：Pr為雷達(dá)接收的信號(hào)功率；Pt為發(fā)射功率；Gr為總接收增益；G(θ,φ)為距離向和方位向的增益；λ為波長(zhǎng)；S為地距；σ為目標(biāo)的雷達(dá)截面積。

雷達(dá)截面積σ是在假設(shè)雷達(dá)所接收的回波信號(hào)功率密度與實(shí)際目標(biāo)在同等條件下所接收的目標(biāo)回波信號(hào)功率密度相同的理想條件下，定義的實(shí)際目標(biāo)的等效雷達(dá)截面積，是表征地物目標(biāo)散射特性的常用參數(shù)。

對(duì)于分布目標(biāo)，雷達(dá)截面積σ可表示為后向散射系數(shù)的函數(shù)，即

σ=σ0Ac,

(2)

式中：σ0為單位面積的雷達(dá)截面積(即后向散射系數(shù))，可用來(lái)表征目標(biāo)對(duì)電磁波的散射能力；Ac為雷達(dá)地面分辨單元的面積。

星載SAR傳感器在發(fā)射之前，通常會(huì)利用角反射器在地面試驗(yàn)場(chǎng)進(jìn)行多次絕對(duì)輻射定標(biāo)實(shí)驗(yàn)； 并且大多數(shù)SAR傳感器會(huì)將輻射定標(biāo)模型所需的參數(shù)簡(jiǎn)化成角度、天線增益、定標(biāo)常量等[5]，并存儲(chǔ)在數(shù)據(jù)產(chǎn)品的頭文件中供用戶使用。

2 典型SAR傳感器絕對(duì)輻射定標(biāo)方法

與傳統(tǒng)SAR相比，新一代SAR傳感器普遍具有高空間分辨率、高靈敏輻射強(qiáng)度、多極化模式和多掃描帶寬等特點(diǎn)，表1列出5種典型SAR傳感器的參數(shù)[6-11]。

表1 5種新型SAR傳感器的參數(shù)

2.1 ENVISAT ASAR數(shù)據(jù)輻射定標(biāo)

ENVISAT ASAR[6]是迄今為止在環(huán)境監(jiān)測(cè)領(lǐng)域應(yīng)用最為廣泛的星載合成孔徑雷達(dá)，為遙感應(yīng)用領(lǐng)域的學(xué)者提供了大量高質(zhì)量的數(shù)據(jù)。ASAR的L-1 DGRP數(shù)據(jù)包括： ASA_IMP_1P,ASA_IMM_1P,ASA_APP_1P, ASA_APM_1P, ASA_WSM_1P, ASA_IMG_1P和 ASA_APG_1P共7種，通常以“.N1”文件格式分發(fā)。圖像中任一像元(i,j)的后向散射系數(shù)σ0為

(3)

對(duì)后向散射系數(shù)求平均，以db為單位表示，其計(jì)算公式為

(4)

式中:U和V分別為平均窗口的行和列；N=UV。

以上定標(biāo)模型涉及的參數(shù)均可從“.N1”格式的文件中直接或間接獲取。DN是數(shù)據(jù)產(chǎn)品的像元值，直接從SAR數(shù)據(jù)產(chǎn)品文件中的雷達(dá)數(shù)據(jù)記錄中讀??； 標(biāo)定常數(shù)K從ASAR數(shù)據(jù)文件的元數(shù)據(jù)中讀取，該參數(shù)與極化方式有關(guān)，具體讀取方式見表2。

表2 ASAR數(shù)據(jù)絕對(duì)定標(biāo)常數(shù)K

入射角αi,j需要通過(guò)計(jì)算得到，在注記數(shù)據(jù)集下的Geolocation Grid ADSRs中記錄了SAR的入射角信息，即影像方位向的行和距離向的列均存在對(duì)應(yīng)的11組入射角，各像元(i,j)的入射角通過(guò)插值得到。由于影像在方位向不超過(guò)100 km時(shí)可認(rèn)為入射角隨影像方位向行號(hào)變化的影響可以忽略不計(jì)，因此入射角只隨列號(hào)而變化。通常選擇一元二次方程作為插值函數(shù)，利用最小二乘法估計(jì)，即可計(jì)算得到列號(hào)j對(duì)應(yīng)的像元入射角αi,j。

2.2 Radarsat-2數(shù)據(jù)輻射定標(biāo)

2007年發(fā)射的Radarsat2[7]是Radarsat1的后續(xù)衛(wèi)星，它引領(lǐng)了星載SAR的多項(xiàng)技術(shù)革命，如第一次出現(xiàn)了全極化模式。Radarsat2的L-1 DGRP數(shù)據(jù)以Product.xml文件組織，包含TIFF格式的影像數(shù)據(jù)和XML格式的增益偏移參數(shù)文件，其絕對(duì)輻射定標(biāo)過(guò)程可以分為獲取雷達(dá)截面積影像和計(jì)算入射角。

1)獲取雷達(dá)截面積影像。Radarsat2影像雷達(dá)截面積σ為

σ=(DN2+A)/G ,

(5)

式中：DN為雷達(dá)影像像元值；A為增益偏移參數(shù)；G為第j像元(即第j列)的增益。其中偏移和增益參數(shù)可以在lutSigma.xml文件中獲取。

2)計(jì)算像元的入射角。雷達(dá)入射角與地球半徑和軌道高度有關(guān)，其幾何關(guān)系為

(6)

式中：h為軌道高度；R為地球半徑；s為斜距；Sj為圖像中第j列的地距。對(duì)于L-1DGRP，在升軌右視和降軌左視情況下(根據(jù)Product.xml文件中的passDirection和antennaPointing值判斷)，Sj的計(jì)算公式為

Sj=a+jLb+(jL)2c+(jL)3d+(jL)4e+(jL)5f ;

(7)

在降軌右視和升軌左視情況下，Sj的計(jì)算公式為

Sj=a+(J-j)Lb+[(J-j)L]2c+[(J-j)L]3d+[(J-j)L]4e+[(J-j)L]5f ,

(8)

式中：L為像元間隔；J為影像總列數(shù)；a,b,c,d,e,f分別為地距到斜距的6個(gè)轉(zhuǎn)換參數(shù)，可從Product.xml文件中的groundToSlantRangeCoefficients獲取。此外，Sj也可以通過(guò)最小入射角和最大入射角插值得到。

綜合以上過(guò)程，可得到雷達(dá)后向散射系數(shù)圖像σ0(以db為單位表示)，即

(9)

2.3 Cosmoskymed數(shù)據(jù)輻射定標(biāo)

Cosmoskymed[8]是意大利研發(fā)的COSMO-skymed高分辨率雷達(dá)衛(wèi)星星座的代稱，該星座共有4顆衛(wèi)星，重訪能力高，具備全球范圍觀測(cè)能力。Cosmoskymed的L-1 DGRP數(shù)據(jù)(有DGM和GEC兩種)通常以HDF5格式(“.H5”)分發(fā)。其絕對(duì)輻射定標(biāo)過(guò)程為

(10)

式中： |imginp(i,j)|2為計(jì)算強(qiáng)度影像表達(dá)式；αref為參考入射角；F為尺度因子；K為定標(biāo)常數(shù)。以db為單位表示，則

σ0(i,j)db=10lg[σ0(i,j)] ,

(11)

以上參數(shù)均可在HDF5文件的元數(shù)據(jù)中找到，其中αref取自Reference Slant Range；F取自Rescaling Factor；K取自Calibration Constant。

為了減少系統(tǒng)噪聲的影響，通常可以采用一個(gè)窗口平均值替代，即

(12)

式中：r,c分別為窗口的行和列；wr,wc分別為窗口的行寬和列寬。

2.4 TerraSAR-X數(shù)據(jù)輻射定標(biāo)

TerraSAR-X[9]是一顆成像分辨率和軌道精度都非常高的SAR衛(wèi)星，在大范圍對(duì)地觀測(cè)和干涉測(cè)量領(lǐng)域具有明顯優(yōu)勢(shì)。TerraSAR-X的L-1 DGRP數(shù)據(jù)也有2種： 多視地距探測(cè)產(chǎn)品(multi look ground range detected,MGD)和地理編碼橢球校正產(chǎn)品(geocoding of terrain correction,GEC)，通常以“.XML”格式組織，包含TIFF影像數(shù)據(jù)、快視數(shù)據(jù)和元數(shù)據(jù)等。TerraSAR-X的L-1 DGRP數(shù)據(jù)為絕對(duì)輻射定標(biāo)處理提供了定標(biāo)常量、等效噪聲分布和本地入射角參數(shù)，具體過(guò)程如式(13)―(16)所示，即

σ0=(ks|DN|2-B)sinθloc,

(13)

(14)

(15)

(16)

式(13)―(16)中：ks為輻射定標(biāo)常量，對(duì)應(yīng)XML格式元數(shù)據(jù)文件中calibrationConstant節(jié)點(diǎn)下calFactor的參數(shù)值；DN為像元的后向散射強(qiáng)度值，影像數(shù)據(jù)位于IMAGEDATA文件夾；B為強(qiáng)度等效噪聲，反映不同噪聲分布模型對(duì)信號(hào)的影響，所涉及的參數(shù)均可在元數(shù)據(jù)文件中獲取(由于B的影響很小，在有些情況下該參數(shù)可以忽略)；θloc為本地入射角，其計(jì)算所需的參數(shù)GGIM可以在georef.xml文件中獲?。籨eg和coffi分別為多項(xiàng)式擬合的階數(shù)和系數(shù)；τ為快時(shí)間；τref為參考快時(shí)間；τmin和τmax分別為快時(shí)間的最小值和最大值。

2.5 Sentinel1數(shù)據(jù)輻射定標(biāo)

歐空局于2014年發(fā)射的新一代SAR衛(wèi)星Sentinel1[11](亦稱“哨兵1號(hào)”)是ENVISAT衛(wèi)星的后繼衛(wèi)星，繼續(xù)擔(dān)負(fù)著對(duì)地環(huán)境觀測(cè)的任務(wù)，其數(shù)據(jù)可以在歐空局官方網(wǎng)站免費(fèi)申請(qǐng)獲得。Sentinel1的L-1 DGRP數(shù)據(jù)代碼為GRD，以“manifest.safe”索引文件組織，包含影像數(shù)據(jù)(measurement文件夾)、快視數(shù)據(jù)(preview文件夾)、軌道及標(biāo)定參數(shù)數(shù)據(jù)(annotation文件夾)以及說(shuō)明數(shù)據(jù)(support文件夾)。對(duì)于Sentinel1的L-1 DGRP數(shù)據(jù)的絕對(duì)輻射定標(biāo)，數(shù)據(jù)產(chǎn)品本身提供了一個(gè)定標(biāo)矢量，將影像的強(qiáng)度值轉(zhuǎn)化為后向散射系數(shù)，具體轉(zhuǎn)換過(guò)程為

(17)

式中：DN為TIFF格式地距影像中對(duì)應(yīng)像元的像素值；Ssigma為定標(biāo)參數(shù)，該參數(shù)可通過(guò)“annotation”文件夾下的“Calibration”文件夾中的XML元數(shù)據(jù)文件sigmaNought域得到。元數(shù)據(jù)中Ssigma是一個(gè)查找表(LUT)，為行方向和列方向具有一定間隔的矢量數(shù)組； 對(duì)于任意像元的Ssigma，需要通過(guò)插值得到，插值方法往往選擇雙線性插值方法。

3 實(shí)驗(yàn)與分析

本文以Sentinel1的L-1 DGRP數(shù)據(jù)為例，采用VC++編程語(yǔ)言和GDAL柵格數(shù)據(jù)讀寫庫(kù)，驗(yàn)證了絕對(duì)輻射定標(biāo)方法。實(shí)驗(yàn)數(shù)據(jù)為2015年3月15日獲取的覆蓋環(huán)渤海地區(qū)Sentinel1的L-1 DGRP，影像大小為19 342行×25 267列，空間分辨率為5 m×20 m，干涉寬幅掃描(IW)成像模式，掃描帶寬為250 km。根據(jù)2.5節(jié)介紹的方法對(duì)Sentinel1的L-1 DGRP數(shù)據(jù)進(jìn)行絕對(duì)輻射定標(biāo)。首先讀取影像的像元值，Sentinel1影像采用16 bit記錄強(qiáng)度信息，因此像元值在0～65 535之間； 然后獲取每個(gè)像元定標(biāo)參數(shù)值S，由于數(shù)據(jù)產(chǎn)品提供了30行×633列的定標(biāo)參數(shù)矩陣以及對(duì)應(yīng)的行列坐標(biāo)矩陣，因此采用雙線性內(nèi)插方法計(jì)算其他行列坐標(biāo)像元的定標(biāo)參數(shù)值； 最后根據(jù)式(17)計(jì)算得到每個(gè)像元的后向散射系數(shù)值，即得到輻射定標(biāo)處理結(jié)果。

圖1(a)為用本文方法的軟件實(shí)現(xiàn)的Sentinel1的L-1 DGRP數(shù)據(jù)輻射定標(biāo)結(jié)果，圖1(b)為用歐空局的S1-ToolBox軟件對(duì)相同數(shù)據(jù)進(jìn)行輻射定標(biāo)得到的結(jié)果，兩者在視覺(jué)上差別甚小。

(a) 本文方法 (b) 歐空局S1-ToolBox軟件

圖1 2種不同方法軟件輻射標(biāo)定結(jié)果

Fig.1 Calibration results by using two kinds of software from different methods

為了進(jìn)一步分析本文方法絕對(duì)輻射定標(biāo)處理結(jié)果的精度，從以上2個(gè)輻射定標(biāo)結(jié)果中隨機(jī)選取了12組樣本點(diǎn)的后向散射系數(shù)值進(jìn)行對(duì)比，誤差分析結(jié)果如表3所示。

表3 2種方法標(biāo)定精度對(duì)比

可以看出，2種方法處理結(jié)果相對(duì)誤差的最大值為0.003，最小值為0?？紤]到計(jì)算機(jī)插值處理精度的問(wèn)題，以及忽略有效數(shù)字的差別，本文方法的軟件處理得到的結(jié)果與歐空局S1-ToolBox軟件處理的結(jié)果影像是基本相同的，從而證實(shí)了本文所介紹方法的正確性。

4 結(jié)論

輻射定標(biāo)是SAR影像地表參數(shù)定量反演中的關(guān)鍵步驟，本文在明確幾種新型SAR一級(jí)地距產(chǎn)品(L-1 DGRP)數(shù)據(jù)組織結(jié)構(gòu)特點(diǎn)的基礎(chǔ)上，結(jié)合輻射定標(biāo)基本原理，闡述了幾種新型SAR的L-1 DGRP輻射定標(biāo)方法，結(jié)論如下：

1) 針對(duì)ENVISAT ASAR,Radarsat2,Cosmoskymed,TerraSAR-X和Sentinel1的L-1 DGRP數(shù)據(jù)，發(fā)展了利用產(chǎn)品的定標(biāo)參數(shù)和入射角插值方法、進(jìn)而利用輻射定標(biāo)公式進(jìn)行后向散射系數(shù)計(jì)算。

2) 以Sentinel1的L-1 DGRP數(shù)據(jù)為例，通過(guò)對(duì)2種方法的軟件處理結(jié)果對(duì)比研究，驗(yàn)證了本文所述算法的正確性。

3)下一步的研究將充分考慮地形因子對(duì)新型SAR的L-1 DGRP數(shù)據(jù)輻射定標(biāo)的影響，并發(fā)展更為嚴(yán)密的輻射定標(biāo)模型。

志謝： 感謝歐洲空間局提供ESA S1 ToolBox軟件。

[1] Desnos Y L,Li Z Y,Gao Z H,et al.The dragon programme-status and achievements[C]//Proceedings of the Envisat Symposium 2007.Montreux,Switzerland:ESA,2007.

[2] Laur H,Bally P,Meadows P,et al.ERS SAR Calibration.Derivation of the Backscattering Coefficientσ0in ESA ERS SAR PRI Products[R].ESA/ESRIN ES-TN-RS-PM-HL09.2009:3-6.

[3] Cumming I G,Wong F H.Digital Processing of Synthetic Aperture Radar Data:Algorithms and Implementation[M].Boston,MA:Artech House Remote Sensing Library,2005:56-123.

[4] 袁孝康.星載合成孔徑雷達(dá)導(dǎo)論[M].北京:國(guó)防工業(yè)出版社,2003. Yuan X K.Introduction to Synthetic Aperture Radar[M].Beijing:National Defense Industry Press,2003.

[5] Schwerdt M,Brautigam B,Bachmann M,et al.TerraSAR-X calibration-first results[C]//Proceedings of the 2007 IEEE International Geoscience and Remote Sensing Symposium.Barcelona:IEEE,2007:3932-3935.

[6] Rosich B,Meadows P.Absolute Calibration of ASAR Level 1 Products Generated with PF-ASAR[R].ESA-ESRIN,ENVI-CLVL-EOPG-TN-03-0010,2004:5-8.

[7] Slade B.RADARSAT-2 Product Description[R].RN-SP-52-1238,2009:2-25.

[8] Agenzia Spaziale Italiana.COSMO-SkyMed SAR Products Handbook[M].Agenzia:ASI,2007:399-405.

[9] Fritz T,Werninghaus R.TerraSAR-X Ground Segment Level 1b Product Format Specification[R].Berlin:Clustert Applied Remote Sensing(CAF),German Aerospace Center(DLR),2007:156-170.

[10]Rosenqvist A,Shimada M,Suzuki S,et al.Operational performance of the ALOS global systematic acquisition strategy and observation plans for ALOS-2 PALSAR-2[J].Remote Sensing of Environment,2014,155:3-12.

[11]Geudtner D,Torres R,Snoeij P,et al.Sentinel1 system capabilities and applications[C]//Proceedings of the 2014 IEEE Geoscience and Remote Sensing Symposium.Quebec City,QC:IEEE,2014:1457-1460.

(責(zé)任編輯： 邢宇)

Absolute radiometric calibration of level-1 detected ground range products of new SAR sensors

DU Weina1， XU Aigong1， SONG Yaoxin2， SUN Huasheng1

(1.LiaoningTechnicalUniversity，F(xiàn)uxin123000，China； 2.TrafficManagementResearchInstituteofMinistryofPublicSecurity，Wuxi214151,China)

For the current situation of the lack of the new SAR sensor data preprocessing software, this paper introduced in detail the methods of absolute radiometric calibration and the parameter acquisition for several new SAR sensor level-1 detected products, such as ENVISAT ASAR，Radarsat2，Cosmoskymed，TerraSAR-X and Sentinel1. In addition, the absolute radiometric calibration process was achieved by programming with the level-1 detected ground range products(L-1 DGRP) data of Sentinel1 sensor, and C++ programming language was used to achieve the absolute radiation of the calibration process. At last, the radiometric calibration results produced by the method developed in this paper and implemented in the authors’ software were compared with those by ESA S1 ToolBox, the freely distributed SAR data processing tool by European Space Agency, and it is shown that the two numerical back scattering systems are basically the same. The radiometric calibration method developed in this paper is proved to be correct by the program implementation.

new SAR sensors； absolute radiometric calibration； level-1 detected ground range products(L-1 DGRP)； back scattering coefficient； Sentinel1

10.6046/gtzyyg.2016.04.05

杜偉娜,徐愛(ài)功,宋耀鑫,等.新型SAR傳感器一級(jí)地距產(chǎn)品絕對(duì)輻射定標(biāo)方法[J].國(guó)土資源遙感,2016,28(4)：30-34.(Du W N,Xu A G,Song Y X,et al.Absolute radiometric calibration of level-1 detected ground range products of new SAR sensors[J].Remote Sensing for Land and Resources,2016,28(4)：30-34.)

2015-05-13；

2015-07-08

國(guó)家自然科學(xué)基金項(xiàng)目“遙感數(shù)據(jù)的空間分辨率和波段數(shù)對(duì)土地覆蓋制圖的影響研究”(編號(hào)： 41201454)資助。

TP 751.1

A

1001-070X(2016)04-0030-05

杜偉娜(1988-),女,碩士研究生,主要研究方向?yàn)镾AR圖像處理與應(yīng)用。Email： RSwendu@163.com。