常雅軒，于穎，崔綠園，侯吉宇，曾德偉

(東北林業(yè)大學(xué) 林學(xué)院，哈爾濱 150040)

基于HJ-1A星HSI數(shù)據(jù)的植被吸收光合有效輻射APAR估算

常雅軒，于穎*，崔綠園，侯吉宇，曾德偉

(東北林業(yè)大學(xué) 林學(xué)院，哈爾濱 150040)

為了更準(zhǔn)確地測(cè)量植被對(duì)太陽(yáng)輻射的有效吸收量，本文對(duì)基于HJ-1A衛(wèi)星的HSI數(shù)據(jù)進(jìn)行植被吸收光合有效輻射APAR的研究，為陰陽(yáng)葉的吸收有效光合輻射量的分離計(jì)算提供基礎(chǔ)，并為森林生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力計(jì)算提供準(zhǔn)確的參數(shù)數(shù)據(jù)。研究的主要步驟有以下三步：首先對(duì)一景包含漠河氣象站的HSI數(shù)據(jù)進(jìn)行預(yù)處理，其中包括HSI數(shù)據(jù)的格式轉(zhuǎn)換，去條帶處理以及大氣校正等；在此基礎(chǔ)上應(yīng)用逐日太陽(yáng)輻射模型計(jì)算日總太陽(yáng)輻射，并通過(guò)NDVI計(jì)算得到光合有效輻射吸收比率FPAR；最后用FPAR與太陽(yáng)輻射共同計(jì)算出APAR。最終得到處理后的圖像，并且驗(yàn)證了逐日太陽(yáng)輻射模型和研究區(qū)域氣象站數(shù)據(jù)的相關(guān)性，通過(guò)NDVI和SR的結(jié)合計(jì)算出FPAR，精度得到了提高；植被吸收光合有效輻射APAR計(jì)算完成后，將灰度圖像顯示為專題圖像。從專題圖中可以得出不同地區(qū)森林的APAR值，并最終推算出全年的APAR。

HSI數(shù)據(jù)處理；太陽(yáng)輻射模型；FPAR；APAR

0 引言

吸收光合有效輻射(Absorbed Photosynthetically Active Radiation，APAR)指植被冠層吸收的參與光合生物量積累的光合有效輻射部分[1]。森林生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力在研究碳循環(huán)和全球氣候變化中具有重要地位，綠色植物冠層吸收光合有效輻射是CASA等碳循環(huán)模型中計(jì)算凈初級(jí)生產(chǎn)力的基礎(chǔ)參數(shù)，因此，準(zhǔn)確估算區(qū)域APAR對(duì)研究區(qū)域森林生態(tài)系統(tǒng)的光合作用、凈初級(jí)生產(chǎn)力和碳循環(huán)具有十分重要的意義，并且研究APAR是進(jìn)一步研究植被陰陽(yáng)葉APAR分離與計(jì)算的基礎(chǔ)[2]。最初，由于輻射觀測(cè)點(diǎn)極少，獲得APAR的觀測(cè)資料較少，人們主要運(yùn)用氣候?qū)W計(jì)算方法，包括理論和經(jīng)驗(yàn)方法，獲得任一地區(qū)的APAR分布情況，但具有依賴于地面的氣候觀測(cè)資料、計(jì)算方法參數(shù)多、計(jì)算過(guò)程復(fù)雜的缺陷，后來(lái)，一部分學(xué)者開(kāi)始嘗試?yán)眠b感信息來(lái)反演大范圍地區(qū)的APAR分布。

