應用月平均SST資料對遙感反演數(shù)據(jù)進行質(zhì)量控制

1、MODIS遙感反演SST的比照測試和分析技術報告概述：利用浮標資料對MODIS的SST海外表溫度反演的結(jié)果進行比照測試，是 SeaDAS衛(wèi)星遙感工程的一個重要目標。目的是利用浮標實測海溫資料對SST反演結(jié)果進行校正，以減少反演結(jié)果的誤差，擴大海溫遙感反演產(chǎn)品的使用范圍， 為更好的利用衛(wèi)星遙感資料效勞于海洋監(jiān)測工作建立一個良好的根底。技術背景MODIS 中等分辨率成像光譜儀是一個擁有 36個通道的可見光和紅外波 段光譜輻射計，波段范圍從0.345卩m 可見光到14.235卩m 熱紅外。各通 道的量化等級為12Bits。MODIS第1和第2通道的空間分辨率為250米，第3到第7通道的空間分 辨率為

2、500米，其余通道的空間分辨率為1000米，掃描觀測寬度達2330公里可 以同時提供反映海洋水色，葉綠素、懸浮物濃度等海洋特征信息。MODIS被兩顆對地觀測衛(wèi)星AQUA和TERRA作為有效負載，其中TERRA 是上午星、AQUA是下午星。我單位目前使用SeaDAS5.0軟件進行MODIS數(shù)據(jù)處理。SeaDAS 5.0反演 海外表溫度的算法是美國航空航天局 NASA開發(fā)的“邁阿密探路人算法，有兩 套，分別利用了 3.7卩m和11卩m附近的兩個紅外大氣窗口。其中晝間使用11卩m算法，即利用11.03卩m,12.03卩m通道，夜間使用3.7卩m和11卩m通道。 由于晝間觀測數(shù)據(jù)不僅可提供海外表溫度，

3、而且可提供葉綠素等海洋水色數(shù)據(jù)， 我單位目前僅接收和處理 MODIS的晝間數(shù)據(jù)。我們的比測試驗也是針對晝間的 11卩m算法。“邁阿密探路人在11卩m的海溫反演算法如下：SST = G +6爲1 +G為日1G心其中，SST是反演的海外表溫度，T31是第31通道11.03卩m的輻射亮溫,T32是第32通道12.03卩m的輻射亮溫。B是衛(wèi)星天頂角。C1C4是四個系數(shù)。 SeaDAS軟件中內(nèi)置的各個系數(shù)估計值如下：兩種大氣條件干和濕T32-T31 < L 7T32 - T31 ；、0. 7C1L 228552L 692521C20. 9576555也 9558419C30. M821960. 0

4、873754L 7746311.199584表1:MPSST算法的系數(shù)圖1 :MPSST算法反演獲得的東海區(qū)域海外表溫度圖(2007.10.24 2:38UTC)為了對遙感反演獲得的海外表溫度進行比照測試， 需要選擇適宜的其它測量 手段，要求這些測量手段有足夠的時間分辨率，同時又有較高的精度。目前，現(xiàn)場測量海溫的方法主要有：岸站、志愿船、資料浮標。其中，岸站受地理條件限制，測量的海水溫度是近岸水溫，受大陸、島嶼的 地溫影響較大。志愿船不但資料較少，時間上不連續(xù)，空間分布有不確定性，難 以和遙感資料匹配。最適于海外表溫度比照測試的測量手段是海洋資料浮標。國家海洋局東海分局在我國近海布置了 3個大

5、型海洋浮標。這些資料浮標直接測量水面下20-30厘米的海水溫度，平時每3小時傳輸一次水溫資料 水溫傳感器的技術參數(shù)如下：參數(shù)測量范圍測量精度備注水溫(C)3 + 35±).5 C通訊方式為VSAT表2:海洋資料浮標的水溫傳感器參數(shù)海洋資料浮標由錨系固定，布放完成后坐標根本不變動，這給進行比照測試 帶來了方便。在了解了遙感反演算法和選擇了適宜的比照測量手段后， 我們就可以設計和 實施海溫比照測試。比照測試的方案設計遙感反演數(shù)據(jù)和浮標測量數(shù)據(jù)進行比照的前提是： 在時間和空間上具有可比 性。最理想的情況是，遙感反演數(shù)據(jù)和浮標測量數(shù)據(jù)的時間差和空間距離都為 0 但事實上，由于衛(wèi)星軌道精度存在

