剖面法與水面之上法測(cè)量水面下遙感反射率的比較.docx

1、第46卷增刊】廈門大學(xué)學(xué)報(bào)(自然科學(xué)版)2007年8月JournalofXiamenUniversity(NaturalScience)Vol.46Sup.1Aug,2007剖面法與水面之上法測(cè)量水面下遙感反射率的比較隋曉飛微，商少平帶，馬曉鑫3,商少凌(1.近海海洋環(huán)境科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室(廈門大學(xué))，2.廈門大學(xué)海洋與環(huán)境學(xué)院福建廈門361005»3.YS】公司.福建廈門361005)摘要：臺(tái)灣海峽、福建沿岸及兩個(gè)淡水湖的表面和剖面光譜測(cè)敏結(jié)果表明，分別由水面之上法和剖面法得到的恰于水面下的遙感反射率"”(s)和r(p»之間存在一定差異.在債值上，波段平均的無偏

2、百分比誤差參數(shù)在1.7%68.3%之間.但50%以上站次<20%$在光譜形狀上兩者有很好的相似性各站次的®都在0.93以上：490nm和555nm波段比無偏百分比誤差參數(shù)勤(M/?！?在0.5%1】.0%之間，且近70%的站次做GU/卷)<5%.在平靜湖泊與、(P)量值的統(tǒng)計(jì)差與美國(guó)NASA水色遙感工作組對(duì)離水輻射率所做的類似統(tǒng)計(jì)相當(dāng)，但在海洋水體的結(jié)果遠(yuǎn)大于此，分析測(cè)量環(huán)境條件變化造成的光學(xué)參數(shù)測(cè)量誤差是主要原因.關(guān)鍵詞：水面之上法；剖面法；遙感反射率中圖分類號(hào):TP722.4文獻(xiàn)標(biāo)識(shí)碼：A恰于水面下的遙感反射率(remotesensingreflectancejust

3、belowthesurface,rn)是主要的水體表觀光學(xué)參數(shù)之一,也是眾多海洋水色遙感算法中的一個(gè)基礎(chǔ)物理量，表達(dá)如下：f=L,(0-)/Ed(0-)(1)其中L“和E分別為恰于水面下的向上輻亮度和向下福照度與離水輻亮度(LQ相比，由于能夠部分消除不同入射光照條件對(duì)光譜數(shù)據(jù)的影響，因此更能反映水體不同成分的吸收和散射信息，從而在反演水體成分濃度的算法中得到更廣泛的應(yīng)用.然而現(xiàn)場(chǎng)光譜測(cè)量很難直接測(cè)錄到恰于水面下的光學(xué)量，只能通過其它手段間接得到.目前，主要有兩種方法可以獲得rn：一是通過海表遙感反射率(RQ推算得到，稱之為水面之上法；二是利用剖面光譜儀測(cè)量水下不同水層的上和E，以此反推得到恰于

4、水面下的匚和E”從而得到小，稱之為剖面法.水面之上法利用下式計(jì)算兇：Dr”=丁+劫"(2)其中T=t-t+/n2"一是從水表面下到水表面上的輻射透過率是從水表面上到水表面下的輻射透過率，為水體折射系數(shù)"是水一氣界面反射率;Q是光場(chǎng)分布參數(shù)，為上行輻亮度與上行輻照度的比值.收稿日期：2007-04-10M金項(xiàng)目：國(guó)家自然科學(xué)基金(40376031,40331004,40521003),863it劃(2002AA639540,2006AA09A302),教育部創(chuàng)新群體計(jì)劃，廈門大學(xué)新世紀(jì)人才計(jì)劃聯(lián)合資助通機(jī)作者：sishang©文章編號(hào):0438-0479(

5、2007)S1-OOO6-O61995年，Mobley通過HYDROLIGHT對(duì)R”和"，進(jìn)行計(jì)算，認(rèn)為對(duì)于光學(xué)深水和位于天底方位觀測(cè)的傳感器而言Tq0.52,7Qq1.7叫剖面法利用剖面儀測(cè)量得到的水下不同水層的L.和E，通過局部線性回歸法進(jìn)行線性擬合，反推得到恰于水面下的L.(O-)和耳(0)，通過式(1)計(jì)算得到疽七水面之上法具有現(xiàn)場(chǎng)操作簡(jiǎn)單、現(xiàn)場(chǎng)實(shí)驗(yàn)成本低等優(yōu)點(diǎn)，但受海面、天氣條件等的影響大，會(huì)給數(shù)據(jù)結(jié)果帶來很多不確定的因素.相對(duì)而言，剖面法則受現(xiàn)場(chǎng)環(huán)境因素影響較小，但也有其局限性，如受儀器自身的自陰影，水下儀器的姿態(tài)等因素影響而對(duì)數(shù)據(jù)結(jié)果帶來一定的誤差.關(guān)于兩種方法所獲得的

6、U的一致性，已見數(shù)例報(bào)道，而對(duì)于遙感半分析反演算法的關(guān)鍵參數(shù)Z，兩種來源數(shù)據(jù)的一致性究竟如何，迄今未見報(bào)道.本研究利用GER1500型地物高光譜儀以及基于GER1500型地物高光譜儀研發(fā)的光纖式剖面光譜測(cè)量系統(tǒng)3】，在臺(tái)灣海峽南部海域、福建沿岸水域和兩個(gè)陸地平靜湖泊同步測(cè)量了表面光譜數(shù)據(jù)和剖面光譜數(shù)據(jù)，對(duì)分別由水面之上法和剖面法得到的結(jié)果進(jìn)行比較，并分析其影響因素.1方法1.1儀器與測(cè)量方法在水面之上法中，光譜測(cè)量采用的儀器為GER1500型便攜式高光譜儀，其參數(shù)如表1所示.表1GER1500參數(shù)表Tab.1TheparametersofGER1500儀器型號(hào)探測(cè)器光譜分辨率/nm光譜范圍波段

7、數(shù)視場(chǎng)角掃描時(shí)間/ms尺寸/cmGER1500512Si1.535010505124°（無光纖）23°（帶光纖），5（可設(shè)定）15X8X26由于目標(biāo)水域基本上都屬于二類水體，所以表面光譜測(cè)量按照唐軍武等（2004年1月）提出的二類水體水面以上光譜的測(cè)量方法.儀器觀測(cè)平面與太陽(yáng)入射平面的夾角為40°<%<135°（背向太陽(yáng)方向），儀器與海面法線方向的夾角30°<%<45°.這種觀測(cè)幾何也是目前國(guó)際水色SIMBIOS計(jì)劃中推薦的方法，其優(yōu)點(diǎn)是：與剖面觀測(cè)的固有差異較小，能更好地避免太陽(yáng)直射反射和船體自身對(duì)水體光場(chǎng)的

8、破壞，缺點(diǎn)是天空光的均勻性較差，船舶的晃動(dòng)會(huì)產(chǎn)生一定的影響.而后通過下式獲得R”（Q：pj一Lw（人）._>4財(cái)（“o=再（"其中，L”為離水輻射率，E（0+，人）為海面入射輻照度，L,為水面向上總輻亮度，Lb為天空福亮度，L.ye為漫反射板輻亮度"=50%（漫反射板反射率）/為氣一水界面對(duì)天空光的反射率，取值0.028“*.剖面光譜測(cè)量采用針對(duì)近海高污水體光譜測(cè)量開發(fā)的光纖式剖面高光譜測(cè)量系統(tǒng)，該系統(tǒng)的設(shè)計(jì)思路是以光纖為媒介，通過光纖將水下光信號(hào)引出至水面上的表面光譜儀內(nèi)，使原本用于測(cè)量表面光譜的GER1500型表面光譜儀擴(kuò)展為可以進(jìn)行水下剖面光譜測(cè)量的剖面光譜儀，

9、其概念圖如圖1所示，詳細(xì)說明參見文獻(xiàn)8.其中，在兩個(gè)海上航次的剖面光譜測(cè)量中，使用的是完全版的剖面光譜測(cè)量系統(tǒng)，其光纖總長(zhǎng)15m,浮子系統(tǒng)靠長(zhǎng)3m的硬質(zhì)連桿推離船體，水下探頭依靠浮子系統(tǒng)維持海面0基準(zhǔn)的控制，用液位計(jì)測(cè)量水下探頭的水深；而在兩個(gè)淺水平靜湖泊進(jìn)行測(cè)量時(shí)，使用的是簡(jiǎn)易版的剖面光譜測(cè)量系統(tǒng)，其光纖總長(zhǎng)5m,利用光纖上標(biāo)志深度的標(biāo)簽觀測(cè)探頭的水深囚1.2數(shù)據(jù)采集與處理采用上述儀器與方法進(jìn)行了4個(gè)批次的測(cè)量，分別是2006年5月到7月廈門大學(xué)校園內(nèi)的兩個(gè)平靜淡水湖一芙蓉湖和情人湖（圖2a）.2004年7月25日至8月6日臺(tái)灣海峽南部（圖2b）.2004年9月1018日福建沿岸（圖2c）,

10、總測(cè)量站次為77次，得到表面與剖面光譜匹配、且通過質(zhì)景控制的數(shù)據(jù)，共21組.按上述兩種不同方法計(jì)算恰于水面下的遙感反射率，其中由水面之上法得到的標(biāo)志為七（s），由剖面法得到的標(biāo)志為腿）.圖1剖面測(cè)量系統(tǒng)概念圖Fig.1Theconceptualchartofprofilingspectroradiometer1.3比較方法為了進(jìn)一步分析比較r（p）和r（S）之間的差異，參照StanfordB.Hooker等人在評(píng)價(jià)水面之上法和剖面法得到的離水輻射率時(shí)采用的方法，從量值和光譜形狀變化的角度去分析比較兩者之間的異同.用400700nm波長(zhǎng)范圍內(nèi)共188個(gè)波段的平均無偏百分比誤差參數(shù)（Unbiase

11、dpercentdifference,簡(jiǎn)稱UPD）表征量值上的差異，用rn（p）和rn（s）的R2以及490.555nm波段比的無偏百分比誤差參數(shù)來說明在光譜形狀上的差異.計(jì)算公式如下：#（*）一200"胡成（義,刃（4）媯（偵）一20。xA（®+xB（m式中X'和XB分別表示由兩種不同途徑得到的光譜平均的UPD由下式計(jì)算：盟（志£所（5）490.555nm波段比的UPD由下式計(jì)算：（人490/人555）=200X|XA（%9o）/XAS55）-XB（Lo）/XB（%G|XA（A49o）/XA（A555）+XBa）/XB（原55）（6）2結(jié)果與討論各站次r

12、n（s）和rn（p）隨波長(zhǎng)變化的光譜曲線圖如圖3所示，二者比較的誤差統(tǒng)計(jì)結(jié)果如表2所示.從光譜曲線圖（圖3）粗略分析，匕（s）和七3）存在一定差異，有的站次rn（s）高于（力），而有的站次相反.圖2各航次站位圖a.情人湖、芙蓉湖位置圖；b.福建沿岸航次站位圖；c.臺(tái)灣海峽航次站位圖Fig.2Samplingstations就f（$）與rn（p）波段平均的量值差異而言，波段平均無偏百分比誤差明（入）變化顯著，范圍在1.7%68.3%之間，平均值為25.6%，但50%以上站次<20%.其中兩個(gè)平靜淡水湖的在1.7%13.2%之間，平均值為7.3%,與StanfordB.Hooker等人對(duì)離水

13、輻射率進(jìn)行比較的結(jié)果（3.2%17%）相似.而在臺(tái)灣海峽水域和福建沿岸（Atll）在5.9%68.3%之間，變化范圍很大且平均值為38.37%,遠(yuǎn)大于Hooker等人的離水輻射率比較結(jié)果.這與現(xiàn)場(chǎng)測(cè)量條件、光場(chǎng)環(huán)境以及目標(biāo)水體隨站位和測(cè)量時(shí)間不同而引起的太陽(yáng)高度角的差異等實(shí)驗(yàn)背景的不同而呈現(xiàn)的各異性有很大的關(guān)系.這些實(shí)驗(yàn)背景的差異不僅影響到直接測(cè)量參數(shù)的準(zhǔn)確性，同時(shí)對(duì)隨后的遙感反射率推算產(chǎn)生影響.具體而言，在計(jì)算rn時(shí)使用的一些參數(shù)，如公式（2）中的光場(chǎng)分布參數(shù)Q以及公式（3）中的氣一水界面反射率/等，是隨著現(xiàn)場(chǎng)測(cè)房條件（如太陽(yáng)角度，海面風(fēng)速等）的不同而變化的，但這種變化，目前很難度量，本文仿

14、通常的做法文獻(xiàn)2等采用固定值，因此我們所觀察到的差異在一定程度上可能是由Q和廣的變化帶來的.表3列舉了有代表性的5個(gè)站次實(shí)驗(yàn)時(shí)的現(xiàn)場(chǎng)條件，其中060514情人湖站次符合準(zhǔn)確測(cè)量光學(xué)參數(shù)的標(biāo)準(zhǔn)條件，現(xiàn)場(chǎng)光學(xué)參數(shù)測(cè)量準(zhǔn)確，無偏百分比誤差很?。?.7%山而其它幾個(gè)站次因?yàn)樘?yáng)角度過大，云移動(dòng)導(dǎo)致光場(chǎng)變化，水流過急或水面有白沫等不同原因而使現(xiàn)場(chǎng)光學(xué)測(cè)量和后期rn的計(jì)算結(jié)果誤差增大，導(dǎo)致<A.h）均在40%以上.就光譜形狀的相似度而言，所有的21個(gè)站次數(shù)據(jù)，R2均在0.93以上，其中有15個(gè)站次在0.99以上.各站次490.555nm波段比值的無偏百分比誤差0j（A49O/A555）范圍在0.5%