Hipps等研究了冬小麥抽穗前后APAR與LAI的函數(shù)關(guān)系[3]；Meeree等研究得到植物的光合作用產(chǎn)物與APAR呈近似線性關(guān)系[4]；Chen(1999)根據(jù)兩葉模型，將冠層植被分為陰葉和陽(yáng)葉，分別計(jì)算陰陽(yáng)葉的太陽(yáng)輻射，可用于陰陽(yáng)葉APAR的分離與計(jì)算[5]。在國(guó)內(nèi)，張喜旺等利用MODISNDVI數(shù)據(jù)研究了2010年河南省APAR分布[6]；李剛等基于MODIS數(shù)據(jù)研究了內(nèi)蒙古地區(qū)APAR及其時(shí)空變化特征[7]；陳新芳等研究對(duì)于大范圍APAR的監(jiān)測(cè)和估算，通過(guò)對(duì)PAR和基于遙感數(shù)據(jù)計(jì)算的FPAR來(lái)實(shí)現(xiàn)[8]。本研究基于高光譜遙感數(shù)據(jù)分三步估算APAR。高光譜遙感數(shù)據(jù)具有高光譜分辨率，被廣泛應(yīng)用于植被長(zhǎng)勢(shì)檢測(cè)及植被生物化學(xué)信息提取，環(huán)境資源減災(zāi)HJ-1A衛(wèi)星搭載超光譜成像儀(Hyper-SpectralImager，簡(jiǎn)稱HSI)，其數(shù)據(jù)具有高光譜分辨率、較短的時(shí)間分辨率和較低的信噪比，且對(duì)廣大用戶免費(fèi)下載，是本文研究大興安嶺漠河地區(qū)很好的遙感數(shù)據(jù)源。

1 研究區(qū)域概況及數(shù)據(jù)

研究區(qū)域是黑龍江省大興安嶺地區(qū)，該地區(qū)位于東經(jīng)121°12′～127°00′、北緯50°00′～53°33′。東部緊挨小興安嶺，西部和內(nèi)蒙古緊密相連，南部與松嫩平原交界，北部與俄羅斯只有一江之隔。大興安嶺地區(qū)平均海拔573 m，年平均氣溫-2.6℃，年平均降水量428.5～526.8 mm，該地區(qū)與全國(guó)其他地區(qū)均不相同，是典型的寒溫帶大陸性季風(fēng)氣候，全年無(wú)霜期80～110 d，冰封期180～200 d。太陽(yáng)輻射數(shù)以百萬(wàn)計(jì)共4 500耳，年日照2 600 h，年累積溫度2 100℃。

HJ-1A星HSI數(shù)據(jù)根據(jù)其校正程度共分為0～5級(jí)產(chǎn)品，分別為原始數(shù)據(jù)產(chǎn)品、輻射校正產(chǎn)品、系統(tǒng)幾何校正產(chǎn)品、幾何精校正產(chǎn)品、高程校正產(chǎn)品和標(biāo)準(zhǔn)鑲嵌圖像產(chǎn)品[9]。普通用戶通過(guò)訪問(wèn)中國(guó)資源衛(wèi)星應(yīng)用中心的網(wǎng)站就可以免費(fèi)獲得HSI的2級(jí)產(chǎn)品，本文采用HSI的2級(jí)產(chǎn)品。研究選取了包含122°22′E，53°28′N(xiāo)漠河氣象站的黑龍江省大興安嶺地區(qū)的一景圖像，該景圖像的產(chǎn)品號(hào)是1049632，圖像的詳細(xì)信息見(jiàn)表1。

表1 實(shí)驗(yàn)數(shù)據(jù)的詳細(xì)參數(shù)

2 研究方法

本研究的主要流程包含三步：首先對(duì)一景包含漠河氣象站的HSI數(shù)據(jù)進(jìn)行預(yù)處理，包括HSI數(shù)據(jù)的格式轉(zhuǎn)換，去條帶處理以及大氣校正等；第二，在此基礎(chǔ)上應(yīng)用逐日太陽(yáng)輻射模型計(jì)算日總太陽(yáng)輻射，并通過(guò)歸一化植被指數(shù)(NDVI)和比值植被指數(shù)(SR)計(jì)算得到光合有效輻射吸收比率FPAR；最后用FPAR與太陽(yáng)輻射共同計(jì)算出APAR。流程如圖1所示。

圖1 技術(shù)路線圖Fig.1 Technology roadmap

2.1 HSI數(shù)據(jù)預(yù)處理

