锡林浩特草原陆表热红外发射率光谱野外测量.pdf

书书书锡林浩特草原陆表热红外发射率光谱野外测量张勇，杨虎，郑照军，张立军，陆其峰，李贵才，杨忠东（中国遥感卫星辐射测量与定标重点开放实验室 国家卫星气象中心，北京 ）摘要：锡林浩特草原陆表发射率光谱是利用陆表场地进行卫星遥感器红外通道绝对辐射定标和产品真实性检验的关键因子之一。基于光谱平滑的温度与发射率分离反演迭代算法，利用高精度的波曼（）傅立叶变换红外光谱仪和红外标准板，对锡林浩特草原陆表发射率光谱进行测量。获得了不同时间和观测条件下测量的陆表发射率光谱数据。利用该发射率光谱测量结果，可以利用锡林浩特草原对目前国内外主流的遥感卫星热红外通道进行在轨场地绝对辐射定标，并对热红外遥感产品进行真实性检验。关键词：草原；热红外遥感；发射率光谱；野外测量；真实性检验中图分类号：；文献标识码：文章编号：（）地表发射率是描述地表红外辐射特性的基本参数，也是利用陆表作为卫星遥感器红外通道场地绝对辐射定标靶区和进行陆表热红外遥感产品真实性检验所必须掌握的关键陆表参数之一。目前，国内外的红外卫星遥感器在轨场地绝对辐射定标主要是利用高海拔、大气干洁、人为扰动少、温度场分布均匀的高原湖泊作为定标靶区，可以达到很好的定标精度；但是，这些目标区域往往水表温度较低，且低于卫星对地观测的绝大部分目标；只能满足辐射定标线性低端的精度要求。如果对定标结果向辐射响应的高端进行线性外推，则无法保证定标的精度要求。因此，必须找到一个辐射响应高端的陆面目标作为定标靶区进行在轨定标。要利用高辐亮度的陆表进行辐射定标，则陆表发射率是必须考虑的关键因子之一，也是关系到定标精度和成败与否的关键参数。近些年来，无论是室内还是野外，人们都在努力探索各种方法来测量地表发射率。但实验室测量的发射率数据往往不能满足实际遥感应用的需求，因为发射率不仅依赖于地表物体的组成成分，而且与物体的表面状态（表面粗糙度等）及物理性质（介电常数，含水量、温度等）有关，并随着所测定的辐射能的波长（）、观测角度（）等条件的变化而变化。在发射率光谱的野外测量方面，国内外开展了大量的研究工作。等、和 、和 、等为了不同的应用目的，开展了土壤热红外发射率光谱的野外测量试验；目前国内在这方面也开展了重要的研究，张仁华和田国良提出黑体筒法来测量土壤的发射率；唐新斋等利用黑体筒封闭法测定发射率，该方法大大简化了发射率的测量及计算过程，并针对这种方法分析了目标物体、环境的出射辐照度的波动对发射率精确测定的影响；肖青和柳钦火利用 傅里叶变换红外光谱仪测量得到了土壤的发射率光谱，并分析了不同组分、粒径及含水量土壤的热红外发射率变化规律；王先兵等利用傅里叶变换红外光谱仪（）进行了地物目标发射率光谱的野外测量方法研究。近年来，利用遥感手段进行草地研究，包括草原产草量遥感监测、草地植被盖度遥感研究，以及草地植被长势遥感监测 等方面也得到了极大进展。本研究基于 红外光谱仪和红外标准板探讨了锡林浩特陆表发射率光谱的野外测量方法，并对不同时间、不同陆表结构测量获得的发射率光谱进行比较分析，表明在垂直测量的模式下，不同的陆表几何结构对发射率光谱的影响主要表现在绝对数值的差异上，对谱形的变化影响很小，获取的发射率具有较好的一致性，发射率光谱测量精度较高，可以利用该组发射率光谱数据应用于国内外不同热红外遥感器的在轨场地绝对辐射定标和热红外遥感产品真实性检验。第 卷第期 ，草业学报 年 月收稿日期：；改回日期：基金项目：公益性行业（气象）科研专项项目（）和国家自然科学基金（）资助。作者简介：张勇（），男，陕西镇巴人，助理研究员，博士。：材料与方法 热红外发射率光谱的野外测量原理在 的热红外光谱范围内，由辐射传输理论可知，假定地表为朗伯表面，大气下行热辐射各向同性，卫星所接收到的光谱辐射是太阳光谱辐射、大气及地物目标三者相互作用的总贡献。主要包括部分：第部分地物目标发射的热辐射，其大小由物体表面温度和发射率及目标与卫星之间的大气透过率决定；第部分为地物目标对大气下行辐射、环境背景辐射和太阳入射中的热辐射部分的总辐射量的反射辐射，通常情况下太阳入射中的热辐射部分忽略不计；第部分为地物目标与卫星之间的大气上行辐射，它与大气中吸收气体的含量和物理状态有关。因此，卫星入瞳处分谱的热辐射可以表达为：犔犜 犗 犃（）（）犔犅 犅（，犜犛）（）犔犪 狋 犿（）犪 狋 犿犔犪 狋 犿（）（）式中，犔犜 犗 犃（）为传感器接收到的光谱辐亮度；犔犅 犅（，犜犛）为地表真实温度犜犛所对应的黑体辐亮度；（）为地表发射率；犪 狋 犿为大气透过率；犔犪 狋 犿（）、犔犪 狋 犿（）分别为大气下行和上行辐射。利用 红外光谱仪在测量地表时探测器距离地面高度在以内；根据 模式的模拟结果，在 的波段内，该距离下大气的上行辐射犔犪 狋 犿（）小于探测器接收的总辐射能量的，大气透过率犪 狋 犿近乎为，因此可以忽略这项；于是利用光谱仪近地表观测时，式（）可以进一步简化为：犔犛（）（）犔犅 犅（，犜犛）（）犔犪 狋 犿（）（）式中，犔犛（）为仪器在近地表测量目标得到的光谱辐亮度；犔犪 狋 犿（）为大气下行辐射。由式（）可以求解出地面观测目标的发射率光谱（）：（）犔犛（）犔犪 狋 犿（）犔犅 犅（，犜犛）犔犪 狋 犿（）（）式中，表明地表发射率光谱测量必需获得被测目标表面的出射辐射、大气下行辐射和目标表面温度个物理量。发射率测量过程中需把握的个关键问题：一是大气下行辐射或环境辐射的精确测量；二是测准物体表面温度。为解决这两方面问题，本研究在分析现有测量方法的基础上，利用高精度的 傅立叶变换红外光谱测量仪，并采用红外标准板对环境辐射进行了测量，同时利用光谱平滑迭代（）算法进行温度和发射率的分离反演。野外测量系统本研究用于野外的整个光谱测量系统主要由红外光谱仪、计算机、红外标准板、黑体、辐射测温仪、接触式温度计、自动气象站，以及 定位设备等组成。傅立叶变换红外光谱辐射仪是加拿大 仪器公司 系列产品中的一种。它的光学系统是个 干涉仪。红外傅立叶变换光谱辐射仪配备有个中长波红外探测器（碲镉汞）（，），个短波红外探测器（锑化铟）（，）和个近红外探测器（硅）（，）。用于锡林浩特草原陆表发射率测量的是中长波红外探测器。的光谱分辨率分别为，或。本研究野外测量时选用的分辨率为。另外，专门配备了台便携式的黑体源 （美国 公司生产），用作红外光谱仪辐射定标。本试验中所采用的红外标准板是由 公司生产。红外标准板在近红外、中红外和热红外波段具有 的反射率。在 波段每 间隔的反射率都进行了精确的测量，可以追溯到美国国家标准研究所的标准。在天顶角 内每隔 给出了其反射率变化的方向特性，如图所示，可以看出该红外标准板在 的波段内具有很高的反射率，并且整个板面具有较好的朗伯性。野外测量步骤）选择地势平坦，地表均匀的区域进行测量；）在近地表无风无云的天气条件下，测量时间应选择在太阳高度角变化缓慢，地表温度较稳定的时间段；或没有太阳，且地表已经处于热平衡状态下的夜间进行观测；）利用红外辐射测温仪获取地表亮温初值，然后在该温度附近选择个温度点，将黑体控温在这几个点上对 红外光谱仪进行定标测量；）利用 测量红外标准板，采集数据；）用接触式点温计测量 （），图红外标准板发射率方向特性示意图犉 犻 犵 犇 犻 狉 犲 犮 狋 犻 狅 狀 犪 犾 犮 犺 犪 狉 犪 犮 狋 犲 狉 犻 狊 狋 犻 犮 狊狅 犳 犻 狀 犳 狉 犪 狉 犲 犱犵 狅 犾 犱 犲 狀犫 狅 犪 狉 犱红外标准板的物理温度；）然后用 测量地面目标区域；）利用 定位仪获取观测区域地理位置数据。结果与分析 热红外发射率光谱的野外测量目标区域本次发射率测量试验于 年月在锡林浩特国家观象台牧业气象试验站进行。在试验牧场区域选择地势平坦，地表植被分布较均匀的区域进行测量，试验区主要分布着克氏针茅（犛 狋 犻 狆 犪犽 狉 狔 犾 狅 狏 犻 犻）、羊草（犔 犲 狔 犿 狌 狊犮 犺 犻 狀 犲 狀 狊 犻 狊）和蒙古葱（犃 犾 犾 犻 狌 犿犿 狅 狀 犵 狅 犾 犻 犮 狌 犿）等草原植被，属于典型的草原型地表覆盖。本次红外辐射特性测量一共进行了次，分别在 年月 日凌晨点、日凌晨点、日下午 点和 日下午进行，对天然草场、草地植被和裸土目标进行了 次重复测量。测点位置如图所示。观测目标区域地表情况如图所示。红外光谱仪辐射定标辐射定标是将红外光谱仪观测到的辐射量与仪器的输出计数值之间建立相关关系的物理过程。观测目标出射的光谱辐亮度由普朗克函数给出，如下：犔（）犆 （犆犜）（）式中，犆 ，单位：（），为第辐射常数；犆 ，为第辐射常数；为波数，单位：；犜为温度，单位：；犔（）为光谱辐亮度，单位：（）。式（）是普朗克函数用波数表达的形式，用波长表达的形式为：犕（）犆 （犆 犜）（）因为，犱 犱，故：第 卷第期草业学报 年图锡林浩特草原陆表发射率测量示意图犉 犻 犵 犔 狅 犮 犪 狋 犻 狅 狀 狊狅 犳 犲 犿 犻 狊 狊 犻 狏 犻 狋 狔犿 犲 犪 狊 狌 狉 犲 犿 犲 狀 狋 狊 犻 狀犇 狌 狀 犺 狌 犪 狀 犵犮 犪 犾 犻 犫 狉 犪 狋 犻 狅 狀狊 犻 狋 犲图观测目标区地表照片犉 犻 犵 犘 犺 狅 狋 狅 犵 狉 犪 狆 犺 狊狅 犳犿 犲 狋 狉 犻 犮 犪 犾 狋 犪 狉 犵 犲 狋 狊 狌 狉 犳 犪 犮 犲犔（）犱 犕（）犱（）理想的红外光谱仪，其辐射响应是线性的。也就是说，每个光谱通道的目标测量信号都与探测器测得的辐射光谱功率成正比。探测器探测到的功率包含两部分：目标的光谱辐射功率和辐射仪自身的热辐射。如果辐射仪自身的热辐射是不能忽略的话，仪器系统特性测量（系统定标）就至少需要个点，如图，且点的连线不通过原点。个待标定的测值可表示为：犛犕 犲 犪 狊 狌 狉 犲 犱（）犓（）犔犛 狅 狌 狉 犮 犲（）犕犛 狋 狉 犪 狔（）（）式中，犛犕 犲 犪 狊 狌 狉 犲 犱（）是测量的复数原始光谱；犓（）是仪器的响应函数（）；犔犛 狅 狌 狉 犮 犲（）是目标光谱辐亮度（）；犕犛 狋 狉 犪 狔（）是杂散辐射的辐亮度（）。（），由于分束器的偏差（波长折射指数的不同）和探图目标光谱辐亮度与探测器接收到的光谱功率线性关系犉 犻 犵 犔 犻 狀 犲 犪 狉狉 犲 犾 犪 狋 犻 狅 狀 狊 犺 犻 狆犫 犲 狋 狑 犲 犲 狀狋 犪 狉 犵 犲 狋 狉 犪 犱 犻 犪 狀 犮 犲犪 狀 犱犱 犲 狋 犲 犮 狋 狅 狉 犻 狀 犮 犲 狆 狋狆 狅 狑 犲 狉测器的各级放大步骤的误差（各频率电子延迟的不同），干涉图总有一定的不对称性。这种不对称使得干涉图（光谱）的傅 立叶变 换 产 生 了 虚 数 项。在 红外光谱仪的设计中，定标算法直接用复数光谱即可进行校正，而无须从复数光谱中分离出光谱的实部。辐射定标分成步。第，确定探测器的个辐射定标系数：增益犓（）和截距犕犛 狋 狉 犪 狔（）；第，将待定标的测值产生个标定过的光谱。要计算获得增益和截距需要次测量，即个高温黑体测量（比待测目标高）和个低温黑体测量（比待测目标低），这就是如图所示的点定标。由式（）可推出测量黑体获得的辐亮度：犔犎（）犛犕 犲 犪 狊 狌 狉 犲 犱犎（）犓（）犕犛 狋 狉 犪 狔（）（）犔犆（）犛犕 犲 犪 狊 狌 狉 犲 犱犆（）犓（）犕犛 狋 狉 犪 狔（）（）式中，犔犎（）和犔犆（）是由黑体发射率与普朗克函数计算出的黑体在某温度处的光谱辐亮度，为：犔狓（）狓（）犆 （犆 犜）（）式中，狓犆对应低温黑体，狓犎对应高温黑体；狓（）是狓黑体的发射率；犜狓是狓黑体的温度（犓）求解（）、（）两式，可得到增益和截距为：犓（）犛犕 犲 犪 狊 狌 狉 犲 犱犎（）犛犕 犲 犪 狊 狌 狉 犲 犱犆（）犔犎（）犔犆（）（）犕犛 狋 狉 犪 狔（）犔犎（）犛犕 犲 犪 狊 狌 狉 犲 犱犆（）犔犆（）犛犕 犲 犪 狊 狌 狉 犲 犱犎（）犛犕 犲 犪 狊 狌 狉 犲 犱犎（）犛犕 犲 犪 狊 狌 狉 犲 犱犆（）（）就可得到被标定的光谱为：犛犆 犪 犾 犻 犫 狉 狋 犲 犱（）犛犕 犲 犪 狊 狌 狉 犲 犱（）犓（）犕犛 狋 狉 犪 狔（）（）以上方法要求测量定标光谱和测量目标光谱时的仪器配置保持一致。所有的光谱信息将均包含在复数光谱的实部，而虚部只包含噪音。在野外测量中，首先用红外辐射测温仪测量目标区域的地表亮温初始值；为了保证定标精度和监测仪器响应的线性情况，采用点法进行辐射定标，即以辐射测温仪的测值为中心，以 为间隔对 黑体进行温控。然后利用 红外光谱仪对黑体进行测量，从而完成对仪器的辐射定标。大气下行辐射的测量大气下行辐射到达地表被反射的部分能量是目标表面整个出射辐射场的一个非常重要的组成部分，在测量地物目标的真实发射率光谱时必须扣除反射的大气下行辐射的影响。在大气自身的发射和散射中，水汽通常是最主要的贡献成分，并且云是最为重要的辐射源。因此，为了减小大气因素的影响，通常选择在大气水气含量低、晴朗少云的天气条件下进行地物目标发射率光谱的野外测量。这样有助于减少大气的热辐射贡献，并且也可以保证在利用 光谱仪测量地表辐射亮度和大气下行辐射的时间段内天空向下的热辐射总量基本保持不变，从而保证发射率测量的准确性。大气下行辐射的测量通常是利用一个朗伯体漫反射标准板来辅助完成的。在本次测量试验中使用的是 公司生产红外标准板，它是在铝制的底板上面镀上漫射金薄膜。大气下行辐射测量通常是在测量目标之前测量红外标准板的辐亮度。为了准确测量向下大气辐射被目标表面反射的辐射分量，金板的位置和方位第 卷第期草业学报 年应当与地物目标的位置和方位一致。将金板放置在地物目标上，就可以使两者的方位和位置基本达到一致。根据 光谱仪对红外标准板测量输出的电压计数值犇犖犌 狅 犾 犱，利用公式（）就可以计算出红外标准板的出射辐射亮度犔犌 狅 犾 犱：犔犌 狅 犾 犱（）犇犖犌 狅 犾 犱（）犓（）犕犛 狋 狉 犪 狔（）（）式中，犓（）和犕犛 狋 狉 犪 狔（）分别为仪器的增益和截距。红外标准板的反射率变化范围为 （图），它并不是一个理想的反射表面（发射率约为 ），因而在精确计算大气下行辐射时需要扣除红外标准板的热辐射贡献量。为了修正在测量大气下行辐射时由红外标准板自身的发射辐射引起的误差，在测量过程中利用接触式温度计来监测红外标准板的温度，接触式温度计的测量精度为。根据接触式温度计的测量值可以计算出红外标准板的热辐射贡献，进而利用公式（）来计算大气下行辐射：犔犪 狋 犿（）犔犌 狅 犾 犱（）犌 狅 犾 犱（）犔犅 犅（，犜犌 狅 犾 犱）犌 狅 犾 犱（）（）式中，犔犪 狋 犿（）是大气下行辐射，犔犌 狅 犾 犱（）是光谱仪测量的红外标准板出射辐亮度，犌 狅 犾 犱（）是红外标准板的光谱发射率，犔犅 犅（，犜犌 狅 犾 犱）是红外标准板温度的普朗克分布函数。由于野外测量时天气条件并不是稳定不变的，所以及时准确地测量大气下行辐射是至关重要的。如果在测量地表辐亮度和红外标准板辐射亮度的间隔时间内，大气条件发生变化，引起大气下行的辐射场发生变化，造成实际测量的大气下行辐射并不等于测量地表时到达目标表面的天空背景热辐射，那么计算出的地表发射率包含的误差将正比于大气辐射场变化的大小。此外，红外标准板镀金表面也不是一个理想的朗伯体漫反射表面，在实际的野外测量过程中，应当避免在太阳的镜面反射角进行测量；并且在测量过程中避免红外标准板受到外部环境的热辐射影响。观测目标辐亮度测量自然界中大多数地物都是固态的，这些物体表面固有的物理特性 发射率是由其物质组成和表面结构决定的。因此，在测量固体目标表面辐射亮度时，应当对测量目标的物质组成、表面结构和热惯量特性有所了解。尤其是热惯量特性对野外地物目标光谱辐射亮度的测量是十分重要的，因为热惯量是物质对改变其表面温度的外部条件的抵抗能力。如果物体具有很高的热惯量，则对环境温度变化具有很强的抵抗能力，那么在一次测量期间，风和强烈的太阳照射并不显著影响这类目标表面的温度。而在同样的风力和太阳照射等外部环境条件下，对于具有低热惯量的目标，在光谱仪的干涉图扫描期间目标表面温度较易发生变化，造成测量到的目标光谱辐射亮度带有较大的误差。因此，在利用 光谱仪测量锡林浩特草原这类较低热惯量目标时，应尽可能将外部环境的变化控制到最小。在野外利用 光谱仪测量地物目标的辐射亮度时，通常需要将光谱仪的光学头部保持在距离地面目标小于的高度，以减少大气上行辐射的影响。光谱仪测量地物目标时需要通过多次扫描来叠加平均条目标光谱曲线，以达到有效的信噪比。为了控制由于扫描次数过多，持续时间过长，引起仪器内部温度的变化，通常是在保证采集的光谱数据有足够高的信噪比的情况用最短的时间对目标进行或 次扫描，来获得条目标光谱曲线。目标温度与发射率的分离反演由辐射测量数据反演地物的温度和发射率，本质上是一个病态问题。无论光谱细分到何种程度，都属于由犖个方程组成方程组求解犖个未知数（个地表温度和犖个波段的发射率值），虽然光谱细分可以产生一些有用的约束条件，但它并不能改变问题的本质。近 年的研究都集中在如何采取一定的假设和近似，构造多余观测（或者减少待反演参数），使方程组完备，形成了一些代表性的温度发射率分离算法 ，但这些假设随测量条件和应用领域的不同，其合理性和适应性也不尽相同。本研究针对红外高光谱测量数据，基于 算法对温度和发射率光谱进行分离反演。算法的基本思想是，对于高光谱数据而言，地物的发射率光谱要比大气的下行辐射平滑得多，精确的地表温度能够很好地消除地表发射率光谱中的大气吸收线，得到个较为平滑的 （），地表发射率光谱，否则地表发射率光谱中就会有残留的大气吸收线，必然会使地表的发射率光谱变得粗糙。给定一个地表温度初值，以初值为基础，以一定的温差间隔产生个温度分布范围，计算每个温度对应的地表发射率曲线的平滑度指数，具有最小平滑度指数的发射率曲线对应的温度作为地表的温度，最后由式（）得到地表的发射率。可见地物温度确定的精度以及发射率提取的精度和平滑度指数的定义是紧密联系在一起的，程洁等 针对已有的平滑度指数以及建立在物理特性基础上可能的平滑度指数定义进行分析，评价它们对 算法精度的影响，给出最优的平滑度指数。本试验中所用的 傅立叶变换红外光谱仪的光谱分辨率为，肖青和柳钦火用此光谱分辨率模拟计算了仪器噪声对反演精度的影响。当噪声小于 （等效于信噪比大于 ）时，算法可以获得优于 的反演精度（实际 光谱仪的定标精度优于 ），发射率的反演精度优于 。发射率光谱野外测量结果本次测量分别在 年月 日凌晨点、日凌晨点、日下午 点和 日下午进行，对天然草场、草地植被和裸土目标进行了 次重复测量。对有效数据进行处理，取平均后获得的发射率光谱如图所示。图锡林浩特草原陆表发射率光谱测量结果犉 犻 犵 犕 犲 犪 狊 狌 狉 犲 犱犲 犿 犻 狊 狊 犻 狏 犻 狋 狔 狊 狆 犲 犮 狋 狉 狌 犿狅 犳犡 犻 犾 犻 狀 犺 犪 狅 狋 犲犿 犲 犪 犱 狅 狑狊 狌 狉 犳 犪 犮 犲不同条件测量的发射率光谱，天然草地、铝箔背景上人工均匀种植的草丛（草丛密度：丛）、铝箔背景上草均匀平铺（仪器视场完全被草填充）和割去草后的裸土种不同测量条件（图）。不同时间和陆表状态测量获得的发射率光谱存在差别，并且测量结果不同程度的有大气的吸收线残留。这主要是因为不同的观测时间天气状况不同导致的结果，无云的晴空特别稳定，可以精确测量大气下行辐射，进行温度发射率分离时可以避免大气影响的残留，使反演的光谱曲线比较平滑；另外不同的观测时间，陆表的含水量也存在差别；第方面就是不同的测量条件，导致观测的混合像元效应不同，从而影响发射率的测量结果。结论由 傅立叶变换红外光谱仪、红外标准板、定标黑体、辐射测温仪、接触式温度计等仪器设备第 卷第期草业学报 年组成的测量系统，可以精确的测量锡林浩特草原陆表发射率光谱。通过不同时次和背景条件下的测量，发现锡林浩特草原陆表发射率光谱具较好的时空稳定性，不同测量间的误差在 以内。对试验数据的处理分析表明，基于平滑度指数的 温度与发射率分离算法，不仅限制条件少、适用性强，而且精度较高，算法易于实现。利用本研究的发射率光谱测量结果，可以在锡林浩特草原对目前国内外主流的遥感卫星热红外通道进行在轨场地绝对辐射定标，并对热红外遥感产品进行真实性检验。参考文献：赵英时遥感应用分析原理与方法北京：科学出版社，（）：，：，：，：，：张仁华，田国良常温物体比辐射率的测量科学通报，：唐新斋，张仁华，苏红波，等比辐射率封闭法测定技术及误差分析遥感学报，（）：肖青，柳钦火热红外比辐射率光谱的野外测量方法与土壤热红外比辐射率特性研究红外与毫米波学报，（）：王先兵，杨世植，乔延利，等地物目标热红外光谱发射率的野外测量方法研究大气与环境光学学报，（）：陶伟国，徐斌，杨秀春草原产草量遥感估算方法发展趋势及影响因素草业学报，（）：池宏康，周广胜，许振柱，等草地植被盖度的近距离遥感测定草业学报，（）：徐斌，陶伟国，杨秀春，等我国退牧还草工程重点县草原植被长势遥感监测草业学报，（）：，：，：，：，：，：，：，犲 狋 犪 犾 （），（）：，：，（）：（），：程洁，肖青，柳钦火，等平滑度指数对光谱迭代平滑温度发射率分离算法精度的影响分析大气与环境光学学报，（）：犉 犻 犲 犾 犱犿 犲 犪 狊 狌 狉 犲 犿 犲 狀 狋 狊狅 犳犿 犲 犪 犱 狅 狑狊 狌 狉 犳 犪 犮 犲 犲 犿 犻 狊 狊 犻 狏 犻 狋 狔狊 狆 犲 犮 狋 狉 犪犪 狋 狋 犺 犲犡 犻 犾 犻 狀 犺 犪 狅 狋 犲犵 狉 犪 狊 狊 犾 犪 狀 犱狅 犳犆 犺 犻 狀 犪 ，（，）犃 犫 狊 狋 狉 犪 犮 狋：（），犓 犲 狔狑 狅 狉 犱 狊：；第 卷第期草业学报 年