基于机载高光谱成像技术的海洋生态监测 - 海洋生态监测仪器

" U+ |' b" w3 ~

本期作者：研发部应用工程师 吴亚蓝

( v. E/ o7 [ {- a3 n

应用关键词：海洋生态、水质监测

1 引言

高光谱成像作为一种新型的光谱成像技术，在海洋生态监测领域中充分展现了其自身的技术优势。由于高光谱成像光谱仪具有光谱覆盖范围广、分辨率高和波段多、图谱合一等优点，因此已成为海洋总氮、总磷等的有效监测手段。它既可用于海水中高锰酸盐浓度、悬浮泥沙含量、某些污染物和表层水温探测，还可用于海冰、海岸带等地表的探测。

- y5 k6 X- w& ]* t% U7 ? t

国内海洋遥感应用基础研究主要是一些数学模型的构建。在海洋水质监测应用方面，只有可见光波段的光谱能够观测水下的状况。另外，陆源污染、海水养殖、滩涂等海岸带典型要素的光谱特性研究工作也在开展，研究人员以航空/卫星的高光谱图像为数据源，选取陆源污染、海水养殖、滩涂为监测要素，进行上述要素的光谱波段敏感性研究，获得其探测的最佳波段，并进一步发展准确、快速识别和探测技术。在海洋表面温度测量、海洋表层悬浮泥沙浓度的定性或半定量的观测、海洋动力现象的研究等方面都开展了相应的研究。通过水质遥感监测分析水体反射光谱特征与水质参数浓度之间的关系，建立水质参数的反演算法。这是一种不仅省时省力，而且还可以利用卫星影像进行全面、及时的水质监测的方法。

6 x+ k0 b0 E; H6 \+ W d' X

2 材料与试验部分

! Y6 @5 C6 X, W) @8 b

2.1 研究区域

野外试验飞行总共为两次，均在中国的某段海面上进行。

2.2 数据采集设备

2 z0 m: x8 N' P D

本次试验采用大疆六旋翼无人机M600 Pro(无人机净重约10 kg，最大载重约 6.5 kg)，在无人机遥感平台上搭载江苏双利合谱科技有限公司自主研发的高光谱成像光谱仪GaiaSky-mini-2，该成像系统采用拥有自主知识产权的内置推扫方式获取地面图像(其主要参数见表1)。无人机的飞行高度均为300m，分别采集了33景和37景图像，无人机遥感平台的实景图如图1所示。

表1 GaiaSky-mini2 无人机载成像高光谱仪系统参数

7 G# I/ o6 b( x- I" n8 G

序号

项目

5 m: h i. T' _/ \- n2 }! J: {- Y

参数

1

, f! u% z) ~+ w7 H, Y! g

光谱扫描范围/nm

6 c& S+ l7 C# }9 i, d

400~1000

: s" Y# a5 |' B, m2 Z1 W

2

光谱分辨率/nm

3.5 nm

