森林是地球上分布范围广的陆地生态系统类型之一，区域乃至全球尺度的森林地上生物量遥感估算一直是生态学研究的重要内容。森林作为陆地生态系统大的碳库和土壤碳库的重要输入源，因其巨大的碳储量和固碳能力而备受关注[1-2]。森林地上生物量的研究有助于对区域乃至全球的碳循环深入地理解和认识[3]。森林地上生物量不仅是森林固碳能力的重要标志，也是森林生态系统服务功能和生态资产评价的重要参数。利用遥感技术可以获取区域乃至全球尺度的陆表参数。光学、微波遥感数据和野外观测数据集成应用是提高区域尺度森林地上生物量估算精度的一种有方法[4-5]。Basuki等分别采用混像元分解、离散小波变换和Brovey变换将Landsat-7 ETM+和ALOS/PALSAR进行图像融来估算热带森林地上生物量双鸭山罐体保温厂家，结果表明采用光学和微波图像融的方法能够提高生物量估算的精度[4]。Treuhaft等认为基于森林结构剖面的生物量估算比微波或光学遥感的方法更精确，通过融高光谱和雷达数据所测得的叶面积指数和叶面积密度来提高生物量估算的精度[6]。
由于地形对卫星接收的反射信号的影响，经验统计模型估算的森林地上生物量精度受到很大影响。利用植被指数和生物量建立的统计关系受土壤等下垫面背景的光谱特征影响比较明显[7]。统计模型只能用于特定的研究区域和数据源，同时植被冠层的反射特也受到外界条件的影响，利用植被的光谱特征和植被结构参数建立的统计关系不适在区域尺度和多传感器之间应用。为了克服下垫面背景的空间异质对经验统计关系的影响，利用遥感物理模型估算森林生态参数是一种可行的方法[8]。遥感物理模型通过输入林分几何结构、地形特征、纯像元的光谱信息来模拟图像像元的光谱特，输出林分的生态参数。Soenen等人采用几何光学模型和查找表方法从遥感图像SPOT-5估算了加拿大卡南纳斯基斯市落基山的森林参数，研究表明几何光学模型适对复杂地形条件下的森林植被进行遥感监测研究[9]。
贺兰山是我国西部温带草原与荒漠的分界线和连接青藏高原、蒙古高原及华北植物区系的枢纽[10-11]。国内研究人员对贺兰山的森林生长状况、群落特征、生物多样在海拔梯度上的变化规律等进行了深入研究[12-15]。穆天民采用生态学统计的方法对贺兰山地区的青海云杉林地上生物量进行了研究，结果表明乔木层的生物量达到127.26 t/hm2，青海云杉林的地上生物量不亚于东北地区的云冷杉和落叶松林[16]。近年来很少见到利用遥感技术研究该地区森林分布状况、水文、生物量、碳循环、气候变化和生态资产评价等内容的相关报道。由于贺兰山受不理土地利用和气候变化的影响，生态系统其脆弱。因此，利用遥感技术对贺兰山的森林生态系统的生态参数和服务功能进行深入系统的研究尤为要。
本文选择位于干旱半干旱地区的贺兰山自然保护区的乔木林(主要为针叶林、阔叶林和混交林)为研究对象，建立基于光学遥感HJ1B和微波遥感ALOS/PALSAR等多源遥感数据的森林地上生物量反演模型，并进行森林地上生物量制图，为进一步对该区域的森林生态系统服务功能和生态资产进行综评价提供理论依据。
1 研究区与数据处理 1.1 研究区概况本文的研究区选择位于宁夏和内蒙古自治区交界的贺兰山自然保护区境内，中心坐标为38°43′47.9″N，105°55′7.4″E(图 1)。贺兰山地处内陆，属于典型的大陆气候，冬季严寒，夏季干燥炎热，降水量稀少，年平均气温8.2—8.6℃，年平均降水量250—438 mm，年平均蒸发量1600—1800 mm之间[15]。主要的乔木树种有青海云杉(Picea crassifolia)、油松(Pinus tabulaeformis)、山杨(Populus davidiana Dode)、灰榆(Ulmus glaucescens Franch)和杜松(Juniperus rigida)。灌木树种主要有小叶金露梅、蒙古扁桃、蒙古绣线菊、华西银露梅和忍冬。
1.2 遥感图像数据的获取与预处理研究中使用的图像数据包括HJ1B-CCD1多光谱图像和ALOS/PALSAR雷达图像数据。光学图像数据HJ1B—CCD1获取时间为2010年7月20日，卫星轨道17/68，图像分辨率为30 m。PALSAR是ALOS卫星上携带的一个L波段的主动式成孔径雷达传感器，工作波段为L波段。原始图像为level 1.1双化数据，化方式为HH、HV化，其中HH表示水平化(horizontal polarization)，HV表示交叉化(cross polarization)。入射角为38.722°，数据类型为单视复数据(SLC)，像元大小为9.368 m × 3.164 m，成像时间分别为2010年6月30日和2010年8月15日。
HJ1B-CCD1图像的预处理包括几何精校正、辐射定标和大气校正。利用HJ1B数据处理工具将灰度值(DN)值转为辐亮度值，利用野外采集的研究区GPS控制点在ENVI4.8中对影像进行几何精校正和大气校正，将辐亮度图像校正为反射率图像。
ALOS/PALSAR图像的预处理采用GAMMA软件进行，利用分辨率为90 m的SRTM(shuttle radar topography mission)将ALOS/PALSAR校正为像元大小为25m分辨率的图像。为了校正地形对雷达信号强度的影响，采用Zhou[17]提出的地形校正算法计算雷达成像斜距平面与地面的夹角，再进行几何平面的地形辐射纠正得到雷达地形辐射校正后的图像。后将HJ1B和ALOS/PALSAR进行了图像配准，选择了20个控制点，配准精度在0.2个像元以内，均方根误差RMSE为0.1。
1.3 遥感图像融分类为了得到较高精度的分类图像，先对配准后的HJ1B图像(图 2a)和ALOS/PALSAR SAR图像进行了图像融处理，然后再进行图像分类。采用离散小波变换将预处理后的2010年6月30日HH化SAR图像(图 2b)、HJ1B反射率图像、主成分变换图像和归一化植被指数(NDVI)图像进行了图像融处理。
决策树分类法已经广泛应用于遥感影像信息提取和土地利用/土地覆被分类研究中[18]。采用决策树分类将融后的图像进行分类处理，分类结果如图 2c所示。图像共分为10类，分别为建筑、农田、针叶林、灌木林、阔叶林、混交林、沙地、草地、水体和道路。分类的总体精度为85.60%，Kappa系数为0.8269。
1.4 野外数据采集与处理2012年7月13日—7月30日在研究区开展了野外实地考察和森林参数测量。在研究区内布设样地，选择地形坡度小、植被分布相对比较均匀，布设了50个有代表的30 m×30 m标准样地(针叶林29个，平均树龄105a，郁闭度40%—80%；阔叶林2个，平均树龄38a，郁闭度30%—60%；混交林19个，平均树龄73a，郁闭度40%—70%)，用手持GPS记录每个样地的四角坐标，选择胸径在5 cm以上的林木实行每木检尺，测量胸径和树高，随机选择几株林木测量枝下高，东西和南北两个方向的冠幅。在室内对测量数据进行处理。
1.5 样方生物量计算森林生物量模型的建立采用Kittredge提出的CAR(Constant Allometric Ratio)算法[19]，其数学表达式为：Y=aXb。穆天民[16]、刘斌[20]、梁咏亮[21]等人分别对贺兰山的青海云杉、天然油松林、灰榆疏林等的生物量模型进行了研究，建立了以胸径(DBH)、树高(H)、胸径平方乘树高(D2H)等为自变量的生物量幂指数回归模型，根据这些数学模型计算得到单株树的实测生物量。样地内所有的林木生物量累加得到各样地生物量总和。生物量总和除以样地的面积得到林分地上生物量(t/hm2)，计算结果见表 1。
2 反演方法 2.1 GOMS-MFM反演算法Li-Strahler几何光学模型(geometric-optical mutual shadowing双鸭山罐体保温厂家，GOMS)[22]主要用于遥感影像像元尺度的森林结构参数的反演。多元前向模式(multiple forward mode，MFM)是利用几何光学模型和查找表(LUT)反演生态结构参数的常用方法，它由遥感物理模型正向模拟计算建立参数查找表数据库[23-24]。利用MFM算法进行参数反演，需要从查找表里搜索和图像像元反射率比较接近的值作为模型的优解。为了得到比较好的反演结果，需要建立比较大的查找表[25]。国外一些学者利用几何光学模型和MFM算法从多光谱图像数据反演了森林结构参数[8-9, 26-27]。
GOMS模型的输入参数包括树的结构、地形和光谱信息(表 2)。树的结构参数包括：树密度(λ，trees/m2)、树冠短半轴(r，m)、树冠长半轴(b，m)、树冠中心高度(h，m)、树高分布(dh，m)。根据样地实测数据统计结果输入树的结构参数的小值、大值和增加值。地形信息包括坡度(°)和坡向(°)，利用分辨率为30 m的ASTER GDEM数据提取研究区的地形坡度和坡向图，统计坡度和坡向的变化范围。模型的光谱信息包括：光照冠层、光照背景和阴影等纯像元组分端元的光谱信息。采用ENVI软件纯像元指数(pixel purity index，PPI)从HJ1B反射率图像提取3种森林类型(阔叶林、针叶林和混交林)的纯像元的组分端元光谱信息。根据GOMS模型的输入格式要求，将3种森林类型输入参数的变化范围和增加值在表 2中分别列出。
Soenen等人以树冠表面积(SA)和测量的单木的生物量(B，t/株)建立了单木生物量幂指数生长方程来计算林分地上生物量[9]。为了利用遥感图像估算像元尺度的森林地上生物量，根据野外测量的3种森林类型(阔叶林、针叶林和混交林)结构参数数据，以实测树高H为自变量建立3种林型单木生物量和树高H之间的幂指数函数回归模型如(1) 式所示：
(1)式中,B为单株树地上生物量(kg/株)，H为树高(H，m)，b0和b1为回归方程的系数。由公式(1) 计算得到单株树的地上生物量。利用GOMS模型模拟计算建立参数查找表时，根据树密度参数λ和公式(1) 建立的树高和生物量回归模型可以得到像元尺度的森林地上生物量。
利用HJ1B图像和GOMS模型反演森林参数的过程包括两步：(1) 光谱模拟和林分地上生物量数据库建立。由公式(1) 建立的生物量模型计算得到单木的地上生物量B(B，kg/株)，像元尺度的森林地上生物量AGB(t/hm2)为：AGB=λB，λ为树密度。GOMS模型正向模拟输出4个波段反射率和生物量信息，然后建立参数查找表文件。(2) 模型解的查找。采用光谱距离的均方根误差(RMSE)优化技术从查找表中搜索图像中每个像元4个波段反射率和模型模拟的对应波段反射率接近的值，作为模型的优解。通过查找模型的优解便得到每个像素的生物量。RMSE的计算公式如式(2) 所示[28-29]：
(2)式中,ρi表示图像的反射率，ρ^i表示模拟的反射率，nb为图像的波段个数。
2.2 线光谱混分解线光谱混分解模型(linear spectral mixing model，LSMM)是混像元分解(spectral mixture analysis，SMA)的常用方法，它假设图像中每个像元的反射率是由该像元中每种地物(纯像元)反射率以及它们的地物类型在该像元内所占像元面积的比例作为权重系数的线组[30-31]。它的数学表达式如式(3，4) ：
(3) (4)式中,Ri,j是i个纯像元、j个波段的反射率，Ki是i个纯像元的面积百分比，vj是j个波段的误差。m为组分的个数，p为图像波段的个数。
2.3 混像元分解SMA模型利用地形校正模型和ASTER GDEM数据对配准后的HJ1B反射率图像进行地形校正，消除了地形变形对遥感图像反射率的影响，将坡地光谱信号校正为水平地面信号[32]。对地形校正后的图像HJ1B进行线光谱混分解，分解后得到植被组分图像veg(图 3a)。分别将植被组分图像(veg)和雷达图像ALOS/PALSAR的HH和HV化的雷达后向散射系数图像进行离散小波融。利用S-Plus软件的线回归将样地实测生物量和不同的融图像(图 3b)对应的图像像元值建立了回归模型，从建立的回归模型中选择精度好的模型，作为后的回归模型如式(5) ：
(5)式中,veg_hh表示HJ1B混像元分解后的植被组分图像veg和2010年6月30日的PALSAR的HH化融后的图像。决定系数R2=0.61，均方根误差RMSE=17.11 t/hm2，F值为20.0752，P＜0.01。
2.4 PALSAR多元逐步回归模型采用S-Plus软件的多元逐步回归算法，铝皮保温将两个时相雷达图像HH和HV化雷达后向散射系数和样地实测生物量进行多元逐步回归分析，建立的地上生物量模型如式(6) ：
(6)式中,σHV0630、σHH0815分别表示6月30日获取图像的HV化雷达后向散射系数和8月15日获取图像的HH化雷达后向散射系数。决定系数R2=0.63，均方根误差RMSE=12.94 t/hm2，F值为21.62，P＜0.01。
3 结果与分析 3.1 3种森林类型地上生物量模型根据3种不同的森林类型，分别从实测数据中选择200株树，采用公式(1) 计算像元尺度地上生物量，拟3种森林(阔叶林、针叶林和混交林)地上生物量模型的回归系数、R2、RMSE(kg/株)和F值如表 3所示(P＜0.01) 。利用该生物量回归模型模拟计算得到不同森林类型的林分地上生物量(t/hm2)。
邮箱：215114768@qq.com 3.2 生物量制图利用本文介绍的3种算法对贺兰山研究区的森林进行生物量制图，研究区的森林地上生物量分布制图分别如图 4所示。根据遥感图像分类图(图 2c)，将分类图上非森林覆盖区域(包括建筑物、道路、草地、灌木)并为一类，定义为非林地。图 4a为利用GOMS模型和MFM方法估算的生物量分布图，其中非林地约占27.57%，乔木林地面积占整个图像区域的72.43%，乔木林大生物量为332.87 t/hm2。生物量在范围0—50、50—100、100—150、150—200、200—250、>250 t/hm2，分别占总图像覆盖区域面积的28.54%、19.49%、12.07%、6.17%、3.03%、3.12%，其中0—150 t/hm2约占60.11%，>150 t/hm2约占12.32%。图 4b为利用SMA模型估算的生物量分布图，其中非林地约占29.39%，乔木林大生物量为241.93 t/hm2。生物量在范围0—50、50—100、100—150、150—200、200—250 t/hm2，分别占总图像覆盖区域面积的12.00%、29.79%、24.97%、3.70%、0.14%，其中0—150 t/hm2约占66.76%，>150 t/hm2约占3.84%。图 4c为利用PALSAR回归模型估算的生物量分布图，其中非林地约占34.25%，乔木林大生物量为177.35t/hm2。生物量在范围0—50、50—100、100—150、150—200 t/hm2，分别占总图像覆盖区域面积的24.35%、32.50%、8.64%、0.25%，其中0—150 t/hm2约占65.49%，>150 t/hm2约占0.25%。
由3种算法的统计结果可以看到：3种算法估算的生物量主要分布在范围0—150 t/hm2，不同方法估算的生物量结果之间存在一定的差异。采用混像元分解和PALSAR雷达后向散射系数图像估算的生物量大于150 t/hm2时达到饱和，建立的回归模型无法正确反演，反演的结果存在较大的误差[6]。SMA方法将PALSAR和光学图像进行融处理减少了生物量估算的误差。
3.3 估算误差分析为了进一步对3种方法估算的生物量进行误差分析，根据野外实测生物量数据，将MFM、混像元分解SMA和PALSAR多元逐步回归等3种方法估算的森林地上生物量进行分析。将模型估算地上生物量和实测生物量的差值作为模型估算误差。选择其余的20个实测样地(plots)生物量数据作为真值，计算它们的平均值和标准差(SD)。把实测值落在区间avg-SD和avg+SD内的值记为plots ＜1 SD，反之则记为plots > 1 SD。根据这个划分标准，这20个实测生物量值分成两组，其中plots ＜1 SD共有11个，plots > 1 SD共有9个数据。然后分别计算实测生物量、模型估算的生物量、模型估算生物量误差、模型估算生物量plots ＜1 SD误差(11个数据)和plots > 1 SD(9个数据)误差的平均值、标准差、小值、小值和变化范围，统计结果如表 4和图 5所示。
根据表 4和图 5的误差统计结果分析可以看到：MFM、SMA和PALSAR多元逐步回归等3种方法的平均误差分别为8.0、8.2、15.1 t/hm2，MFM和SMA估算误差比较接近，MFM误差小；SMA法估算生物量的小值、大值和变化范围比实测生物量略宽，而MFM法估算生物量变化范围比实测值变化范围窄，PALSAR回归模型估算生物量变化范围宽；MFM法和SMA法估算生物量的平均值和实测结果比较接近，而PALSAR回归法估算生物量平均值比实测结果小；MFM法和SMA法估算生物量标准差比实测结果略低，而PALSAR回归模型估算生物量标准差要大于实测结果；11个数据误差(plots ＜1 SD)，MFM法和SMA法的平均误差和标准差都接近总误差(20个数据)，而PALSAR回归比总误差(20个数据)大；9个数据误差(plots >1 SD)，与＜plots 1 SD情况正好相反。
由以上的分析可以看出，采用经验统计和基于遥感物理模型的反演方法从遥感图像数据估算森林地上生物量，不同的反演算法估算的森林地上生物量的精度有一定的差异，原因主要有以下几点：(1) 森林下垫面背景对卫星图像信息有一定的影响。数学模型很难真实的描模拟森林的三维结构及下垫面背景，会造成一定的反演误差。(2) 地形对反演精度的影响。地形改变了太阳、地表和卫星传感器三者所构成的相对观测几何变化，从而造成方向反射的变化[33]，地形对遥感图像的质量和反演精度的影响比较大。(3) 野外实测数据误差的影响。由于地形和林下植被的影响，在野外很难保证实测样地大小为标准样地，同时卫星接收的图像数据受地形坡度和坡向的影响，野外样地的位置很难与图像完全匹配，这对反演的验证精度有一定的影响。
3.4 估算结果验证为了检验参数反演过程的过拟现象，有要结野外实测数据对模型进行验证。选择20个样地实测生物量和3种方法估算的生物量进行了对比，用实测生物量值验证模型，结果见图 6。图 6分别验证了MFM和SMA法估算生物量结果，R2分别为0.61和0.60，实测值与估算值之间的RMSE分别为8.33、9.417 t/hm2，图中各点均分布在1∶1线附近。图 6还验证了PALSAR回归法估算生物量结果，R2=0.39，RMSE=14.89 t/hm2。
3种模型估算的生物量和实测生物量的误差变化范围和模型的验证结果进一步表明：MFM和SMA两种方法估算的生物量结果比较接近实测值，而PALSAR多元逐步回归模型估算的生物量精度差。由此可见，采用PALSAR多元回归模型估算的山地森林生物量与实测的生物量相差比较大，精度比较低，决定系数小于0.5，该方法估算的生物量精度明显低于其它两种方法估算的生物量。这表明采用几何光学模型和混像元分解SMA更适估算山地森林地上生物量。
3.5 讨论植被的生态结构参数不能通过卫星图像数据直接获取，需要利用经验统计或遥感物理模型反演的方法从卫星图像估算得到。利用Li-Strahler几何光学模型反演森林生物量不需要较多的地面验证数据，模型的输入参数不需要每棵树的详细信息。利用几何光学模型模拟建立查找表时，根据野外直接测量的数据，如树高(H)、胸径(DBH)、胸径平方乘树高(D2H)等数据与生物量建立回归模型，然后通过中间变量(胸径或树高)从遥感图像上估算森林地上生物量。
由于地形引起雷达成像几何关系的变化，进而导致雷达后向散射系数的变化。利用雷达图像PALSAR估算森林结构参数时，地形对估算精度影响较大。因此，采用雷达图像估算森林参数精度较低，反演的结果存在较大的误差。结雷达图像和光学图像能够减少生物量估算的误差，改善地形对森林参数估算的精度的影响。采用离散小波变换将混像元分解后的组分图像和雷达图像进行融，然后建立生物量多元逐步回归模型，该模型比利用雷达后向散射系数建立的生物量模型估算精度更高，这是由于组分图像改善了地形对雷达后向散射系数的影响。因此，利用混像元分解和离散小波融可以有的提高森林生物量的估算精度。
采用光谱距离的均方根误差(RMSE)小从查找表中搜索模型的优解，由于测量误差和模型模拟数据的问题，方程的解不唯一或无解，从而导致方程的病态问题[34-35]，在实际应用中可以结一些其它的植被指数如归一化植被指数NDVI或比值植被指数(SR)来查找方程解。采用查找表进行查找时，如果图像的覆盖范围比较大，同时建立的查找表数据库比较大时，在查找过程中查找速度会比较慢，需要对查找表结构和查找算法进行优化。
4 结论本文探讨了利用GOMS几何光学模型和经验统计模型估算山地森林地上生物量的方法。利用几何光学模型模拟时将可直接测量的参数，如胸径(DBH)或树高(H)作为中间变量，建立树高和森林地上生物量回归模型，计算像元尺度的林分地上生物量，通过MFM方法从卫星图像HJ1B数据反演研究区的森林地上生物量，同时将反演结果和经验统计模型方法进行了对比。研究结果表明，MFM方法比经验统计模型估算森林地上生物量精度更高。采用几何光学模型反演的生物量精度好，利用混像元分解和图像融可以改善地形对森林参数估算的精度的影响，利用PALSAR图像估算的森林生物量的精度差。因此，几何光学模型和混像元分解模型适估算森林地上生物量。
利用几何光学模型进行前向模拟的输入参数不需要较多的先验信息和野外样地调查数据，该方法更适于通过多时相多传感器的遥感图像反演区域尺度的植被生态结构参数。因此，本文介绍的方法进行森林地上生物量遥感监测研究具有一定的应用潜力。然而，本研究所用的野外观测样方的数量相对偏少，采用统计方法所建立的生物量模型的适用需要进行一步验证。在今后的研究中，可以考虑通过增加野外观测的数量和小样本回归的方法，比如机器学习算法建立回归模型双鸭山罐体保温厂家，并对模型的适用进行进一步的验证。
致谢: 中国林业科学研究院资源信息研究所谭炳香老师帮助写作，特此致谢。 相关词条:离心玻璃棉隔热条PA66生产设备
钢绞线厂家