海洋学报

传感器	波段名称	波长/nm	分辨率/m
HY-1C卫星CZI	蓝（波段1）	420～500	50
绿（波段2）	520～600	50
红（波段3）	610～690	50
近红外（波段4）	760～890	50

传感器

波段名称

波长/nm

分辨率/m

HY-1C卫星CZI

蓝（波段1）

420～500

绿（波段2）

520～600

红（波段3）

610～690

近红外（波段4）

760～890

传感器	波段名称	波长/nm	分辨率/m
HY-1C卫星CZI	蓝（波段1）	420～500	50
绿（波段2）	520～600	50
红（波段3）	610～690	50
近红外（波段4）	760～890	50

传感器

波段名称

波长/nm

分辨率/m

HY-1C卫星CZI

蓝（波段1）

420～500

绿（波段2）

520～600

红（波段3）

610～690

近红外（波段4）

760～890

海冰类型	特征描述
固定冰	固定冰是指与海岸、岛屿或海底冻结在一起，不能做水平运动，而能随海面升降做垂直运动的海冰。固定冰的封冻程度和范围对沿岸水产养殖业等有较大影响，固定冰的挤压对岸基和堤坝等建筑可能造成破坏。
白冰	当灰（白）冰厚度继续增加时，冰块颜色由灰白色过渡到白色，表面凹凸不平，形状不规则，成为白冰。厚度大于30 cm。
灰（白）冰	由于灰冰和灰白冰边界模糊，将两类作为一类进行区分。灰冰厚度为 10～15 cm ，灰白冰厚度为15～30 cm，因此，两类合为一类后，灰（白）冰厚度为10～30 cm，灰（白）冰表面比较粗糙，呈灰色或灰白色。
初生冰	当海上气温下降到海水冰点，海水开始结冰形成初生冰，多呈针状、薄片状、浆糊状或绵状。厚度在10 cm以内。

海冰类型

特征描述

固定冰

固定冰是指与海岸、岛屿或海底冻结在一起，不能做水平运动，而能随海面升降做垂直运动的海冰。固定冰的封冻程度和范围对沿岸水产养殖业等有较大影响，固定冰的挤压对岸基和堤坝等建筑可能造成破坏。

白冰

当灰（白）冰厚度继续增加时，冰块颜色由灰白色过渡到白色，表面凹凸不平，形状不规则，成为白冰。厚度大于30 cm。

灰（白）冰

由于灰冰和灰白冰边界模糊，将两类作为一类进行区分。灰冰厚度为 10～15 cm ，灰白冰厚度为15～30 cm，因此，两类合为一类后，灰（白）冰厚度为10～30 cm，灰（白）冰表面比较粗糙，呈灰色或灰白色。

初生冰

当海上气温下降到海水冰点，海水开始结冰形成初生冰，多呈针状、薄片状、浆糊状或绵状。厚度在10 cm以内。

海冰类型	特征描述
固定冰	固定冰是指与海岸、岛屿或海底冻结在一起，不能做水平运动，而能随海面升降做垂直运动的海冰。固定冰的封冻程度和范围对沿岸水产养殖业等有较大影响，固定冰的挤压对岸基和堤坝等建筑可能造成破坏。
白冰	当灰（白）冰厚度继续增加时，冰块颜色由灰白色过渡到白色，表面凹凸不平，形状不规则，成为白冰。厚度大于30 cm。
灰（白）冰	由于灰冰和灰白冰边界模糊，将两类作为一类进行区分。灰冰厚度为 10～15 cm ，灰白冰厚度为15～30 cm，因此，两类合为一类后，灰（白）冰厚度为10～30 cm，灰（白）冰表面比较粗糙，呈灰色或灰白色。
初生冰	当海上气温下降到海水冰点，海水开始结冰形成初生冰，多呈针状、薄片状、浆糊状或绵状。厚度在10 cm以内。

海冰类型

特征描述

固定冰

白冰

当灰（白）冰厚度继续增加时，冰块颜色由灰白色过渡到白色，表面凹凸不平，形状不规则，成为白冰。厚度大于30 cm。

灰（白）冰

初生冰

当海上气温下降到海水冰点，海水开始结冰形成初生冰，多呈针状、薄片状、浆糊状或绵状。厚度在10 cm以内。

特征标识	定义公式
注：BLUE、GREEN、RED、NIR分别表示HY-1C卫星CZI数据的蓝、绿、红、近红外波段。
NDWI_H	$\rm \left( {NIR - BLUE} \right)/(NIR + BLUE) $
B_G	$\rm (BLUE - GREEN)/(BLUE + GREEN) $
B_R	$ \rm (BLUE - RED)/(BLUE + RED) $
G_R	$\rm (GREEN - RED)/(GREEN + RED) $

特征标识

定义公式

注：BLUE、GREEN、RED、NIR分别表示HY-1C卫星CZI数据的蓝、绿、红、近红外波段。

NDWI_H

$\rm \left( {NIR - BLUE} \right)/(NIR + BLUE) $

B_G

$\rm (BLUE - GREEN)/(BLUE + GREEN) $

B_R

$ \rm (BLUE - RED)/(BLUE + RED) $

G_R

$\rm (GREEN - RED)/(GREEN + RED) $

特征标识	定义公式
注：BLUE、GREEN、RED、NIR分别表示HY-1C卫星CZI数据的蓝、绿、红、近红外波段。
NDWI_H	$\rm \left( {NIR - BLUE} \right)/(NIR + BLUE) $
B_G	$\rm (BLUE - GREEN)/(BLUE + GREEN) $
B_R	$ \rm (BLUE - RED)/(BLUE + RED) $
G_R	$\rm (GREEN - RED)/(GREEN + RED) $

特征标识

定义公式

注：BLUE、GREEN、RED、NIR分别表示HY-1C卫星CZI数据的蓝、绿、红、近红外波段。

NDWI_H

$\rm \left( {NIR - BLUE} \right)/(NIR + BLUE) $

B_G

$\rm (BLUE - GREEN)/(BLUE + GREEN) $

B_R

$ \rm (BLUE - RED)/(BLUE + RED) $

G_R

$\rm (GREEN - RED)/(GREEN + RED) $

	波段1	波段2	波段3	波段4
