基于多次散射的植被_土壤二向反射模型
第8卷第3期2004年5月
遥 感 学 报
Vol. 8, No. 3May, 2004
-4619(2004) 03-0193-08 文章编号:1007
基于多次散射的植被-土壤二向反射模型
邓孺孺
1, 2
, 田国良, 柳钦火
11
(11中国科学院遥感应用研究所, 北京 100101; 21中山大学地理学院, 广州 510275)
摘 要: 运算过于复杂是目前多数植被二向反射模型存在的一个突出问题。对冠层的各次散射、透射和地面土壤反射分别进行分析和考虑, 建立各次反射光谱分量与叶面积指数的函数关系。在Verstraete (1986) 的物理模型的基础上加上多次散射和土壤反射的分量, 建立了基于多次散射的植被-土壤二向反射模型, 并证明仅考虑前三次散射即可满足一般精度要求。该模型克服了原模型的不足, 并可根据叶面积指数等参量直接计算二向反射率, 避免了运算量巨大的叠代运算, 故便于应用。实验证明, 该模型具有较高的精度。关键词: 二向反射模型; 植物冠层模型; 多次散射; 叶面积指数; 土壤含水量中图分类号: TP701 文献标识码: A
迭代运算求整个冠层的方向反射率。较之Suits 模
1 引 言
植被二向反射模型的建立对定量遥感有重要意义, 是进行土壤含水量调查、作物长势调查、作物估产及地表反照率反演等方面定量遥感必须解决的问题。在过去的30年中, 国际上针对植物的光谱特征
进行了大量的研究, 建立了为数众多的植被冠层模型。目前占主流地位的是辐射传输模型和几何光学模型以及近几年发展起来的核驱动模型。辐射传输模型以辐射传输方程为基础, 假定植被冠层为一层或数层水平无限延伸的介质, 每层中的植被元素可以看作吸收和散射微粒, 这些微粒的几何光学特性参数可以给定。然后用辐射传输方程的近似来研究冠层的二向辐射传输特性或反演植被元素特征值。辐射传输模型中最典型的是Suits 模型和SAIL 模型。Suits 模型[1]将叶片、树杆等植被分解成水平和垂直两个分量, 然后根据辐射传输方程对向下、向上的直射光与散射光进行迭代运算, 至数值稳定时得到冠层方向反射率。Suits 模型对植被生长状态等方面有较好反映。SAIL 模型[2]假定植被冠层中叶片倾向为随机分布, 而倾角服从某种特定的分布, 然后根据入射方向和观测方向上叶片分布的统计几何关系计算冠层中9个辐射传输分量的系数。再通过
收稿日期:2002-11-01; 修订日期:2003-01-21
型, SAI L 模型反映更一般的情况, Suits 模型相当于SAIL 模型的一个特例。辐射传输模型适合于对农
作物、草地及密林的遥感, 但对稀疏的林地, 其效果不好。辐射传输模型无论是Suits 模型还是SAIL 模型, 即使在所有参数已知的情况下对地表二向反射率的求解仍要通过迭代运算, 不但运算量巨大, 而且稳定性不好。
几何光学模型以几何光学理论为基础, 从植株等宏观个体入手, 将地面目标假定为一定形状的几何体, 引入了光照植被、阴影植被、光照地面和阴影地面4个分量的概念, 根据该4个参数在不同光照和观测条件下的几何光学关系建立二向反射分布模型。该模型的代表为L-i Strahler 几何光学模型[3]和Jupp 等人提出的模型。几何光学模型能够解释热点现象, 且能够较好地描述稀疏林地和粗糙地表的二向性反射特征。但初期的模型没有考虑多次散射的影响, 对浓密的林地效果不如辐射传输模型。为了能够使适用于浓密的林地, 李小文等进一步提出了考虑树冠相互阴影效应的GOMS 模型[5]。其后为进一步更好考虑冠层多次散射的影响, 李小文等人在间隙率模型的基础上, 进一步提出了几何光学-辐射传输(GO -RT) 混合模型, 从而实现了几何光学模型向辐射传输模型的逼近。
[6][4]
基金项目:国家自然科学基金项目(49899271和49971062) 和国家基础研究计划项目(G2000077900)资助。
作者简介:邓孺孺(1963) ) , 男, 副教授, 1984年本科毕业于武汉地质学院, 1989年于中国地质大学获理学硕士学位, 现在中山大学任教并在中国科学院遥感应用研究所攻读博士学位。目前主要从事辐射传输模型和环境遥感方面的研究, 已发表论文20余篇。
194 遥 感 学 报第8卷
近年来又出现了将物理模型与统计模型相结合的/半经验模型0。/半经验模型0将BRDF 等价为体散射/核0、表面散射/核0和一个表征各向同性常量的加权和, 于是称为/核驱动模参数0, 操作起来较为容易。
核驱动模型虽然有效地考虑反射率角度的变化, 但是这3个参数无法反演地表组分。而前两类模型由于需要复杂运算或对地表已知参数要求过高, 不便于实际运用。提高遥感的定量水平, 需要一个能根据遥感数据定量地求解地面参数, 又便于实际应用的植被遥感模型。本文对冠层的各次散射、透射和土壤反射分别进行分析, 建立各次反射光谱分量与叶面积指数的函数关系。在Michel M 1Ver -straete 等[9]的物理模型的基础上增加了多次散射和土壤反射的分量, 建立一个基于各次散射的植被-土壤二向反射模型。模型的求解无需迭代运算, 便于实际运用。
冠层的第2, 3和n 次散射率。211 一次散射的二向反射率
设植被冠层的厚度为h , 叶面积指数为L , 垂向坐标用z 表示。植被层顶点z 为零, 冠层底部z =h ; 太阳天顶角和仪器观测天顶角分别为H s 和H v 。入射光进入冠层后受到叶子的阻截而削弱, 当通过的冠层厚度为z 时, 其透过率为:
S 1(z )
T 1=exp -=exp -1
Q
z 0
J 1(z ) +(z )
d z 1
(2)
其中:S 1(z ) 是高度在z 以上部分冠层的光学厚度; +(z ) =d L /d h 叶面指数密度, 即单位面积之上单位厚度的叶面积, 单位为m 2m -3, J 1(z ) 是冠层在入射方向上的消光系数:
J 1(z ) =3cos (14z
式中:(1为阳光与叶面法线之间的夹角, 3cos (14z
为垂向坐标z 之上对所有叶片计算cos (1的平均值, 可以根椐植物叶倾角分布函数计算出; L 1=cos H s 。如入射光强为I 0, 那么进入水平冠层的光强在冠层顶部为:I =I 0L 1, 进入厚度为z 的冠层后, 未被叶片拦截的光强为IT (z ) , T (z ) 为透过率。根据Hapke . s [8]等的推导, 阳光直接照射到冠层高度为z +d z 部分的光强为: IT 1(z+d z) =I 0L 1T 1(z) exp -d z
z 1-d S 1(z )
=I 0L 1T 1(z ) exp 1
2 基于多次散射的植被二向反射模型
植被反射光主要为许多叶子反射光的贡献之和, 由于叶子的位置、倾角和倾向不同, 相当一部分
光经叶子反射或透射后又被其它叶子所阻挡, 再次发生反射或透射。这一过程可重复多次, 即发生多次散射作用, 我们分别将经历一次散射, 两次散射和n 次散射后的入射光称为一次散射光, 二次散射光和n 次散射光。散射光的能量主要集中在第1次散射。随散射次数的增加, 散射强度近呈等比级数锐减。故一些辐射传输模型中, 只考虑一次散射, 忽略多次散射, 如Michel M 1Verstraete
[9]
z+d z
其中:d S 1(z ) =
的物理模型。
该模型当叶面指数较大时, 会产生明显的误差。尤其在近红外波段, 由于叶子是透射率和反射率都很高, 多次散射光是不应忽略的。
研究表明, 植物叶在可见光波段反射率和透射率极低, 在红外波段无论是反射或透射都很高, 而且均接近于朗伯体。假设植被叶在所有波段是朗伯体, 植被冠层由无数倾角、倾向均随机分布的叶子组成, 植被层是水平均匀(但垂向上可能是不均匀的) , 且是无限延伸的。植被覆盖区被传感器探测到的地面光谱由植被冠层散射光和土壤散射光组成。植被冠层和土壤散射光由一至多次散射组成。植被冠层的反射率R p 为:
R p =Q 1+Q 2+Q 3+, +Q n
(1)
[10) 12]
Q
z
z+d z
3cos (14z +(z ) d z =J 1(z ) +(z ) d z
冠层中z 到z +d z 的水平薄层所拦截的太阳直射光
为:
IT 1(z+d z) -IT 1(z) =I 0L 1T 1(z ) e xp 1-e xp -U I 0L 1T 1(z)
d S 1(z)
1
d S 1(z)
=I 0T 1(z) d S 1(z) 1
设X 为单张叶片的一次散射反照率, L 2=cos H v , P (g) 散射相函数, 散射角g 为入射光与散射光的夹角, 有以下关系:
cos g =cos H os (U s cos H v +sin H s sin H v c s -U v ) 出冠层的一次散射光为:
d E 1=I 0T 1(z ) T 2(z ) d S 1(z )
4P
(3)
则冠层的z 到z +d z 的水平薄层在H v 方向反射并透
其中:Q 1为叶丛的一次散射率, Q 2, Q 3和Q n 分别为
第3期邓孺孺等:基于多次散射的植被-土壤二向反射模型
[2]
195
整个冠层的一次散射光E 1是上式的积分:
E 1=I 0
4P
其中:
T 2(z ) =exp -S 2(z )
=exp -2
T (z ) T (z ) d S (z )
1
2
1
h
叶倾角无关。将某次散射分为向上半球散射E n +
和向下半球散射E n -两部分(n 为散射次数) , 这时对于一水平薄层的冠层, 有3种情况的散射:其一是半球入射光来自上方, 散射方向为上半球方向(上向
) , 相当于水平冠层薄层的上表面对上球面入射光的球面反射, 我们用X C +表示其半球散射率; 其二是半球入射光来自下方, 散射方向为下半球方向散射, 相当于水平冠层薄层的下表面对下球面入射光的球面反射, 用X C -表示其半球散射率; 其三是半球入射光来自上方, 散射方向为下半球方向(下向) , 或半球入射光来自下方, 散射为上半方向, 相当于球面入射光穿过水平冠层薄层后产生的球面透射光, 用X t 表示其半球散射率。
Q
z
k 2(z ) +(z )
d z 2
式中:S 2(z ) , k 2(z ) 分别为沿H v 方向的光学厚度和消光系数:
k 2(z ) =3cos (24z
(2为观测方向与叶面法线之间的夹角。故得冠层一次散射的二向反射率为:
E 1
I 0L 1L 2h =T 1(z )T 2(z ) d S 1(z ) 4PL 1L 20
Michel M 1Verstraete [9]等给出了上式积分结果的Q s , U s ; H v , U v =1H
一般式, 我们对其进行了简化。实验表明, 当入射天顶角不是很大时, 叶片的正面反射率、背面反射率和透射率在可见光波段反射率均很低, 而在红外波段很接近[10, 12]并接近朗伯反射。为方便起见, 设叶子为朗伯体, 且其正、反面的反射率和透射率相等, 叶子的方向为均匀分布, 则有
Q l l =2PQ =
4P 4P 2
Q l
Q ; J , J 分别为l 为叶子上面的反射率。令R l =212J 1(z ) 和J 2(z ) 的平均值; +(z ) =L /h, 有:
-R l P(g)
H , U ; H , U Q s s v v =1-e 1
1L 2+J 2L 1
J 1J 2
+12
图1 水平冠层二次散射几何关系示意图
Fi g. 1 Sketch map for geometric relationship of scattering of the
second time
(4)
J 1, J 2和P (g ) 取决于叶子角度分布, 可根据叶倾角分布函数导出。当叶倾角为随机(均匀) 分布时, 有:
J 1=J 2=
P(g )
_
_
_
如图1, 对于冠层高度为z 的水平薄层, 其二次散射是将来自各高度的冠层一次反射光进行散射, 可分为上半球入射和下半球入射两部分。上半球入射光E 1+是z 以上的冠层的下向散射光(相应的散射率为X t ) 的积分, 可表示为:
E 1+=I X t
Q Q cos
2P 2
s , n l cos v n
__
_d H l d U l (5)
Q T (z ) T(z
1
z
1
-z ) d S (z )
I =I 0cos H s 。下半球入射光E 1-是z 以下的冠层的上向散射光(相应的散射率为X r +) 的积分, 可表示为:
E 1-=I X r+
s , v 和n l 分别为入射方向, 观测方向和叶面法线方向的单位失量。
212 二次和多次散射的散射率
因叶片反射与透射率接近, 且反射与透射均近于各向均匀, 故植物冠层反射与透射方向纷纭杂乱,
接近于朗伯体, 散射的次数越多, 其各向均匀的程度越高。设二次以上的散射为朗伯体散射。对于半球散射, 冠层消光系直接取决于叶面积指数密度, 而与
Q T (z c ) T (z c -z ) d S (z c )
h z
1
1
高度为z 的水平薄层对上半球入射光再次散射(二次散射) , 对上向散射光, 其散射率为X r +, 而对于下半球入射光, 其再次散射(二次散射) 的上向散射光所对应的散射率为X t 。该冠层薄层总的二次散射强度为:
196 遥 感 学 报第8卷
3
d E 2=
Q X X T (z c ) T(z c -z ) d S (z c ) +
X X Q T (z ) T (z -z ) d S (z ) d S (z )
I
h z
r+
1
t 1
1
z
l -2L -45-e (4+12L +8L 2) -e (13) Q 3+=8
如只考虑三次散射, 观测方向冠层反射率为:R P H 1, H 2, U 1-U 2=Q 1H 1, U 1; H 2, U 2+Q 2++Q 3+
(14)
213 植被冠层的透射21311 透射光的组成
即使在密林区, 也有一部分阳光透过树冠层射到地面。透射出冠层的光有两种, 一种有直接透射光, 即直接从叶丛间隙射到地面的阳光, 另一种为以叶片散射后再透射到地面的散射光。后者包括了多次散射的透射光。故冠层透射率T p 可由下式表示:
T p (H (15) 1) =T 1(H 1) +T 2+T 3+, +T n 其中:T 1为入射光的直接透射率, T 2, T 3和T n 分别为第2, 3和n 次散射光透射到地面的光的比率。
21312 一次透射光和散射透射光的求解
第一次透射率为:
S 1(h)
T 1=exp -=exp -1
-L
t r+
11
等式右边第一项为下半球入射光的二次散射强度, 第二项为上半球入射光的二次散射强度。再次
散射的散射光传播出冠层还要受冠层叶片的阻截削减, 即乘以T (z 1) , 故向上透出冠层的二次散射光为:
d E 2+=IT(z 1)
t r+
X X Q T (z ) T(z
z 1z 10
h
X r+X t T (z c ) T (z c -z 1) d S (z c ) +
1
-z ) d S (z ) d S (z 1)
式中:S (z ) , T (z ) 分别为高度在z 以上冠层部分垂向上的光学厚度和透射率。因对于半球散射, 冠层消光系直接取决于叶面积指数密度, 而与叶倾角无关
[1]
, 故:
d S (z ) =+(z ) d z
(6)
T(z ) =exp -S (z =exp -整层冠层的二次散射率为:Q 2+=
h
h
Q +(z ) d z
z 0
1
(7)
Q k (z ) +(z ) d z
h 0
1
Q T (z ) X X T(z c ) T(z c -z ) d S (z c ) +
X X Q T (z) T(z -z ) d S (z ) d S (z ) (8)
1z
z
r+
t
1
1
(16)
T 1(H 1) =e
第二次散射透率:
T 2=R l
即:
T 2=R l L e -第三次散射透率:
T
3=R 2T(h -z 1) d S (z 1) l
L
h
J 1
1
(17)
r-t
11
因设叶片的正面反射率、背面反射率和透射率三者相等, 其值为Q l , 则有
Q l
X r +=X r -=X t =X ==R l ,
2
上式可写成: Q 2+U R l
Q T(h -z ) T (z ) d S (z )
(18)
Q
2h 03
R l
T(z1) d S (z1)
Q
h 02
T(z )T z -z 1d S (z) (9)
Q
h Q T(z)T
h 0
z 1-z
d S (z)
积分得:
T 3=R 2l
-e
4
3L
同理可得三次散射光的反射率贡献为:
Q 3
+
=
-
Q T(z ) d S (z) Q T
@Q T (z ) T z -z
2
h 01R
n+
l
h 0
n
n h 0
h h
-L -L e +L e +L 2e -42
L
z 1-z
2
1
d S (z 1)
第n +1次散射透射率:
(10)
T n +
1
=
1
R n+l
h 0
h
(19)
d S (z )
h 0
n +1次散射光的反射率贡献为:
Q n+1=
Q T(z ) d S (z ) Q T z -z
@d S (z ) , Q T z -
z d S (z ) @Q T (z ) T z -z d S (z ) (11)
n-
1
n
n-1
1
2
1
h 0
1
Q T(h-z ) d S (z ) Q T z -z
@d S (z ) , Q T z -z d S (z ) @Q T (z ) T z -z d S (z ) (20)
n
n
n-1
n
h
n-1
1
2
1
h 0
1
214 植被冠层-土壤的综合反射
植被覆盖区的反射光谱是由植被冠散射光和透出冠层的下伏土壤反射光谱之和, 即:
R =R p +Q s
(21)
将(6) , (7) 式分别带入(9) , (10) 式积分得:
R l -Q 1-e 2+=2
2
2L
-2L e
-2L
(12)
第3期邓孺孺等:基于多次散射的植被-土壤二向反射模型 197
土壤的入射光部分为冠层的直接透射光, 这部
分光有方向性; 另一部分为多次散射的透射光, 这部分光为半球方向入射。直接透射光经土壤反射后仍会与冠层发生作用:一方面因叶阻截而削弱, 另一方面也有一部分被冠层散射回地面, 而这种反射反复进行。散射下地面的光为半球方向的, 故总的土壤表面入射光为:
I s =I T p (H 1) +I T p (H 1) R s R p (H 0) + I T p (H 0) 1) R s R p (H I T p (H 0) 1) R s R p (H
n y ]时:
I s =I T p (H 1)
1-R s R p (H 0其中:H 0=0, 为垂向上的入射天顶角。向上透射出冠层的土壤反射光为:
E s =I s T p (H 0) R s =I T p (H 0) T p (H 1)
R s
1-R s R p (H 02
+, +
Q b , 观测天顶角0b , Q l =016, R s =012, 太阳天顶角451) Q 4分别为1) 4散射光; R p ) 冠层散射反射率, Q s ) 土壤反射率贡献, R p +
Q =850nm s ) 总反射率, K
图2 植被覆盖区各反射组分模拟反射率曲线图
n
Fig. 2 Simulate reflectance curves of different component
of reflectance
(1) 总的反射率和各次散射反射率在LAI
时均随LAI 的增加而快速增加; 在13后, 总的反射率和各次散射反射率趋于稳定, 随LAI 的增加极其缓慢。
(2) 一次散射对总的反射率的贡献最大, 随散射次数的增加, 各次散射的贡献成近等比级数减小, 在LAI 1时, 前三次散射与土壤反射之和也大于总散射的98%。所以一般情况下, 仅考虑植被的前三次散射即可满足精度要求。
(3) 土壤反射光对总的反射率贡献随LAI 的增大而减小, 当LAI >3时, 土壤反射光趋于0。
21512 模型二向反射率分布特征
设土壤含水量为0, 叶面积指数为2的情况下, 用公式(24) 进行模拟计算。主平面上, 不同太阳天顶角下, 用(24) 式计算的各观测方向的反射率如图3。
之所乘以T p (H 0) 是因为经土壤反射后, 入射光已变为半球散射光, 可不考虑其方向性。又因T p 已考虑了多次散射的透射, 故E s 已包含了土壤反射光经冠层散射并透射出冠层的那一部分光。写成反射率形式为:
Q s =
T p (H 0) T p (H 1) R s
1-R s R p (H 0)
(22)
其中:H 0=0。R s 为土壤反射率, 当考虑土壤湿度时可将下式[13]代入:
R s =R s 0e
-A V
w s ws
(23)
R s 0) 干土反射率, A ws ) 土壤水吸收率, V ws ) 土壤含水量。
将(14) , (22) 和(23) 式代入(21) 式可求植被冠层-土壤的综合反射率R 。 R H 1, U 1; H 2, U 2=Q 1, U 1; H 2, U 2+Q 1H 2++Q 3++
p 0p 1s
1-R s R p (H 0)
(24)
215 植被冠层-土壤反射率的组成特征分析21511 冠层直接反射光与二、三次散射及土壤反射光贡献的比较
以红外波段为例, 假设植物叶反射率为016(为便于与其它文献对比, 比实际反射率稍大) , 叶倾角
为均匀分布, 入射天顶角为45b , 观测天顶角为0b , 并分别代入(4) , (
12) , (13) , (11) 和(4) 式得冠层Q 1) Q 4和冠层总的反射率R p 随面积指数的变化曲线如图2。
从图2中可以得到以下几方面特征:
Q l =0. 6, R s =0. 2, v sw =0, LAI =2, ze s 为观测天顶角
图3 不同太阳和观测天顶角下, 太阳主平面上的计算反
射率曲线
Fig. 3 Calculated reflectance curves in principle plane under
different solar zenith
198 遥 感 学 报第8卷
图3中太阳一侧观测天顶角为负, 相对一侧观测天顶为正, 下同。
当太阳天顶角为50b 不同观测方位角反射率曲线如图4所示
:
表1 3种植被覆盖情况下土壤含水量计算误差对比Table 1 Calculating errors of water content under three
different vegetative coverings
测点号LAI 土壤湿度平均绝对误差平均相对误差最大绝对误差(a) (b) (c)
[1**********]2
3132311715
[***********]007519
[***********]
[**************]1
03313
Q l =0. 6, R s =0. 2, v sw =0, LAI =2, ze
s 为观测天顶角, H b , s =50
az v 为观测方位与太阳方位角的夹角
图4 不同观测方位角和天顶角下的计算反射率曲线Fig. 4 Calculated reflectance curves under different observing
azimuth and zenith
3 土壤-植被反射模型的验证实验
实验地点为北京市顺义县赵全营镇的冬小麦麦地, 时间为2001-04-12, 所用光谱测量仪器为ASD 连续光谱仪。共进行了两项实验, 其一是检验垂直观
测条件下模型对LAI 和土壤含水量计算精度。另一项实验检验多角度观测情况下模型的二向反射精度。首先测量了干土壤和叶子反射率曲线作为计算的基础数据。
311 对LAI 、土壤含水量计算精度检验
该项实验是在垂直观测条件下测量LAI 、土壤含水量不同的3个样点的光谱曲线, 并测量LAI 、土壤含水量。然后将LAI 、土壤含水量代入公式(23) , 公式(24) 计算出光谱曲线, 将计算光谱曲线与测量光谱曲线进行对比, 检验模型对LAI 、土壤含水量计算精度。
图5为3个点的计算光谱曲线和测量光谱曲线, 表1为该3个点测量光谱与计算光谱反射率在350) 2500nm 范围内的平均误差对比。除了水吸收带附近外, 二者仍很接近。表1为该3个点全波段的平均误差。反射率平均绝对误差均在0101内, 接近TM 等仪器的探测精度。312 模型二向反射特征的检验
该项实验的方法是在已量测叶面积指数和土壤
图5 不同土壤含水量和LAI 的情况下麦地计算反射率曲线与测量反射率曲线对比((a) 至(c) 分别为3个土壤含
水量和LAI 不同的观测点)
Fig. 5 Calculated and measured reflectance curves of wheat fields of di fferent LAI and water con ten t((a) to (c) are three
observed points)
含水量的麦地进行主平面和垂直太阳主平面上多角度光谱测量, 得到不同观测角下各波段的光谱曲线; 然后用以上模型根据叶面积指数和土壤含水量计算
第3期邓孺孺等:基于多次散射的植被-土壤二向反射模型 199
各波段在主平面和垂直太阳主平面上、不同观测角
下的反射率, 将之与测量值进行对比。为简单起见, 本模型未考虑热点的作用, 为避免热点的影响, 采用与主平面夹角为10b 的平面代替主平面。图6为叶面积指数LAI 为113, 土壤含水量为6%, 太阳天顶角为45b 时, 太阳主平面和垂直太阳主平面上850nm 波段在不同观测天顶角下的测量反射率曲线与计算反射率曲线, 表2为计算反射率和测量反射率的误差分析。从图6可以看出计算曲线和测量曲线均有相同的趋势, 在垂直太阳主平面上, 二者拟合得很好, 平均绝对误差仅为0101; 受热点和观测条件等因素的影响, 在主平面上, 误差略大, 平均绝对误差为01015。两个平面上的相关系数均大于0197。事实上由于所测麦地的小麦并不茂密, 仪器视域范围有限, 在不同观测角度下落入仪器观测范围内的植物-土壤面积比不尽相同, 这就不可避免地造成一定的误差。考虑到该方面的影响因素, 以上误差是可
以接受的。
K =850nm, Q l =0. 56, R s =0. 24, v s w =6%,LAI =1. 3。ze 天顶角 , H b 。az s =50
s 为观测
v 为观测方位与太阳方位角的夹角
图6 太阳主平面和垂直太阳主平面上, 测量反射率和计
算反射率曲线对比
Fig. 6 Comparing of calculated and measured reflectance curves in principle plane and the plane perpendicular to principle plane. a z
v is the angle between solar azi mu th and observing azimuth
表2 太阳主平面和垂直太阳主平面上, 850nm 波长测量反射率和计算反射率误差对比
Table 2 Error comparing of calculated and measured reflectance on the spectrum of 850nm in principle plane and the
plane perpendicular to principle plane
太阳主平面
观测天顶角/(b )
-60-40-200204060平均误差相关系数
计算值[***********][**************]00
测量值[***********][***********]101501983
绝对误差[**************]-01016-[**************]
计算值[***********][**************]00
垂直太阳主平面
测量值[***********][***********]101001973
绝对误差-01013-[***********]09-01011-01016
于应用, 且可以根据要求的精度决定考虑的散射次
4 结论与讨论
在现有众多被普遍接受的植被二向反射模型中, 普遍存在不便于应用的问题。要么需要进行运算量巨大, 且导致不稳定性增加的迭代运算; 要么要求已知很多实际上很难获得的参数, 使得植被二向反射模型与应用存在相当的距离。本文提出的模型对冠层散射的各分量分别进行考虑, 所要求的参数均可通过实验获得, 无需进行迭代运算, 在实际中易
数, 以减少模型不必要的复杂性。从以上的模型计算和实验检验中还可得到以下结论:(1) 冠层各次散射分量随散射次数的增加而快速减少, 一般情况下, 只考虑前3次散射即可满足精度要求。(2) 反射光谱的方向性由一次散射决定, 二次以上的散射对总体反射光的方向性影响很小。参考文献(References)
[1] Gwynn H. Sui ts. The Calculation of the Di recti onal Reflectance of a
200 遥 感 学 报
Vege tation Canopy [J ]. Re mote Sensing o f Environme nt . 1972, 2(2) :117) 125.
第8卷
and Kerne-l driven Models of B-i di recti onal Reflectance[J]. Journal o f Ge ophysical Re se arch . 1995, 100(D10) :21077) 21089. [8] Hapke, B. W. B-i di recti onal Reflec tance Spectroscopy, 4. The Ex -ti nction Coefficient and The Opposition Effec t [C ]. Icarus, 1986, 67:264) 280.
[9] M ic hel M Verstraete, Bernard Pi nty, Robert E Dickinson. A Phys -i
cal Model of B-i directional Reflectance of Vegetation Canopys Theory [J]. Journal o f Geo physical Researc h . 1995, (D8):11755) 11765.
[10] David M Gates, Harry J Keegan et al . Spectral Properties of Plants
[J]. Applied Optic s . 1965, 4(1) :11) 20.
[11] Joseph T Woolley. Reflectance and Trans mi ttance of Light by Leaves
[J]. Plant Physiology . 1971, 47(5) :656) 662.
[12] H. T. Breece Ó, R A Holmes. B-i directi onal Scattering Charac ter -i stics of Heal thy Green Soybean and Corn leaves i n Vi vo[J]. Applied Optics . , 1971, 10(1) :119) 127.
[13] Deng R uru, Tian Guliang e t al . Research on Remote Sensing Model
for Soil Water on Rough Surface[J]. Journal o f Re mote Sensing [J], (waiti ng for publish) , 邓孺孺, 田国良等. 粗糙表面土壤含水量遥感模型研究[J]. 遥感学报(待刊) .
[2] Verhoef W. Light Scattering by Leaf Layers wi th Applicati on to Canopy
Mideli ng :The SAIL Model [J]. Remote Sensing o f Environment , 1984, 16(2) :125) 141.
[3] Li Xiaowen, A H Strahler. Geometric -Optical Modeli ng of a Conifer
Forest Canopy [J ]. IEEE Transactions on Geosc ienc e and Remote Sensing . 1985, G E -23(5) :705) 721.
[4] Jupp D L B, Walker, J, Penridge, L, K. Interpre tation of Vegeta -tion Structure in Landsat MSS Imagery:A Case Study i n Dis turbed Sem-i aric Eucalypt Woodland. Part2. Mode-l based Anal ysis[J]. J. Envi ronmental Manage me nt , 1986, 23:35) 57.
[5] Li, Xiaowen, A H Strahler. Geo metric -optical B-i directional Re -flectance Modeling of the Di screte Crown Vege tation Canopy:Effect of Crown Shape and Mutual Shadowing [J ]. IEEE Trans. Geosic Re m. Se ns . 1992, GE -30, 276) 292.
[6] Li Xiao wen, AH Strahler. Geome tric -optical B-i direc tional Relectance
Modeli ng of a Xonifer Forest Canopy[J]. IEEE Trans. Geosc i. Re m. Sens . 1986, G E -24, 906) 919.
[7] Wanner W, Li Xiaowen, A H Strahler. On the Derivation of Kernels
B-i directional Reflectance Model of Canopy and S oil Based on Mult-i Scatterings
DENG Ru -ru 1, 2, TI AN Guo -liang 1, LI U Qin -huo 1
(1. Institute o f Remote Se nsing Applications, Chinese Acade my o f Sciences, Be ijing 100101, China ; 2. Center f o r Remote Se nsing Applic ation , Zhong Shan U nive rsity , Guangzhou 510275, China)
Abstract: It . s a primary question for many B-i directional models of Canopy presently that it . s difficult for them to be used in practice due to their complexity. In this paper, mult-i scatterings and transmission taking place in canopy of ve ge -tation, and reflectance of soil are analyzed respectively, and then the functions between every component of scattering and LAI are put forward. Improving Verstraete . s(1986) Physical Model of B -i directional Reflectance of Vegetation Canopy by adding mult-i scattering and soil reflectance in the model, a B-i directional reflectance model of canopy -soil as proposed and it was proved that the model . s precision can be satisfied only when the first three scattering being taken into account. The shortcoming of the original model has been overcome in this model, and B -i directional reflectance can be calculated
directly according to LAI and some other factors, avoiding overlap calculation which results tre mendous amount of compu -tation, so the model is easier applied in practice. The precision of the model was verified by experiments. Key words: b-i direc tional reflectance model; canopy model; muti-l scattering; LAI , soil water content