王明常，牛雪峰，陳圣波，王亞楠，汪自軍

1.吉林大學(xué)地球探測(cè)科學(xué)與技術(shù)學(xué)院，長(zhǎng)春 130026

2.千葉工業(yè)大學(xué)建筑與城市環(huán)境學(xué)院，日本 千葉 275－0016

3.長(zhǎng)春市測(cè)繪院，長(zhǎng)春 130021

4.上海工程衛(wèi)星研究所，上海 200240

0 引言

葉面積指數(shù)（leaf area index，LAI）是一個(gè)關(guān)鍵的植被特征參量。LAI是植被冠層結(jié)構(gòu)的重要量化指標(biāo)，反映植物葉面數(shù)量、冠層結(jié)構(gòu)、植物群落及其環(huán)境效應(yīng)，為植物冠層表面的物質(zhì)和能量交換提供定量信息，與蒸騰作用、光合作用、呼吸作物等密切相關(guān)。隨著遙感成為生態(tài)學(xué)研究的一種重要工具，人們?cè)絹?lái)越重視通過(guò)衛(wèi)星遙感來(lái)定量估算LAI。

多角度遙感數(shù)據(jù)通過(guò)更多的觀測(cè)角度來(lái)獲得更高幾何和物理精度，為植被的研究提供了可靠信息。針對(duì)如何結(jié)合二向性反射模型，利用多角度信息進(jìn)行植被參數(shù)的反演并提高反演精度，國(guó)內(nèi)外學(xué)者做了很多不同的研究。Miina等［1］根據(jù)不同類(lèi)型森林的各向異性輻射傳輸模型特征，利用CHRIS（compact high resolution imaging spectrometer）數(shù)據(jù)反演了冠層結(jié)構(gòu)，建立了森林和地面光譜分離的方法；Verrelst等［2］利用CHRIS數(shù)據(jù)進(jìn)行植被指數(shù)分析，并分析其變化特征；Kamalesh等［3］結(jié)合冠層反射模型，反演了葉面積指數(shù)；Quaife等［4］利用CHRIS數(shù)據(jù)在進(jìn)行了大氣校正和幾何校正后反演了地表BRDF（bidirectional reflectance distribution function）；Barnsley等［5］根據(jù)冠層結(jié)構(gòu)的光學(xué)性質(zhì)，對(duì)不同的土地覆蓋和表面類(lèi)型進(jìn)行分類(lèi)；Vuolo等［6］根據(jù)冠層反射模型，利用CHRIS數(shù)據(jù)反演了葉面積指數(shù)；Zbyněk等［7］利用DART（discrete anisotropic radiative transfer）模型模擬了植被近紅外的冠層輻射特征；蓋利亞等［8］、申茜等［9］、董廣香等［10－11］、柳彩霞等［12］對(duì) CHRIS數(shù)據(jù)的大氣校正、幾何校正等預(yù)處理方法進(jìn)行了研究，對(duì)水體葉綠素濃度進(jìn)行了反演，獲得了理想的研究結(jié)果。在利用CHRIS數(shù)據(jù)進(jìn)行植被結(jié)構(gòu)參數(shù)反演方面，國(guó)內(nèi)相應(yīng)的研究很少。目前利用遙感估算葉面積指數(shù)主要有2種方法：經(jīng)驗(yàn)統(tǒng)計(jì)方法和物理模型方法。筆者利用各向異性輻射傳輸（DART）模型建立不同覆蓋情況下的針葉林反射率查找表，通過(guò)對(duì)CHRIS數(shù)據(jù)進(jìn)行條帶去除、大氣校正、幾何校正和地形校正，根據(jù)建立的查找表匹配反演葉面積指數(shù)，并對(duì)反演的敏感性進(jìn)行分析。

1 DART模型原理與查找表的建立