15、11.0%之間，且近70%的站次5%,遠(yuǎn)小于上述量值上的差異，與Hooker等人1.4%7.2%的結(jié)果在可比范圍之內(nèi).這表明每個(gè)站次的"（P）和i（s）之間雖然在量值上存在著差異，但其差異在400700nm波段范圍內(nèi)幾乎不存在波長(zhǎng)依從特征.綜合衡ft分析各站次譜型圖、線性相關(guān)系數(shù)、W0Q以及02<A49O/A555）,-方面表明恰于水面下表2各站次誤差統(tǒng)計(jì)結(jié)果Tab.2Statisticsofrn(/»)andr(.s')站次R2回（如）/%(490/S5$)/%060512-芙蓉湖0.998912.52.8060514-情人湖0.99971.71.00606

16、11-芙蓉湖0.99892.30.4060611-情人湖0.995313.29.6060706-芙蓉湖0.99517.010.7040801-F30.945861.16.4040802-F30.998317.64.5040804-F20.999520.83.2040804-F30.994249.04.8040804-F40.978529.611.0040805-G10.969017.53.9040805-F00.994217.67.6040805-F10.998747.85.4040805-F20.999825.60.87040805-F30.988914.00.8040805-F40.977

17、668.34.6040806-G10.991739.97.8040806-F10.995427.51.1040914-F30.93935.90.5040914-F120.986127.11.4040917-F630.993632.58.1的遙感反射率4的光譜曲線形態(tài)，光譜波段比值的測(cè)量可以達(dá)到較高的質(zhì)量水平，但七的量值測(cè)酸還存在較大的誤差.另一方面，兩個(gè)陸地淡水湖的實(shí)驗(yàn)結(jié)果優(yōu)于臺(tái)灣海峽南部水域和福建沿岸海域的實(shí)驗(yàn)結(jié)果，分析其原因主要有兩方面，一是相對(duì)于海上測(cè)量而言在平靜的陸地淡水湖測(cè)量時(shí)，可以比較精確的把握現(xiàn)場(chǎng)測(cè)量幾何，而且對(duì)于測(cè)量結(jié)果影響很大的船的晃動(dòng)、自陰影和海面波浪起伏等客觀因素的影響會(huì)