均值	5.687 4	3.56	5.234 4	2.27
方差	0.049 2	0.047 5	0.045 9	0.044 9
同质性	2.927 5	2.152 5	2.559	1.66
对比度	0.056 5	0.052	0.051	0.048
相异性	0.807 9	0.776	0.785 4	0.67
熵	4.003	2.354 1	3.293 5	1.69
二阶矩	3.540 8	1.897 6	2.76	1.41
相关性	0.570 3	0.297 1	0.456 2	0.261 5

波段1

波段2

波段3

波段4

均值

5.687 4

3.56

5.234 4

2.27

方差

0.049 2

0.047 5

0.045 9

0.044 9

同质性

2.927 5

2.152 5

2.559

1.66

对比度

0.056 5

0.052

0.051

0.048

相异性

0.807 9

0.776

0.785 4

0.67

熵

4.003

2.354 1

3.293 5

1.69

二阶矩

3.540 8

1.897 6

2.76

1.41

	波段1	波段2	波段3	波段4
均值	5.687 4	3.56	5.234 4	2.27
方差	0.049 2	0.047 5	0.045 9	0.044 9
同质性	2.927 5	2.152 5	2.559	1.66
对比度	0.056 5	0.052	0.051	0.048
相异性	0.807 9	0.776	0.785 4	0.67
熵	4.003	2.354 1	3.293 5	1.69
二阶矩	3.540 8	1.897 6	2.76	1.41
相关性	0.570 3	0.297 1	0.456 2	0.261 5

相关性

0.570 3

0.297 1

0.456 2

0.261 5

分类方法	海冰类型	生产者精度/%	用户精度/%	总体精度/%	kappa系数
基于光谱特征的SVM海冰分类	固定冰	66.84	39.06	74.58	0.66
白冰	33.05	76.02
灰（白）冰	98.70	74.81
初生冰	69.08	92.62
海水	99.89	99.35
基于纹理特征的SVM海冰分类	固定冰	66.84	35.72	73.41	0.65
白冰	34.87	79.74
灰（白）冰	98.70	73.93
初生冰	69.08	91.37
海水	99	99.35
基于最优特征集的SVM海冰分类	固定冰	85.43	64.44	84.89	0.80
白冰	60.69	85.14
灰（白）冰	96.14	82.14
初生冰	84.22	96.65
海水	99.91	99.40

分类方法

海冰类型

生产者精度/%

用户精度/%

总体精度/%

kappa系数

基于光谱特征的SVM海冰分类

固定冰

66.84

39.06

74.58

0.66

白冰

33.05

76.02

灰（白）冰

98.70

74.81

初生冰

69.08

92.62

海水

99.89

99.35

基于纹理特征的SVM海冰分类

固定冰

66.84

35.72

73.41

0.65

白冰

34.87

79.74

灰（白）冰

98.70

73.93

初生冰

69.08

91.37

海水

99.35

基于最优特征集的SVM海冰分类

固定冰

85.43

64.44

84.89

0.80

白冰

60.69

85.14

灰（白）冰

96.14

82.14

初生冰

84.22

96.65

海水

99.91

99.40

分类方法	海冰类型	生产者精度/%	用户精度/%	总体精度/%	kappa系数
基于光谱特征的SVM海冰分类	固定冰	66.84	39.06	74.58	0.66
白冰	33.05	76.02
灰（白）冰	98.70	74.81
初生冰	69.08	92.62
海水	99.89	99.35
基于纹理特征的SVM海冰分类	固定冰	66.84	35.72	73.41	0.65
白冰	34.87	79.74
灰（白）冰	98.70	73.93
初生冰	69.08	91.37
海水	99	99.35
基于最优特征集的SVM海冰分类	固定冰	85.43	64.44	84.89	0.80
白冰	60.69	85.14
灰（白）冰	96.14	82.14
初生冰	84.22	96.65
海水	99.91	99.40

分类方法

海冰类型

生产者精度/%

用户精度/%

总体精度/%

kappa系数

基于光谱特征的SVM海冰分类

固定冰

66.84

39.06

74.58

0.66

白冰

33.05

76.02

灰（白）冰

98.70

74.81

初生冰

69.08

92.62

海水

99.89

99.35

基于纹理特征的SVM海冰分类

固定冰

66.84

35.72

73.41

0.65

白冰

34.87

79.74

灰（白）冰

98.70

73.93

初生冰

69.08

91.37

海水

99.35

基于最优特征集的SVM海冰分类

固定冰

85.43

64.44

84.89

0.80

白冰

60.69

85.14

灰（白）冰

96.14

82.14

初生冰

84.22

96.65

海水

99.91

99.40

影像获取时间	海冰类型	生产者精度/%	用户精度/%	总体精度/%	kappa系数
2020年12月19日	固定冰	93.41	81.28	93.67	0.89
灰（白）冰	89.76	91.94
海水	94.91	99.78
初生冰	92.93	79.35
2021年1月10日	固定冰	70.54	98.59	91.75	0.88
灰（白）冰	99.03	76.31
海水	97.85	99.85
初生冰	96.10	87.43

影像获取时间

海冰类型

生产者精度/%

用户精度/%

总体精度/%

kappa系数

2020年12月19日

固定冰

93.41

81.28

93.67

0.89

灰（白）冰

89.76

91.94

海水

94.91

99.78

初生冰

92.93

79.35

2021年1月10日

固定冰

70.54

98.59

91.75

0.88

灰（白）冰

99.03

76.31

海水

97.85

99.85

初生冰

96.10

87.43

影像获取时间	海冰类型	生产者精度/%	用户精度/%	总体精度/%	kappa系数
2020年12月19日	固定冰	93.41	81.28	93.67	0.89
灰（白）冰	89.76	91.94
海水	94.91	99.78
初生冰	92.93	79.35
2021年1月10日	固定冰	70.54	98.59	91.75	0.88
灰（白）冰	99.03	76.31
海水	97.85	99.85
初生冰	96.10	87.43

影像获取时间

海冰类型

生产者精度/%

用户精度/%

总体精度/%

kappa系数

2020年12月19日

固定冰

93.41

81.28

93.67

0.89

灰（白）冰

89.76

91.94

海水

94.91

99.78

初生冰

92.93

79.35

2021年1月10日

固定冰

70.54

98.59

91.75

0.88

灰（白）冰

