冰川冻土, 2022, 44(1): 316-326 doi: 10.7522/j.issn.1000-0240.2022.0039

冰冻圈技术

基于AVHRR影像的北半球积雪识别算法

王京达,1,2, 郝晓华,2, 和栋材1, 王建2,3, 李弘毅2, 赵琴2

1.太原理工大学 矿业工程学院, 山西 太原 030024

2.中国科学院 西北生态环境资源研究院, 甘肃 兰州 730000

3.江苏省地理信息资源开发与利用协同创新中心, 江苏 南京 210023

Snow discrimination algorithm in the Northern Hemisphere based on AVHRR image

WANG Jingda,1,2, HAO Xiaohua,2, HE Dongcai1, WANG Jian2,3, LI Hongyi2, ZHAO Qin2

1.College of Mining Engineering，Taiyuan University of Technology，Taiyuan 030024，China

2.Northwest Institute of Eco-Environment and Resources，Chinese Academy of Sciences，Lanzhou 730000，China

3.Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application，Nanjing 210023，China

通讯作者: 郝晓华，研究员，主要从事积雪遥感、积雪与气候变化研究. E-mail： haoxh@lzb.ac.cn

收稿日期: 2021-11-21   修回日期: 2022-01-18  

基金资助: 国家自然科学基金项目.  41971325.  42171391.  41971399
高分辨率对地观测系统国家重大专项.  21-Y20B01-9001-19/22

Received: 2021-11-21   Revised: 2022-01-18  

作者简介 About authors

王京达，硕士研究生，主要从事积雪遥感研究.E-mail：wangjingda0591@link.tyut.edu.cn , E-mail：wangjingda0591@link.tyut.edu.cn

摘要

针对2000年前北半球较高时空分辨率和高精度的历史积雪范围数据缺失问题，利用NOAA-AVHRR地表反射率数据，以Landsat-5 TM生成的积雪范围影像作为参考真值，优化基于多指标的多级决策树积雪识别算法的阈值，并结合云雪混淆区分技术，生成了北半球AVHRR 1981—1999年L1级逐日积雪范围数据集。此外，针对AVHRR在高纬度地区数据完全缺失和低纬度地区数据部分缺失问题，利用微波雪深数据集进行填充，生成了北半球L2级逐日积雪范围数据集。最后，利用北半球1981—1999年间2 546个气象台站记录的雪深数据和939景Landsat-5 TM参考积雪范围影像作为验证数据，对AVHRR积雪范围数据集进行了精度验证。结果表明：L1级和L2级数据集的总体精度分别为81.8%和82.2%，用户精度分别为83.7%和83.8%，生产者精度分别为81.7%和84.2%，说明算法精度较高，错分误差和漏分误差均比较均衡。进一步利用Landsat-5 TM参考积雪范围影像对L2级数据集进行面上精度评估，发现L2级数据集的总体精度为90.3%，用户精度为90.2%，生产者精度为99.1%，L2级数据集精度较高。生成的北半球历史数据集可为全球积雪变化研究提供有效数据补充。

关键词： AVHRR ; 北半球 ; 积雪识别 ; 积雪范围 ; NDSI

Abstract

Due to the lack of snow cover extent dataset with slightly high spatiotemporal resolution and high accuracy in the Northern Hemisphere from 1981 to 1999， we used NOAA-AVHRR version 4 surface reflectance data as the basic input data and Landsat-5 TM snow maps as reference maps， to obtain the optimal threshold of snow discriminant algorithm based on multi-indicator multi-level decision tree. Then， combined with the cloud discrimination algorithm， we produced the Northern Hemisphere AVHRR L1 （Level 1） daily snow cover extent dataset. In addition， aimed at the complete lack of AVHRR dataset in high latitudes and partial lack of data in low latitudes， we filled the gaps of AVHRR L1 dataset by the Northern Hemisphere 0.25° snow depth dataset， and generated AVHRR L2 （Level 2） daily snow cover extent dataset. Finally， taking 2 546 ground snow depth stations measurements from 1981 to 1999 and 939 Landsat-5 TM snow maps as validation data， accuracy of the Northern Hemisphere AVHRR daily snow cover dataset was accessed. The results showed that overall accuracy （OA） of AVHRR L1 and L2 dataset is 81.8% and 82.2%， user’s accuracy （UA） is 83.7% and 83.8%， and producer’s accuracy （PA） is 81.7% and 84.2%， respectively. These accuracies were relatively high， and commission error and omission error were relatively balanced. Furthermore， we used Landsat-5 TM snow maps for L2 dataset to perform accuracy assessment. The results showed that OA of L2 dataset is 90.3%， UA is 90.2% and PA is 99.1%. Therefore， our product has virtually provided more reliable snow knowledge over the Northern Hemisphere， and thereby can better serve for study on global snow cover change.

Keywords： AVHRR ; Northern Hemisphere ; snow discrimination ; snow cover extent ; NDSI