6、誤差，衛(wèi)星遙感傳感器分辨率的問題， 遙感反 演數(shù)據(jù)不可能完全吻合浮標測量數(shù)據(jù)。因此，必須容許衛(wèi)星遙感數(shù)據(jù)和浮標測量 的數(shù)據(jù)在時間和空間范圍內(nèi)有一定的誤差， 我們總的要求是，用于進行遙感比照 的數(shù)據(jù)本身和真實海外表溫度之間的誤差限不能大于 1度。我們設計符合比測要求的數(shù)據(jù)標準:工程允許誤差時間15分鐘空間約2千米，0.018緯距表3:選擇遙感數(shù)據(jù)比測的時間和空間標準通常，海外表溫度的時間變化率不超過1攝氏度/小時，空間變化率不超過0.4 度/千米。因此，符合上述標準的數(shù)據(jù)，其誤差限可以表達為E = sqrt(e1A2+e2A2+emA2)(其中E是總誤差限，el是時間誤差限，在我們這里約為 &#

7、177;).25,e2是空間誤差限，約為±).8,em是儀器測量精度的誤差限，為±).5)，因此， 估計由于時間、空間和浮標測量儀器引入的誤差，其誤差限為±1.01度。根本符合要求。未經(jīng)質(zhì)控的遙感反演SST數(shù)據(jù)的比照測試結(jié)果從2007年3月3日到2007年7月27 日，147天內(nèi)，共選取符合篩選標準的數(shù)據(jù)71組，進行了比照測試，全部數(shù)據(jù)細節(jié)詳見附件1結(jié)果的摘要如下：工程結(jié)果樣本數(shù)量71取小誤差0.04取大誤差-25.45誤差的標準差8.59表4：未經(jīng)質(zhì)控的遙感反演數(shù)據(jù)比照測試結(jié)果摘要從上表可以看到，如果對遙感反演獲得的SST不做任何處理，其誤差標準差 可達8C多，

8、其應用價值很低。造成這種情況的原因是多方面的，東海海域是陸緣海，由于海水的深度較淺， 太陽加熱的作用比擬明顯、沿岸上空的大氣狀況比擬復雜，尤其云量較多時，對 紅外輻射吸收嚴重。因此，未經(jīng)質(zhì)控的數(shù)據(jù)不宜直接作為海溫遙感產(chǎn)品使用。應用月平均SST資料對遙感反演數(shù)據(jù)進行質(zhì)量控制為了對遙感反演數(shù)據(jù)進行質(zhì)量控制， 一條合理的思路是，通過經(jīng)驗判定合理 的數(shù)據(jù)區(qū)間，將明顯不合理的數(shù)據(jù)丟棄，這樣能夠大大改善遙感數(shù)據(jù)的質(zhì)量。海水的熱容量巨大，日溫差波動范圍小，每月的波動范圍也十分有限。我 們通過統(tǒng)計F09浮標2004-2006三年的資料其中缺少2004年4、5月份資料的發(fā)現(xiàn)，浮標測量的月度海溫極值 偏離當月海溫

9、平均值一般不超過 4C,從不超過5C,而月度海溫極值距多年的月平均海溫值的偏差一般不超過 3C,從不超過4C。分析發(fā)現(xiàn)F17和F18浮標的資料也滿足上述條件。2月份9多年平均月極值最大偏差1月16.27-1.72月13.98-2.153月13.92-2.84月14.863.935月20.793.076月25.092.817月27.842.058月28.762.849月27.73.010月25.6-2.611月22.902.212月19.742.36表5： F9浮標的多年平均值和月度極值相對平均值的偏差由此，我們可以采取這樣的質(zhì)量控制措施：遙感反演的海溫，在當月海溫 多年平均值芳C范圍內(nèi)的，是有

