一种基于高光谱遥感影像的近海养殖筏提取指数的提取方法与流程

1.本发明属于遥感影像的分类提取技术领域，尤其涉及一种基于高光谱遥感影像的近海养殖筏提取指数。


背景技术：

2.根据2016年全球水产养殖产量统计，中国的水产养殖产量占世界总产量的一半以上，是世界上最大的水产养殖国家，而筏式养殖是海水养殖中最重要的方式之一。但是近些年来为追求经济利益，沿海水产养殖产业迅速扩大导致沿海环境遭受巨大压力，养殖过程中营养物质和抗生素的过度供应导致富营养化、水污染、赤潮等环境问题。因此，对沿海养殖渔业进行遥感监测，系统而深入地研究浮筏养殖用海面积及分布，了解养殖区的空间分布及发展，为渔业管理部门合理规划养殖用海提供辅助决策信息，对于海岸带的科学管理和可持续发展至关重要。
3.近年来随着对近海养殖区观测资料的日益稳定和丰富，对养殖区的养殖范围和养殖密度，甚至于对养殖类型进行识别和提取等日益成为海洋遥感研究中的热点问题。目前传统浮筏养殖的识别和提取常用的方法大致可分为三类：人工目视解译，面向对象分类，机器学习分类等。人工目视解译工作量大，成本高，操作繁琐；面向对象分类需要同时利用光谱信息和形状信息且对数据的分辨率要求高，普适性低，不易实现；机器学习分类的方法需对特定研究区的养殖筏选取大量样本进行训练，实现难度大，缺乏可迁移性。
4.研究常用的数据源为中高分辨率的多光谱遥感数据，其中最普遍的为landsat数据，除了中分辨率卫星遥感影像，spot、ikonos以及我国近年发射的高分系列卫星等高分辨率影像也是该研究方向重要的数据来源。但是近海养殖浮筏与背景水体光谱曲线极为接近，仅在可见光的红黄光波段存在明显的差异。目前研究常用的多光谱数据波段数量很少，依靠宽波段信息较难区分养殖筏与背景水体。高光谱数据光谱分辨率高，能够区分地物的细微差异特征，可以充分体现出养殖筏与背景海水的光谱差异。


技术实现要素：