DART模型［13］全面考慮了地物目標(biāo)空間分布不同引起的光譜輻射各向異性問(wèn)題。在可見(jiàn)光和近紅外波段，可以同時(shí)設(shè)置若干個(gè)波段，模擬地球任意位置、任意太陽(yáng)方位角、任意觀測(cè)角度傳感器接收的輻射傳輸。DART有2種模式：R模式，只模擬反射率，太陽(yáng)或大氣作為唯一輻射源；T模式，模擬亮度溫度，輻射源分為地球和大氣兩者結(jié)合或地球、大氣、太陽(yáng)三者結(jié)合2種。DART模型在場(chǎng)景劃分時(shí)水平方向與垂直方向是相互獨(dú)立的，并且垂直方向上的變化大于水平方向。場(chǎng)景內(nèi)所有成分被劃分成由長(zhǎng)方體單元所構(gòu)成的單元集合，根據(jù)每個(gè)單元的中心在三維場(chǎng)景中用坐標(biāo)來(lái)定位。單元格總數(shù)為

其中：ΔX、ΔY、ΔZ是研究場(chǎng)景的長(zhǎng)、寬、高；Δx、Δy、Δz為每個(gè)單元的長(zhǎng)、寬、高。每個(gè)單元被用來(lái)模擬不同類(lèi)型的場(chǎng)景元素，如土壤、樹(shù)干、樹(shù)枝、樹(shù)葉、建筑物、水體等。

1.1 輻射方向上的能量傳輸

利用Hapke傳輸方程的基本思想，在方向Ω上、位置為r處的反光強(qiáng)度I（r，Ω）為

其中：μ，η，ξ分別表示r與z，y，x方向軸之間夾角的余弦；α（r，Ω）是消光系數(shù)；αd（r，Ω′→Ω）則表示從Ω′方向到Ω方向的散射系數(shù)的數(shù)學(xué)微分值；Ω和Ω′分別是光子的入射方向和散射方向。

為區(qū)分初次碰撞與多次碰撞，將方向Ω離散為N個(gè)量Ωn。在三維笛卡爾坐標(biāo)系中，沿Ωij方向的離散角度傳輸方程為

其中：u和v分別是天頂角μ和方位角φ的離散度，U和V是離散度的最大值；Ωuv代表某一特定方向；Cuv是與方向Ωuv有關(guān)的所有方向上的散射權(quán)重；Ω（r，Ωij）為第一次碰撞源量；I（r，Ωuv）是方向Ωuv上的光強(qiáng)度

光強(qiáng)度I（r，Ωuv）是由微分散射系數(shù)αd（r，Ω′→Ω）決定的，它代表平均光強(qiáng)度，描述的是從ΔΩuv至ΔΩij的散射過(guò)程。把離散輻射傳輸方程簡(jiǎn)化為

其中：W（r，Ωij）代表沿Ωij方向上任一圓錐角ΔΩij在位置r處的太陽(yáng)輻射能量通量，W（r，Ωuv）代表沿Ωuv方向上任一圓錐角ΔΩuv在位置r處的太陽(yáng)輻射能量通量。W（r，Ωuv）和W（r，Ωij）分 別 與I（r，Ωuv）·ΔΩuv和I（r，Ωij）·ΔΩij成比例。

1.2 DART模擬樹(shù)的空間分布特征

DART模型用3種方式模擬樹(shù)的分布：一是隨機(jī)的位置和尺寸，即對(duì)于每一樹(shù)種、樹(shù)干和樹(shù)冠的尺寸是滿(mǎn)足高斯分布的，通過(guò)在場(chǎng)景格網(wǎng)內(nèi)節(jié)點(diǎn)處的分布半徑，根據(jù)樹(shù)位置存在的概率來(lái)模擬樹(shù)的分布；二是確切的位置和隨機(jī)的尺寸，設(shè)定每棵樹(shù)的類(lèi)型和位置，但樹(shù)干和樹(shù)冠的尺寸滿(mǎn)足高斯隨機(jī)分布；三是確切的位置和尺寸，每棵樹(shù)的位置和尺寸大小都是給定的，但需要場(chǎng)景單元的大小小于樹(shù)冠的尺寸，否則不會(huì)生成樹(shù)。

1.3 DART模型中樹(shù)的描述

每棵樹(shù)都是由樹(shù)干和樹(shù)冠兩大部分組成。樹(shù)冠是由枝干、樹(shù)枝、樹(shù)葉以及空隙組成。樹(shù)干是由多個(gè)垂直的矩形和一個(gè)水平的正方形六面體重疊而成。

樹(shù)的模擬由如下因子性質(zhì)決定。