99.03

76.31

海水

97.85

99.85

初生冰

96.10

87.43

基于最优特征集的HY-1C卫星海岸带成像仪影像海冰分类方法研究

PDF下载

臧金霞 ¹ , 刘建强 ²^,^* , 殷晓斌 ¹ , 曾韬 ² , 周磊 ¹

海洋学报 | 海洋一号C/D卫星数据处理与典型应用 2022,44(5): 35-46

收起

海洋学报 | 海洋一号C/D卫星数据处理与典型应用 2022, 44(5): 35-46

基于最优特征集的HY-1C卫星海岸带成像仪影像海冰分类方法研究

全屏

臧金霞¹, 刘建强²^,^*, 殷晓斌¹, 曾韬², 周磊¹

作者信息

1.航天宏图信息技术股份有限公司，北京 100195

2.国家卫星海洋应用中心，北京 100081

臧金霞（1992-），女，山东省德州市人，算法工程师，主要从事海洋遥感应用等方面研究。E-mail: zangjinxia@piesat.cn

通讯作者:

刘建强（1964-），男，湖南省益阳市人，研究员，主要从事海洋遥感应用研究与地面应用系统建设。E-mail: jqliu@mail.nsoas.org.cn

Study on sea ice classification of HY-1C satellite coastal zone imager images based on the optimal feature set

Jinxia Zang¹, Jianqiang Liu²^,^*, Xiaobin Yin¹, Tao Zeng², Lei Zhou¹

Affiliations

1. PIESAT Information Technology Co., Ltd., Beijing 100195, China

2. National Satellite Ocean Application Service, Beijing 100081, China

出版时间: 2022-05-01 doi: 10.12284/hyxb2022021

文章导航

摘要

收起

基于海洋一号C（HY-1C）卫星海岸带成像仪（CZI）遥感影像，提出了一种基于最优特征集的支持向量机海冰分类方法。分别提取CZI影像的光谱特征和纹理特征，采用基于距离可分性的判据进行特征选择，得到最优特征集，以最优特征集作为支持向量机分类器输入，分别对3期辽东湾海域CZI影像开展海冰分类实验和结果分析。结果表明：本文方法得到的海冰分类结果精度优于仅利用光谱特征或纹理特征的海冰分类精度；基于本文方法的3期影像的海冰分类精度均较高，2020年12月19日、2021年1月10日与2021年1月16日的海冰分类总体精度分别为93.67%、91.75%、84.89%，均在80%以上；利用海冰分类结果图估算海冰面积，发现3期辽东湾海冰面积依次增大，最大约为11 998.98 km²。

关键词

海洋一号C卫星 / 光谱特征 / 纹理特征 / 最优特征集 / 海冰分类

Abstract

收起

A support vector machine (SVM) sea ice classification method of Haiyang-1C (HY-1C) satellite coastal zone imager (CZI) images based on the optimal feature set is proposed in this paper. The spectral features and the texture features of CZI images are extracted, and then distance separability criterion is used for feature selection to obtain the optimal feature set. The sea ice classification experiment and analysis of the three CZI images of Liaodong Bay are carried out based on SVM classification method with the optimal feature set as the input of the classifier. The results show that the sea ice classification accuracy obtained by the proposed method is better than that of only using the spectral features or the texture features. The sea ice classification accuracy of December 19, 2020, January 10, 2021 and January 16, 2021 are 93.67%, 91.75% and 84.89%, respectively, all above 80%. The sea ice area of Liaodong Bay is estimated according to the sea ice classification map. It is found that the sea ice area of Liaodong Bay in the three images increased successively, and the maximum area is about 11 998.98 km².

Key words

Haiyang-1C satellite / spectral features / texture features / optimal feature set / sea ice classification

引用本文

臧金霞, 刘建强, 殷晓斌, 曾韬, 周磊. 基于最优特征集的HY-1C卫星海岸带成像仪影像海冰分类方法研究. 海洋学报, 2022 , 44 (5) : 35 -46 . DOI: 10.12284/hyxb2022021

Jinxia Zang, Jianqiang Liu, Xiaobin Yin, Tao Zeng, Lei Zhou. Study on sea ice classification of HY-1C satellite coastal zone imager images based on the optimal feature set[J]. Haiyang Xuebao, 2022 , 44 (5) : 35 -46 . DOI: 10.12284/hyxb2022021

正文

收起

1 引言

收起

海冰一般是指由海水直接冻结而成的咸水冰，是极地海域和某些高纬度地区最突出的海洋灾害^[1]。我国渤海北部海域纬度较高，每年的11月末至翌年3月初出现海冰，主要为一年冰。海冰的冻结和漂移对海上生产、航运及油气勘探等都有不同程度的影响。近几十年来的几次严重冰情甚至造成了海冰灾害^[2]，尤其在辽东湾海域，由于纬度高且受到冬季西北风的影响，最易结冰。海冰监测成为海情监测的重点业务。海冰分类是海冰监测的前提，海冰密集度、海冰面积、海冰最大边缘线等信息都是由海冰分类识别衍生而来^[3]。

目前，卫星遥感技术是海冰监测的主要技术手段，相比目测法、器测法等传统监测方法，利用卫星遥感影像可实现实时、快速、大面积地海冰监测^[4]。合成孔径雷达（Synthetic Aperture Radar, SAR）卫星数据是常见的海冰监测遥感数据，国内外学者基于SAR卫星数据开展了大量的海冰分类方法研究^[4-9]。光学影像虽然易受云雾影响，但其价格低廉，大部分能够免费获取，时间分辨率高，也成为海冰监测主要技术手段。目前，国内外光学卫星海冰监测研究多集中在低空间分辨率卫星遥感数据，如MODIS数据、AVHRR数据等^[10-12]。这些数据覆盖范围广，重访周期短，可实现大范围实时海冰监测，但其空间分辨率较低，因此海冰分类精度不高。中高分辨率光学卫星数据具有较高的空间分辨率，在海冰类型判读等方面具有独特优势，可进行精细化海冰分类。王姝力等^[13]基于Landsat-8和Sentinel-1A数据对渤海湾进行海冰分类，对比分析了两种数据分别利用不同分类算法分类得到的海冰分类精度。张晰等^[14]对Landsat-8和GF-1卫星的海冰探测能力进行对比分析，发现在冰水区分能力方面GF-1数据优于Landsat-8数据，在海冰类型识别方面，Landsat-8数据优于GF-1数据。Barbieux等^[15]对Landsat-8数据的红、近红外和短波红外波段进行波段组合，构建了冰湖指数（Icy Lakes Index），有效实现了冰水分离。王志勇等^[16]利用Sentinel-2、Landsat-8、Landsat-7 中高分辨率影像对辽东湾海冰进行监测，利用归一化水指数（NDWI）实现冰水分离，并进行海冰厚度反演和海冰分类研究。以上研究工作虽利用中高分辨率光学数据开展海冰监测研究，但仅利用其光谱信息，对中高分辨率光学数据的纹理信息未有效利用。