HSI數(shù)據(jù)的預(yù)處理包括圖像旋轉(zhuǎn)、去條紋、大氣校正和MNF噪聲分離變換等過(guò)程。HJ-1A星 HSI數(shù)據(jù) 2 級(jí)產(chǎn)品存在像元灰度值比較小、與周?chē)忻黠@差異的垂直條帶。條帶噪聲的存在嚴(yán)重影響了圖像的質(zhì)量和應(yīng)用，使有用波段信息不能完全被人們獲得和利用。因此，本文選取了全局去條紋法對(duì)條紋進(jìn)行去除(圖2)[10]。HSI數(shù)據(jù)去條紋處理時(shí)認(rèn)為圖像是按列垂直排列的，所以在進(jìn)行去條紋處理前要對(duì)圖像先進(jìn)行旋轉(zhuǎn)。之后，對(duì)進(jìn)行條帶處理后的HSI數(shù)據(jù)采用ENVI的FLAASH模塊進(jìn)行大氣校正[11]。

圖2 去條紋處理前后圖像對(duì)比Fig.2 Image comparison before and after stripe dealing

對(duì)比大氣校正前后的光譜曲線可知(圖3)，去條紋的噪聲處理和大氣校正后原始圖像在近紅外波段的反射峰值有所提高，突出了植被特征，在藍(lán)光處出現(xiàn)了吸收谷，消除了大氣的多次散射影響。最后，對(duì)處理前后的圖像做MNF噪聲變換，圖像處理后圖像信息主要集中的波段噪聲較少。

2.2 逐日太陽(yáng)輻射的模擬

太陽(yáng)輻射是地球表層上的物理、生物和化學(xué)過(guò)程的主要能源，也是生態(tài)系統(tǒng)過(guò)程模型、水文模擬模型和生物物理模型研究中的必要參數(shù)[12]。本文選取大興安嶺研究區(qū)域中的漠河氣象站，將該氣象站的氣象數(shù)據(jù)和輻射數(shù)據(jù)作為逐日太陽(yáng)輻射模型中緯度、日照時(shí)數(shù)和實(shí)測(cè)日太陽(yáng)輻射的輸入數(shù)據(jù)，對(duì)逐日太陽(yáng)輻射進(jìn)行模擬。

太陽(yáng)輻射是太陽(yáng)直接輻射和散射輻射之和，時(shí)間尺度不同太陽(yáng)總輻射值不同。利用逐日太陽(yáng)輻射模型可以計(jì)算出一天之中的太陽(yáng)總輻射，一天當(dāng)中太陽(yáng)總輻射與當(dāng)天日長(zhǎng)和日照時(shí)數(shù)有關(guān)，逐日太陽(yáng)輻射利用Angtrom-Prescott方程計(jì)算[13]：

H=HL×(a+b×S/SL)。

(1)

式中：H為日實(shí)測(cè)太陽(yáng)總輻射；HL為晴天狀態(tài)下的地面總輻射；S為日照時(shí)數(shù)；SL為日長(zhǎng)；a和b為兩個(gè)經(jīng)驗(yàn)參數(shù)，本文引用童成立等利用實(shí)測(cè)總輻射和日照百分率月平均值等資料計(jì)算得到的適合我國(guó)大部分地區(qū)的經(jīng)驗(yàn)參數(shù)a為0.248，b為0.752[13]。

圖3 圖像處理前后植被光譜曲線對(duì)比Fig.3 Vegetation spectral curve comparison beforeand after image processing

理想狀態(tài)下(晴天狀態(tài))大氣上空太陽(yáng)輻射通過(guò)大氣層到達(dá)地面的輻射量可以表示為：

HL=0.8×HO。

(2)

式中：HL為理想狀態(tài)下(晴天狀態(tài))的地面總輻射；HO為大氣上空的太陽(yáng)總輻射。

大氣上空的太陽(yáng)總輻射HO計(jì)算公式：

