基于alospalsar数据的红树林hh、hv极化后向散射特征提取.docxVIP

基于alospalsar数据的红树林hh、hv极化后向散射特征提取 1 林木资源监测的重要性 红树林分布在位于亚热带海岸和陆地交汇处的河口区。它是世界上最生产力的生态系统之一。这是湿地和国际生物多样性保护的重要组成部分。红树林湿地不仅具有促淤沉积、改善环境、护堤防波等生态功能,而且还为许多海洋动物提供重要的栖息地和食物。随着经济的发展,由于围海造地、围海养殖、砍伐等人为因素,再加上工农业废水对近海污染的日益严重,红树林资源遭受严重破坏,面积日趋减少,生态功能退化。因此,及时获取红树林生长面积的变化成为加强红树林保护的迫切需求。 由于红树林生长于近海滩涂上,红树林内泥、沙松软,且植被郁密度高,很难进入其中进行实地监测。在全球范围内,人们广泛应用遥感技术进行红树林资源监测,大部分的研究是利用光学遥感数据开展红树林资源监测。随着雷达技术的发展,人们日益认识到雷达相较于光学遥感在红树林监测上所具有的优势:雷达的穿云透雾特点使其具有全天时、全天候的观测能力,这对于生长于多云多雨热带亚热带的红树林,监测上不受时间和气象的限制,特别是更有利于获得较低潮位的数据;L波段雷达波长较长,能反映冠层以及冠层下的枝干信息;雷达的多极化方式,能获取红树林更多的信息。 人们最早应用机载雷达系统进行红树林监测,通过机载的多波段多极化的SAR数据,对红树林的生物量等进行定量分析。星载雷达系统的发展,为定期监测红树林资源提供了更好的数据源,Hervé Trebosse等利用ERS 数据(C波段)更新快的特点,监测红树林的海岸演变情况,并利用监测结果更新沿海海图。随着ENVISAT ASAR和ALOS PALSAR 雷达传感器的相续投入使用,其所具有的多极化、多入射角的重复观测能力,为更好地研究红树林提供了可能。 2006年1月成功发射的先进对地观测卫星ALOS是JERS-1的后继星,其上搭载的高性能PALSAR雷达传感器,具有更高的空间分辨率和多种极化方式。本文利用多时相的ALOS PALSAR数据,进行漳江口红树林资源的提取研究。 2 生态保护区概况 福建漳江口红树林国家级自然保护区位于福建省漳州市云霄县漳江入海口(图1),地理位置为117°24′07″～117°30′00″E,23°53′45″～23°56′00″N,保护区总面积2 360 hm2,为福建省最重要的湿地生态系统类型的国家级自然保护区。该保护区的红树林植被类型主要有秋茄林和白骨壤林两个群系,滨海盐沼为天然互花米草盐沼群系,红树林高2～3 m,郁密度80%～90%。 3 多视数据处理 研究所用的PALSAR FDB双极化Level 1A数据,包括2007年7月13日、8月28日、10月13日3个时相的HH、HV双极化共6个通道数据。包括原始前视、星下点和后视图像文件,分别附带独立的辐射定标和几何定标参数文件。数据处理流程如图2所示。 从Level 1A数据中提取出的为单视复数(SLC)数据,它包含幅度信息A和相位信息?。本研究应用的是后向散射能量信息,需将幅度和相位信息转换为强度信息。用幅度信息A和相位信息?记录的微波信息服从方程(1)的分布。 Pz1,z2(z1,z2)=1πσexp(?z21+z22σ)(1)Ρz1,z2(z1,z2)=1πσexp(-z12+z22σ)(1) 其中:z1=Acos?,z2=Asin?,Pz1,z2为随机变量为z1、z2的概率密度函数。估计每个微波分辨单元的后向散射强度信息,在同相信息z1=Acos?和正交信息z2=Asin?已知条件下,代入方程(2)便得到每个分辨单元的最大似然估计平均强度I。 I=z2112+z2222=(Acos?)2+(Asin?)2=A2(2) 数据转为强度信息后,对图像进行多视处理。进行多视处理,会增强辐射分辨率,降低斑点噪声水平,同时也会降低空间分辨率,但有利于识别出更多的空间分布特征及地物特征信息。为保持较高的图像分辨率,设置距离向和方位相视数分别为1和3,距离向和方位向分辨率约为9.4 m*9.5 m。 多视处理后,对图像进行辐射定标,使不同景数据以及同一景数据不同区域的辐射特征具有可比性,辐射定标运用方程(3)和方程(4)进行。 I0=z2112+z2222×10-11.5(3) σ0=10.0×log10z2112+z2222-115.0 (4) 其中:I0是辐射定标后向散射强度,σ0是以dB表示的后向散射系数,代表空间平均。 将由距离和方位表示的雷达坐标图像转换为UTM坐标图像,单景雷达图像的处理就完成了。要使多时相图像具有可比性从而进行分析,就必须对不同时相的图像进行高精度配准,以一景图像作为基准图像,运用相关匹配技术,将其它图像以基准图像为参考图像进行匹配。配准后的图像误差控制在1个像元大小以内。 噪声