为发挥中高分辨率光学遥感影像纹理信息在海冰分类中的作用，本文以辽东湾海域为研究区，基于海洋一号C（HY-1C）卫星海岸带成像仪（Coastal Zone Imager，CZI）数据开展海冰分类研究。首先提取其光谱信息和纹理信息，然后进行特征提取与筛选，构建最优特征集，最后将最优特征集输入到支持向量机（Support Vector Machine，SVM）分类器中进行海冰分类，比较分析多特征海冰分类与单一特征海冰分类结果精度的差异。利用海冰分类结果进一步估算海冰总面积与各海冰类型面积。

2 数据和方法

收起

2.1 研究区域和数据

2.1.1 研究区

选择辽东湾海域作为研究区，研究区域如图1所示。渤海是北半球纬度最低的大面积结冰海域，由于辽河、黄河等沿海大陆河流的淡水注入，使渤海含盐量较低；另一方面随河流带来的泥沙不断沉积，改变海底和海岸的地貌，造成渤海变浅，使渤海水温极易受陆地气温影响^[17]。这些情况都有利于冬季渤海海冰的生成。渤海是典型的一年冰海区，每年11月末到翌年3月初总会出现不同程度的结冰现象^[18]。辽东湾位于渤海北部，纬度最高，是渤海海冰灾害最严重的区域。

2.1.2 数据

本文所用数据为海洋一号C卫星海岸带成像仪数据，数据下载地址为https://osdds.nsoas.org.cn/#/。其空间分辨率为50 m，幅宽为1 000 km，设置有蓝、绿、红、近红外4个波段，具体信息见表1。选取2020年12月19日、2021年1月10日和2021年1月16日共3期数据，如图2所示。这3个时期的海冰面积随着气温下降不断增加，由初冰期向盛冰期过渡，海冰类型有所变化。

2.2 研究方法

本文的总体技术流程如图3所示。首先进行海冰类型划分与数据预处理，对预处理后的CZI影像提取光谱特征和纹理特征；然后进行特征选择，得到最优特征集；最后基于最优特征集进行SVM海冰分类，对分类结果进行对比验证分析。

2.2.1 海冰类型划分

国际气象组织（World Meteorological Organization，WMO）对海冰类型的划分标准是根据两极和高纬度区域的海冰情况定义的，包括多种方式，如海冰生长过程、海冰形态、海冰表面特征、海冰运动状态等^[17]。结合海冰划分标准，根据遥感影像上的海冰表征和海冰生长过程，经专家解译，本文将实验海域冰水覆盖类型分为5类：固定冰、白冰、灰（白）冰、初生冰、海水。各海冰类型的特征见表2 。图4展示了HY-1C CZI遥感影像中4种海冰类型的表征。由于海冰形成初期不会形成白冰类型，在进行分类时，将初冰期影像中白冰类型去除，分为固定冰、灰（白）冰、初生冰、海水4类。

2.2.2 数据预处理

本文利用PIE-Basic软件进行预处理，主要包括图像裁剪、图像配准、海陆掩膜等。为减少数据运算量，对影像进行裁剪，裁剪出辽东湾海域影像；对多时相影像进行图像配准，以2021年1月16日CZI影像为基准，对2020年12月19日和2021年1月10日的两景影像进行配准，各选取约500个配准点，XY误差小于1 m；利用海陆掩膜矢量文件对影像进行海陆掩膜处理，去除陆地对海冰提取的影响。

2.2.3 特征提取与选择

图像分类的关键是图像特征提取，常用的特征包括光谱特征、纹理特征、形状特征等。研究表明多特征组合后的分类精度高于单一特征，但特征过多可能会造成信息冗余，降低分类精度^[19]。因此需要进行特征选择，去除冗余信息。本文提取了HY-1C卫星CZI影像的光谱特征和纹理特征，对多特征进行特征选择，得到最优特征集。

光学影像光谱特征丰富，在冰水识别方面，由于海冰反射率比水体反射率高，本文利用水体指数区分海水与海冰。2012年，Haq等^[20]提出了NDWI_H指数，通过实验得出该指数比NDWI更适用于冰水识别，因此选择NDWI_H作为光谱特征之一。张晰等^[14]分别对GF1和Landsat8影像波段计算得到多种光谱特征，其中GF1影像的B_G、B_R、G_R等光谱特征对海冰分类的效果较好。考虑到HY-1C CZI波段设置与GF1卫星波段设置相似，本文提取了CZI影像的B_G、B_R、G_R 3个光谱特征以及张晰等^[14]提取的其他光谱特征。通过对比分析发现B_G、B_R、G_R 3个光谱特征海冰分类效果较好。本文选用的4个光谱特征的计算公式见表3。

纹理特征分析是目标识别、图像分割和图像分类中的重要方法^[21]。在遥感影像中，海冰类型是粗糙、不规则的，海水相对光滑均匀，因此纹理特征特别适用于海冰类型识别。本文利用灰度共生矩阵方法（GLCM）进行纹理特征提取。GLCM是一种通过研究灰度的空间相关特性来描述纹理的常用方法，它描述了成对像素的灰度组合分布^[22]。利用GLCM方法对HY-1C影像的蓝、绿、红、近红外4个波段分别进行纹理特征提取。经过多次实验，纹理特征计算参数设置为：滑动窗口为3×3，步长为1，方向为45°，灰度级为64。每个波段提取的纹理特征包括：均值、方差、熵、对比度、二阶矩、同质性、相关性、相异性，共得到32个纹理特征。过多的特征会引起冗余，影响计算效率，降低分类精度^[23]，因此需要进行特征选择以降低特征维数，提高计算效率，改善分类效果。本文利用基于距离可分性的特征选择判据^[24]，特征选择的原则是，使各类地物类间离散度较大而类内离散度较小。