10、效數(shù)據(jù)，通過質(zhì)量控制，否那么，判定為是無效數(shù) 據(jù)，不做處理。在浮標定點的海域，使用浮標多年的月度平均數(shù)據(jù)作為平均月海 溫F17和F18布置時間短，只取得了 這些數(shù)據(jù)已經(jīng)求岀，見附件 2,在沒有浮標資料覆蓋的海域，可采用大范圍海溫等溫線圖資料，同時質(zhì)量 控制門限可放寬到多年平均值 站C范圍。經(jīng)質(zhì)量控制的遙感SST數(shù)據(jù)的比照測試結(jié)果和回歸分析在符合篩選標準的全部71條數(shù)據(jù)中，通過質(zhì)量控制的共計 49條，占69%, 這說明選取的質(zhì)量控制標準寬嚴適度， 既考慮了數(shù)據(jù)的可靠性，又兼顧了遙感數(shù) 據(jù)的可獲取性。我們對49組經(jīng)質(zhì)控數(shù)據(jù)的比測結(jié)果詳見附件 3,摘要如下：工程結(jié)果樣本數(shù)量49最小誤差0.04最大誤

11、差-6.48誤差的標準差2.49表6:未經(jīng)質(zhì)控的遙感反演數(shù)據(jù)比照測試結(jié)果摘要可見，通過質(zhì)量的數(shù)據(jù)減少了誤差，提高了數(shù)據(jù)的可靠性。為了進一步提高海溫遙感的精度，我們用經(jīng)過質(zhì)控的遙感反演的 SST和浮 標測得的海溫分別作為自變量和函數(shù)進行了線性回歸分析 ，得到了根據(jù)反演SST 結(jié)果預測實際海外表溫度的線性回歸方程。1年的資料，沒有多年的月海溫平均值。反演SST和浮標實測SST回歸分析測-測量結(jié)果回歸方程圖2 MODIS反演SST和實測SST回歸分析回歸統(tǒng)計Multiple R0.950309344R Square0.90308785Adjusted R Square0.90102589標準誤差1.

12、583839791觀測值49表7 MODIS 反演SST和實測SST回歸統(tǒng)計回歸分析的結(jié)果說明，建立的線性回歸方程能夠有效改良反演 SST的精度, 回歸方程的預測值的誤差標準差降低到 1.584C,相關性系數(shù)到達0.90,作為大范 圍海溫監(jiān)測產(chǎn)品，具備了實用價值。結(jié)論SeaDAS 5.0提供的SST反演算法是基于全球海域得出的，對于東海近海海 域，本身存在著一定誤差。又由于 SeaDAS 5.0內(nèi)置的云檢測算法不完善，不能 有效過濾由于水汽吸收而造成的海溫反演異常值，造成的誤差更大。這些問題制約了 MODIS遙感海溫產(chǎn)品的使用。我們通過用浮標數(shù)據(jù)和MODIS遙感反演SST數(shù)據(jù)的比照測試，尋找并確定 了排除反演SST異常值的方法，并建立了根據(jù)多年海溫月平均值的數(shù)據(jù)質(zhì)控措 施。為了進一步提升數(shù)據(jù)精度，我們對經(jīng)質(zhì)控的反演SST和實測值作了線性回歸 分析，求出了根據(jù)反演SST預測實際海溫的線性回歸方程，進一步降低了遙感數(shù) 據(jù)的誤差。雖然這個初步的SST比照測試和改良的工作已經(jīng)到達了工程任務書的要求。 但是限于時間、經(jīng)費等因素，在提高遙感海溫產(chǎn)品的精度，拓展遙感產(chǎn)品的業(yè)務 化應用方面， 仍有很大的提升余地。 例如，可以將數(shù)據(jù)質(zhì)控措施和回歸方程集成 到SeaDAS 5.0軟件內(nèi)部，實現(xiàn)高精度的SST產(chǎn)品的穩(wěn)定輸出，以及啟動水汽吸 收更小的夜間4卩m海溫資料等。我們希望在后續(xù)的工程中能夠