HO=(1/π)·Gsc·EO·(cosΦ·cosδ·sinWs+(π/180)·sinΦ·sinδ·Ws。

(3)

式中：Gsc為太陽(yáng)常數(shù)，其值一般為1 367 W/m2(相當(dāng)于118.108 MJ/(m2·d))；EO為地球軌道偏心率校正因子；Φ為緯度；δ為太陽(yáng)赤緯；Ws為時(shí)角(角度)。

2.3 光合有效輻射

到達(dá)地球表面的全部太陽(yáng)輻射只有一部分特定波段的太陽(yáng)輻射能被綠色植物用來(lái)進(jìn)行光合作用，這部分能量被稱作光合有效輻射，其中光合有效輻射量與到達(dá)地表的太陽(yáng)總輻射有關(guān)，利用公式計(jì)算：

PAR=H×ηPAR。

(4)

式中：PAR為光合有效輻射；H為日實(shí)測(cè)太陽(yáng)總輻射；ηPAR為光合有效輻射占太陽(yáng)輻射的比值，不同方法得到的ηPAR有差異，國(guó)內(nèi)學(xué)者估算PAR時(shí)多采用ηPAR為0.49～0.50，本文采用Szeice測(cè)定的 為0.5[14]。

2.4 植物吸收光合有效輻射

植物吸收光合有效輻射與植被類型和植被覆蓋狀況有關(guān)，用植物光合有效輻射吸收率來(lái)表示植物對(duì)光合輻射的吸收效率，研究表明光合有效輻射吸收率與NDVI和植被類型有關(guān)，本文APAR的估算是基于CASA模型的，計(jì)算APAR是以像元為單位逐一計(jì)算的，植物吸收光合有效輻射可以表示為[15]：

APAR=PAR×FPAR。

(5)

式中：APAR為植物吸收光合有效輻射；PAR為光合有效輻射；FPAR為光合有效輻射吸收率，利用植物吸收光合有效輻射計(jì)算公式就可以計(jì)算出每一個(gè)像元的植物吸收光合有效輻射，從而確定整個(gè)研究區(qū)域的總的植物吸收光合有效輻射。

2.5 光合有效輻射吸收比例FPAR的計(jì)算

2.5.1 歸一化植被指數(shù)NDVI計(jì)算FPAR

植被對(duì)太陽(yáng)有效輻射的吸收比例FPAR取決于植被類型和植被覆蓋情況，由遙感數(shù)據(jù)獲得的歸一化植被指數(shù)FPAR能很好地反映植物覆蓋狀況[16]。在一定范圍內(nèi)，F(xiàn)PAR與NDVI之間存在著線性關(guān)系[17]，這一個(gè)關(guān)系可以根據(jù)FPAR的最大值、最小值以及該種植物類型對(duì)應(yīng)的NDVI的最大值、最小值來(lái)確定。由于FPAR的最大值最小值與植被類型無(wú)關(guān)，分別取0.95和0.001[18]，所以具體公式如下：

(6)

2.5.2 比值植被指數(shù)SR計(jì)算FPAR

進(jìn)一步研究表明，F(xiàn)PAR與比值植被指數(shù)SR也存在著較好的線性關(guān)系[19]，可由以下公式表示：

(7)

其中，F(xiàn)PAR的最大值和最小值分別取0.95和0.001，SR(x,t)由公式(8)確定：

(8)

由于本研究區(qū)域?yàn)榇笈d安嶺區(qū)域，該區(qū)域的植被類型主要為落葉針葉林，所以根據(jù)朱文泉[20]對(duì)中國(guó)各陸地植被類型NDVI和SR的最大值最小值列出的分布表(表2)查得本式中NDVImax為0.738，NDVImin為0.023，SRmax為6.63，SRmin為1.05。

表2 中國(guó)陸地植被類型NDVI和SR的最大值與最小值

3 結(jié)果與分析

3.1 逐日太陽(yáng)輻射的結(jié)果驗(yàn)證

選取了大興安嶺地區(qū)的漠河氣象站進(jìn)行比較分析，其輻射觀測(cè)臺(tái)站地理信息見(jiàn)表3。

表3 輻射觀測(cè)臺(tái)站地理信息

