erramodis雪被产品的分类精度对比研究.docxVIP

erramodis雪被产品的分类精度对比研究 雪是一个重要的地球表面覆盖，具有很高的不稳定性。在我国青藏高原和新疆北部,冬春季积雪经常引发灾害,严重制约草地畜牧业的可持续发展。由于积雪的重要性,准确监测其覆盖范围就成为地球系统科学的一项重要研究目标,尤其对频繁发生雪灾的牧区积雪进行深入研究,对防灾减灾具有特别重要的意义。 遥感技术是监测积雪覆盖范围最有效的方法。自20世纪80年代开始,主要使用NOAA/AVHHR卫星资料,开展积雪遥感监测。但是,NOAA/AVHRR资料在积雪研究中存在如下问题:① 广泛使用的NOAA/AVHRR数据,仅对瑞利散射和臭氧吸收作了校正,对诸如水汽、气溶胶等大气干扰所作的校正极为有限。② NOAA/AVHRR数据受天气状况的严重影响,极大地限制了在积雪监测中的应用。③ 当地表积雪厚度超过20 cm以上时,难以对雪灾评价具有重要意义的雪深这一指标的空间变化进行准确监测。1999年美国发射的Terra卫星和2002年发射的Aqua卫星携带的中分辨率成像光谱仪MODIS(Moderate Resolution Imaging Spectroradiometer)使得上述问题有了极大的改进。 自2000年2月美国航空航天局开始正式接收和分发MODIS/Terra卫星数据产品以来,对MODIS积雪产品的应用研究便成为全球广泛关注的热点问题[12,13,14,15,16,17,18,19]。但是,迄今为止针对MODIS雪被制图算法分类精度方面的研究还较少。本文结合GIS技术,通过对新疆北部地区(以下简称“北疆地区”)MOD10A1数据的处理与分析,探索了MODIS积雪制图算法的分类精度,分析了雪深和土地覆盖等因素对MODIS积雪制图算法分类精度的影响,以期为北疆牧区积雪遥感动态监测提供科学依据。 1 西部土地典型立地类型 北疆地区位于42°~50°N,79°~92°E。该区夏季干热少雨,冬季严寒多雪,1月份极端最低气温为-40 ℃,积雪期长达120天左右,积雪最大厚度60 cm,山区达1.2 m。北疆地区是我国三大积雪分布中心之一,也是新疆主要的畜牧业基地。由于区内地形复杂,气候寒冷,冬、春季雪灾频繁,严重影响着草地畜牧业的可持续发展。 2 材料和方法 2.1 资料与资料使用 利用① 遥感图像:通过美国国家冰雪数据中心网站,订购下载了覆盖北疆地区的2001年11月1日至2005年3月31日4个积雪季期间的Terra/MODIS每日积雪合成产品MOD10A1,共计585天的数字图像。覆盖北疆地区的Terra/MOD10A1雪被图像在正弦曲线地图投影(SIN: Sinusoidal Projection)坐标系统中,有H23和H24两条轨道,南北方向的行号为V04,每日涉及的MOD10A1图像有2幅,编码为h23v04和h24v04。研究区总计使用了1 170幅MOD10A1数字图像。② 气象资料:同时期的北疆地区20个地面气象观测站测量的日降雪量、积雪深度、最低和最高气温及日平均气温等资料。③ 空间数据库:主要有北疆地区行政分区、积雪观测站空间分布、数字高程模型、土地利用、草地资源类型和季节放牧利用草地分区等GIS空间数据库。 2.2 学习方法 2.2.1 图像质量评价 MODIS积雪产品采用了对地观测系统的一种分层数据格式(HDF:Hierachical Data Format)进行保存,除包含积雪分类图像外,还包含大量的图像质量评价数据及产品描述的辅助信息。本项研究利用MODIS数据处理工具软件MRT(MODIS Reprojection Tool),对MOD10A1积雪产品进行了接边处理,将正弦曲线投影转换为Albers地图投影坐标系统,椭球体选为Krasovsky,重采样方法选用最邻近法,图像文件转换为Geotiff格式。采用ArcGIS工作站版的宏语言工具,将Geotiff格式的图像文件转换为GRID格式,格网分辨率设置为500 m。 2.2.2 地面观测数据 MODIS积雪制图算法使用的是可见光和红外通道的光谱反射率数据,因此,在进行精度评价时必须选择晴空条件下的观测数据。为了自动化地判断并筛选一个地面观测台站所在地卫星经过时是否处于晴空状态,该项研究使用了如下判别规则:当地面观测台站所对应的MOD10A1图像的一个像素及其周边24个像素均无云时,该台站在卫星经过时定义为晴空状态,相应的地面观测数据可用于精度分析。据此,在585天的20个地面气象台站观测的数据中,共筛选出3 245个晴空条件下的地面积雪观测数据。在Arc-GIS工作站版GRID模块下,利用定点采样和宏语言方法,提取地面观测台站对应的585个时相的MOD10A1积雪分类编码数据,并从中选出晴空条件下的观测数据用于比较分析。 2.2.3