18、徹底消除或大大減弱'二是由于兩個(gè)湖都是高污渾濁水體，其490nm的水體漫衰減系數(shù)K（490）基本都在2m'1以上，光學(xué)深度在0.30.5m內(nèi)變動(dòng)，因此可以忽略底部和水邊界的影響.雖然光學(xué)深度很淺，但是通過提高深度分辨率（湖里測(cè)量為0.05m,海上測(cè)量為1m）,依然可以推導(dǎo)出相對(duì)可靠的L”和E的垂直變化,從而得到更為準(zhǔn)確的m3）.至此我們知道S（s）與七（p）之間主要存在量值上的差異，且此差異基本不隨波長(zhǎng)而變.我們認(rèn)為產(chǎn)生差異的主要原因，一方面是與波段無關(guān)的測(cè)量幾何波動(dòng)和環(huán)境條件變化導(dǎo)致表面法測(cè)量的R”和剖面法直接測(cè)量的L.S）和E（A,z）的測(cè)量誤差引起，另一方面由測(cè)最環(huán)境條件

19、不同導(dǎo)致的Rr到轉(zhuǎn)換參數(shù)的變化也可能產(chǎn)生一定影響.為此我們嘗試做一個(gè)估算，公式（2）中，丁和相對(duì)來說較為固定，我們?nèi)匀《ㄖ?Q則隨不同水體、太陽(yáng)角度以及儀器觀測(cè)角度在1.77之間變動(dòng)，以Rn值較大的O4O914-F12站次為例，Q分別取1.7、3.75（0R1.7）、7計(jì)算f（s）,并與該站次的仁3）做比較，如圖4所示.可見Q取3.75.7時(shí)，計(jì)算得到的rn（5）差異較大，兩者的也（）為4.2%,而媯（扁,/A炊）為0.3%.表2所示3）與rn（p）之間差異的統(tǒng)計(jì)結(jié)果,Q的取值是3.75,對(duì)于該站位rn（s）與rn（p）的為27.1%,網(wǎng)CUo/%55）為1.4,可見由Q的變化引起的誤差對(duì)04