1）樹(shù)冠形狀，可分為橢球形、半橢球形、截錐形和不規(guī)則四邊形等，每種樹(shù)冠的尺寸對(duì)BRDF均有影響。

2）樹(shù)種j的葉片數(shù)量Bio［j］：

Bio［j］＞0時(shí)，LAI［j］為

式中：Dx和Dy為坐標(biāo)系中沿x軸和y軸的長(zhǎng)度：

l＝l（z），u［j，l］是l層j樹(shù)種的體密度，u［j，l，x，y］是葉片水平體密度，uf［j］＝

則DART模型三維場(chǎng)景單元格內(nèi)葉面積指數(shù)LAIcell［j，x，y，z］為

式中：cell（j，x，y，z）為單元格位置。

當(dāng)Bio［j］＜0時(shí)，LAI［j］為

場(chǎng)景模擬中可設(shè)定是否模擬樹(shù)枝。所模擬樹(shù)干反射率ρtrunk（i，l）和樹(shù)枝反射率ρbranch（i，l）的類(lèi)型可以是朗伯體、朗伯體加鏡面反射或Hapke模型中的一種。冠下樹(shù)干和冠內(nèi)樹(shù)干的幾何參數(shù)包括冠下樹(shù)干的高和直徑，冠內(nèi)樹(shù)干的高和直徑。

1.4 利用DART模型建立查找表

據(jù)查找表的原理與方法，模擬了不同參數(shù)下的植被分布情況。根據(jù)需要利用CHRIS數(shù)據(jù)圖像與查找表進(jìn)行匹配?；贑HRIS數(shù)據(jù)的波段模式，將通過(guò)DART模型建立的查找表設(shè)置為18個(gè)波段。DART模型對(duì)場(chǎng)景共離散了61個(gè)方向［14］，是天頂角在0°～90°，方位角在0°～360°范圍內(nèi)角度的不同組合。設(shè)定模型模擬的枝干和土壤為反射類(lèi)型朗伯體，葉片為針葉林，輸入各組成波譜特性數(shù)據(jù)。CHRIS數(shù)據(jù)的空間分辨率為17m，設(shè)定基本模擬場(chǎng)景大小為17m×17m×17m。考慮了冠高、冠底半徑、冠下桿高、冠下胸徑、葉面積指數(shù)等植被結(jié)構(gòu)參數(shù)作為輸入。通過(guò)DART模型得到植被BRF（bidirectional reflectance factor），將結(jié)構(gòu)參數(shù)與模擬所得的BRF依次對(duì)應(yīng)成表，得到針葉林植被結(jié)構(gòu)參數(shù)的查找表（表1）。

表1 查找表結(jié)構(gòu)Table 1 Look－up－table（LUT）structures

2 研究區(qū)選擇與遙感數(shù)據(jù)處理

2.1 研究區(qū)選擇

研究區(qū)位于吉林省東南部，東經(jīng)127°40′－128°16′、北緯41°35′－42°25′。選擇針葉林區(qū)域?yàn)檠芯繉?duì)象，進(jìn)行針葉林葉面積指數(shù)反演研究。長(zhǎng)白山地區(qū)生物物種資源豐富、生態(tài)系統(tǒng)完整，是歐亞大陸北半部最具有代表性的自然綜合體，不僅是東北地區(qū)生態(tài)系統(tǒng)的重要屏障，也對(duì)東北地區(qū)乃至東北亞地區(qū)的物質(zhì)能量交換、氣候變化起著重要的作用。

2.2 遙感數(shù)據(jù)選擇

歐洲空間局2001年10月發(fā)射的PROBA（project for on board autonomy）實(shí)驗(yàn)衛(wèi)星上搭載的CHRIS多角度傳感器，與太陽(yáng)同步軌道衛(wèi)星，CCD推掃式成像，空間分辨率17m，光譜分辨率10nm，成像像面積13km×13km，軌道高度為615km，共5個(gè)掃描角度（0°、±36°、±55°），為大氣、陸地和海洋的二向性反射研究提供科學(xué)數(shù)據(jù)。這些優(yōu)點(diǎn)有利于生物量評(píng)估以及生物健康狀況的監(jiān)測(cè)，同時(shí)對(duì)植被冠層結(jié)構(gòu)、密度、植被識(shí)別或林木種類(lèi)方面的研究起著重要作用。

