热带人工林生物物理参数及生物量对RADARSAT SAR 信号响应研究.docVIP

热带人工林生物物理参数及生物量 对RADARSAT SAR信号响应研究 王臣立1，2，郭治兴3，牛 铮1，林文鹏1 中国科学院遥感应用研究所遥感科学国家重点实验室，北京 100101；2. 中国科学院研究生院，北京 100039； 3. 广东省生态环境与土壤研究所，广东 广州 510650 摘要：合成孔径侧视雷达穿透冠层而获得树干信息的特点使得其具有在南方NDVI易饱和地区的生物量估测方面优于光学遥感。文章利用地面实测调查数据，首次系统地研究了热带人工林生物物理参数及生物量对RADARSAT-SAR信号响应。结果表明：RADARSAT-SAR后向散射系数与森林生物量、树高、胸径可以用对数模型模拟。对于本实验区人工速生林来说，在森林生物量小于40 t·hm-2时RADARSAT-SAR后向散射系数与森林生物量具有较强相关性。桉树树高与后向散射系数的相关系数大于桉树胸径与后向散射系数关系。松树胸径与后向散射系数相关系数高于树高与后向散射系数的相关系数。研究结果为SAR在热带森林研究中的应用提供了一定的理论基础。 关键词：RADARSAT-SAR；生物量；生物物理参数；热带人工林；后向散射系数 中图分类号：TP79∶S771.8 文献标识码：A 文章编号：1672-2175（2006）01-0115-05地球陆地表面上空的生物量的总数及其分布是我们了解全球碳循环以及它在将来如何变化不确定因素之一。合成孔径侧视雷达能够穿透冠层而获得树干信息的特点使得其具有在南方多云雨地区的生物量估测方面优于光学遥感。利用雷达的穿透特性，可以准确实时的获取全球森林变化,并与光学遥感相互补充。为我们理解全球变化提供可靠的数据。 表1 RADASAT-SAR数据产品参数 Table 1 Radarsat-sar parameters 波段 极化 象元大小 /m 覆盖范围 入射角中心点经度中心点纬度接收日期35.12 110.81 21.34 2002-09-17 RADARSAT-SAR是加拿大于95年11月4日发射的适用于全球环境和土地利用、自然资源监测等。研究区概况 试验林地在广东雷州林业局，位于东经l11°，北纬21，年均温度235。，1月气温最低，均平152。。终年无霜。年降雨量1884毫米，集中于5—9月。6一8月为台风季节。equation of biomass for Eucalyprus urophylla and Pinus caribaea 树种及器官 回归方程 r 树种及器官 回归方程 r 加勒比松 尾叶桉 树干 logW=0.9797 log(D2H)－1.9556 0.997 干材 W=0D2H）0.9768058 0.9972 树枝 logW=1.0919 log(D2H)－2.7325 0.9871 干皮 W=0.0063525（D2H）0.8738162 0.9902 树叶 logW=0.7454 log(D2H)－2.0235 0.9995 树枝 W=0D2H）0.8449044 0.8973 树叶 W=0.0D2H）0.933545 0.9334 将样地林木各径阶的平均胸径和树高代入表2各生物量方程，将各器官生物量累加，得到各径阶 及林分生物量。 1.2 数据处理 1.2.1 图像的斑点躁声去除和几何纠正 选定小波变换做为最佳滤波器进行斑点躁声去除。此工作在PCI中进行。 首先采用EDARSE 软件的雷达模块对RADAR- SAT-SAR进行几何纠正，需要DEM数据，然后再通过控制点拟合3次多项式对雷达图像作进一步的几何校正，以消除误差影响。 1.2.2 后向散射系数计算 公式如下：首先计算雷达象元亮度值： ?oj（dB）= 10*log10[(DNj2 + A3)/A2j] （1） 其中，?oj是第J象元的亮度值，DNj是雷达图象j象元的的灰度值。A2j是j象元的增益，A3是偏移。A3，A2j是与雷达产品的查找表有关的常数，可以在RADARSAT-SAR 数据产品的查找表中找到。再将亮度值转为后向散射系数。 ?oj (dB) = ?oj + 10*log10(sinIj) （2） 其中Ij是J象元的瞬时角，?oj j为后向散射系数。 1.2.3 每个小班平均后向散射系数的计算 在GIS支持下，林场小班分布图与RADARSAT-SAR图像投影转换为统一的UTM投影，林场矢量图与影象图叠加，计算每个小班的象元数和平均后向散射系数，此工作在PCI中进行。通过GIS将空间数据和地面属性数据进行连接。然后通过统计模型对每个小班的平均后向散射系数与林场各个小班的树高（H）、胸径（D）、