20、0914-F12站次的r（s）與rn（p）之間的總體差異的貢獻(xiàn)相對(duì)較小.3結(jié)論本文利用GER1500型便攜式高光譜儀和一套主機(jī)為GER1500型便攜式高光譜儀的光纖式剖面光譜測(cè)量系統(tǒng)，分別在臺(tái)灣海峽南部，福建沿岸海域以及兩個(gè)陸地淡水湖中進(jìn)行了表面和剖面的光譜測(cè)量，積累了21個(gè)站次的同步數(shù)據(jù)，分別計(jì)算出了“3）和f（s），并利用波段平均的無偏百分比誤差或（心）和波段比值的無偏百分比誤差（網(wǎng)（9。/%55）等統(tǒng)計(jì)分析方法，初步分析比較了兩者之間的異同.與美國(guó)NASA水色遙感工作組對(duì)離水輻射率所做的類似統(tǒng)計(jì)相比，海上的網(wǎng)（9。/%55）遠(yuǎn)表3幾個(gè)站次現(xiàn)場(chǎng)實(shí)驗(yàn)記錄Tab.3Recordofthesam

21、plingconditionsforacoupleofrepresentativestations站位測(cè)量時(shí)間海況云信風(fēng)速060514-情人湖13：45水面平靜，水流速很小晴.無云,光照很強(qiáng)近于無風(fēng)040801-F316,10有波紋無波頂晴，無云光照較強(qiáng)微風(fēng)040805-F416：30有波紋無波頂晴，少fit薄云光照很強(qiáng)微風(fēng)040805-F111,25細(xì)小波紋.流速較急，儀器有傾角晴.大做厚云，測(cè)鼠時(shí)光場(chǎng)有變化光照很強(qiáng)微風(fēng)040806-G110：06水面平靜無波，拍剖面時(shí).水面有些許白泳晴，無云光照強(qiáng)烈微風(fēng)0.020040914-F3040914-F120606U-FR0.0150.0100.