根據(jù)逐日太陽(yáng)輻射模型模擬計(jì)算了大氣上界和地面晴空條件下的太陽(yáng)輻射日變化，逐日太陽(yáng)輻射值選取了2013年漠河氣象站全年記錄的日太陽(yáng)實(shí)測(cè)輻射值。結(jié)果表明模擬計(jì)算得到的晴空條件下太陽(yáng)輻射值與實(shí)際觀測(cè)到的逐日太陽(yáng)輻射最大值基本一致(圖4)，少數(shù)奇異點(diǎn)(327.66 MJ/m2)為缺測(cè)數(shù)據(jù)，由于大氣上界的太陽(yáng)輻射穿過(guò)大氣層時(shí)，大氣中的各種成分會(huì)對(duì)太陽(yáng)輻射有吸收、散射作用，最終導(dǎo)致太陽(yáng)輻射的衰減，使得到達(dá)地球表面的太陽(yáng)輻射值小于大氣上界的太陽(yáng)輻射值。

圖4 漠河氣象站地區(qū)大氣上界、地面晴空條件下和實(shí)際觀測(cè)太陽(yáng)輻射日變化Fig.4 Observed diurnal variation of solar radiation under atmospheric boundary and ground clear-sky conditions in Mohe Weather Station area.

利用逐日太陽(yáng)輻射模型計(jì)算得到漠河氣象站一年中每日的太陽(yáng)輻射模擬值與動(dòng)態(tài)監(jiān)測(cè)的太陽(yáng)輻射真實(shí)值，進(jìn)行回歸分析(圖5)，模型模擬結(jié)果與漠河觀測(cè)站點(diǎn)結(jié)果有很高的相關(guān)性(表4)，決定系數(shù)(R2)達(dá)到了0.95，平均偏差MBE為-1.06 MJ/(m2·d)，表明逐日太陽(yáng)輻射模擬值偏小，并且大都小于實(shí)際觀測(cè)的太陽(yáng)輻射值，均方根誤差RMSE反映了模擬值與實(shí)際觀測(cè)值間的偏離程度，為1.94 MJ/(m2·d)，其結(jié)果小于童成立等對(duì)黑龍江漠河的研究結(jié)果[21]，以太陽(yáng)輻射預(yù)測(cè)值為橫坐標(biāo)，殘差為縱坐標(biāo)，得到殘差圖(圖6)，圖中可以看出描繪的點(diǎn)圍繞殘差等于0的直線上下隨機(jī)散布，表明變量實(shí)測(cè)太陽(yáng)輻射與模擬太陽(yáng)輻射有顯著的線性相關(guān)關(guān)系，進(jìn)而表明逐日太陽(yáng)輻射模型可用于模擬漠河氣象站地區(qū)的日太陽(yáng)輻射值，實(shí)現(xiàn)了由漠河氣象站的點(diǎn)輻射值，驗(yàn)證太陽(yáng)輻射模型并應(yīng)用到一景遙感圖像的面輻射值，可用于APAR的計(jì)算，但是日照時(shí)數(shù)和日序數(shù)作為逐日太陽(yáng)輻射模型的輸入?yún)?shù)，在時(shí)間步長(zhǎng)為旬、月和年的情況下，模擬太陽(yáng)輻射值計(jì)算的工作量勢(shì)必增大。

圖5 太陽(yáng)日輻射模擬值與真實(shí)值比較Fig.5 Comparison between simulated and actual values of solar radiation

圖6 模擬太陽(yáng)輻射的殘差圖Fig.6 Simulated residual plots of solar radiation

站點(diǎn)R2MBE(MJ/m-2·d-1)RMSE(MJ/m-2·d-1)黑龍江漠河095282-10568193812

3.2 基于NDVI和SR估算的FPAR結(jié)果

以往研究表明，通過(guò)NDVI估算出來(lái)的FPAR值高于實(shí)測(cè)值，通過(guò)SR估算出來(lái)的FPAR值低于實(shí)測(cè)值，但是SR估算出來(lái)的誤差小于直接由NDVI估算出來(lái)的結(jié)果，所以基于這種情況，Los SO[22]將這兩種方法結(jié)合起來(lái)，取兩者的平均值作為FPAR的估算值，此時(shí)估算的FPAR與實(shí)測(cè)值之間的誤差達(dá)到最小。本研究最終將基于NDVI估算的FPAR和基于SR估算的FPAR的兩個(gè)結(jié)果加和取平均值作為FPAR的最終估算值[23]。