5.本发明的目的是克服现有技术中的不足，提供一种基于高光谱遥感影像的近海养殖筏提取指数的提取方法。
6.这种基于高光谱遥感影像的近海养殖筏提取指数的提取方法，包括以下步骤：
7.步骤1、对高光谱数据进行预处理，预处理包括正射校正、辐射定标、大气校正和影像裁剪；
8.步骤2、分析地物光谱特征，构建高光谱遥感影像的近海养殖筏提取指数，根据高光谱遥感影像的近海养殖筏提取指数的阈值来区分养殖筏、背景水体、植被和其他地物；
9.步骤2.1、分析高光谱遥感影像中养殖筏、背景水体、植被、其他地物的均值光谱曲线，得到高光谱遥感影像的高光谱信息；根据高光谱信息构建养殖筏指数ari，通过等光谱间波段的差值运算得到背景地物信息；其中养殖筏指数ari为：
10.ari＝(ρ
593
‑
ρ
955
) (ρ
1040
‑
ρ
1240
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
11.上式中，ari为养殖筏指数，ρ
593
表示中心波长为593nm的波段值，ρ
955
表示中心波长为955nm的波段值，ρ
1040
表示中心波长为1040nm的波段值，ρ
1240
表示中心波长为1240nm的波段值；ρ
593
‑
ρ
955
得到的植被波段值为负值，表示ρ
593
处植被的反射率低于ρ
955
处植被的反射率，ρ
593
处养殖筏的反射率与背景水体的反射率均高于ρ
955
处养殖筏的反射率与背景水体的反射率，且ρ
593
‑
ρ
955
得到的植被波段值总是小于背景水体波段值；将植被从图像中分离开，分离背景水体和养殖筏；
12.步骤2.2、通过envi中band math工具来计算步骤2.1构建的养殖筏指数；
13.步骤2.3、统计各类地物样本养殖筏指数的计算结果，生成箱型图；根据养殖筏提取指数的阈值来提取养殖筏：
14.min
m
≤ari(m,n)≤max
m
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
15.上式中，ari(m,n)为位置(m,n)的养殖筏指数的计算结果，min
m
、max
m
分别为养殖筏提取指数的阈值下限和阈值上限；若位置(m,n)处的养殖筏指数ari(m,n)满足上式(2)，则执行步骤3，若不满足上式(2)，则判定位于(m,n)的地物为非养殖筏。
16.步骤3、计算养殖筏提取结果的面积信息，过滤面积小于设定值的误分像素。
17.作为优选，步骤3具体包括以下步骤：
18.步骤3.1、将提取到的养殖筏对象转为矢量数据，在arcmap软件中计算提取到的养殖筏对象的面积信息，通过过滤掉面积小于设定值的提取要素来去除误分像素；
19.步骤3.2、若提取要素i的面积满足下式(3)，则该提取要素为养殖筏，若不满足下式(3)，则判定该提取要素为误分像素；
20.s
i
>min
m1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
21.上式中，s
i
为提取要素i的面积；min
m1
为养殖筏面积的最小下限。
22.作为优选，步骤2.3中养殖筏提取指数的阈值下限和阈值上限分别为0.8和1.1。
23.作为优选，步骤3.2中养殖筏面积的最小下限为0.9平方米。
24.本发明的有益效果是：
25.本发明针对近海浮筏养殖遥感信息自动提取应用难度较大的问题，简化现有养殖筏提取的流程与难度，分析高光谱影像上养殖筏的光谱信息，通过增大其与其他地物的光谱差距，构建一种基于高光谱数据的近海水产养殖筏提取指数，增大养殖筏与其他地物的差异根据养殖筏指数结果使用阈值分割提取养殖筏，最后结合面积信息去除误分像素；实现了近海养殖筏的简单、快速、精确提取。
26.本发明构建了养殖筏提取指数来增大养殖筏与其他地物的差异，并根据养殖筏指数结果使用阈值分割提取养殖筏，最后结合面积信息去除误分像素，实现近海养殖筏的简单、快速、精确提取。本发明流程容易实现，是对现有养殖筏提取方法的有效补充，有利于提高养殖筏提取的准确性、反映近海养殖筏的真实空间分布情况，能为养殖筏的管理与保护提供科学依据，具有十分重要的实际应用意义。
27.本发明提出的一种基于高光谱遥感影像的近海养殖筏提取指数，充分利用了高光谱数据“图谱合一”的优点，简化了近海养殖筏提取的难度与流程，较传统分类方法具有明显的优越性。ρ
1040
‑
ρ
1240
的作用是抑制与水和养殖筏光谱曲线相似的潮滩、湿地的指数值，。
附图说明
28.图1为本发明实施例的流程图；
29.图2是资源一号02d高光谱影像中养殖筏、背景水体、植被、其他地物的均值光谱曲线图；
30.图3为各类样本进行养殖筏提取指数(ari)的计算结果统计生成的箱线图；
31.图4
‑
1为本发明实施例中温州市苍南县近海养殖区的真彩色合成原始高光谱影像，图4
‑
2为本发明实施例中养殖筏最终提取结果图。
具体实施方式
32.下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
33.实施例一
34.本技术实施例一提供了一种基于高光谱遥感影像的近海养殖筏提取指数的提取方法，包括以下步骤：
35.步骤1、对高光谱数据进行预处理，预处理包括正射校正、辐射定标、大气校正和影像裁剪；
36.步骤2、分析地物光谱特征，构建高光谱遥感影像的近海养殖筏提取指数，根据高光谱遥感影像的近海养殖筏提取指数的阈值来区分养殖筏、背景水体、植被和其他地物；
37.步骤2.1、分析高光谱遥感影像中养殖筏、背景水体、植被、其他地物的均值光谱曲线，得到高光谱遥感影像的高光谱信息；根据高光谱信息构建养殖筏指数(aquaculture raft index，以下简称ari)，通过等光谱间波段的差值运算得到抑制水体等背景地物信息，增强养殖筏与背景地物间的差异，从而突出筏式养殖区的信息特征；其中养殖筏指数ari为：
38.ari＝(ρ
593
‑
ρ
955
) (ρ
1040
‑
ρ
1240
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
39.上式中，ari为养殖筏指数，ρ
593
表示中心波长为593nm的波段值，ρ
955
表示中心波长为955nm的波段值，ρ
1040
表示中心波长为1040nm的波段值，ρ
1240
表示中心波长为1240nm的波段值；ρ
593
‑
ρ
955
得到的植被波段值为负值，表示ρ
593
处植被的反射率低于ρ
955
处植被的反射率，因为近海养殖悬浮泥沙的存在，ρ
593
处养殖筏的反射率与背景水体的反射率均高于ρ
955
处养殖筏的反射率与背景水体的反射率，且ρ
593
‑
ρ
955
得到的植被波段值总是小于背景水体波段值；将植被从图像中分离开，分离背景水体和养殖筏；
40.步骤2.2、通过envi中band math工具来计算步骤2.1构建的养殖筏指数；
41.步骤2.3、统计各类地物样本养殖筏指数的计算结果，生成箱型图；ρ
1040
‑
ρ
1240
的作用是抑制与水和养殖筏光谱曲线相似的潮滩、湿地的指数值，减小ρ
1040
‑
ρ
1240
指数计算值；通过构建养殖筏指数(ari)增强了养殖筏与其他地物的区别，这样即可根据ari
植被
<ari
其他地物
<air
养殖筏
<ari水体的关系，然后确定养殖筏提取指数的阈值，其中ari
植被
为高光谱遥感影像中植被对应的养殖筏指数，ari
其他地物
为高光谱遥感影像中其他地物对应的养殖筏指数，air
养殖筏
为高光谱遥感影像中养殖筏对应的养殖筏指数，ari
水体
为高光谱遥感影像中背景水体对应的养殖筏指数；根据养殖筏提取指数的阈值来简单快速地提取养殖筏：
42.min
m
≤ari(m,n)≤max
m
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
43.上式中，ari(m,n)为位置(m,n)的养殖筏指数的计算结果，min
m
、max
m
分别为养殖筏提取指数的阈值下限和阈值上限；若位置(m,n)处的养殖筏指数ari(m,n)满足上式(2)，则执行步骤3，若不满足上式(2)，则判定位于(m,n)的地物为非养殖筏。
44.步骤3、计算养殖筏提取结果的面积信息，过滤少量面积小于设定值的误分像素。
45.步骤3.1、将提取到的养殖筏对象转为矢量数据，在arcmap软件中计算提取到的养殖筏对象的面积信息，通过过滤掉面积小于设定值的提取要素来去除误分像素；
46.步骤3.2、若提取要素i的面积满足下式(3)，则该提取要素为养殖筏，若不满足下式(3)，则判定该提取要素为误分像素；
47.s
i
>min
m1
ꢀꢀꢀꢀꢀꢀ
(3)
48.上式中，s
i
为提取要素i的面积；min
m1
为养殖筏面积的最小下限。
49.实施例二
50.随着我国高光谱卫星的发射，资源1号02d等高光谱卫星可以同时采集地物的光谱和图像信息，实现了“图谱合一”，大大提升了数据应用潜力；本技术实施例二在实施例一的基础上针对目前使用多光谱数据进行近海养殖筏提取时存在的“异物同谱”问题，提出一种基于高光谱数据的近海养殖筏提取指数，实现近海养殖筏简单快速的提取。提供了一种基于高光谱数据的近海养殖筏提取指数，以资源一号02d高光谱数据为数据源，通过envi5.3、arcmap10.2软件实现，其提取流程如图1所示：
51.步骤1、对资源一号02d的高光谱数据进行预处理，该预处理主要包括辐射定标、正射校正、大气校正和影像裁剪；
52.步骤2、分析地物光谱特征构建特征指数(高光谱遥感影像的近海养殖筏提取指数)，根据养殖筏提取指数的阈值来区分养殖筏、背景水体、植被和其他地物；
53.步骤2.1、分析如图2所示高光谱影像中养殖筏、背景水体、植被、其他地物的均值光谱曲线，分析各地物光谱关系；针对影像的高光谱信息，构建养殖筏指数(aquaculture raft index，以下简称ari)，通过等光谱间波段的差值运算抑制水体等背景地物信息，增强养殖筏与背景地物间的差异，从而突出筏式养殖区的信息特征。
54.ari＝(ρ
593
‑
ρ
955
) (ρ
1040
‑
ρ
1240
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
55.上式中，ari为本实施例所构建的养殖筏指数，ρ
593
表示中心波长为593nm的波段值(即zy1
‑
02d影像可见光近红外数据中band 23的波段值)，ρ
955
表示中心波长为955nm的波段值(即zy1
‑
02d影像可见光近红外数据中band 72的波段值)，ρ
955
表示中心波长为955nm的波段值，ρ
1040
表示中心波长为1040nm的波段值(即zy1
‑
02d影像可见光近红外数据中band 76的波段值)，ρ
1240
表示中心波长为1240nm的波段值(即zy1
‑
02d影像短波红外数据中band15的波段值)。ρ
593
处植被的反射率低于ρ
955
处的反射率，ρ
593
‑
ρ
955
得到的植被的值即为负值，从而将植被从图像中分离开，因为近海养殖悬浮泥沙的存在ρ
593
处养殖筏与背景水体的反射率高于ρ
955
处的反射率，且ρ
593
‑
ρ
955
的结果值养殖筏总是小于背景水体，以此分离背景水体和养殖筏。ρ
1040
‑
ρ
1240
的作用是抑制与水和养殖筏光谱曲线相似的潮滩、湿地的指数值，以此将其指数计算值减小。通过构建养殖筏指数(ari)增强了养殖筏与其他地物的区别；这样即可根据ari
植被
<ari
其他地物
<air
养殖筏
<ari
水体
的关系，设定阈值来提取养殖筏。
56.步骤2.2、通过envi中band math工具计算所构建的养殖筏指数。
57.步骤2.3、根据各类地物样本的养殖筏指数(ari)计算结果来统计生成如图3所示箱型图，结果显示出了本发明指数对养殖筏、水体、植被、其他地物具有良好的可分离性；然后确定养殖筏提取阈值，根据阈值分割来简单快速地提取养殖筏：
58.0.8≤ari(m,n)≤1.1
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
59.上式中，ari(m,n)为位置(m,n)的养殖筏指数的计算结果，经计算统计，养殖筏的指数阈值下限为0.8和养殖筏的指数阈值上限为1.1；若位置(m,n)处的养殖筏指数ari(m,n)满足上式(2)，则执行步骤3，若不满足上式(2)，则判定位置(m,n)所在的区域的植被位于(m,n)的地物为非养殖筏。
60.步骤3、图4
‑
1为本发明实施例中温州市苍南县近海养殖区的真彩色合成原始高光谱影像，图4
‑
2为本发明实施例中养殖筏最终提取结果图；计算如图4
‑
2提取结果面积信息，过滤掉少量小面积误分像素。
61.步骤3.1、将提取到的养殖筏对象转为矢量数据，在arcmap软件中计算其面积信息，通过过滤掉面积过小的提取要素来去除误提像素。
62.s
i
>0.9(m2)
ꢀꢀꢀꢀꢀ
(3)
63.上式中，s
i
为要素i的面积；0.9平方米为养殖筏面积的最小下限。
64.步骤3.2、若要素i的面积大于0.9平方米，该要素为养殖筏，若不满足，则判定该要素为误提像素。
65.图4
‑
2为本实施例的近海养殖筏提取结果图，通过本发明提出的指数简单、快速、精确地提取出了温州市苍南县的近海养殖浮筏，并且去除了绝大部分其他地物的像素被误分的情况，很好地反映了近海养殖筏的真实空间分布情况，有利于近海养殖浮筏的监测与管理。