黃波 鄧冉

摘? 要：長(zhǎng)江口是上海市重要的水源地，是東海帶魚(yú)的主要產(chǎn)卵地，對(duì)周邊農(nóng)業(yè)、漁業(yè)生產(chǎn)及居民生活都有著重大影響。該文基于HJ-1B/IRS數(shù)據(jù)，采用普適性單通道算法進(jìn)行長(zhǎng)江口溫度反演，利用MODIS溫度產(chǎn)品MOD28和對(duì)比分析。結(jié)果顯示：該方法反演的長(zhǎng)江口溫度和MOD28產(chǎn)品溫度誤差為+0.23℃，相關(guān)性為0.94，反演結(jié)果良好。

關(guān)鍵詞：長(zhǎng)江口;溫度;普適性單通道算法;HJ-1B/IRS;MOD28

中圖分類號(hào) X87文獻(xiàn)標(biāo)識(shí)碼 A文章編號(hào) 1007-7731（2019）12-0145-2

1 引言

長(zhǎng)江口是東海帶魚(yú)的主要產(chǎn)卵地。長(zhǎng)江口的地表及其海表溫度，對(duì)周邊氣候變化以及農(nóng)業(yè)、漁業(yè)生產(chǎn)和人們的生活都有著重大的影響。長(zhǎng)江口的主要地表覆蓋特征決定了其氣候具有周期性和不規(guī)則的變化特征，而其主要取決于熱收支狀況及其時(shí)間變化。研究長(zhǎng)江口溫度，既是研究大氣、長(zhǎng)江口海氣交換的1個(gè)基本參數(shù)，也是了解該區(qū)域氣候變化的一個(gè)重要標(biāo)志。

HJ-1B衛(wèi)星是中國(guó)環(huán)境小衛(wèi)星中重要的環(huán)境和災(zāi)害遙感監(jiān)測(cè)工具。對(duì)于溫度反演，有學(xué)者基于單通道算法提出了適用于所有熱紅外波段德普適性單通道算法計(jì)算溫度。段四波等在2008年修改了算法中的大氣函數(shù)，得到適用于我國(guó)環(huán)境一號(hào)衛(wèi)星紅外相機(jī)中單通道反演溫度的算法。為此，本次研究采用該算法，對(duì)HJ-1B/IRS數(shù)據(jù)進(jìn)行長(zhǎng)江河口海表溫度反演，并且利用MOD28數(shù)據(jù)驗(yàn)證反演結(jié)果[1]。

2 研究區(qū)及數(shù)據(jù)

2.1 研究區(qū) 研究區(qū)為長(zhǎng)江口，為長(zhǎng)江在東海入海口的一段區(qū)域，是上海市重要的水源地，兼具工農(nóng)業(yè)取水、圍墾、通航、生態(tài)、納污等多種功能。處于我國(guó)重要的經(jīng)濟(jì)區(qū)長(zhǎng)三角地區(qū)的核心區(qū)域，地理位置特殊，是中國(guó)最重要的港口之一。

2.2 研究數(shù)據(jù) 研究數(shù)據(jù)是環(huán)境與災(zāi)害監(jiān)測(cè)預(yù)報(bào)小衛(wèi)星星座的1顆衛(wèi)星HJ-1B，該衛(wèi)星搭載了CCD相機(jī)和紅外相機(jī)（IRS），采用的數(shù)據(jù)時(shí)間為冬天無(wú)云的晴天，時(shí)間為2011年1月12日上午10∶26。采用Terra衛(wèi)星上搭載的MODIS 1B數(shù)據(jù)用作研究中的水汽反演數(shù)據(jù)，同時(shí)采用MODIS的標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品MOD28用于最后的反演結(jié)果驗(yàn)證。

3 HJ-1B衛(wèi)星數(shù)據(jù)長(zhǎng)江口溫度反演方法