3.3 植物吸收光合有效輻射APAR的估算結(jié)果

通過(guò)ENVI軟件計(jì)算研究區(qū)遙感圖像的APAR，并利用密度分割工具將APAR的大小分割成13等級(jí)，分不同的顏色顯示，運(yùn)用ArcMap制作成一景遙感圖像中APAR分布圖，具體結(jié)果如圖7所示。

從大興安嶺漠河地區(qū)APAR分布圖可以看出，APAR值為5～6 MJ/(m2·d)和6～7 MJ/(m2·d)的分布最多，占據(jù)了整個(gè)區(qū)域的65%以上，通過(guò)與HSI數(shù)據(jù)的標(biāo)準(zhǔn)假彩色影像圖進(jìn)行對(duì)比，APAR值為5～6 MJ/(m2·d)所分布的地方為大片平地上的森林，APAR值為4～5 MJ/(m2·d)所分布的地方為坡地上的森林，而APAR值為3～4 MJ/(m2·d)所分布 的地方為平地上的草地和低矮的植被，裸地上沒(méi)有APAR值。由于本研究的是大興安嶺漠河地區(qū)一天的APAR值，用一天的APAR值估計(jì)年APAR值約為1 800 MJ/(m2·a)，在對(duì)我國(guó)植被凈第一生產(chǎn)力的估算中[24-26]，東北長(zhǎng)白山、小興安嶺、大興安嶺等林區(qū)的APAR值一般在1 300～2 200 MJ/(m2·a)，所以本研究基本符合要求。

圖7 大興安嶺漠河地區(qū)一景遙感圖像的APAR分布圖Fig.7 APAR distribution of remote sensing images of a scene in Daxinganling Mohe area

4 結(jié)論

植被吸收光合有效輻射是計(jì)算森林生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力的一個(gè)非常重要的參數(shù)，本研究基于HJ-1A星HSI數(shù)據(jù)的APAR研究方法，以一景遙感圖像作為數(shù)據(jù)源，運(yùn)用ENVI對(duì)HSI數(shù)據(jù)進(jìn)行預(yù)處理，根據(jù)漠河氣象站的氣象和輻射數(shù)據(jù)，應(yīng)用逐日太陽(yáng)模型模擬出太陽(yáng)總輻射，利用NDVI和SR估算FPAR，最后基于CASA模型以像元為單位逐一計(jì)算出APAR。主要研究結(jié)果如下：

(1)本文選取了全局去條紋法對(duì)HSI數(shù)據(jù)旋轉(zhuǎn)后垂直分布的條紋進(jìn)行去除，可以在保證信息質(zhì)量的前提下，最大程度的去除條紋的影響，提高圖像質(zhì)量。

(2)本文對(duì)逐日太陽(yáng)輻射模擬結(jié)果進(jìn)行統(tǒng)計(jì)分析，得出模擬值與實(shí)測(cè)值相比普遍偏低，但偏差不大，RMSE為1.94 MJ/m2·d，以日照時(shí)數(shù)和日序數(shù)作為模型輸入?yún)?shù)時(shí)，如果在時(shí)間步長(zhǎng)為旬、月、年時(shí)，模擬太陽(yáng)輻射值的工作量增大，模型方法將在之后研究中進(jìn)一步改進(jìn)。與已有的太陽(yáng)輻射模型相比，該模型具有易操作的特點(diǎn)，所需的輸入資料簡(jiǎn)單(緯度，日照時(shí)數(shù)，日序數(shù))，實(shí)測(cè)太陽(yáng)輻射與模擬太陽(yáng)輻射有顯著的線性相關(guān)關(guān)系，模擬值與實(shí)測(cè)值之間達(dá)到了顯著相關(guān)水平，決定系數(shù)(R2)為0.95，模擬漠河氣象站地區(qū)的每日太陽(yáng)輻射精度高，實(shí)現(xiàn)了由漠河氣象站的點(diǎn)輻射值，到一景遙感圖像的面輻射值，可作為APAR計(jì)算的重要參數(shù)。

