全 文 ：书第 !" 卷第 " 期
!##" 年 " 月
生 态 学 报
$%&$ ’%()(*+%$,+-+%$
./01 !"，-/1 "
2301，!##"
!""#： %%%& ’()*)+,(-& (.
基金项目：国家重点基础研究发展规划资助项目（!##4%56##7#!，!##8%59#:6#6）；中国科学院“百人计划”择优支持资助项目（;<#7"6:9）；中国
科学院大气物理研究所东亚区域气候<环境重点实验室开放基金资助项目
收稿日期：!##4<#7<9#；修订日期：!##"<#9<9#
作者简介：陈艳华（8:;9 =），男，福建省福州市人，硕士生，主要从事植被遥感与遥感生态应用研究1 ’<>?@0：AB?CDEFG HI?1 ?F1 FC
!通讯作者 %/JJIKL/CM@CD ?3HB/J1 ’<>?@0：AB?CEFG HI?1 ?F1 FC
致谢：美国马里兰大学地理系方洪亮博士和荷兰国家航空实验室的 N/3H .IJB/IO教授在研究过程中给予细心指导，在此表示衷心感谢！
/)0.1-",). ,"’2：&BI LJ/PIFH E?K O@C?CF@?00Q K3LL/JHIM RQ HBI 8## &?0ICH S/3CDHB TJ/PIFH /O %B@CIKI $F?MI>Q /O ,F@ICFIK （-/1 ;< #7"6:9），HBI
-?H@/C?0 UIQ VIWI0/L@CD TJ/DJ?> O/J 5?K@F ,F@ICFIK /O %B@C? （*J?CH -/1 !##4%56##7#!，!##8%59#:6#6）
3’(’,4’1 1-"’：!##4<#7<9#；5((’#"’1 1-"’：!##"<#9<9#
6,)+7-#!8：%X’- S?C?@C0Q ICD?DIM @C JI>/HI KICK@CD O/J WIDIH?H@/C IF/0/DQ1 ’<>?@0：FBICQBZCP3G 8491 F/>
基于分类知识利用神经网络反演叶面积指数
陈艳华8，张万昌!，!，雍[ 斌8
（81南京大学国际地球系统科学研究所，南京[ !8##:9；
!1中国科学院大气物理研究所东亚区域气候<环境重点实验室，全球变化东亚区域研究中心，北京[ 8###!:）
摘要：叶面积指数（!"#，)I?O $JI? +CMI\）是陆面过程中一个十分重要的输入参数，其遥感反演方法研究一直是国内外遥感应用 研究的热点问题。基于统计的遥感反演方法由于缺乏物理基础，其可靠性和普适性差。基于物理的冠层反射模型的 !"# 反演 方法克服了上述弊端，但是由于反演过程是病态的，模型反演结果一般不唯一。神经网络算法的介入可在一定程度上改善这一 问题，但是模型反演的病态问题至今仍无法很好地解决。在 T](,$+)模型敏感性分析的基础上提出了一种基于影像分类的神
经网络反演方法，引进了土壤反射指数用于替代原模型中难以确定的土壤背景反射参数，分别针对不同植被类型建立各自的神
经网络，对经过大气纠正后的 )?CMK?H ’&Y ^影像进行了模拟实验并同野外实测 !"#数据进行比较。结果表明，对于 !"#小于 9
的植被区该方法的反演精度比较可靠，而 !"#大于 9 的植被区，反演的 !"#偏小，原因归结为密植被的冠层反射在 !"#大于 9 以
后趋于饱和而无法敏感地表征 !"#的变化所导致的。
关键词：叶面积指数；神经网络；冠层反射模型；T](,$+)
文章编号：8###<#:99（!##"）#"3’"7,’4,.+ *’-9 -7’- ,.1’: 0;,.+ - .’07-* .’"%)7< =-;’1 ). (*-;;,9,(-",). <.)%*’1+’
%X’- S?C8 #$%&’$(%)*$(+ #$,%)%-%& .*’ /(’%0 12,%&3 14)&$4& （/11#），5($6)$7 8$)9&’,)%2，5($6)$7 !8##:9，:0)$( ! ;&7)*$(+ :+)3(%&&’(%& /(,% ",)(，#$,%)%-%& *. "%3*,>0&’)4 ?02,)4,，:"1，@&)6)$7 8###!:，:0)$( !"#$ %"&’&()"$*)+)"$，>??@，>@（@）：>@AB C >@DE,
5=;"7-("：$K ?C @>L/JH?CH IF/0/D@F?0 L?J?>IHIJ @C 0?CM K3JO?FI LJ/FIKKIK，)I?O$JI? +CMI\ （!"#）?CM @HK @CWIJK@/C E@HB
JI>/HI0Q KICKIM M?H? ?JI B/H H/L@FK @C a3?CH@H?H@WI JI>/HI KICK@CD，R/HB M/>IKH@F?00Q ?CM @CHIJC?H@/C?00Q1 ’>L@J@F?0<
JI0?H@/CKB@LK，/E@CD H/ HBI@J KB/JHF/>@CDK @C LBQK@F?0 >IFB?C@K> MIKFJ@LH@/CK，0?Fb JI0@?R@0@HQ
?CM OI?K@R@0@HQ @C ?LL0@F?H@/C1 TBQK@F?00QK，K3FB ?K HB/KI MIWI0/LIM E@HB HBI F?C/LQ JIO0IFH?CFI >/MI0，
/WIJF/>I HBI ?R/WI>ICH@/CIM KB/JHF/>@CDK；B/EIWIJ，!"# @CWIJK@/C E@HB F?C/LQ HBI JIO0IFH?CFI >/MI0 @K 3K3?00Q @00EB@FB >?bIK HBI @CWIJK@/C C/H 3C@a3I1 ’>L0/Q>ICH /O ? CI3J?0 CIHE/Jb @C !"# JIWIJK@/CK F?C @>LJ/WI K3FB @KK3IK H/ ?
FIJH?@C I\HICH，R3H HBI @00/MI0 @CWIJK@/C @K QIH H/ RI JIK/0WIM1 (C HBI R?K@K /O
KICK@H@W@HQ ?C?0QKIK 3K@CD HBI T](,$+) >/MI0，HBI LJIKICH KH3MQ MI>/CKHJ?HIK ?C ?LLJ/?FB HB?H 3KIK ? CI3J?0 CIHE/Jb R?KIM /C @>?DI F0?KK@O@F?H@/C @CF/JL/J?HIM @C T](,$+) O/J ?FF3J?HI JIHJ@IW?0 /O !"#1 5Q @CF03M@CD HBI K/@0 JIO0IFH?CFI @CMI\ @C HBI
!""#： %%%& ’()*)+,(-& (.
)/,+,.-* 0123456 7)8’* ") "-9’ "!’ #*-(’ ): ;),* <-(9+/)=.8 /’:*’("-.(’ #-/-7’"’/; "!-" -/’ 8,::,(=*" ") 8’"’/7,.’，;#’(,:,(
.’=/-* .’"%)/9; -/’ ().;"/=("’8 ()//’;#).8,.+ ") ,.8,>,8=-* "?#’; ): >’+’"-",). ()>’/& @A#’/,7’."; %,"! 6-.8;-" @BC D
8-"- ,.8,(-"’ "!-" "!’ /’"/,’>-* -((=/-(? ,; !,+!’/ :)/ >’+’"-",). %,"! - !"# *’;; "!-. "!/’’，-.8 -; !"# ,.(/’-;’;，/’"/,’>-*
-((=/-(? 8’(/’-;’; -(()/8,.+*?& B!’ #/,7-/? /’-;). ,; -""/,<="’8 ") (-.)#? /’:*’(",). .) *).+’/ <’,.+ ;’.;,",>’ ") !"# %!’.
"!’ >’+’"-",). ,; ")) 8’.;’*? #)#=*-"’8 （!"# E F）&
!"# $%&’(：6’-: 4/’- 5.8’A；.’=/-* .’"%)/9；(-.)#? /’:*’("-.(’ 7)8’*；0123456 叶面积指数（!"#）是指单位水平地表面积上所有绿色叶子表面积的一半［G］，它控制着植被许多生物、物 理过程，如光合、呼吸、蒸腾、碳循环和降水截获等，是表征植被冠层结构的最基本的参量之一，是陆面过程模 型的重要输入参数［G］。叶面积指数的精确估算在碳循环、能量循环和环境影响评价研究中起着非常重要的 作用［H I F］。 近 FJ-来，众多遥感技术 !"#反演方法被提了出来［K］，从方法学角度主要分为两大类：统计模型法和光学 模型法［L］。由于，传统的 !"#反演算法是基于植被指数的统计回归，易受土壤背景等多种外在因素的影响，缺 少可移植性［M］。因此，科学家们转向反演用于描述冠层内电磁辐射传输过程的辐射传输模型，而冠层辐射传 射模型的反演是一个反问题，具有不适定性［N］，即不满足解的存在性、解的惟一性和解对数据的连续依赖性 中的任何一条。这主要有两个方面的原因造成的［O］：（G）不同的模型参数可能产生几乎一样的反射光谱。比 如，叶子水平分布的稀疏植被冠层的光谱反射与叶子垂直分布的浓密植被冠层的光谱反射相似；（H）模型建 模与观测的反射率存在不确定性。模型的不确定性来自叶子散射行为的简化和叶子朗伯特散射假设等，而观 测的反射率数据的不确定性主要由传感器噪声与数据预处理产生。由于反问题的不适定性，将会导致模型反 演的解在参数空间产生一种跳跃，意味着求出来的解可能在分布在整个参数空间，而不收敛于真实解附近。 P-(Q=’7)=8和 R-/’"［S］利用 0123456模型模拟高光谱反射，假定 !"#保持在一个初始值附近来进行冠层反射 模型数值反演，结果表明 !"# 越小，叶子结构参数$ 越大，平均叶倾角 !"%（C’-. *’-: ,.(*,.’8 -.+*’
8,;"/,<=",).;）也越小，说明 !"#和 !"%与叶子结构参数 $对冠层反射率的影响作用是相反的，不同的参数组合 可能对应几乎相似的光谱。因此，模型反演的病态问题仍然存在，至今仍是遥感模型反演研究中的瓶颈问题。 本研究在基于辐射传输模型的叶面积指数反演方法研究!的基础上，采用 0123456 辐射传输模型，针对 不同的植被状况，采取不同波段组合对模型的各植被参数进行敏感性分析，提出一种基于影像分类的神经网 络反演方法，运用先验知识来减少模型反演的不确定性来提高反演的可靠性。应用 HJJH 年 S 月 H 日到 GF 日 在中国西北地区黑河流域张掖绿洲区获取的 !"#野外观测数据，及近乎同时相 @BC D影像（成像于 HJJH 年 S 月 GN 日）对提出的反演算法进行了验证和评估。 )* 研究区概况与数据预处理 试验区选择在我国西北甘肃省黑河流域的张掖绿洲，位于河西走廊中部，我国第二大内陆河黑河贯穿全 境，形成了特有的荒漠绿洲景象。该区年平均气温为 NT，年均降水量 GKJ77左右，为典型的干旱、半干旱生 态类型区，并且该地势平坦，土壤肥沃，物产丰饶，盛产小麦、玉米、水稻、油菜、胡麻等农作物，为我国西北重要 农业区，也是国家西部生态建设的重点地区。 研究所用的数据主要包括 HJJH 年 S 月 H 日到 GF 日间在中国西北地区黑河流域张掖绿洲区由 B14U 仪 器（跟踪辐射与冠层结构测量仪，B/-(,.+ 1-8,-",). -.8 4/(!,"’("=/’ ): U-.)#,’;）获取的 LM 个 !"# 野外观测数 据，及近乎同时相 @BC D影像（成像时间 HJJH 年 S 月 GN 日）。遥感影像前处理工作包括几何精校正、大气校 正和影像分类。影像的精校正使用 GL 个均匀分布的地面控制点进行以保证校正后影像的平均几何误差小于 MONH V 生V 态V 学V 报V V V HN 卷V ! 陈艳华，张万昌，雍斌& 基于 BC的辐射传输模型反演叶面积指数可行性研究& 国土资源遥感，投稿 !""#：$ $%%%& ’()*)+,(-& (. /& 0 个像元。大气校正是根据 1,*-2’3"［4/］和 5!-.+［44］提出的方法。这种方法是利用暗像元推导出相关的大气 参数来进行大气订正的，是完全基于遥感影像信息的，不需要实测成像时复杂的大气参数，因而简单有效。由 于本研究的重点是确定不同植被类型的敏感性参数及其取值范围，需要分类后的各个植被类型的光谱比较集 中，因此影像分类方法是利用 6789:;:（即重复自组织数据分析技术，6"’3-",<’ 7’*=>83+-.,?,.+ 9-"- :.-*@A,A ;’(!.,BC’）非监督分类算法，并且结合野外调查资料与研究区的土地利用资料，将研究区的影像分为非植被 区、密农作物（玉米与小麦等）、农作物、阔叶林、草地和灌木。野外观测 !"#时，尽量考虑到植被种类的不同， 选择地表覆盖较均匀的地区进行观测，并严格按照 ;D:E 仪器操作规范的要求进行。观测场一般以相隔约 40F的 G 根长 0/F左右的测线布设，取这 G 条断面的平均 !"#为该测点的实际 !"#值。 !" 模型反演方法研究 !& #" 模型选择与简介 目前已经发表了众多的冠层辐射传输模型，一般情况下越精确的模型越复杂，并且需要的输入参数也越 多，计算效率也越低，不容易进行模型反演。一个适用于反演的理想模型应该是既有一定的模拟精度又有满 意的运算效率［4H］。本研究所用的辐射传输模型 ID87:6J 模型是由 7:6J 冠层反射模型［4G］和 ID87IKE; 叶 子光学特性模型［4L］耦合而成的。7:6J模型是个模拟精度与运算效率都非常好的一维模型。7:6J 模型是把 植被当作混浊介质，假设叶片方位角分布均匀，考虑任意的叶片倾角，冠层双向反射率作为观测角度的一个函 数。7:6J模型广泛应用于农业植被类型与阔叶植被，而且模拟精度较高［40，4M］。NCCAO 等［4P］对 7:6J 模型又 进行了改进，考虑进了热点影响，计算单次散射对二向反射率的贡献时，考虑进了叶片以及相应的阴影影响。 ID87IKE;模型是一个基于 :**’.等［4Q］的平板模型改进的用于模拟叶子反射与透射的叶片辐射传输模型，表 述了新鲜植被叶片从 L//.F到 H0//.F的光学特性。平板模型假设叶片是一个致密的、均匀地充满了吸收和 散射物质的、表面为朗伯的平板。入射的漫射辐射光在叶子表面发生反射与透射，在叶子内部由叶绿素、水分 子和木质素等微粒吸收。叶子表面的光学属性使用折射指数和吸收系数来表征。在平板模型的基础上， ID87IKE;模型假设叶片是被$>4 层空气间隔分隔的由 $层均质层叠合起来的平板。ID87:6J冠层反射耦 合模型所涉及到的参数及其符号如表 4 所示。本研究是针对 K;R S光谱数据，所以传感器方位角 %与传感器 天顶角 !都取为 /T。 表 #"$%&’()*模型的参数表
+,-./ #" +0/ .123 45 6,7,8/3/72 547 $%&’()* 849/. 项目 6"’F 参数 I-3-F’"’3A 符号 7@F2)*A 外部参数 KU"’3.-* #-3-F’"’3A 太阳天顶角 7)*-3 ?’.,"! -.+*’（T） %ACF 太阳方位角 7)*-3 -?,FC"! -.+*’（T） !AC. 传感器天顶角 7’.A)3 ?’.,"! -.+*’（T） % 传感器方位角 7’.A)3 -?,FC"! -.+*’（T） ! 水平能见度 V)3,?)."-* <,A,2,*,"@（OF） &#’ 冠层结构和叶子物理参数 叶面积指数 J’-= -3’- ,.W’U !"# E-.)#@ A"3C("C3’ -.W *’-= #!@A,(-* #-3-F’"’3A 平均叶倾角 R’-. *’-= ,.(*,.’W -.+*’ W,A"3,2C",).A （T） !"( 叶子内部结构参数 J’-= ,."’3.-* A"3C("C3’ #-3-F’"’3$
叶绿素 :X含量 J’-= (!*)3)#!@** - S 2 ()."’." （!+ $(FH） )- S 2 叶子干物质含量 J’-= W3@ F-""’3 ()."’." （+$ (FH） )F
叶子等效水厚度 J’-= ’BC,<-*’." %-"’3 "!,(O.’AA （(F） )%
叶大小与冠层高度比值 J’-= A,?’ $(3)# !’,+!" 3-",) ’4 土壤反射率 7),* 3’=*’("-.(’ 土壤反射率 7),* 3’=*’("-.(’ "#* PQPHY P 期 Y Y Y 陈艳华Y 等：基于分类知识利用神经网络反演叶面积指数 Y !""#：$ $%%%& ’()*)+,(-& (. !& !" 土壤反射指数的引进 由于 /012345模型使用的土壤背景反射参数是一个土壤反射光谱，不利于模型的有效反演。因此，需要利 用一个简单的参数替代。研究表明土壤反射率的光谱形状具有相似性，可以用一种土壤反射指数来代表特定的 土壤光谱［67］。对于一种土壤类型，影响土壤反射光谱的一个主要因素是土壤湿度。其对于土壤反射率的影响大 于不同土壤类型造成的影响，并且光谱反射曲线的形状不受土壤湿度的影响。土壤反射函数如式 6所示： !!" 8 !!6 9（!!: ; !!6）< #$% （6）
式中 !!6、!!:分别是裸土在波段 !的最小（湿土）和最大（干土）反射率，本研究采用 32=>0 土壤光谱库的最大
最小反射率来替代。#$% （2),* ?’@*’("-.(’ ,.A’B）是土壤反射指数，代表土壤光谱反射位于裸土最大反射和最 小反射之间的比例。某一植被点 #$%的定义及其量值的确定如图 6 所示。
图 6C &’%等值线与 #$%等值线 D,+& 6C 4E)*,.’ )@ &’% -.A #$%
图 6 所示为 /012345 模型模拟的 &’% 等值线与
#$%等值线图。&’%等值线是指在一定的 &’%和 &’(及 一定的外部条件（太阳角、视角、大气条件）下，改变背 景的光学性质时所观测到的叶冠反射率的点对所组成 的线。如图中实线所示，它既不汇聚原点，也不与土壤 线（即 &’% 8 F 的线）平行。#$% 等值线是指在一定的
#$%和 &’(及一定的外部条件下，改变 &’% 时所观测到 的叶冠反射率的点对组成的线，图中用虚线表示。通过 模型模拟可知，&’% 等值线是一条直线，并且线上的点 分别对应一个 #$% 值。#$% 越大，对应的冠层反射率越 大，相应地冠层反射点对越偏离原点。图中假设点 ) 表示某一植被点，点 * 为 &’% 等值线与土壤线的交点， 点 ’是通过点 ) 的 &’% 等值线与 #$% 8 F 等值线的交
点，点 +是通过点 ) 的 &’% 等值线与 #$% 8 6 等值线的 交点，则定义点 ) 的土壤反射指数为点 ) 和点 ’ 之间 图 :C 不同植被类型的平均光谱反射率 D,+& :C 3G’?-+’ E#’("?HI ?’@*’("-.(’ @)? A,@@’?’." G’+’"-",).E 的距离 J)’ J与点 ’ 和点 + 之间的距离 J ’+ J的比值，即 #$% 8 J)’ JJ’+ J。因此，应用这个公式，可以直接从大气纠
正后的遥感影像计算植被像元的 #$%值，然后运用公式 （6）就可计算每个植被像元的土壤背景反射率。 !& #" 不同植被类型的敏感参数确定 为了减少模型反演的不适定性，需要一些先验知识 对参数的解空间进行限制，并且反演时需要选择一些敏 感的参数进行反演［:F］。首先，自由变化的参数个数越 多，生成的神经网络训练数据集大小就成指数增长，并 且训练时需要占用更多的计算机资源与更长的计算时 间。其次，冠层反射率对不同参数的敏感性不相同，冠 层反射率对 &’%、&’( 和 , 参数敏感，而对 -./、-0、-1 和 #2 参数较不敏感。同时，不同的植被类型对应的敏 感参数及其取值范围是有很大差异，若在反演时给出这 些先验知识对解空间进行限定，可以提高反演的精度。因此，首先根据影像分类结果，统计各植被类型的光谱 范围与平均光谱值（如图 :），然后对不同植被类型进行敏感性分析，确定不同植被类型的敏感参数。根据陈 KKL: C 生C 态C 学C 报C C C :L 卷C !""#：$ $%%%& ’()*)+,(-& (. 艳华和张万昌!的各植被参数敏感性分析的研究结果，本研究针对不同植被类型选择的敏感参数及反演波段 组合如表 / 所示。 表 !" 不同植被类型的敏感参数及反演波段组合 #$%&’ !" (’)*+,+-’ -$.+$%&’* /0 1+00’.’), -’2’,$,+/)$)1 %$)1 3/4%+)$,+/) 0/. 4/1’& +)-’.*+/)
植被类型 0’+’"-",). 敏感参数 1’.2,",3’ 3-4,-5*’2 反演波段组合 6-.7 ()85,.-",).
密农作物 9’.2’ (4)# !"#，!"$，%&#，’ :，; 农作物 <4)# !"#，!"$，%&#，’ :，;
阔叶林 9’(,7=)=2 >)4’2" !"#，!"$，%&#，’ :，; 灌丛 1!4=5 !"#，!"$，%&# :，;，?
草地 @4-22 !"#，!"$，%&# :，;，? !& 5" 数据库生成与神经网络模拟 根据不同植被类型的光谱范围，确定不同植被类型敏感参数的数值范围及固定参数的值，如表 : 所示。 !"#的数值范围主要是根据统计的光谱范围确定，而 !"$ 主要是根据文献资料确定［A，B; C BD，/B］，其它固定参数
是根据平均光谱确定（如图 /）。叶子等效水厚度 () 由 EFG H的第 I 波段和第 ? 波段平均反射率确定，第 I
波段和第 ? 波段越小，() 相对较大。由于研究区为夏季的半干旱区，假定植被类型中最大的 () 为 J& J:(8，
如密农作物与农作物两种植被类型，最小的 () 为 J& JBI(8，如草地，而其它植被类型的 () 根据第 I 和第 ? 波
段的吸收状况内插得到。叶绿素含量 (*+主要由 EFG H的前 : 个波段平均反射率确定，前 : 波段平均反射率
越小，(*+值越大，其计算原理与 () 相似。首先，假定植被类型最大的 (*+为 IJ"+ $(8 /，如密农作物与农作物， 最小的 (*+为 /I"+$ (8
/，如草地，而其它植被类型的 (*+根据前 : 个波段的吸收状况内插得到。叶子干物质含
量 (,是假设鲜叶子中叶水含量占鲜叶总重的 ?/K
［//］。然后，利用冠层反射模型分别模拟各植被类型的光谱
反射数据库，将各个光谱反射数据库随机分成 : 份：一份用于神经网络学习来模拟冠层反射与参数之间的关
系，称为训练数据；一份是用于检查以避免在学习过程中过拟合，称为校正数据；另一份是用于评价模拟性能，
称为检验数据。
表 6" 不同植被类型的敏感参数的数值范围及固定参数取值
#$%&’ 6" 7$)2’* /0 *’)*+,+-’ $)1 ,8’ 0+9’1 -$.+$%&’* 0/. 1+00’.’), -’2’,$,+/)*
类别 <*-22 !"# !"$’ (*+ () (, %- %&# 密农作物 9’.2’ (4)# G,. B& A II B& : IJ J& J: J& JB/ J& JI J G-L D& I DI B& A IJ J& J: J& JB/ J& JI B 农作物 <4)# G,. J II B& : IJ J& J: J& JB/ J& JI J G-L : DI B& A IJ J& J: J& JB/ J& JI B 阔叶林 9’(,7=)=2 >)4’2" G,. J& I ;I B& : ;I J& J/D J& JBJ J& JB J G-L :& I II B& A ;I J& J/D J& JBJ J& JB B 灌丛 1!4=5 G,. J IJ B& I :A J& J/: J& JJM J& J/ J G-L /& I DJ B& I :A J& J/: J& JJM J& J/ B 草地 @4-22 G,. J I? B& I /I J& JBI J& JJ? J& B J G-L / D? B& I /I J& JBI J& JJ? J& B B N N !表中参数单位同表 B O-4-8’"’4 =.,"2 ,. "-5*’ : -4’ ,7’.",(-* ") "!)2’ ,. F-5*’ B 神经网络是一种简化的人脑数学模型，它由大量简单的信息处理单元（神经元）广泛连接而成，用以模拟 人脑神经系统的结构和功能的复杂网络系统。现有的研究成果表明神经网络已具有人脑的一些基本特征，如 学习、记忆和归纳等功能，神经网络最大的特征是大规模并行处理，信息的分布式存储，良好的自适应性和自 MA?/N ? 期 N N N 陈艳华N 等：基于分类知识利用神经网络反演叶面积指数 N ! 陈艳华，张万昌，雍斌& 基于 FG的辐射传输模型反演叶面积指数可行性研究& 国土资源遥感，投稿 !""#：$ $%%%& ’()*)+,(-& (. 组织性，高度的容错性和鲁棒性。神经网络通过“学习”的方式来实现某项特定的任务，根据事先定义好的 “学习规则”，对所提供的特定样本进行学习，在此过程中调整神经系统的内部结构，从而达到记忆、归纳和联 想的功能。人工神经网络发展至今，已经建立了多种神经网络模型，其中应用最广泛的是误差反向传递神经 网络，简称 /0神经网络。/0神经网络作为一种有效的非线性逼近方法，近年来在遥感数据参数反演中已有 应用［12］。常用的 /0神经网络由输入层、隐含层（一层或多层）和输出层组成；相邻之间的神经元由权重系数 相互连接；同一层内的神经元之间是平行的，无连接关系。神经网络的这种拓扑结构决定了各处理单元、各层 之间信息传递方向和途径。/0算法的中心思想是根据模拟的输出值与真实输出值之间的误差，反向调整神 经网络各层的权值，使网络总误差最小，从而达到学习的目的。 图 23 /0神经网络结构示意图 4,+& 23 5"67("76’ 8,-+6-9 ): "!’ /0 .’76-* .’"%)6; 本研究使用的是带有一个隐藏层的后向传播神经 网络（简称 /0神经网络），以光谱反射率及土壤反射指 数 !"#为网络的输入，以叶面积指数$%# 为网络的输
出，以双曲正切函数为神经元的激励函数。图 2 为建立
的神经网络结构示意图。每种植被类型都建立了各自
的神经网络，各个神经网络的结构相同。第 < 层 = 个神
经元，激励函数为双曲正切函数；第 1 层 2 个神经元，激
励函数也为双曲正切函数；而第 2 层为 < 个神经元，与
输出个数一致，激励函数为线性函数，因为线性函数的
输出数值范围可以是任意的。各个网络之间的区别主
要是输入不同。然后，模型模拟的各种植被类型的光谱
反射数据对数变换后用标准化方法作正归化处理，采用 >’?’.@’6+AB-6C7-68" 快速学习算法［1=］训练网络。为
避免过拟合，将 B-(D’E的 /-E’6F正则化方法应用于 >’?’.@’6+AB-6C7-68" 学习算法［1G］，从而改善网络的推广
性能。本研究的多次数值试验显示，与其他训练策略和正则化方法相比，该网络在运算速度与推广能力方面
具有一定的优越性。最后，利用训练数据与校正数据对网络进行训练，并且在训练好之后，利用检验数据对训
练的网络进行评价，计算模拟的 $%# 与真实的$%# 之间的相关系数。结果显示，草地类型的相关系数为
HI JKG，其它 2 种类型的相关系数都为 <。
!& "# $%#反演及精度评价与讨论 图 =-为研究区影像分类图，图 =@为利用公式（<）计算得到的土壤反射指数 !"#图。首先利用 LMB N反 射率影像的第 2 和第 = 波段作散点图，然后用曲线方程拟合散点图的边线，即图 < 中的 !"# O H 和 !"# O < 两 条曲线，最后利用公式 < 计算影像的 !"#值。将经过订正后得到的反射率图像与计算得到的 !"#图像作为输 入，得出$%#值如图 =(所示。其中白色区域为非植被区，主要包括水域、道路、城市用地、裸土、沙漠地和盐碱
地等区域。研究区南边的植被主要以农作物为主，植被茂密，大部分 $%# 值都在 1 和 = 之间。左上角是植被 稀少的区域，$%#值在 H 和 < 之间。黑河流经的区域的植被都比远离黑河的植被茂密，这种空间分布格局与实
际情况一致。
图 G 表示反演的 $%#与实测的$%#的比较，均方根误差 "&!’为 H& 2H<，"1为 H& JHG。从图中可以看出，在
$%#小于 2 时，反演的$%#与实测的 $%# 比较一致，相差不大。而当$%# 大于 2 时，反演的 $%# 普遍比实测的$%#偏小，主要原因可能是 $%#大于 2 时在冠层反射率不再对$%#非常敏感，这个结果同陈艳华和张万昌!的
研究结果一致。要想反演 $%# P 2 的密植被区的$%#值，精确确定植被像元的平均叶倾角值 $%( 是非常重要 的，但是实际条件下，大范围$%(的精确观测或反演是相当困难的，需要使用更精确的三维冠层反射模型和
含有更多植被反射信息的多角度或高光谱遥感数据进行 $%#反演，以便提高密植被区的$%#值反演精度。因
HJQ1 3 生3 态3 学3 报3 3 3 1Q 卷3
! 陈艳华，张万昌，雍斌& 基于 MB的辐射传输模型反演叶面积指数可行性研究& 国土资源遥感，已投稿
!""#： %%%& ’()*)+,(-& (.
接彩图 /
01234 2 期 4 4 4 陈艳华4 等：基于分类知识利用神经网络反演叶面积指数 4
!""#： %%%& ’()*)+,(-& (.
图 /0 反演的 !"#值与实测的 !"#的比较
1,+& /0 !"# ’2",3-"’4 %,"! "!’ 3)4’* 5’6272 !"# )82’65’4
此，一维的 9:;<=>? 模型反演 !"# 具有一定的适用范
围，适用于 !"# @ A 的植被较稀少地区，像干旱与半干旱
区域，并不适用于 !"# B A 的密植被区域，对于密植被区
的反演结果具有较大的不确定性。
本方法反演 !"#的影响因素，主要包括植被类型的
分类精度、模型反演时设置的敏感参数的取值范围和固
定参数取值的精度。由于植被类型分类存在误分类现
象，并且模型反演时设置的敏感参数的取值范围及固定
参数取值与真实情况不相符，带有一定的主观性，会给
反演精度带来一定的影响。若 !"$的设定范围大于实 际范围，将造成反演的 !"# 值偏大；若 !"$ 的设置范围
小于实际范围，将造成反演的 !"# 值偏小。同时，固定
参数的取值只能代表一个植被类型的均值，并不可能完
全代表该植被类型中的所有植被状况。
!" 总结与展望
本文提出了一种基于影像分类的神经网络反演 !"# 方法，分别针对不同植被类型建立各自的神经网络，
运用不同植被类别参数的一些先验知识，对模型反演的参数的解空间进行限制，一定程度上减少了模型反演
的不适定性，提高了模型反演的可靠性。实验结果表明，一维的冠层反射模型 9:;<=>?并不适用于浓密植被
区域的反演，主要原因是在密模型区域模型对 !"#的敏感性降低，降低了 !"#的反演精度。
目前每种模型都是基于一定的假设，不同的模型都有其适应的植被类型，比如，针叶林和阔叶林具有不同
的冠层结构与叶片反射特性，建立的冠层反射模型是不同的，因此，在进行大范围 !"# 遥感反演时，需要综合
利用适用于不同植被类型的冠层反射模型来进行 !"# 反演。同时，模型反演时，不同植被类型的敏感性参数
的取值范围的确定有待进一步研究，准确确定不同植被类型的敏感性参数范围对于减少模型反演的不适定性
具有一定的意义。对于遥感影像数据，CDE F图像是属于多光谱数据，光谱包含的信息量较少，目前许多科学
已经把辐射传输模型反演 !"#的目光转向能够提供更多信息的高光谱数据与多角度遥感数据，可以更大程度
地减少模型反演的不适定性，提高反演的可靠性。
#$%$&$’($)：
［G ］0 H!’. I E，9-5*,( J，K6)%. ?，%& ’(& L’6,5-",). -.4 5-*,4-",). )M H-.-4- %,4’ ()-62’ 6’2)*7",). *’-M -6’- ,.4’N 3-#2 72,.+ !,+! 6’2)*7",). 2-"’**,"’
,3-+’6O -.4 +6)7.4 3’-276’3’."2& :’3)"’ <’.2,.+ )M C.5,6).3’."，PQQG，RQ：GS/ GRT&
［P ］0 >JK9& D!’ >."’6.-",).-* J’)2#!’6’UK,)2#!’6’ 96)+6-3：= 2"74O ,. +*)8-* (!-.+’& D!’ ,.,",-* H)6’ 96)V’("2，:’#)6" W) GP，GXXQ&
［A ］0 <’**’62 9 I，<(!,3’* L <& :’3)"’ 2’.2,.+ )M "!’ *-.4 8,)2#!’6’ -.4 8,)+’)(!’3,2"6O ,. "!’ C;< ’6-：<(,’.(’ #6,)6,",’2，3’"!)42 -.4
,3#*’3’."-",).UC;< *-.4 8,)2#!’6’ -.4 8,)+’)(!’3,(-* (O(*’2 #-.’*2& J*)8-* 9*-.’"-6O H!-.+’，GXXA，Y：PYX PXY&
［T ］0 ?,-.+ < ?& Z7-.","-",5’ :’3)"’ <’.2,.+ )M ?-.4 <76M-(’& W’% [)6\：I)!. ],*’O -.4 <).2，>.(&，PQQA&
［/ ］0 ^,3’2 L <，^ .O-_,\!,. [，96,5’""’ I ?，%& ’(& >.5’62,). 3’"!)42 M)6 #!O2,(-**OU8-2’4 3)4’*2& :’3)"’ <’.2,.+ :’5,’%2，PQQQ，GR：ARG TAX&
［S ］0 1-.+ ‘ Z，a!-.+ ] H& D!’ -##*,(-",). )M 6’3)"’*O 2’.2’4 4-"- ") "!’ ’2",3-",). )M "!’ *’-M -6’- ,.4’N& :’3)"’ <’.2,.+ M)6 ?-.4 b :’2)76(’2，
PQQA，（A）：/R SP&
［Y ］0 ?,7 I Z& H*-22,M,(-",). )M ,.5’62’ #6)8*’3 ,. 3-"!’3-",(2 #!O2,(2 ’c7-",). -.4 "!’ 2)*7",). ") "!’ ,**U#)2’4& H)33& ;. =##*& E-"!& =.4
H)3#7"&，GXRA，（T）：TA SA&
［R ］0 H)38-* K，K-6’" 1，%& ’(& :’"6,’5-* )M (-.)#O 8,)#!O2,(-* 5-6,-8*’2 M6)3 8,4,6’(",).-* 6’M*’("-.(’ 72,.+ #6,)6 ,.M)63-",). ") 2)*5’ "!’ ,**U#)2’4
,.5’62’ #6)8*’3& :’3)"’ <’.2,.+ )M C.5,6).3’."，PQQP，RT：G G/&
［X ］0 I-(c7’3)74 <，K-6’" 1& C2",3-",.+ 5’+’"-",). 8,)#!O2,(-* #-6-3’"’62 8O ,.5’62,). )M - 6’M*’("-.(’ 3)4’* ). !,+! 2#’("6-* 6’2)*7",). 4-"-& H6)#
2"67("76’ -.4 *,+!" 3,(6)(*,3-"’：H!-6-("’6,_-",). -.4 =##*,(-",).2& 9-6,2，16-.(’：>W:=，GXXA& AAX A/Q&
［GQ］0 J,*-8’6" E =& =. -"3)2#!’6,( ()66’(",). 3’"!)4 M)6 "!’ -7")3-",( 6’"6,’5-* )M 276M-(’ 6’M*’("-.(’ M6)3 DE ,3-+’2 & >."& I& :’3)"’ <’.2，GXXT，G/
PXYP 0 生0 态0 学0 报0 0 0 PY 卷0
!""#： %%%& ’()*)+,(-& (.
（/0）：1023 1042&
［//］5 6789: ; <，=->-+?,(!, =，@+-% A& BC-*?-",). )D "!’ BDD’(" )D EF’#F)(’GG,.+ )D "!’ H’>)"’*I J’.G’K L-"- ). "!’ 8("?-* BC-#)F-",).，J?FD-(’
J),* M),G"?F’ M-##,.+ NI -. 8##F)-(! OG,.+ P-.KG-"，LBM -.K M’"’)F)*)+,(-* L-"-& :’)(-F") Q."’F，1000，/3（R）：3S 2S&
［/1］5 T-(U?’>)?K J，V-()?F <，!" #$& <)>#-F,G). )D D)?F F-K,-",C’ "F-.GD’F >)K’*G ") G,>?*-"’ #*-." (-.)#,’G F’D*’("-.(’：L,F’(" -.K ,.C’FG’ >)K’& H’>)"’ J’.G,.+ )D B.C,F).>’."，/W43，/4，103 11X& ［/X］5 Y’F!)’D ;& P,+!" G(-""’F,.+ NI *’-D *-I’FG %,"! -##*,(-",). ") (-.)#I F’D*’("-.(’ >)K’*,.+："!’ J8QP >)K’*& H’>)"’ J’.G,.+ )D B.C,F).>’."， /W4R，/2：/13 /R/& ［/R］5 T-(U?’>)?K J，V-F’" Z& EH@JEB<[：8 >)K’* )D *’-D )#",(-* #F)#’F",’G G#’("F- & H’>)"’ J’.G,.+ )D B.C,F).>’."，/WW0，XR：S3 W/& ［/3］5 :)’* 9 J，[!)>#G). H P& Q.C’FG,). )D C’+’"-",). (-.)#I F’D*’("-.(’ >)K’*G D)F ’G",>-",.+ -+F).)>,( C-F,-N*’G5 !& [)"-* ,.C’FG,). )D "!’ J8QP >)K’*& H’>)"’ J’.G,.+ )D B.C,F).>’."，/W4R，/3：1XS 13X& ［/2］5 M-\)F L，J(!--*\’ : V，;,’+-.K <，V*-K V P& 8((?F-(I -.K G’.G,",C,"I -.-*IG’G )D J8QP >)K’* #F’K,("’K F’D*’("-.(’ )D >-,]’& H’>)"’ J’.G,.+ )D B.C,F).>’."，/WW1，R/：2/ S0& ［/S］5 A??G^ 8& [!’ !)"_G#)" ’DD’(" ,. #*-." (-.)#I F’D*’("-.(’，E!)")._Y’+’"-",). Q."’F-(",).G，8##*,(-",). ,. @#",(-* H’>)"’ J’.G,.+ -.K E*-." B()*)+I& 9’% =)F^：OJ8 J#F,.+’F Y’F*-+，/WW/& /XW /3W& ［/4］5 8**’. ; 8，:-?G>-. 7 ;，H,(!-FKG). 8 T -.K [!)>-G T H& Q."’F-(",). )D QG)"F)#,( P,+!" %,"! - <)>#-(" E*-." P’-D& T& @#"& J)(& 8>，/W2W， 3W（/0）：/XS2 /XSW& ［/W］5 P,-.+ J P，J"F-!*’F 8 7& 8. -.-*I",( VHLZ >)K’* )D (-.)#I F-K,-",C’ "F-.GD’F -.K ,"G ,.C’FG,).& QBBB [F-.G-(",). ). :’)G(,’.(’ -.K H’>)"’ J’.G,.+，/WWX，X/：/04/ /0W1& ［10］5 P, ‘ ;，;-.+ T L，!" #$& [!’ 8##*,(-",). )D EF,). A.)%*’K+’ ,. H’>)"’ J’.G,.+ H’"F,’C-& J(,’.(’ ,. ［1/］5 Z-.+ 7 P，P,-.+ J P& H’"F,’C,.+ *’-D -F’- ,.K’a ?G,.+ - +’.’",( -*+)F,"!> %,"! - (-.)#I F-K,-",C’ "F-.GD’F >)K’*& H’>)"’ J’.G,.+ )D B.C,F).>’."G，
100X，43：13S 1S0&
［11］5 E,’"F) <’((-")，G"b#!’.’ Z*-GG’，!" #$& L’"’(",.+ C’+’"-",). *’-D %-"’F ()."’." ?G,.+ F’D*’("-.(’ ,. "!’ )#",(-* K)>-,.& H’>)"’ J’.G,.+ )D
B.C,F).>’."，100/，SS：11 XX&
［1X］5 Z-.+ 7 P，P,-.+ J P& H’"F,’C,.+ *’-D -F’- ,.K’a %,"! - .’?F-* .’"%)F^ >’"!)K：G,>?*-",). -.K C-*,K-",).& QBBB [F-.G-(",).G ). :’)G(,’.(’ -.K
H’>)"’ J’.G,.+，100X，R/（2）：1031 1021&
［1R］5 7-+-. M [，M’.!-\ M& [F-,.,.+ D’’K_D)F%-FK .’"%)F^G %,"! "!’ M-FU?-FK" -*+)F,"!>& QBBB [F-.G-(",).G ). 9’?F-* 9’"%)F^G，/WWR，3（2）：W4W
WWX&
［13］5 L’>?"! 7，V’-*’ M& 9’?F-* 9’"%)F^ [))*N)a OG’FcG :?,K’& 9-",(^：[!’ M-"! ;)F^G，/WWS& /XS /WR&
参考文献：
［2 ］5 方秀琴，张万昌& 叶面积指数（%&’）的遥感定量方法综述&国土资源遥感，100X，（X）：34 d 21&
［S ］5 刘佳琦&数学物理方程反问题的分类及不适定性问题求解&应用数学与计算数学，/W4X，（R）：RX d 2X&
［10］5 李小文，王锦地，等& 先验知识在遥感反演中的作用& 中国科学（L辑），/WW4，14（/）：2S d S1&
XWS15 S 期 5 5 5 陈艳华5 等：基于分类知识利用神经网络反演叶面积指数 5