3 j& W- o( E3 F

3

( i# q! G- J! F: U1 A

成像镜头/mm

9 B% e4 j, }9 l

18.5

4

5 C4 F2 S; C& n" ~

光谱通道数

/ V( ^1 [+ l/ I- F: _ v

360

5

: j4 V2 m4 h8 u9 o2 ^

全幅像素

1936×1456

+ C, }5 M; ^. B1 P- o

6

传感器

CCD Sony ICX 674

% P& W7 Q1 m, y

8 B+ Y8 O; o1 m+ V: l1 K: d9 H0 T2 k0 z

& J& k# X" M% Z6 q# F8 `2 w% Y. S" t

图1 无人机高光谱成像系统实景图

0 |$ R5 n3 m8 {3 g9 a5 r4 e+ f: B! G

& `( z( I9 I# k E7 ` ^, T

2.3 无人机高光谱数据的预处理与分析

" \5 e4 p! G. F4 G% c

无人机高光谱图像的预处理在江苏双利合谱科技有限公司自主研发的SpecView软件中进行，包括镜像变换、黑白帧校准（如公式1所示）。

a% I( n8 Z% W7 C' M# J

7 }* ^' v) M# o( t9 A% u+ l" c

（1）

. p( H, _; r1 o9 E) J

式中，Rref是黑白校正过的图像的反射率值，DNraw 是原始图像的DN值，DNwhite为白板的白帧数据，DNdark 是相机的系统误差DN值。

. k+ W k! E, Q' B V- Z, C9 ~

考虑到无人机飞到一定高度后，高光谱成像仪获取的高光谱影像数据可能会受到大气、水汽等因素的影响。为了消除这些因素的影响，我们（系统）在无人机起飞之前，在拍摄区域放置一块经过国家计量院标定过的2m*2m灰布，在高光谱影像获取的时候，只需要在其中的一景高光谱影像中覆盖到灰布即可。消除大气、水汽等因素影响的方法如公式2所示。

（2）

+ z5 p% K3 P' t4 p: V# T

! f0 K, u1 }. L: p& W( [

式中，Rfixed是消除大气、水汽等因素后的图像光谱反射率，Rref是经过黑白校正后的图像反射率，Rstandard是经过国家计量院标定的灰布的光谱反射率，Rgrayref是经过黑白校正后图像中灰布的光谱反射率。

- L7 y5 i) k; r1 j9 `; s

2.4 无人机高光谱影像拼接

* f+ P5 W' _+ a3 ^

无人机高光谱影像的拼接采用江苏双利合谱科技有限公司自主研发的无人机高光谱拼接软件SpecStitcher进行拼接，该拼接软件有图像筛选、拼接预览、投影方式选择、拼接方法选择、重采样方法选择、是否匀色、拼接结果格式输出选择等功能。拼接软件界面如图2所示。

图2 无人机高光谱拼接软件SpecStitcher

. L$ i0 `$ I2 o) d

3 结论与分析

3.1 拼接结果预览

图3.1（左）及图3.2（左）为利用无人机高光谱拼接软件SpecStitcher对海面及其海岸的无人机高光谱影像的三波段拼接效果预览图（RGB分别代表640 nm/550 nm/460 nm最邻近波长）。从拼接结果来看，第一次与第二次采集到的海面影像数据拼接效果较好。提取出陆地和海水的光谱曲线如图3所示。

A4 M! O; `0 l( z/ t

图3.1第一次数据无人机高光谱拼接预览图

% C7 X7 ^0 t4 i

# ]0 G) ~, d. D! x1 u1 M: O9 v

5 U/ @" |# U) m! N1 u

图3.1第一次海水数据提取

! F4 o0 c# @* u. S

图3.2第二次数据无人机高光谱拼接预览图

3 ]* S5 r+ y( N- Q7 L4 J

图3.2第二次海水数据提取

: F3 T1 ~7 D1 [& |

a6 u( p5 n% i/ _: W; Z9 u. h

图4 海水与陆地光谱对比

9 r9 w0 e' Y; d, o0 @% [1 N9 C6 J

3.2 无人机高光谱影像河流提取方法

无人机高光谱影像不仅包含有河流，还有其他的树木、杂草、土壤、道路、建筑物等，从影像中提取感兴趣的目标物，通常的方法有非监督分类、监督分类、决策树、面向对象分类等，本研究采用监督分类的马氏距离法对拼接好的无人机高光谱影像进行分类，从而提取了两次数据中的海水信息，并对第二次海水进行可见总氧、总磷和高锰酸盐等水质参数的反演。两次飞行提取出的海水信息如图3.1（右）及图3.2（右）所示。

+ L! x1 [, f8 l( `

3.3 总氮、总磷和高锰酸盐的反演

水体中总氮、总磷和高锰酸盐含量是衡量水质的重要指标。常规氮磷和高锰酸盐的测量方法需要长时间的高温、高压消解，且消解的温度、时间和试剂对测定的结果均有较大的影响，整个操作烦琐、费时和耗力；因此本研究试图利用高光谱遥感技术，通过对水体中氮、磷和高锰酸盐光谱的测定，探索水体氮、磷和高锰酸盐与反射光谱特征的关系，建立氮、磷浓度和高锰酸盐的反演模型，为湖泊、水库和河流等大型内陆水体氮磷及高锰酸盐遥感定量监测提供理论依据。目前的一些研究仅根据总氮、总磷、高锰酸盐与海水水质之间具有的密切相关关系，建立总氮、总磷和高锰酸盐的遥感信息模型。

以第二次无人机飞行数据为例，利用机载高光谱影像反演海水的总磷浓度、总氮和高锰酸盐浓度的分布图，如图5所示。

: p( @$ i) c; {4 c e1 p$ y

' y" e2 B j n* C7 ]

图5 海水水质衡量反演图—(a) 总磷

) d: m' v4 U6 l" Y

图5 海水水质衡量反演图—(b) 总氮

4 L' ~( x# |- Y- K4 `( N2 V/ U

图5 海水水质衡量反演图—(c)高锰酸盐

4 X1 Z0 L# e1 t9 ]9 [% }+ r* }

4 结论与讨论

7 T l3 h7 y0 B! D: l3 p

本文以中国某段海面作为研究对象，利用无人机搭载高光谱成像仪获取的海面光谱图像信息，采用拼接软件对获取的影像进行拼接以得到一段完整的海面及海滩的高光谱影像数据；基于海水的固有特性，提取了两次无人机高光谱数据中的海水信息；进而根据现有报道的地面尺度研究海洋、湖泊富营养以及水污染的高光谱数据模型，选择了稳定性较好的经典模型，对两段海面进行富营养（总氮、总磷和高锰酸盐）的反演研究。结合现场调研以及高光谱影像反演效果来看，本实验可以很好地反映出海面的富营养和受污染程度。该研究为水利、水电部门和环保部门提供了一种监测海水富营养及受污染的分布情况的有效方法，为相关技术人员判断污染源以及寻找水体富营养的原因提供了有效的技术手段。

% C) m: Y, |3 Z+ L

如您对高光谱遥感在农业、水体、自然资源等多领域的应用感兴趣，请关注双利合谱公众微信号！或者拨打13520222630电话咨询！

* K8 r3 _, N7 D% C$ Y ?