對(duì)CASA模型的算法進(jìn)行了改進(jìn)和簡(jiǎn)化，改進(jìn)了其中的FPAR算法，由于由NDVI計(jì)算出的FPAR值偏高，由SR估計(jì)的FPAR偏小，其中SR由NDVI計(jì)算得到，所以將NDVI和SR的影響結(jié)合起來(lái)考慮，使得最后得到的FPAR結(jié)果更接近于真實(shí)值。

本研究獲得的NDVI數(shù)據(jù)在無(wú)植被地區(qū)有時(shí)會(huì)出現(xiàn)較大的值，因此這需要將遙感數(shù)據(jù)和地面植被數(shù)據(jù)在計(jì)算前后進(jìn)行比對(duì)，盡量消除這種由于數(shù)據(jù)本身原因而導(dǎo)致結(jié)果誤差的特殊情況。

[1]王保林，王晶杰，楊勇，等.植被光合有效輻射吸收分量及最大光能利用率算法的改進(jìn)[J].草業(yè)學(xué)報(bào)，2013，22(5)：220-228.

[2]何明珠.遙感數(shù)據(jù)與兩葉光能利用率模型同化的總初級(jí)生產(chǎn)力模擬[D].南京：南京大學(xué)，2013.

[3]Hipps L E，Asrar G，Kanemasu E T.Assessing the interception of photosynthetically active radiation in winter wheat[J].Agricultural Meteorology，1983，28：253-259.

[4]Meere K，McCree J.A solarineter for measuring photosynthetically active radiation[J].Agric Meteorol，1966，3(5-6)：353-366.

[5]Chen J M，Liu J，Cihlar J，et al.Daily canopy photosynthesis model through temporal and spatial scalingforremote sensing applications[J].Ecological Modelling，1999，124(2-3)：99-119.

[6]張喜旺，于寧，秦奮，等.基于遙感的植被吸收光合有效輻射估算[J].生態(tài)科學(xué)，2013，32(5)：604-608.

[7]李剛，辛?xí)云剑醯例?，?基于MODIS數(shù)據(jù)的草地植被光合過(guò)程參數(shù)估算——以內(nèi)蒙古自治區(qū)為例[J].中國(guó)草地學(xué)報(bào)，2008，33(3)：1-7.

[8]陳新芳，安樹(shù)青，陳鏡明，等.森林生態(tài)系統(tǒng)生物物理參數(shù)遙感反演研究進(jìn)展[J].生態(tài)學(xué)雜志，2005，24(9)：1074-1079.

[9]鈕立明，蒙繼華，吳炳方，等.HJ-1A星HSI數(shù)據(jù)2級(jí)產(chǎn)品處理流程研究[J].國(guó)土資源遙感，2011(1)：77-82.

[10]高海亮，顧行發(fā)，余濤，等.Hyperion遙感影像噪聲去除方法研究[J].遙感信息，2014，29(3)：3-7.

[11]Gao B，Montes M J，Davis C O.Atmospheric correction algorithms for hyperspectral remote sensing data of land and ocean[J].Remote Sensing of Environment，2009，113(1)：17-24.

[12]齊月，房世波，周文佐.近50年來(lái)中國(guó)地面太陽(yáng)輻射變化及其空間分布[J].生態(tài)學(xué)報(bào)，2014，34(24)：7444-7453.

[13]左大康，王懿賢，陳建綏.中國(guó)地區(qū)太陽(yáng)總輻射的空間分布特征[J].氣象學(xué)報(bào)，1963(1)：78-96.

[14]Szeicz G.Solar radiation for plant growth[J].Journal of Applied Ecology，1974，11(11)：617-636.

[15]董丹，倪健.利用CASA模型模擬西南喀斯特植被凈第一性生產(chǎn)力[J].生態(tài)學(xué)報(bào)，2011，31(7)：1855-1866.

[16]董泰鋒，蒙繼華，吳炳方.基于遙感的光合有效輻射吸收比率(FPAR)估算方法綜述[J].生態(tài)學(xué)報(bào)，2012，32(22)：7190-7201.