PDF (4463KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

王京达, 郝晓华, 和栋材, 王建, 李弘毅, 赵琴. 基于AVHRR影像的北半球积雪识别算法[J]. 冰川冻土, 2022, 44(1): 316-326 doi:10.7522/j.issn.1000-0240.2022.0039

WANG Jingda, HAO Xiaohua, HE Dongcai, WANG Jian, LI Hongyi, ZHAO Qin. Snow discrimination algorithm in the Northern Hemisphere based on AVHRR image[J]. Journal of Glaciology and Geocryology, 2022, 44(1): 316-326 doi:10.7522/j.issn.1000-0240.2022.0039

0 引言

积雪作为冰冻圈的重要组成部分，在地球表面分布较广，大约98%的季节性积雪位于北半球，面积最大可达4.7×107 km2［1］，且变化较为活跃，是影响气候变化的重要因子［2-4］。在全球水循环中，积雪的积累和消融过程调节着水的年内再分配，是干旱、半干旱地区春季最重要的淡水资源。积雪是气候变化的指示器［5-6］，积雪范围、积雪深度及雪水当量等信息是水文模型和气候模式的重要输入参数［7-8］，全球尺度下的积雪参数监测对于气象、水循环、气候变化以及能量平衡领域的研究具有重要的意义［9］。因此，准确获取北半球的积雪信息可为寒区水文循环、水资源管理、雪灾预警等重大研究提供必要可靠的参考资料［10］。

北半球积雪范围和气候变化研究需要长时间序列、高时空分辨率、高精度的积雪范围数据。地面观测积雪范围资料是目前研究积雪长期变化特征最广泛、最可靠的资料，具有时间尺度长、连续性好的优点，但地面观测资料无法准确获取积雪范围［11-12］。可用于北半球研究的积雪范围产品有微波产品和光学产品，微波积雪范围产品主要来源于SMMR、SSM/I［13-14］和AMSR-E［15-16］传感器。微波系列产品时间序列较长，但空间分辨率较低，目前可用的长时间序列微波积雪范围产品的空间分辨率为25 km。光学积雪范围产品主要来源于AVHRR［17］传感器、MODIS［18-19］传感器，以及Landsat系列卫星［20］和哨兵（Sentinel）系列卫星［21］。MODIS积雪范围产品（2000年至今）时空分辨率较高，精度也较高，广泛用于全球的水文和气候变化研究，但对于气候变化研究来说，时间序列仅有约20年，较短的时间序列无法满足研究；Landsat时间序列长，空间分辨率高，但时间分辨率低（16 d），无法准确监测季节性积雪的变化；Sentinel系列卫星虽然空间分辨率高，但时间序列短，时间分辨率较低，也无法满足气候变化研究需求。AVHRR产品具有时间序列长（1981年至今）、时空分辨率较高的特点（时间分辨率为1 d，空间分辨率为5 km），是研究北半球气候变化的理想数据。

一些学者利用AVHRR发展了积雪范围产品。Hüsler等［22］利用1 km分辨率的AVHRR原始数据集，开发了一种改进的积雪范围检索算法来识别欧洲阿尔卑斯山的积雪范围；Zhou等［23］利用决策树算法基于1 km分辨率的AVHRR LAC数据和HRPT数据，以中亚地区为例开发了一种长时间序列积雪范围的识别算法。但此类研究主要集中在小区域的积雪识别，算法不具有全球普适性。日本宇航局（Japan Aerospace Exploration Agency， JAXA）利用AVHRR和MODIS结合发展的北半球1978—2015年5 km分辨率的JASMES（JAXA Satellite Monitoring for Environmental Studies）积雪范围产品，弥补了大尺度长时间序列逐日积雪范围产品数据的缺失，但是该产品经验证发现在青藏高原精度较低［24］。Hao等［25］在中国区域发展了多级决策树积雪识别算法，并生成了一套长时间序列积雪范围产品，但算法中的阈值对于全球产品是否适用仍需要探索。

本研究的目的是利用AVHRR数据，获取多级决策树积雪识别算法中适合北半球的阈值，生成1981—2000年北半球历史积雪范围数据集，从而为气候变化研究提供更长时间序列、更高精度的逐日积雪范围数据集。研究主要解决两方面问题：①如何通过有限的6个波段，结合研究区给定合理阈值；②如何填补AVHRR在高纬度地区数据完全缺失和低纬度地区数据部分缺失问题，即如何填补数据空缺值。主要研究内容包括：①以高空间分辨率的无云Landsat-5 TM参考积雪范围影像作为参考真值，获取北半球积雪识别算法合理阈值；②填充积雪范围数据中的空缺值，制作北半球AVHRR逐日积雪范围数据集；③利用北半球气象台站的雪深数据和Landsat-5 TM作为参考真值，验证数据集精度。该数据集将为北半球积雪变化研究提供重要支撑。

1 数据来源

1.1　AVHRR地表反射率数据

NOAA-AVHRR地表反射率数据第四版（Surface Reflectance data Version 4， SR V4）被用作本研究的主要数据，来源于美国国家环境信息中心（https：//www.ncdc.noaa.gov/），它是由NOAA GAC（Global Area Coverage）L1数据经过一系列的定标、辐射校正、大气校正及几何校正后得到的地表反射率数据和亮温数据。该数据时间范围为1981年6月24日至2019年5月16日，时间分辨率为1 d，空间分辨率为5 km。本研究使用1981年6月24日至2000年2月27日时间范围内的数据，数据被转换投影坐标为椭球体 WGS84 Albers 投影坐标，包含各个波段的地表反射率、亮温和质量控制标记（quality control bit flags assessment， QA）等，具体如表1所示。

表1   NOAA-AVHRR地表反射率数据

Table 1  NOAA-AVHRR surface reflectivity data

数据缩写波长/μm中心波段/μm描述类型
SREFL_CH1Ref10.58~0.680.64地表反射率
SREFL_CH2Ref20.725~1.0000.86地表反射率
SREFL_CH3Ref31.58~1.643.75地表反射率
BT_CH3BT373.55~3.933.75亮温
BT_CH4BT1110.3~11.311.0亮温
BT_CH5BT1211.5~12.512.0亮温
quality control bit flags assessmentQA质量控制标记

新窗口打开| 下载CSV


1.2　Landsat-5 TM数据

Landsat-5 TM专题成像仪（thematic mapper， TM）是搭载在Landsat-5陆地观测卫星的传感器，它于1984年3月发射成功直到2013年6月退役。本研究通过SNOWMAP算法［18，26］获取Landsat-5 TM二值积雪范围，即雪和非雪。SNOWMAP算法的核心是通过归一化雪被指数（normalized difference snow index， NDSI）区分积雪和非雪，计算公式为

NDSIL=Band2-Band5Band2+Band5

式中：NDSIL为Landsat-5 TM的归一化雪被指数；Band2为Landsat-5 TM影像的第二波段（0.52~0.60 μm）；Band5为第五波段（1.55~1.75 μm）。算法设定NDSIL≥0.4来区分雪和非雪，同时为了排除水体的干扰，设定第四波段（0.76~0.90 μm）的阈值大于0.11。利用该算法获取Landsat-5 TM空间分辨率为30 m的二值积雪范围后，通过最近邻法将空间分辨率30 m的Landsat-5 TM二值积雪范围影像重采样为5 km分辨率，转换过程中设定5 km像元中30 m积雪像元占比≥50%作为积雪像元，否则为非雪像元，投影与AVHRR保持一致，最终将获取的5 km分辨率Landsat-5 TM二值积雪范围影像作为参考真值。

Landsat-5 TM参考积雪范围影像分为两组：第一组参考真值用作训练AVHRR积雪识别算法阈值，选取1986年、1987年、1996年、1997年和1998年积雪季（11月1日开始至次年3月31日）的无云Landsat-5 TM作为参考积雪范围影像，共1 094景，训练数据共包含积雪像元数为156 661个，非雪像元数为696 296个；第二组参考真值用于评估AVHRR积雪范围数据集，筛选1988年、1989年和1990年积雪季Landsat-5 TM无云且积雪占比大于50%的影像作为验证数据集，共939景，验证数据分布如图1所示。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   北半球气象台站和Landsat-5 TM验证数据集的空间分布

Fig. 1   Spatial distributions of meteorological stations and Landsat-5 TM validation data in the Northern Hemisphere


1.3　全球长时间序列逐日雪深数据集

本研究利用微波雪深数据集来填补AVHRR中的空缺值，使用Che等［27］采用SMMR和SSM/I2生成的全球长时间序列逐日雪深数据集。该微波雪深数据集，采用动态亮温梯度算法，通过建立不同频率亮温差与实测雪深在空间和季节上的动态关系［28-29］获取雪深。该数据集来源于时空三极环境大数据平台（http：//poles.tpdc.ac.cn/zh-hans/）。本研究使用的时间范围为1981年6月24日至2000年2月27日，数据集空间分辨率为25 km，利用最近邻法将该数据集重采样为5 km，并转换为与AVHRR相一致的投影。

1.4　地面雪深观测数据

雪深数据用于对AVHRR积雪识别算法和数据进行评估，本研究使用的北半球雪深数据来自于NOAA的国家环境信息中心网（http：//www.ncdc.noaa.gov/），是由世界气象组织（World Meteorological Organization， WMO）提供的气象台站数据。本研究使用1981—1999年2 546个气象台站积雪季的实测雪深数据，使用的气象台站分布如图1所示。

1.5　辅助数据

本研究共使用两种辅助数据，一种是数字高程模型，另一种是气候再分析数据-地表温度产品。数字高程模型数据（Digital Elevation Model， DEM）为SRTM（Shuttle Radar Topography Mission）90 m数字高程数据库第4版本数据，是由美国国家航天局联合美国国防部国家地理情报局通过航天飞机雷达测绘的数据产品（http：//srtm.csi.cgiar.org/），本研究中将STRM数据用作云雪区分和积雪识别算法辅助数据，也用来作为消除积雪伪像元的辅助数据。

ERA5气候再分析数据-地表温度逐日产品［30］是由欧盟资助欧洲中距离天气预报中心（European Center for Medium-Range Weather Forecast， ECMWF）执行哥白尼气候变化服务推出的适用于气候研究的数据集。本研究中将ERA5数据集也作为辅助数据用来消除可能的积雪伪像元，空间分辨率为10 km，同样采用最近邻法将上述数据重采样为5 km，并转换为与AVHRR相一致的投影。

2 方法

本研究AVHRR北半球积雪识别算法流程如图2所示，主要分为以下四个步骤：

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   北半球AVHRR积雪识别算法流程

Fig. 2   Flowchart of AVHRR snow discrimination algorithm in the Northern Hemisphere


（1） 数据预处理

根据AVHRR SR V4中自带的质量控制标记剔除反射率无效值。

（2） 云识别算法

针对冰云和积雪容易混淆的问题，结合云雪混淆区分技术，提高了积雪识别精度。

（3） 积雪识别算法

利用大量的样本进行训练，针对北半球积雪研究区的特性，获取多级决策树积雪识别算法的最优阈值，建立适合北半球的积雪识别算法，通过云、雪识别算法生成AVHRR L1级积雪范围数据集。

（4） 空缺值插补

针对AVHRR高纬度地区数据完全缺失和低纬度地区部分数据缺失问题，利用微波雪深数据集填充空缺数据。从而生成了一套北半球AVHRR 1981—1999年逐日积雪范围数据集，即AVHRR L2级积雪范围数据集。

2.1　数据预处理

根据AVHRR SR V4中自带的质量控制标记，如表2所示，剔除AVHRR各个通道的反射率无效值。有效值被作为下一步云、雪、非雪的区分，无效值被认为是空缺值。

表2   AVHRR SR V4的质量控制标记描述

Table 2  Descriptions of quality control bit flags assessment of AVHRR SR V4

比特值描述
12通道5值无效
11通道4值无效
10通道3值无效
9通道2值无效
8通道1值无效
7通道1~5值有效

新窗口打开| 下载CSV


2.2　云识别算法

关于云的识别算法，很多科学家也做了研究［10，31］，Hao等［25］针对AVHRR SR V4产品发展了一套云检测算法，该算法可以有效解决云雪混淆问题。因此本研究结合Hao等［25］生成的AVHRR云雪混淆区分技术进行云检测，算法如表3所示。

表3   云检测及其阈值

Table 3  Cloud detection tests and their threshold values

土地类型序号是/否高程/mRef1Ref2Ref3Ref1-Ref2NDVIBT11/K（BT37-BT11）/K（BT11-BT12）/K
高寒地（高程>300 m，BT11<260 K）A1是<3 000≥240>14.5
A2是≥3 000≥240>15.5
A3是<240>21.0
A4是>0.1>0.02>25.5
其他土地B1是<260>14.0
B2是>-0.02<310>10.5
B3是>0.3>-0.02<293>11.5
B4是>0.4>-0.03<293>11.5>-1
B5是>0.4<278>11.5>-1
B6是>0.3>0.02>11.5
B7否>0.5>288
B8否>310
B9否>1 000<0.4<-0.04>275
B10否<-0.04>300

注：“是/否”判别中，“是”表示若符合这行阈值条件，判别为云；“否”表示若符合这行阈值条件，判别为非云。

新窗口打开| 下载CSV


2.3　积雪识别算法

由于Hao等［25］开发的多级决策树积雪识别算法是面对中国地区来提取积雪，在北半球算法的阈值还有待改进。因此本研究采用该算法的方案（图3），结合Landsat-5 TM参考积雪范围影像作为参考真值，来获取北半球积雪识别算法最优阈值。根据前人对AVHRR积雪识别的分类研究［10，32］，涉及雪阈值算法的变量主要有Ref1、（Ref3/Ref2）和BT11三个一级阈值变量，NDVI和（Ref3-Ref2）两个二级阈值变量，NDSI一个三级阈值变量，其中NDSI的区分度最好。NDSI通常使用绿色（大约0.50 μm波长）和短波红外波段（大约1.60 μm波长）计算。由于在AVHRR SR V4中没有1.60 μm附近的短波红外观测，参考Hori等［10］使用3.75 μm波长的Ref3用于NDSI的计算，NDSI、NDVI可由下式计算得到。

NDSI=Ref1-Ref3Ref1+Ref3
NDVI=Ref2-Ref1Ref1+Ref2

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   多级决策树积雪识别算法流程

Fig. 3   Flowchart of snow discriminant algorithm based on multi-level decision tree


（1） 一级指标阈值

如图3所示，一级指标为Ref1、（Ref3/Ref2）、DEM与BT11。一级指标确立原则是排除大量的非雪数据，使尽可能多的雪被识别进而再用下一级指标进行区分。一级指标由于仅仅是初级筛选，因此一级指标阈值采用Hao等［25］算法的阈值。表4展示了一级指标和对应的阈值，即当Ref1>0.14，（Ref3/Ref2）<0.32、DEM<1 300 m并且BT11<280 K或者DEM≥1 300 m并且BT11<281.5 K时，则认为可能为积雪，否则为非雪，进行下一级别区分。

表4   北半球积雪识别算法最优阈值选取

Table 4  Selection of optimal threshold for snow discriminant algorithm in the Northern Hemisphere

积雪分类一级指标二级指标三级指标
Ref1BT11/KDEM/mRef3/Ref2NDVIRef3-Ref2NDSI
雪1>0.14<280.0<1 300<0.32<-0.11
>0.14<281.5≥1 300<0.32<-0.11
雪2>0.14<280.0<1 300<0.32<-0.96
>0.14<281.5≥1 300<0.32<-0.96
雪3>0.14<280.0<1 300<0.32≥-0.11≥-0.96≥0.76
>0.14<281.5≥1 300<0.32≥-0.11≥-0.96≥0.76

新窗口打开| 下载CSV


（2） 二级指标阈值

如图3所示，二级指标是NDVI和（Ref3-Ref2）。将一级指标筛选后可能为积雪的数据进行分段统计，将非雪样本1%各指标对应的反射率值作为二级指标阈值，通过计算：获取NDVI最优阈值-0.11，即满足NDVI<-0.11设定为第一类积雪（雪1）；获取（Ref3-Ref2）最优阈值-0.96，即满足（Ref3-Ref2）<-0.96设定为第二类积雪（雪2），不符合条件的被认为可能为积雪，进行下一级别区分。

（3） 三级指标阈值

如图3所示，三级指标是NDSI。利用二级指标筛选后可能为积雪的数据，通过计算最佳总体精度的方法得出三级指标，即逐步长方式计算最佳总体精度（步长设为0.01），最佳总体精度对应的NDSI值即为三级指标。最优阈值的确定方法如图4所示。图4（a）展示了积雪和非雪的NDSI频率分布，选择具有最高总体精度OA（88%）时对应的NDSI值作为最优阈值（0.76）。计算结果显示NDSI>0.76，符合条件的设定为第三类积雪（雪3），不符合条件的被判别为非雪。北半球积雪识别算法的最优阈值如表4所示。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   三级指标NDSI积雪和非雪频率分布（a）和总体精度（b）

Fig. 4   Frequency distribution of third-level indicator NDSI with and without snow （a） and overall accuracy （b）


2.4　空缺值插补

微波数据具有不受云干扰的特点，可对AVHRR的空缺值进行插补。利用Che等［27］、Dai等［28］生成的全球被动微波数据逐日雪深数据集进行插补。插补分为两个步骤：

（1）将25 km空间分辨率的雪深数据利用重采样方法（最近邻法）降尺度到5 km，投影和AVHRR一致。

（2）根据Che等［27］研究结果，插补规则如下：①纬度处于0°~60° N时，若AVHRR的云像元对应当日雪深≥2 cm，则该云像元被分类为积雪像元；若雪深<2 cm，则该云像元被分类为非雪像元。②纬度处于60°~90° N时，若所有像元对应当日雪深≥2 cm，则该像元被分类为积雪像元；若雪深<2 cm，则该像元被分类为非雪像元。

3 结果与精度评估

3.1　结果

AVHRR L1级和L2级逐日积雪范围数据集示意图如图5所示。图5（a）展示了通过云、雪识别后的北半球AVHRR L1级逐日积雪范围数据集示意图，其中雪像元赋值为1，云像元赋值为255，陆地赋值为0。图5（b）是北半球AVHRR L2级逐日积雪范围数据集示意图，其中雪像元赋值为1，云像元赋值为250，陆地赋值为0。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   L1级（a）和L2级（b）数据集示意图（选取日期为1998-01-13）

Fig. 5   Schematic diagram of datasets at Level 1（a） and Level 2（b） （The selected date is 1998-01-13）


3.2　精度评估方法

3.2.1　基于气象台站雪深

利用1981—1999年间北半球2 546个气象台站的雪深数据分别对AVHRR L1级和L2级积雪范围数据集进行精度评估。设定气象台站雪深大于等于1 cm时认为该气象台站所在的像元为积雪像元，雪深小于1 cm时认为该气象台站所在的像元为非雪像元［33］。将气象台站积雪信息与AVHRR各级积雪范围数据集相应位置的积雪信息进行比对验证。利用气象台站雪深数据转化为雪和非雪的二值数据，并作为参考真值，然后分别构建混淆矩阵分析总体精度及误差。具体的SS、SN、NS和NN的定义参见表5，其中SS指AVHRR积雪范围数据集和参考真值两者都判别为雪的像元个数，NS指AVHRR积雪范围数据集判别为有雪且参考真值判别为非雪的像元个数，SN指AVHRR积雪范围数据集判别为非雪且参考真值判别为雪的像元个数，NN指两者都判别为非雪的像元个数。

表5   精度验证混淆矩阵的指标含义

Table 5  Index implication of confusion matrix for accuracy validation

雪非雪
雪SSSN
非雪NSNN

新窗口打开| 下载CSV


分析数据集的总体精度（OA）、用户精度（UA）、制图精度（PA）和F-Score（FS）［34］。其中，OA为雪和非雪像元被正确分类的总体精度；UA为气象台站雪深数据被分类为积雪的概率，作为用户精度；PA为AVHRR积雪范围数据集被分类为积雪的概率，作为制图精度；FS为UA和PA的调和平均值，是综合考虑了UA和PA的结果，当FS较高时说明UA和PA较为均衡。

计算公式为

OA=SS+NNSS+SN+NS+NN
UA=SSSS+NS
PA=SSSS+SN
FS=2⋅UA⋅PAUA+PA

3.2.2　基于Landsat-5 TM

基于Landsat-5 TM的验证同样使用3.2.1节中构建混淆矩阵的方法，进一步评估AVHRR积雪范围数据集（L2级）的精度。将AVHRR积雪范围数据集分别在北半球和中国、北美洲、欧亚大陆、西伯利亚、北极地区这五个北半球重要积雪区［35］进行精度评估，分析数据集的OA、UA、PA和FS。（Landsat-5 TM参考积雪范围被当作参考真值）

3.3　验证结果与分析
3.3.1　基于气象台站雪深的精度评估

利用北半球1981—1999年的共2 546个实测雪深气象台站，分别对AVHRR L1级和L2级积雪范围数据集进行精度评估，计算结果如表6所示，AVHRR L1级数据集的总体精度为81.8%，AVHRR L2级数据集的总体精度为82.2%，总体验证精度较高。如图6所示，总体精度年际变化比较稳定，都在80%左右，其中1996年精度最高达到85.9%。AVHRR L1级数据集的UA为83.7%，PA为81.7%，AVHRR L2级数据集的UA为83.8%，PA为84.2%，表明AVHRR L1级数据集漏分误差略大于错分误差，尤其是在1993年，由于AVHRR大面积的数据空缺，导致漏分误差较大，而AVHRR L2级数据集，制图精度有所提高，表明雪深插值在保证精度的情况下，不仅可以填补空缺数据，而且可以有效降低漏分误差，预测正确的积雪像元个数在AVHRR识别的积雪像元中和气象台站记录的积雪信息中的比例较高。FS是平衡UA和PA的指数，AVHRR L1级和L2级数据集的FS分别为82.9%和85.1%，FS较高，表明漏分误差和错分误差都比较均衡，算法可有效解决积雪识别问题。

表6   气象台站验证评估结果

Table 6  Verification and evaluation results of meteorological stations

AVHRR数据集SSSNNSNNOA/%UA/%PA/%FS/%
L1级328 23967 72840 694179 92681.883.781.782.9
L2级870 093110 51274 793189 36482.283.884.285.1

注：SS、SN、NS和NN对应的数字为像元数。

新窗口打开| 下载CSV


图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   L1级（a）和L2级（b）数据集的OA、UA和PA年际变化

Fig. 6   Interannual variations of OA， UA and PA of datasets at Level 1 （a） and Level 2 （b）


3.3.2　基于Landsat-5 TM的精度评估

在空间上，利用1.2节中的筛选方法，获取了939景Landsat-5 TM参考积雪范围影像作为参考真值，其中中国155景、北美洲382景、欧亚大陆455景、西伯利亚174景及北极地区63景。使用3.2.2节中的验证方法，对待检验的AVHRR积雪范围数据集（L2级）进行评估。

表7展示了北半球和五个北半球重要积雪区（中国、北美洲、欧亚大陆、西伯利亚及北极地区）的验证评估结果。结果表明，AVHRR积雪范围数据集的总体验证结果OA、UA和PA分别为90.3%、90.2%和99.1%，该北半球积雪范围数据集的总体精度、制图精度和用户精度都较高，数据具有一定可靠性。其中PA最高达到99.1%，说明漏分误差较低，预测正确的积雪像元个数在AVHRR识别的积雪像元中和Landsat-5 TM识别的积雪像元中的比例较高。FS为93.4%，是UA和PA的平衡，较高则说明漏分误差和错分误差较为均衡，积雪识别算法比较有效。对中国、欧亚大陆和西伯利亚地区的总体验证结果OA、UA、PA和FS都在90.0%以上，北美洲地区总体精度OA为85.2%，总体精度在中纬度地区较高，而在北极地区，OA为73.7%，UA为73.7%，PA为99.9%，说明错分误差大于漏分误差，且漏分误差为0.1%，这是由于AVHRR在高纬度地区数据的完全缺失，使用2.4节中的数据和方法对其进行插补，该数据集是通过微波数据和气象台站数据融合的数据集，说明该数据集错分误差较高，漏分误差较低。

表7   Landsat-5 TM验证评估结果

Table 7  Verification and evaluation results of Landsat-5 TM

积雪区SSSNNSNNOA/%UA/%PA/%FS/%
北半球1172 5107 206126 32529 78290.390.299.193.4
中国117 5771 41410 7183 30191.591.498.593.6
北美洲456 87995669 8996 85485.284.899.790.0
欧亚大陆532 4614 95935 75912 62292.693.098.895.2
西伯利亚187 8391 12611 5591 18393.093.199.395.6
北极地区96 4661240 8112073.773.799.980.7

注：SS、SN、NS和NN对应的数字为像元数。

新窗口打开| 下载CSV


5景获取的Landsat-5 TM验证样本的详细信息如表8所示，OA、UA、PA和FS的平均值分别为93.4%、93.7%、99.2%和96.2%，数据总体验证精度较高，具有一定可靠性。每景验证样本对应的Landsat-5 TM参考积雪范围、AVHRR积雪范围，以及AVHRR和TM的两者对比如图7所示。说明在大范围覆盖的积雪区域，AVHRR积雪范围与Landsat-5 TM参考积雪范围是较吻合的，而在积雪边界区，AVHRR积雪识别算法有未能识别出的Landsat-5 TM参考积雪范围影像的积雪像元，也有更多识别出的积雪像元，这可能是因为AVHRR积雪范围数据集对较低积雪范围比例的像元识别能力较不准确。总体而言，在积雪范围识别上算法存在一定的漏分和错分情况，但算法在保持漏分误差和错分误差较为均衡的前提下，总体精度较高，数据具有可靠性。

表8   Landsat-5 TM验证样本信息

Table 8  Verification sample information of Landsat-5 TM

行列号过境日期云量雪量百分比/%OA/%UA/%PA/%FS/%积雪区
精度指标平均值93.493.799.296.2
133/0341989-11-08无云84.995.897.198.097.5中国
038/0321989-02-15无云90.296.196.499.497.8北美洲
154/0281990-12-16无云77.696.595.899.797.6欧亚大陆
139/0231990-03-26无云96.591.292.098.995.3西伯利亚
159/0101989-03-19无云98.687.387.3100.093.2北极地区

新窗口打开| 下载CSV


图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   Landsat-5 TM参考积雪范围影像与AVHRR积雪范围数据集对比

Fig. 7   Comparison between Landsat-5 TM snow map and AVHRR snow cover extent datasets： China （a）， North America （b）， Eurasia （c）， Siberia （d） and Arctic region （e）


4 结论

本研究基于AVHRR地表反射率数据，获取了基于多级决策树积雪识别算法的各指标最优阈值，并结合云雪混淆区分技术，生成了北半球AVHRR 1981—1999年L1级逐日积雪范围数据集。利用微波雪深数据插值填补了AVHRR的空缺数据，生成了北半球AVHRR 1981—1999年L2级逐日积雪范围数据集。最后利用气象台站实测雪深数据对L1级和L2级数据集分别进行精度验证，利用高空间分辨率遥感影像对L2级数据集进一步进行精度评估。

通过气象台站雪深数据对L1级积雪范围数据集验证，评估结果表明：数据集总体精度为81.8%，用户精度为83.7%，生产者精度为81.7%，数据集总体精度较高，漏分误差比错分误差略高。FS为82.9%，表明漏分误差和错分误差比较均衡，整体上，北半球积雪识别算法精度较高。

通过气象台站雪深数据对L2级积雪范围数据集验证，评估结果表明：数据集总体精度为82.2%，用户精度为83.8%，生产者精度为84.2%，FS为85.1%，L2级数据集总体精度较高，漏分误差和错分误差相近。同时结合Landsat-5 TM参考积雪范围影像，在面尺度上对L2级数据集进行进一步验证，结果表明：L2级数据集总体精度为90.3%，用户精度为90.2%，生产者精度为99.1%，在面尺度上，AVHRR L2级积雪范围数据集展现了更高的精度，表明插值方法可以有效解决AVHRR在高纬度地区数据的完全缺失和低纬度地区数据的部分缺失问题。该算法将为北半球积雪范围提取提供新的方案，L1级和L2级数据集也为北半球气候变化研究提供了重要数据支撑。

参考文献

Armstrong R L， Brodzik M J.

Recent Northern Hemisphere snow extent： a comparison of data derived from visible and microwave satellite sensors

［J］. Geophysical Research Letters， 2001， 28（19）： 3673-3676.

[本文引用: 1]

Warren S G.

Optical properties of snow

［J］. Reviews of Geophysics， 1982， 20（1）： 67-89.

[本文引用: 1]

Brown R D， Goodison B E.

Interannual variability in reconstructed Canadian snow cover， 1915—1992

［J］. Journal of Climate， 1996， 9（6）： 1299-1318.

Wang Jian， Li Hongxing， Hao Xiaohua， et al.

Remote sensing for snow hydrology in China： challenges and perspectives

［J］. Journal of Applied Remote Sensing， 2014， 8（1）： 084687.

[本文引用: 1]

Estilow T W， Young A H， Robinson D A.

A long-term Northern Hemisphere snow cover extent data record for climate studies and monitoring

［J］. Earth System Science Data， 2015， 7（1）： 137-142.

[本文引用: 1]

Che Tao， Hao Xiaohua， Dai Liyun， et al.

Snow cover variation and its impacts over the Qinghai-Tibet Plateau

［J］. Bulletin of Chinese Academy of Sciences， 2019， 34（11）： 1247-1253.

[本文引用: 1]

车涛， 郝晓华， 戴礼云， 等.

青藏高原积雪变化及其影响

［J］. 中国科学院院刊， 2019， 34（11）： 1247-1253.

[本文引用: 1]

Douville H， Royer J F.

Sensitivity of the Asian summer monsoon to an anomalous Eurasian snow cover within the Météo-France GCM

［J］. Climate Dynamics， 1996， 12（7）： 449-466.

[本文引用: 1]

McGowan H， Callow J N， Soderholm J， et al.

Global warming in the context of 2000 years of Australian alpine temperature and snow cover

［J］. Scientific Reports， 2018， 8： 4394.

[本文引用: 1]

Madan N J.

Snow ecology： an interdisciplinary examination of snow-covered ecosystems

［J］. Journal of Ecology， 2001， 89： 1096-1099.

[本文引用: 1]

Hori M， Sugiura K， Kobayashi K， et al.

A 38-year （1978—2015） Northern Hemisphere daily snow cover extent product derived using consistent objective criteria from satellite-borne optical sensors

［J］. Remote Sensing of Environment， 2017， 191： 402-418.

[本文引用: 4]

Huang Xiaodong， Hao Xiaohua， Yang Yongshun， et al.

Advances in snow-cover monitoring using optical remote sensing

［J］. Pratacultural Science， 2012， 29（1）： 35-43.

[本文引用: 1]

黄晓东， 郝晓华， 杨永顺， 等.

光学积雪遥感研究进展

［J］. 草业科学， 2012， 29（1）： 35-43.

[本文引用: 1]

Brown R， Derksen C， Wang Libo.

A multi-data set analysis of variability and change in Arctic spring snow cover extent， 1967—2008

［J］. Journal of Geophysical Research： Atmospheres， 2010， 115： D16111.

[本文引用: 1]

Che Tao， Li Xin.

Spatial distribution and temporal variation of snow water resources in China during

1993—2002［J］. Journal of Glaciology and Geocryology， 2005， 27（1）： 64-67.

[本文引用: 1]

车涛， 李新.

1993—2002年中国积雪水资源时空分布与变化特征

［J］. 冰川冻土， 2005， 27（1）： 64-67.

[本文引用: 1]

Che Tao， Li Xin.

Retrieval of snow depth in China by passive microwave remote sensing data and its accuracy assessment

［J］. Remote Sensing Technology and Application， 2004， 19（5）： 301-306.

[本文引用: 1]

车涛， 李新.

利用被动微波遥感数据反演我国积雪深度及其精度评价

［J］. 遥感技术与应用， 2004， 19（5）： 301-306.

[本文引用: 1]

Feng Qisheng， Zhang Xuetong， Liang Tiangang.

Dynamic monitoring of snow cover based on MOD

10A1 and AMSR-E in the north of Xinjiang Province， China［J］. Acta Prataculturae Sinica， 2009， 18（1）： 125-133.

[本文引用: 1]

冯琦胜， 张学通， 梁天刚.

基于MOD10A1和AMSR-E的北疆牧区积雪动态监测研究

［J］. 草业学报， 2009， 18（1）： 125-133.

[本文引用: 1]

Yu Hui， Feng Qisheng， Zhang Xuetong， et al.

An approach for monitoring snow depth based on AMSR-E data in the pastoral area of Northern Xinjiang

［J］. Acta Prataculturae Sinica， 2009， 18（4）： 210-216.

[本文引用: 1]

于惠， 冯琦胜， 张学通， 等.

基于AMSR-E信息的北疆牧区雪深遥感监测模型方法初探

［J］. 草业学报， 2009， 18（4）： 210-216.

[本文引用: 1]

Hartman R K， Rost A A， Anderson D M.

Operational processing of multi-source snow data

［C］// Proceedings of the 63rd Annual Western Snow Conference. Brush Prairie， WA， USA： Western Snow Conference， 1995： 147-151.

[本文引用: 1]

Hall D K， Riggs G A， Salomonson V V， et al.

MODIS snow-cover products

［J］. Remote Sensing of Environment， 2002， 83（1/2）： 181-194.

[本文引用: 2]

Gao Yang， Hao Xiaohua， He Dongcai， et al.

Snow cover mapping algorithm in the Tibetan Plateau based on NDSI threshold optimization of different land cover types

［J］. Journal of Glaciology and Geocryology， 2019， 41（5）： 1162-1172.

[本文引用: 1]

高扬， 郝晓华， 和栋材， 等.

基于不同土地覆盖类型NDSI阈值优化下的青藏高原积雪判别

［J］. 冰川冻土， 2019， 41（5）： 1162-1172.

[本文引用: 1]

Chokmani K， Dever K， Bernier M， et al.

Adaptation of the SNOWMAP algorithm for snow mapping over eastern Canada using Landsat-TM imagery

［J］. Hydrological Sciences Journal， 2010， 55（4）： 649-660.

[本文引用: 1]

Beijing Eastdawn Information Technology Co.， Ltd.

ESA’s next generation satellite Sentinel series

［J］. Remote Sensing Information， 2011（6）： 119-120.

[本文引用: 1]

北京东方道迩信息技术股份有限公司.

欧空局下一代卫星Sentinel系列

［J］. 遥感信息， 2011（6）： 119-120.

[本文引用: 1]

Hüsler F， Jonas T， Wunderle S， et al.

Validation of a modified snow cover retrieval algorithm from historical 1-km AVHRR data over the European Alps

［J］. Remote Sensing of Environment， 2012， 121： 497-515.

[本文引用: 1]

Zhou Hang， Aizen E， Aizen V.

Deriving long term snow cover extent dataset from AVHRR and MODIS data： Central Asia case study

［J］. Remote Sensing of Environment， 2013， 136： 146-162.

[本文引用: 1]

Wang Xuan， Hao Xiaohua， Wang Jian， et al.

Accuracy evaluation of long time series AVHRR snow cover area products in China

［J］. Remote Sensing Technology and Application， 2018， 33（6）： 994-1003.

[本文引用: 1]

王轩， 郝晓华， 王建， 等.

中国地区AVHRR长时间序列积雪范围产品精度评估

［J］. 遥感技术与应用， 2018， 33（6）： 994-1003.

[本文引用: 1]

Hao Xiaohua， Huang Guanghui， Che Tao， et al.

The NIEER AVHRR snow cover extent product over China： a long-term daily snow record for regional climate research

［J］. Earth System Science Data， 2021， 13（10）： 4711-4726.

[本文引用: 5]

Hall D K， Salomonson V V， Riggs G A.

Development of methods for mapping global snow cover using moderate resolution imaging spectroradiometer data

［J］. Remote Sensing of Environment， 1995， 54（2）： 127-140.

[本文引用: 1]

Che Tao， Li Xin， Jin Rui， et al.

Snow depth derived from passive microwave remote-sensing data in China

［J］. Annals of Glaciology， 2008， 49（1）： 145-154.

[本文引用: 3]

Dai Liyun， Che Tao， Ding Yongjian.

Inter-calibrating SMMR， SSM/I and SSMI/S data to improve the consistency of snow-depth products in China

［J］. Remote Sensing， 2015， 7（6）： 7212-7230.

[本文引用: 2]

Dai Liyun， Che Tao， Ding Yongjian， et al.

Evaluation of snow cover and snow depth on the Qinghai-Tibetan Plateau derived from passive microwave remote sensing

［J］. The Cryosphere， 2017， 11（4）： 1933-1948.

[本文引用: 1]

Di Napoli C， Barnard C， Prudhomme C， et al.

ERA5-HEAT： a global gridded historical dataset of human thermal comfort indices from climate reanalysis

［J］. Geoscience Data Journal， 2021， 8（1）： 2-10.

[本文引用: 1]

Stamnes K， Li Wei， Eide H， et al.

ADEOS-II/GLI snow/ice products： Part I： scientific basis

［J］. Remote Sensing of Environment， 2007， 111（2/3）： 258-273.

[本文引用: 1]

Hori M， Aoki T， Stamnes K， et al.

ADEOS-II/GLI snow/ice products： Part III： retrieved results

［J］. Remote Sensing of Environment， 2007， 111（2/3）： 291-336.

[本文引用: 1]

Zhao Qin.

Research on the characteristics of the temporal and spatial changes of snow phenology in China and its response to vegetation phenology

［D］. Taiyuan： Taiyuan University of Technology， 2021.

[本文引用: 1]

赵琴.

中国积雪物候时空变化特征及对植被物候的响应研究

［D］. 太原： 太原理工大学， 2021.

[本文引用: 1]

Olofsson P， Foody G M， Stehman S V， et al.

Making better use of accuracy data in land change studies： estimating accuracy and area and quantifying uncertainty using stratified estimation

［J］. Remote Sensing of Environment， 2013， 129： 122-131.

[本文引用: 1]

Li Zhen，

Liu Jiuliang

Tian Bangsen. Spatial and temporal series analysis of snow cover extent and snow water equivalent for satellite passive microwave data in the Northern Hemisphere （ 1978—2010）［C］// 2012 IEEE International Geoscience and Remote Sensing Symposium. New York： IEEE， 2012： 4871-4874.

[本文引用: 1]

/

〈 〉

陇ICP备2021001824号-10
　 甘公网安备 62010202000676号
网站版权 © 《冰川冻土》编辑部
地址：兰州市天水中路8号　
邮编：730000　电话：0931-8260767　
E-mail：edjgg@lzb.ac.cn
总访问量 　 今日访问 　 在线人数
技术支持：北京玛格泰克科技发展有限公司