2.3 CHRIS數(shù)據(jù)預(yù)處理

植被遙感的目的是通過(guò)傳感器觀測(cè)的表觀光學(xué)參數(shù)來(lái)反演地表植被參數(shù)，遙感數(shù)據(jù)預(yù)處理是植被定量遙感的前提。預(yù)處理主要針對(duì)CHRIS特征進(jìn)行條帶噪聲去除、大氣校正［15］、地形校正。

CHRIS數(shù)據(jù)主要有2種噪聲：隨機(jī)噪聲和條帶噪聲。CHRIS數(shù)據(jù)主要在推掃式傳感器CCD獲得圖像的過(guò)程中形成，條帶噪聲主要來(lái)源于傳感器CCD和衛(wèi)星平臺(tái)的在軌運(yùn)動(dòng)。CHRIS數(shù)據(jù)植被模式條帶噪聲的去除，為大氣校正和定量反演提供基礎(chǔ)。根據(jù)大氣輻射傳輸原理，去除大氣對(duì)輻射傳輸?shù)挠绊?，使衛(wèi)星數(shù)據(jù)能夠真實(shí)地反映地表植被的反射特征；利用MORTAN模型原理，對(duì)CHRIS去條帶數(shù)據(jù)進(jìn)行了大氣校正［15］。

由于地表的自然起伏，不同區(qū)域接收到的有效光照有很大差異。在遙感圖像上，陰坡上的像元接收到較弱的照度而具有較小的亮度值，陽(yáng)坡上的像元卻接收到較強(qiáng)的照度而具有較大的亮度值。這種光譜信息的變化，嚴(yán)重影響了山區(qū)遙感數(shù)據(jù)的反演精度。研究區(qū)地形起伏較大，為了提高反演精度，通過(guò)地形校正來(lái)消除地形變化對(duì)光譜的影響，恢復(fù)地物在水平條件下的真實(shí)反射率或輻亮度，這是定量遙感反演的前提。地形校正采用余弦校正模型［16］：

其中：LH和LT分別為某像元在水平地面和傾斜坡面上的輻亮度；β為太陽(yáng)天頂角；α為太陽(yáng)入射角。利用余弦校正模型對(duì)CHRIS數(shù)據(jù)進(jìn)行了輻射校正，圖1為校正前后的反射率對(duì)比圖，根據(jù)地形校正原理，校正后的數(shù)據(jù)更能真實(shí)地反映輻射傳輸?shù)奶攸c(diǎn)。

圖1 地形校正前后CHRIS的反射率值Fig.1 Reflectivity values corresponding to CHRIS data before and after terrain correction

3 葉面積指數(shù)反演

利用葉面積指數(shù)查找表，將預(yù)處理后的CHRIS數(shù)據(jù)與之匹配?；谧钚《朔ㄔ?，將CHRIS數(shù)據(jù)與查找表進(jìn)行匹配，即

其中：Rch代表CHRIS數(shù)據(jù)18個(gè)波段的BRF值；RLUT為經(jīng)過(guò)角度插值篩選后的與CHRIS數(shù)據(jù)角度相符的18個(gè)波段的BRF值；Nch為CHRIS數(shù)據(jù)的角度數(shù)。求出R（r＊）后，找出所在的查找表文件索引，便能找到該文件對(duì)應(yīng)的植被參數(shù)值。根據(jù)匹配算法，選擇CHRIS圖像進(jìn)行植被結(jié)構(gòu)參數(shù)反演。由于研究區(qū)是針葉林區(qū)，葉面積指數(shù)偏低，LAI值主要為1～4，反演結(jié)果（圖2）能夠充分地說(shuō)明反演模型的準(zhǔn)確。LAI分布和變化呈一定規(guī)律性，高大樹(shù)木LAI值較矮小樹(shù)木大；某一參數(shù)值變大其整體參數(shù)值也是變大的，滿(mǎn)足植被的現(xiàn)實(shí)生長(zhǎng)狀況。

4 葉面積指數(shù)對(duì)BRDF的敏感性分析