定义类间离散度S_b、类内离散度S_w和判别数J，一般J值越大，携带信息量越丰富。其中：

(1)

$ {S _{\rm{w}}} = \sum\limits_{i = 1}^n {(\sigma )} \text{，} $

(2)

$ {S _{\rm{b}}} = {\sum\limits_{i = 1}^n {(\overline {ui} - \overline u )} ^2} \text{，} $

(3)

$ J = {S _{\rm{b}}}/{S _{\rm{w}}} \text{，} $

式中，$ \sigma $为样本方差；$\overline{u i}$为第i类的均值；$ \overline{u} $为总体均值。

为进一步验证特征选择结果，利用LDA（Linear Discriminant Analysis）算法判断4个波段纹理特征的分类准确度^[25]。LDA算法是在计算出类内离散度S_w和类间离散度S_b的基础上，进一步计算出最佳投影方向${\boldsymbol{ \omega}}$，将多维特征投影到一个方向上，可直观看出分类的效果。

当投影为最佳投影时，满足：

(4)

$ {J_{{\rm{FFC}}}} = \frac{{{{\boldsymbol{ \omega}} ^{\rm{T}}} {S _{\rm{b}}}{\boldsymbol{ \omega}} }}{{{{\boldsymbol{ \omega}} ^{\rm{T}}} {S _{\rm{w}}}{\boldsymbol{ \omega}} }} \text{，} $

定义拉格朗日函数$L({\boldsymbol{ \omega}},\lambda )$，$ \lambda $为拉格朗日乘子，

(5)

$ L({\boldsymbol{ \omega}},\lambda ) = {{\boldsymbol{ \omega}}^{\rm{T}}}{S _{\rm{b}}}{\boldsymbol{ \omega}} - \lambda ({{\boldsymbol{ \omega}} ^{\rm{T}}}{S _{\rm{w}}}{\boldsymbol{ \omega}} ) \text{，} $

对${\boldsymbol{ \omega}}$求导并令偏导数为0可得到，最佳投影方向由对应特征值的特征向量$ ({\phi _1},{\phi _2},\cdots,{\phi _n}) $组成，可表示为

(6)

${\boldsymbol{ \omega}} = ({\phi _1},{\phi _2},\cdots,{\phi _n}) . $

2.2.4 海冰分类

选择SVM分类方法进行海冰分类，SVM分类是定义在特征空间上的能够使间隔最大的线性分类器，通过核函数转化为求解凸二次规划问题，其通过类别边缘的训练样本，即支持向量，寻找两类间最优的分类超平面^[26]。在非线性变换空间，通过最小化分类器的泛化误差和采用基于核函数的公式化，来计算两类间的线性判别函数。相比神经网络等分类方法，其具有样本小、泛化能力强、计算复杂度适中、结构设计简单等优点，是海冰分类识别领域广泛使用的一种分类方法。SVM分类中，核函数选择至关重要，常用的核函数有线性核函数、多项式核函数、Sigmoid核函数、径向基核函数。

线性核函数为

(7)

$ K({f_s},{f_t}) = \left\langle {{f_s},{f_t}} \right\rangle, $

多项式核函数为

(8)

$ K({f_s},{f_t}) = {\left[ {\left\langle {{f_s},{f_t}} \right\rangle + b} \right]^d},d = 1,2,\cdots, $

Sigmoid核函数为

(9)

$ K({f_s},{f_t}) = \tanh \left[ {v\left\langle {{f_s},{f_t}} \right\rangle + c} \right], $

径向基核函数为

(10)

$ K({f_s},{f_t}) = \exp ( - \gamma {\left\| {{f_s} - {f_t}} \right\|^2}), $

式中，$ b,v,c,\gamma $为核函数参数。

本文选择径向基核函数（式（10）），径向基核函数是SVM海冰分类中的常用核函数。核函数$ \gamma $参数设置为0.25，惩罚参数设置为100。

SVM属于监督分类，选取训练样本是监督分类必不可少的环节。本文利用PIE-Basic软件，以专家解译图为参考，在待分类的每景影像中各均匀选取约120个训练样本，50个验证样本。

3 实验结果与分析

收起

3.1 特征选择结果分析

由于2021年1月16日处于盛冰期，已出现白冰类型，特征选择以2021年1月16日数据为例进行展示。经过波段运算提取出的NDWI_H、B_G、G_R、B_R共4个光谱特征如图5所示。为直观判断4种光谱特征对5种海冰类型的区分度，基于选取的训练样本，以光谱特征值为横轴，分类类型为竖轴作散点图，如图6所示。从图6可看出，4种光谱特征中，海水和固定冰光谱特征值差异最大，初生冰与固定冰、白冰之间的光谱特征值差异也较大，能够较好区分。但海水与初生冰光谱特征值重叠较多，固定冰、白冰与灰（白）冰之间光谱特征值重叠也较多，因此仅依据光谱特征值难以有效区分5种分类类型。

利用GLCM方法提取了CZI影像4个波段的共32个纹理特征，为减少信息冗余，需进行特征选择。经过计算，得到4个波段共32个纹理特征的J值，如表4所示。

整体来看，波段1的J值普遍较高，其中均值、同质性、熵、二阶矩4个纹理特征J值最高，按J值大小排列依次为均值、熵、二阶矩、同质性。利用LDA算法将4个波段中J值最高的4个纹理特征（均值、同质性、熵、二阶矩）分别进行分类投影，如图7所示。

从图7中可看出，4个波段的纹理特征均能将5种类型区分开。5种类型中，4个波段的白冰分类效果均较差，波段2、波段3和波段4的固定冰分类效果较差，波段1中的固定冰分类效果较好。总体来看，J值最高的4个纹理特征可实现海冰分类。综合分析5种类型的分类效果，波段1的分类结果优于其他波段分类结果，与J值分析结果相符。因此，选择波段1的4个纹理特征（图8）。

将NDWI_H、B_G、B_R、G_R 4个光谱特征与均值、同质性、熵、二阶矩4个纹理特征进行组合，得到最优特征集。