本次HJ-1B衛(wèi)星的溫度反演，主要采用HJ-1B衛(wèi)星的IRS傳感器中的第4波段，通過(guò)前人的研究的普適性單通道算法參數(shù)，得到長(zhǎng)江口溫度，此算法中考慮到大氣的影響。首先，在資源環(huán)境應(yīng)用中心下載HJ-1A1B（環(huán)境小衛(wèi)星）數(shù)據(jù)處理程序，ENVI功能擴(kuò)展中加入程序?qū)?shù)據(jù)進(jìn)行預(yù)處理，在預(yù)處理階段主要是涉及到輻射校正和幾何校正，輻射校正。對(duì)處理的結(jié)果用普適性單窗算法可以得到SST，該算法的主要公式如下：

在水汽含量參數(shù)的確定方面，通過(guò)資料文獻(xiàn)的閱讀，對(duì)前人的研究做出分析，發(fā)現(xiàn)HJ-1B數(shù)據(jù)做水汽含量的估算比較復(fù)雜，本次采用在MODIS同時(shí)期的數(shù)據(jù)得到水汽含量，MODIS的水汽產(chǎn)品和通過(guò)MODIS 1B數(shù)據(jù)做出來(lái)的水汽含量對(duì)比，發(fā)現(xiàn)用MODIS 1B數(shù)據(jù)反演出來(lái)的同時(shí)期的水汽含量更加適合HJ-1B/IRS數(shù)據(jù)海表溫度反演研究。因此，可以用MODIS 第2和19波段來(lái)反演大氣水汽含量，對(duì)于MODIS圖像中的任何一個(gè)像元，其大氣水汽含量用下式估計(jì)：

4 結(jié)論與討論

對(duì)于HJ-1B數(shù)據(jù)，采用普適性單通道算法反演2011年1月12日長(zhǎng)江口溫度。為了得到的結(jié)果方便和MOD28數(shù)據(jù)對(duì)比分析，將反演結(jié)果分辨率從150m降到1km，得到結(jié)果如圖2所示（左）。對(duì)于MODIS溫度產(chǎn)品MOD28，對(duì)數(shù)據(jù)采用SeaDAS軟件進(jìn)行格式轉(zhuǎn)換，轉(zhuǎn)換格式后用ENVI對(duì)圖像進(jìn)行幾何校正和裁剪，利用波段運(yùn)算對(duì)像元的灰度值進(jìn)行修改，使其成為溫度值，得到溫度結(jié)果如圖2所示（右）。

根據(jù)反演結(jié)果可以看出，HJ-1B/IRS和溫度產(chǎn)品MOD28具有一定的一致性，兩者長(zhǎng)江口的溫度都基本集中在6～8℃。對(duì)于HJ-1B/IRS數(shù)據(jù)，反演結(jié)果溫度最高為10.75℃，平均溫度為2.54℃;對(duì)于溫度產(chǎn)品，反演結(jié)果溫度最高為10.82℃，平均溫度為2.31℃。HJ-1B/IRS反演結(jié)果相對(duì)產(chǎn)品溫度誤差為+0.23℃，利用ENVI相關(guān)性分析結(jié)果顯示相關(guān)性為0.94。因此，HJ-1B衛(wèi)星數(shù)據(jù)采用普適性單通道算法在冬季反演長(zhǎng)江口溫度可行。

參考文獻(xiàn)

[1]段四波，閻廣建，錢永剛，等.利用HJ-1B模擬數(shù)據(jù)反演地表溫度的兩種單通道算法[J].自然科學(xué)進(jìn)展，2008（09）：1001-1008.

[2]黃妙芬，毛志華，邢旭峰，等.HJ-1B/IRS水溫反演模型及監(jiān)測(cè)示范[J].國(guó)土資源遙感，2011（02）：81-86.

[3]龔紹琦，陸品廷，孫德勇，等.HJ-1B/IRS熱紅外數(shù)據(jù)反演太湖水溫的方法比較[J].湖泊科學(xué)，2016（03）：645-653.

（責(zé)編：張宏民）