選擇490nm、700nm 2個(gè)波段葉面積指數(shù)對(duì)植被BRDF值的敏感性進(jìn)行分析。利用DART模型模擬單棵樹(shù)葉面積指數(shù)對(duì)植被BRDF值的敏感性分析時(shí)，模擬LAI值設(shè)定為 2、3、4、5、6，植被BRDF值的變化情況如圖3所示。相同背景條件下，在近紅外波段和可見(jiàn)光波段，植被BRDF值對(duì)葉面積指數(shù)變化的敏感程度是不同的。在490nm時(shí)，隨著LAI值的增大，BRDF值的變化先增大后減?。辉?00nm時(shí)，BRDF值隨LAI值增大而減小。LAI值為4～6時(shí)，BRDF值變化較小。沿著天頂角增大的方向上BRDF值也是呈下降趨勢(shì)的，在接近天頂附近時(shí)值為最小。由此可見(jiàn)，葉面積指數(shù)對(duì)植被BRDF值存在一定影響，不是十分敏感。

5 結(jié)論與討論

利用多角度高光譜CHRIS數(shù)據(jù)反演葉面積指數(shù)，能夠真實(shí)、快速地獲取研究區(qū)葉面積指數(shù)，從物理機(jī)制反演森林LAI，取得較好的效果。

葉面積指數(shù)對(duì)植被BRDF值存在一定影響。植被BRDF隨著LAI值的增大，在490nm波長(zhǎng)時(shí)BRDF值先增大后減小，在700nm波長(zhǎng)時(shí)BRDF值減小；沿著天頂角增大的方向上BRDF值下降，在天頂附近時(shí)值為最小。

DART模型建立的查找表也存在一定的偏差：1）在場(chǎng)景設(shè)置模擬中由于尺度的選擇會(huì)造成一定的影響；2）在模擬過(guò)程中選擇的植被理化參數(shù)都來(lái)自于國(guó)外數(shù)據(jù)庫(kù)，存在一定的偏差。因此，需要通過(guò)實(shí)驗(yàn)室和野外采集參數(shù)進(jìn)行迭代的方法，顯示出多角度高光譜CHRIS數(shù)據(jù)反演葉面積具有的優(yōu)勢(shì)。

圖2 LAI反演結(jié)果Fig.2 Retrieval results of LAI

圖3 不同觀測(cè)方向葉面積指數(shù)對(duì)BRDF的影響Fig.3 Sensitive influences of LAI to BRDFvalues in the different discrete direction

（References）：

［1］Miina R，Mait L，Matti M，et al.Multi－Angular Reflectance Properties of a Hemiboreal Forest：An Analysis Using CHRIS PROBA Data［J］.Remote Sensing of Environment，2008，112：2627－2642.

［2］Verrelst J，Koetz B，Kneubühler M，et al.Directional Sensitivity Analysis of Vegetation Indices from Multiangular CHRIS/PROBA Data［C］//Proc.ISPRS Commission VII Mid－Term Symposium “Remote Sensing：from Pixels to Processes”.Enschede：［s.n.］，2006：677－683.

［3］Kamalesh V L，Luckman A，Grey W M F，et al.Estimation of Leaf Area Index from PROBA/CHRIS Hyperspectral Multi－Angular Data［C］//Measuring Change in the Earth System.［S.l.］：University of Exeter，2008：1－4.

［4］Quaife T，Barnsley M，Thackrah G，et al.Mission Status of CHRIS/PROBA and Its Utility for Sampling the EarthSurface BRDF［C］//RSPS Proceedings.Londen：［s.n.］，2001：128－134.

［5］Barnsley M J，Settle J J，Cutter M A，et al.The PROBA/CHRIS Mission：A Low－Cost Smallsatfor Hyperspectral Multi－Angle Observations of the Earth Surface and Atmosphere［J］.IEEE Trans on Geoscience and Remote Sensing ，2004，42（7）：1512－1520.

［6］Vuolo F，Dini L，Durso G.Retrieval of Leaf Area Index from CHRIS/PROBA Data：An Analysis of the Directional and Spectral Information Content［J］.International Journal of Remote Sensing，2008，17（29）：5063－5072.