3.2 海冰分类结果分析

基于提取的最优特征集进行SVM海冰分类实验，将2021年1月16日的CZI影像分为固定冰、白冰、灰（白）冰、初生冰、海水共5种类型。为验证不同特征对分类精度的影响，分别开展基于光谱特征和纹理特征的海冰分类实验，与基于最优特征集的海冰分类结果进行对比分析，分类结果如图9所示。

采用主观和客观评价标准进行结果评价。主观评价为目视对分类结果进行评价；客观评价为利用分类评价指标评价分类精度。图9d为海冰类型解译图，与图9d对比来看，3种分类结果均能将海冰类型大致区分开，但具体到某一类海冰类型，3种分类结果的精度不同。在基于光谱特征的SVM海冰分类图像中（图9a），沿岸区域的初生冰被错分类为灰（白）冰，且固定冰和白冰之间出现混淆，部分白冰错分类为固定冰；在基于纹理特征的SVM海冰分类图像中（图9b），海水和初生冰出现一定程度混淆，固定冰和白冰也不能较好区分，部分白冰误分类为固定冰。在基于最优特征集的SVM分类结果中，5种分类类型均得到较好区分；相比基于光谱特征或者纹理特征的分类结果，白冰和固定冰的区分效果有明显改善，说明光谱特征和纹理特征之间具有较好的互补性，改善了海冰类型提取精度。以图9d海冰类型解译图为参考，选取验证样本，采用混淆矩阵方法进行客观评价，分类结果精度见表5。

由表5可知，本文提出的基于最优特征集的SVM海冰分类总体精度和kappa系数最高，总体精度为84.89%，kappa系数为0.80。其次为基于光谱特征的SVM海冰分类，总体精度为74.58%，kappa系数为0.66。基于纹理特征的SVM海冰分类总体精度和kappa系数最低，总体精度为73.41%，kappa系数为0.65。通过特征选择得到的最优特征集包含光谱和纹理两种特征，最大程度保留了分类有效信息，同时降低了特征冗余，特征之间实现互补，提高了分类精度。在海冰类型方面，3种分类方法中，白冰的生产者精度均较低，固定冰的用户精度较低，主要是由于白冰与固定冰的特征相近，导致部分白冰被分类为固定冰，这与主观评价的结论一致。基于光谱特征或纹理特征的分类方法中，初生冰的生产者精度较低，结合主观评价结果，发现在基于光谱特征的分类结果中，沿岸部分的初生冰被误分为灰（白）冰，基于纹理特征的分类结果中，部分初生冰被误分为海水。

为验证本文方法的有效性，排除单景影像的偶然因素影响，分别对2020年12月19日和2021年1月10日的HY-1C CZI遥感影像进行海冰分类。结合两个时期的影像表征，考虑到形成白冰的几率较小，因此海冰分类类型分为固定冰、灰（白）冰、初生冰、海水4种类型。分类结果如图10所示。

总体来看，2021年1月10日的海冰面积相较2020年12月19日大幅增加。2020年12月19日的海冰类型中初生冰面积较大，2021年1月10日海冰类型中灰（白）冰海冰类型占比最大。由于1月份气温降低，海水大面积冻结，形成大面积的灰（白）冰类型。从分类图中可看出，海冰类型中灰（白）冰与固定冰之间存在一定程度地混淆。利用混淆矩阵客观评价分类精度，精度评价见表6。

由表6可知，2020年12月19日遥感影像分类总体精度为93.67%，kappa系数为0.89，2021年1月10日遥感影像海冰分类总体精度为91.75%，kappa系数为0.88。分类精度均较高，能够达到一般海冰分类精度要求，进一步验证了本文方法的有效性。对比3个时期的海冰分类精度，发现精度从高至低依次为：2020年12月19日、2021年1月10日、2021年1月16日。随着气温降低，海冰类型增加白冰类型，易误分类海冰类型所占比例增大，分类精度逐渐降低。

利用海冰分类结果估算3个时期的海冰总面积与各海冰类型面积，计算公式为：海冰面积=海冰像元数×2 500/1 000 000，海冰面积单位为平方千米。为直观对比3个时期的海冰面积和类型变化，做海冰面积变化折线图与海冰类型柱状图，如图11所示。从图11a中可得到3个时期的海冰面积逐渐增大，分别为2 635.14 km²、10 674.27 km²、11 998.98 km²，国家卫星海洋应用中心通报的3个日期的海冰面积分别为2 365 km²、10 987 km²、12 078 km²，对比发现本文提取的海冰面积较为准确。从海冰类型来看，2021年1月16日的固定冰、白冰、灰（白）冰3种类型面积大于另外两个日期的面积；2021年1月10日的初生冰面积最大，说明随着气温降低，此阶段海水结冰迅速，出现大量初生冰；到2021年1月16日，大量初生冰转为灰（白）冰、白冰或者固定冰，从而初生冰面积减少，其他海冰类型面积增加。

4 结论

收起

本文基于HY-1C CZI影像，提出了一种有效的海冰分类方法。该方法通过提取HY-1C CZI影像的光谱特征和纹理特征，对多特征进行特征选择，得到最优特征集，基于最优特征集开展SVM海冰分类研究。主要结论如下：

（1）在特征提取与选择方面，选择NDWI_H、B_G、B_R、G_R 4个光谱特征与均值、同质性、熵、二阶矩4个纹理特征组成最优特征集，既保留了分类有效信息，实现信息互补，又降低了信息冗余。

（2）在海冰分类方面，本文提出的基于最优特征集的分类方法优于仅基于光谱特征或纹理特征的分类方法。利用本文方法的3期HY-1C CZI遥感影像分类精度均在80%以上。对海冰分类图进一步估算得到海冰总面积，与国家卫星海洋应用中心通报的海冰面积相差较小，进一步验证了本文方法的有效性。

综上所述，本文提出的基于最优特征集的SVM海冰分类方法是有效的，同时证明了HY-1C CZI数据在海冰分类应用中的可行性，但在特征提取方面，本文仅应用了光谱和纹理特征，更多的特征提取方法有待进一步研究。

基金

收起

国家重点研发计划(2018YFB054900)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

崔艳荣, 邹斌, 韩震, 等. 卷积神经网络在SAR遥感海冰分类中的应用可行性分析[J]. 海洋预报, 2019, 36(5): 77−85.