[17]Ruimy A，Saugier B.Methodology for the estimation of terrestrial net primary production from remotely sensed data[J].Journal of Geophysical Research，1994，97：18515-18521.

[18]高艷妮，于貴瑞，張黎，等.中國(guó)陸地生態(tài)系統(tǒng)凈初級(jí)生產(chǎn)力變化特征——基于過(guò)程模型和遙感模型的評(píng)估結(jié)果[J].地理科學(xué)進(jìn)展，2012，31(1)：109-117.

[19]Los S O，Justice C O，Tucker C J.A global 1by 1NDVI dataset for climate studies derived from the GIMMS continental NDVI data[J].International Journal of Remote Sensing，1994，15：3493-3518.

[20]朱文泉，潘耀忠，張錦水.中國(guó)陸地植被凈初級(jí)生產(chǎn)力遙感估算[J].植物生態(tài)學(xué)報(bào)，2007，31(3)：413-424.

[21]童成立，張文菊，湯陽(yáng)，等.逐日太陽(yáng)輻射的模擬計(jì)算[J].中國(guó)農(nóng)業(yè)氣象，2005，26(3)：165-169.

[22]Los S O.linkages between global vegetation and climate：an analysis based on NOAA advanced very high resolution radiometer data[D].Oregon：National Oceanic and Atmospheric Administration，1998.

[23]Peng D，Zhang B，Liu L.Comparing spatiotemporal patterns in Eurasian FPAR derived from two NDVI-based methods[J].International Journal of Digital Earth，2012，5(4)：1-16.

[24]陳正華，麻清源，王建，等.利用CASA模型估算黑河流域凈第一性生產(chǎn)力[J].自然資源學(xué)報(bào)，2008，23(2)：263-273.

[25]呂丹桔，施心陵，吳云鵬.基于Labview的植物紅光吸收動(dòng)態(tài)監(jiān)測(cè)[J].林業(yè)機(jī)械與木工設(shè)備，2015，43(1)：35-36+42.

[26]孫存舉，雷建華，唐果山，等.基于SPOT NDVI影像的阿壩州植被變化遙感監(jiān)測(cè)[J].四川林業(yè)科技，2014，35(6)：56-59.

Vegetation Absorbed Photosynthetically Active Radiation Estimates Based on HJ-1A Satellite HSI Data

Chang Yaxuan1，Yu Ying*，Cui Lvyuan2，Hou Jiyu2，Zeng Dewei2

(College of Forest，Northeast Forestry University，Harbin 150040)

In order to measure the effective absorption of solar radiation by vegetation，HSI data obtained by HJ-1A satellite was used to estimate the absorbed photosynthetic active radiation.The basis was provided for the separate calculation of sunlit and shadedAPAR.The accurate parameter data was also provided for the net primary productivity calculation of forest ecosystem.This study includes three main steps：Firstly，the HSI data containing Mohe weather station were pre-processed，which includes format conversion，stripe-removing and atmospheric correction and so on.Secondly，daily solar radiation model was used to calculate the daily total solar radiation，and then theFPARwas calculated byNDVI.FPARand solar radiation were used to obtain theAPAR.Finally the processed image was obtained and the correlation between daily solar radiation model and regional meteorological station data was verified.The precision of FPAR calculated by combination ofNDVIandSRwas improved.The grayscale image was displayed as a thematic image when theAPARcalculation of vegetation was completed.TheAPARvalues of different regions was obtained from the thematic maps，and the annualAPARwas finally calculated.

HSI data processing；solar radiation model；FPAR；APAR

2016-08-30

大學(xué)生創(chuàng)新項(xiàng)目(20150225068)

常雅軒，本科生。研究方向：遙感與地理信息系統(tǒng)。

*通信作者：于穎，博士，講師。研究方向：定量遙感與地理信息系統(tǒng)。E-mail：yuying4458@163.com

常雅軒，于穎，崔綠園，等.基于HJ-1A星HSI數(shù)據(jù)的植被吸收光合有效輻射APAR估算[J].森林工程，2017，33(2)：22-27.

S 771

A

1001-005X(2017)02-0022-06