22、0050.0200.0150.0100.005060611-QR4000.0140.0120.0100.0080.0060.0040.0020040804-F3040917-F63014012010008006004002o.o.o.oo.ao.Wavelength/nm圖3各站次和r(j)諧形圖(p)IFig.3Ther(p)andrn(s)ofeachstationWavelength/nm圖4040914-F12站r(p)和不同Q值下的、($)Fig.4Thern(p)andrn(s)fordifferentQof040914-F12大與其統(tǒng)計(jì)結(jié)果，但海上的R2、婢GU/幅)以及平靜湖泊

23、的3個(gè)誤差統(tǒng)計(jì)都與之相當(dāng).說明恰于水面下的遙感反射率光譜M的曲線形態(tài)、光譜波段比值的測(cè)量質(zhì)量可以達(dá)到較高的水平，而量值的測(cè)量:質(zhì)量不容易控制，其原因主要是由與波段無關(guān)的測(cè)景幾何、環(huán)境條件變動(dòng)等原因而造成的光學(xué)量測(cè)成誤差和f計(jì)算公式中光學(xué)參數(shù)的變化引起的.進(jìn)一步根據(jù)可能的光學(xué)參數(shù)變化范圍估算的結(jié)果表明，光學(xué)參數(shù)變化對(duì)總體誤差的貢獻(xiàn)相對(duì)較小，所以盡管衛(wèi)星遙感類似于采用水面之上法，但對(duì)于主流的多波段比值水色遙感算法而言，其與剖面法的差異對(duì)反演結(jié)果造成的影響是可以忽略的.致謝：現(xiàn)場(chǎng)數(shù)據(jù)采集中得到董強(qiáng)和朱任等同學(xué)的大力協(xié)助與支持.謹(jǐn)此一并致謝.參考文獻(xiàn)：1 周良明,劉玉光，郭佩芳.渤海和北黃海水色光語(yǔ)的