Cui Yanrong, Zou Bin, Han Zhen, et al. Feasibility analysis of convolutional neural networks in remote sensing sea ice classification[J]. Marine Forecasts, 2019, 36(5): 77−85.

刘眉洁, 戴永寿, 张杰, 等. 高分辨率全极化合成孔径雷达数据海冰二次分类方法研究[J]. 海洋学报, 2013, 35(4): 80−87.

Liu Meijie, Dai Yongshou, Zhang Jie, et al. The research on the object-based method of sea ice classification of high-resolution quad-polarization SAR data[J]. Haiyang Xuebao, 2013, 35(4): 80−87.

张晰. 极化SAR渤海海冰厚度探测研究[D]. 青岛: 中国海洋大学, 2011.

Zhang Xi. Research on sea ice thickness detection by polarimetric SAR in Bohai Sea[D]. Qingdao: Ocean University of China, 2011.

Tan Wenxia, LeDrew E. Monitoring Arctic sea ice phenology change using hypertemporal remotely sensed data: 1989–2010[J]. Theoretical and Applied Climatology, 2016, 125(1/2): 353−363.

Geldsetzer T, Arkett M, Zagon T, et al. All-season compact-polarimetry C-band SAR observations of sea ice[J]. Canadian Journal of Remote Sensing, 2015, 41(5): 485−504.

Scheuchl B, Caves R, Cumming I, et al. Automated sea ice classification using spaceborne polarimetric SAR data[C]//IGARSS 2001. Scanning the Present and Resolving the Future. Proceedings. IEEE 2001 International Geoscience and Remote Sensing Symposium. Sydney: IEEE, 2001: 3117−3119.

Singha S, Johansson M, Hughes N, et al. Arctic sea ice characterization using spaceborne fully polarimetric L-, C-, and X-Band SAR with validation by airborne measurements[J]. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(7): 3715−3734.

刘惠颖. 宽幅多极化SAR海冰信息提取方法与类型识别研究[D]. 北京: 中国科学院遥感与数字地球研究所, 2017.

Liu Huiying. Study on parameter retrieval and type recognition of sea ice using wide-swath and multi-polarization SAR data[D]. Beijing: Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, 2017.

赵泉华, 郭世波, 李晓丽, 等. 利用目标分解特征的全极化SAR海冰分类[J]. 测绘学报, 2018, 47(12): 1609−1620.