［7］Zbyněk M，Emmanuel M，Lucie H，et al.Influence of Woody Elements of a Norway Spruce Canopy on Nadir Reflectance Simulated by the DART Model at Very High Spatial Resolution［J］.Remote Sensing of Environment，2008，112：1－18.

［8］蓋利亞，劉正軍，張繼賢.CHRIS/PROBA高光譜數(shù)據(jù)的預(yù)處理［J］.測(cè)繪工程，2008，17（1）：40－43.Gai Liya，Liu Zhengjun，Zhang Jixian.Preprocessing of CHRIS/PROBA Hyperspectral Data［J］.Engineering of Surveying and Mapping，2008，17（1）：40－43.

［9］申茜，張兵，李俊生，等.航天高光譜遙感器CHRIS的水體圖像大氣校正［J］.測(cè)繪學(xué)報(bào)，2008，37（4）：476－481.Shen Qian， Zhang Bing， Li Junsheng，et al.Atmospheric Correction for Water Body Images Acquired by Spaceborne Hyperspectral Sensor CHRIS［J］.Acta Geodaetica et Cartographica Sinica，2008，37（4）：476－481.

［10］董廣香，劉正軍，張繼賢，等.基于CHRIS/PROBA數(shù)據(jù)和地面實(shí)測(cè)光譜的三峽庫(kù)區(qū)水體葉綠素濃度反演研究［J］.遙感信息，2008（1）：60－65.Dong Guangxiang，Liu Zhengjun，Zhang Jixian，et al.The Preliminary Result of Chlorophyll： A Retrieval in Three Gorges Dam Using CHRIS/PROBA Mode 2Data Combining with in Situ Measured Data［J］.Remote Sensing Information，2008（1）：60－65.

［11］董廣香，張繼賢，劉正軍.CHRIS/PROBA數(shù)據(jù)條帶噪聲去除方法比較［J］.遙感信息，2006（6）：36－39.Dong Guangxiang，Zhang Jixian，Liu Zhengjun.A Comparison of Several Destriping Methods for CHRIS/PROBA Data［J］.Remote Sensing Information，2006（6）：36－39.

［12］柳彩霞，傅南翔，郭子祺，等.基于人工神經(jīng)網(wǎng)絡(luò)的CHRIS數(shù)據(jù)內(nèi)陸水體葉綠素濃度反演研究［J］.安徽農(nóng)業(yè)科學(xué)，2009，37（26）：12654－12656.Liu Caixia，F(xiàn)u Nanxiang，Guo Ziqi.et al.Inversion Study on Chlorophyll Concentration of Inland Water from CHRIS Data Based on Artificial Neural Networks［J］.Journal of Anhui Agricultural Sciences，2009，37（26）：12654－12656.

［13］Gastellu－Etchegorry J P，Demarez V，Pinel V，et al.Modeling Radiative Transfer in Heterogeneous 3－D Vegetation Canopies［J］.Remote Sensing of Environment，1996，58：131－156.

［14］王亞楠.長(zhǎng)白山地區(qū)植被結(jié)構(gòu)參數(shù)多角度遙感反演［D］.長(zhǎng)春：吉林大學(xué)，2010.Wang Yanan.Inversion of Vegetation Structure Parameters from Multi－Angle CHRIS Data in Changbai Mt［D］.Changchun：Jilin University，2010.

［15］王明常，王亞楠，陳圣波，等.多角度高光譜CHRIS/Proba植被模式數(shù)據(jù)大氣校正［J］.吉林大學(xué)學(xué)報(bào)：地球科學(xué)版，2011，41（2）：609－614.Wang Mingchang，Wang Ya’nan，Chen Shengbo，et al.Atmospheric Correction Research on Vegetation Patterns of Multi－Angle Hyperspectral CHRIS/PROBA Data［J］.Journal of Jilin University：Earth Science Edition，2011，41（2）：609－614.

［16］段四波，閻廣建.山區(qū)遙感圖像地形校正模型研究綜述［J］.北京師范大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，43（3）：362－366.Duan Sibo，Yan Guangjian.A Review of Models for Topographic Correction of Remotely Sensed Images in Mountainous Area［J］.Journal of Beijing Normal University：Natural Science，2007，43（3）：362－366.