24、特征分析J.海洋通報(bào).2005.24(2)：1319.2 GordonHR.BrownOB.EvansRH.etal.Asemi-analyticradiancemodelofoceancolorJ.JournalofGeophysicalResearch.1988,93»1090910924.3 LeeZP.CarderKL.ArnoneRA.Derivinginherentopticalpropertiesfromwatercolor：amultibandquasi-analyticalalgorithmforopticallydeepwatersJ.AppliedOptics.

25、2002.41(27),5755-5772.4 楊安安.李銅基.陳清蓮，等.二類水體表觀光學(xué)特性的測(cè)量與分析一剖面法方法研究J.海洋技術(shù)，2005,26(3):111-115.5 唐軍武，田國(guó)良，汪小勇.等.水體光譜測(cè)重與分析I；水面以上測(cè)量法J.遙感學(xué)報(bào).2004.8(1)：3744.6 李銅基,陳清蓮.韓冰.海面之上法中海面反射率和觀測(cè)幾何模擬研究口.海洋技術(shù),2005,24(4):112-117.7 汪小勇，李銅基.唐軍武，等.二類水體表觀光學(xué)特性的測(cè)量與分析一水面之上法方法研究口.海洋技術(shù)，2004.23(2),1-6.8 馬曉鑫.光纖式高光譜剖面測(cè)量系統(tǒng)D.廈門|廈門大學(xué).2006.

26、9 HookerSB.LazinG,ZibordiG,etal.Anevaluationofa-boveandin-watermethodsfordeterminingwater-leaving.radiancesJ.JournalofAtmosphericandOceanicTechnology,2002,19：486-513.TheComparisonBetweenAbove-andIn-WaterMethodsforDeterminingRemoteSensingReflectanceJustBelowSurfaceSUIXiao-fei】'2,SHANGShao-ping1-2

27、,MAXiao-xin3,SHANGShao-ling1-(1.StateKeyLaboratoryofMarineEnvironmentalScience,XiamenUniversity,2. CollegeofOceanographyandEnvironmentalScience,XiamenUniversity.Xiamen361005.China；TheYSIChinaInc,Xiamen361005.China)Abstract：WeusedGER1500spectrometerandahome-madeprofilertomeasureremotesensingreflectanceintwopools.the