Zhao Quanhua, Guo Shibo, Li Xiaoli, et al. Polarimetric SAR sea ice classification based on target decompositional features[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(12): 1609−1620.

Riggs G A, Hall D K, Ackerman S A. Sea ice extent and classification mapping with the moderate resolution imaging spectroradiometer airborne simulator[J]. Remote Sensing of Environment, 1999, 68(2): 152−163.

史培军, 范一大, 哈斯, 等. 利用AVHRR和MODIS数据测算海冰资源量——以渤海海冰资源测算为例[J]. 自然资源学报, 2002, 17(2): 138−143.

Shi Peijun, Fan Yida, Ha Si, et al. Calculating gross sea ice resource using AVHRR and MODIS data[J]. Journal of Natural Resources, 2002, 17(2): 138−143.

吴龙涛, 吴辉碇, 孙兰涛, 等. MODIS渤海海冰遥感资料反演[J]. 中国海洋大学学报(自然科学版), 2006, 36(2): 173−179.

Wu Longtao, Wu Huiding, Sun Lantao, et al. Retrieval of sea ice in the Bohai Sea from MODIS data[J]. Periodical of Ocean University of China, 2006, 36(2): 173−179.

王姝力, 王志勇, 王磊. 基于Landsat-8和Sentinel-1A辽东湾海冰分类研究[J]. 北京测绘, 2019, 33(12): 1486−1492.

Wang Shuli, Wang Zhiyong, Wang Lei. Study of sea ice classification of Landsat-8 and Sentinel-1A in Liaodong Bay[J]. Beijing Surveying and Mapping, 2019, 33(12): 1486−1492.

张晰, 张杰, 孟俊敏. Landsat-8与GF-1卫星渤海海冰探测能力对比研究[J]. 海洋科学, 2015, 39(2): 50−56.

Zhang Xi, Zhang Jie, Meng Junmin. Comparison of sea ice detection ability of Landsat-8 and GF-1 in the Bohai Sea[J]. Marine Sciences, 2015, 39(2): 50−56.

Barbieux K, Charitsi A, Merminod B. Icy lakes extraction and water-ice classification using Landsat-8 OLI multispectral data[J]. International Journal of Remote Sensing, 2018, 39(11): 3646−3678.

王志勇, 王丽华, 刘健, 等. 基于多源中高分辨率遥感数据提取渤海辽东湾海冰要素信息[J]. 自然灾害学报, 2021, 30(1): 174−182.

Wang Zhiyong, Wang Lihua, Liu Jian, et al. Extraction of sea ice element information in Liaodong Bay of the Bohai Sea based on multi-source medium and high resolution remote sensing data[J]. Journal of Natural Disasters, 2021, 30(1): 174−182.

刘眉洁. 基于高分辨率极化SAR的海冰分类和厚度探测方法研究[D]. 青岛: 中国石油大学(华东), 2016.

Liu Meijie. Research on the sea ice classification and thickness detection with high-resolution and polarimetric SAR data[D]. Qingdao: China University of Petroleum (East China), 2016.

庞海洋, 孔祥生, 孙志伟, 等. 基于遥感和气象数据对辽东湾海冰变化预测研究[J]. 海洋与湖沼, 2018, 49(4): 725−733.

Pang Haiyang, Kong Xiangsheng, Sun Zhiwei, et al. The forecast model of sea ice changes in Liaodong Bay using remote sensing and meteorological data[J]. Oceanologia et Limnologia Sinica, 2018, 49(4): 725−733.

罗会兰, 郭敏杰, 孔繁胜. 集成多特征与稀疏编码的图像分类方法[J]. 模式识别与人工智能, 2014, 27(4): 345−355.

Luo Huilan, Guo Minjie, Kong Fansheng. Image classification method by combining multi-features and sparse coding[J]. Pattern Recognition and Artificial Intelligence, 2014, 27(4): 345−355.

Haq M A, Jain K, Menon K P R. Monitoring glacial lakes using remote sensing Tech[C]//Proceedings of the 14th Annual International Conference and Exhibition on Geospatial Information Technology and Applications. Gurgaon: India Geospatial Forum, 2012.

张明, 吕晓琪, 张晓峰, 等. 结合纹理特征的SVM海冰分类方法研究[J]. 海洋学报, 2018, 40(11): 149−156.

Zhang Ming, Lü Xiaoqi, Zhang Xiaofeng, et al. Research on SVM sea ice classification based on texture features[J]. Haiyang Xuebao, 2018, 40(11): 149−156.

逯跃锋, 和鑫, 陆黎娟, 等. 基于纹理分析的SAR海冰图像分类方法[J]. 山东理工大学学报(自然科学版), 2019, 33(1): 51−55.

Lu Yuefeng, He Xin, Lu Lijuan, et al. Research on classification method of the SAR sea ice image based on texture analysis[J]. Journal of Shandong University of Technology (Natural Science Edition), 2019, 33(1): 51−55.

Kiala Z, Mutanga O, Odindi J, et al. Feature selection on Sentinel-2 multispectral imagery for mapping a landscape infested by parthenium weed[J]. Remote Sensing, 2019, 11(16): 1892.

边肇祺, 张学工. 模式识别[M]. 北京: 清华大学出版社, 2004: 9-42.

Bian Zhaoqi, Zhang Xuegong. Pattern Recognition[M]. Beijing: Tsinghua University Press, 2004: 9−42.

Bogdanov A V, Sandven S, Johannessen O M, et al. Multisensor approach to automated classification of sea ice image data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(7): 1648−1664.

沈杨. 结合MRF与ν-SVM的SAR海冰图像分类[D]. 合肥: 合肥工业大学, 2015.

Shen Yang. Combining MRF and v-SVM for SAR sea ice image classification[D]. Hefei: Hefei University of Technology, 2015.

2022年第44卷第5期

PDF下载

135

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2022021

接收时间：2021-03-03
首发时间：2026-02-01
出版时间：2022-05-01

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2021-03-03
修回日期：2021-07-30

基金

国家重点研发计划(2018YFB054900)

作者信息

1.航天宏图信息技术股份有限公司，北京 100195

2.国家卫星海洋应用中心，北京 100081

通讯作者:

刘建强（1964-），男，湖南省益阳市人，研究员，主要从事海洋遥感应用研究与地面应用系统建设。E-mail: jqliu@mail.nsoas.org.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2022021

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

传感器	波段名称	波长/nm	分辨率/m
HY-1C卫星CZI	蓝（波段1）	420～500	50
绿（波段2）	520～600	50
红（波段3）	610～690	50
近红外（波段4）	760～890	50

传感器

波段名称

波长/nm

分辨率/m

HY-1C卫星CZI

蓝（波段1）

420～500

绿（波段2）

520～600

红（波段3）

610～690

近红外（波段4）

760～890

海冰类型	特征描述
固定冰	固定冰是指与海岸、岛屿或海底冻结在一起，不能做水平运动，而能随海面升降做垂直运动的海冰。固定冰的封冻程度和范围对沿岸水产养殖业等有较大影响，固定冰的挤压对岸基和堤坝等建筑可能造成破坏。
白冰	当灰（白）冰厚度继续增加时，冰块颜色由灰白色过渡到白色，表面凹凸不平，形状不规则，成为白冰。厚度大于30 cm。
灰（白）冰	由于灰冰和灰白冰边界模糊，将两类作为一类进行区分。灰冰厚度为 10～15 cm ，灰白冰厚度为15～30 cm，因此，两类合为一类后，灰（白）冰厚度为10～30 cm，灰（白）冰表面比较粗糙，呈灰色或灰白色。
初生冰	当海上气温下降到海水冰点，海水开始结冰形成初生冰，多呈针状、薄片状、浆糊状或绵状。厚度在10 cm以内。

海冰类型

特征描述

固定冰

白冰

当灰（白）冰厚度继续增加时，冰块颜色由灰白色过渡到白色，表面凹凸不平，形状不规则，成为白冰。厚度大于30 cm。

灰（白）冰

初生冰

当海上气温下降到海水冰点，海水开始结冰形成初生冰，多呈针状、薄片状、浆糊状或绵状。厚度在10 cm以内。

特征标识	定义公式
注：BLUE、GREEN、RED、NIR分别表示HY-1C卫星CZI数据的蓝、绿、红、近红外波段。
NDWI_H	$\rm \left( {NIR - BLUE} \right)/(NIR + BLUE) $
B_G	$\rm (BLUE - GREEN)/(BLUE + GREEN) $
B_R	$ \rm (BLUE - RED)/(BLUE + RED) $
G_R	$\rm (GREEN - RED)/(GREEN + RED) $

特征标识

定义公式

注：BLUE、GREEN、RED、NIR分别表示HY-1C卫星CZI数据的蓝、绿、红、近红外波段。

NDWI_H

$\rm \left( {NIR - BLUE} \right)/(NIR + BLUE) $

B_G

$\rm (BLUE - GREEN)/(BLUE + GREEN) $

B_R

$ \rm (BLUE - RED)/(BLUE + RED) $

G_R

$\rm (GREEN - RED)/(GREEN + RED) $

	波段1	波段2	波段3	波段4
均值	5.687 4	3.56	5.234 4	2.27
方差	0.049 2	0.047 5	0.045 9	0.044 9
同质性	2.927 5	2.152 5	2.559	1.66
对比度	0.056 5	0.052	0.051	0.048
相异性	0.807 9	0.776	0.785 4	0.67
熵	4.003	2.354 1	3.293 5	1.69
二阶矩	3.540 8	1.897 6	2.76	1.41
相关性	0.570 3	0.297 1	0.456 2	0.261 5

波段1

波段2

波段3

波段4

均值

5.687 4

3.56

5.234 4

2.27

方差

0.049 2

0.047 5

0.045 9

0.044 9

同质性

2.927 5

2.152 5

2.559

1.66

对比度

0.056 5

0.052

0.051

0.048

相异性

0.807 9

0.776

0.785 4

0.67

熵

4.003

2.354 1

3.293 5

1.69

二阶矩

3.540 8

1.897 6

2.76

1.41