一种基于珞珈一号夜光遥感数据的夜间海洋船舶检测方法

摘要

针对目前夜间海洋船舶检测方法的不足，提出了一种基于珞珈一号夜光遥感数据的夜间海洋船舶检测方法。利用珞珈一号夜光遥感数据，通过采集海上油气平台燃烧源和陆地掩膜数据剔除影像中的杂光，并利用中值滤波技术对影像进行去噪处理，最后利用双参数恒虚警率检测器检测影像中的船舶。结果表明，该方法可以有效降低虚警目标，能够实现大范围夜间海上船舶高精度检测。

说明书

技术领域

本发明属于夜光遥感影像处理与应用技术领域，特别涉及一种基于珞珈一号夜光遥感数据的夜间海洋船舶检测方法。

背景技术

海上空间及资源在人类发展中发挥着举足轻重的作用，船舶的海上贸易航行以及渔业的远洋捕捞活动也因此日渐频繁。尽管目前对大面积海上的船舶监测方法众多，例如利用光学遥感影像和合成孔径雷达(Synthetic Aperture Radar, SAR)影像。但依旧存在众多局限性，其中光学遥感影像无法对夜间的船舶进行检测，而SAR方法数据量大并且成本相对较高，对于大面积海域的船舶检测缺乏通用性。因此开发更加多样化的海上船舶提取方法对于海上船舶的安全监测具有重要价值。

早期，由Suomi国家极轨道伙伴关系卫星(Suomi National Polar-orbitingPartnership, SNPP)搭载的可见光红外成像辐射仪(Visible Infrared ImagingRadiometer Suite, VIIRS)获取的夜光遥感数据能够实现渔业资源的管理以及夜间海上船舶的识别。Elvidge等人（参考文献：Elvidge C, Zhizhin M, Baugh K, et al.Automatic boat identification system for VIIRS low light imaging data[J].Remote Sensing, 2015, 7(3): 3020-3036.）利用该数据进行海上渔场的分析以及船舶的检测，其通过阀值分割方法识别海上渔场范围以及船舶灯光，并剔除了海上油气平台废气燃烧光源以及其他杂光的影响，可以快速的实现大范围的渔场确定以及部分船舶的识别。但由于数据和方法的限制，导致其对于渔场有很好的识别效果，而对大范围复杂海况的单个船舶难以实现高精度检测。此外，专利“一种用于夜光遥感数据的灯光渔船自动识别方法，公开号：CN109086701A”公开了一种利用夜光遥感数据识别渔船灯光的方法，该方法能够准确识别渔船灯光的范围并实现部分船舶的提取，但是该方法使用SNPP/VIIRS夜光数据分辨率有限，并且没有剔除海上油气平台废气燃烧光源及其他杂光的影响，因此对于大范围独立船舶的检测精度较低。

发明内容

针对现有技术的不足或改进需求，本发明提出了一种基于珞珈一号夜光遥感数据的夜间海洋船舶检测方法，其目的在于克服现有技术对于大范围夜间海洋船舶提取时，数据量大及精度低的问题。

基于珞珈一号夜光遥感数据的夜间海洋船舶检测方法，包括以下步骤：

(1) 采集包含海洋区域的珞珈一号夜光遥感数据；

(2) 采集海上油气平台燃烧源的位置数据，构建掩膜对夜光遥感数据中的海上油气平台废气燃烧光源进行遮蔽；

(3) 根据陆地掩膜范围，构建缓冲区遮盖陆地灯光以及近海岸灯光；

(4) 将步骤(3)所得的影像输入中值滤波器进行去噪处理；

(5) 将步骤(4)获得的影像数据输入双参数恒虚警率检测器中进行船舶检测。

所述步骤(2)中遮蔽海上油气平台灯光的条件为：以海上燃烧源位置点为圆心，构建半径为0-2km的圆形掩膜遮盖油气平台灯光，掩膜的辐射值为0。

优选地，所述圆形掩膜半径为1km。

所述步骤(3)中遮盖陆地灯光以及近海岸灯光的条件为：根据陆地掩膜范围，构建宽度为0-4km的缓冲区遮盖陆地灯光以及近海岸灯光，掩膜的辐射值为0。

优选地，所述缓冲区宽度为2km。

进一步地，所述步骤(4)中，采用宽度为3×3的中值滤波器对影像进行去噪处理。

进一步，所述步骤(5)中，将步骤(4)获得的影像输入双参数恒虚警率检测器中进行船舶检测包含以下子步骤：

(51) 设置双参数恒虚警率检测器参数，其中目标窗口大小为1×1至3×3像素，优选2×2像素；保护窗口大小为8×8至12×12像素，优选9×9像素；背景窗口大小为16×16至24×24像素，优选18×18像素；虚警率为1e-3至1e-7，优选1e-5；所述目标窗口、保护窗口和背景窗口均为正方形；

(52) 对双参数恒虚警率检测的结果进行二值化，将大于检测阀值的目标的DN值赋值为1，小于检测阀值的目标的DN值赋值为0；

(53) 将DN值赋值为1的联通的栅格区域标记为一个船舶目标，之后输出所有船舶目标的坐标信息。

采用上述技术方案所产生的有益效果在于：

本发明基于珞珈一号夜光遥感图像特点，通过夜间海上燃烧源数据剔除油气平台，并将中值滤波与双参数恒虚警率方法结合，可以有效降低虚警目标，实现大范围夜间海上船舶高精度检测。在北部湾区域的应用实例中，体现出了巨大的实用价值。对于海洋船舶的安全监测，休渔期探测非法捕捞活动以及监控海上非法入侵等都具有重要意义。

附图说明

图1为本发明的实施流程图；

图2为实例1船舶检测结果；其中2(a)为珞珈一号原始夜光影像图，图2(b)为预处理后的珞珈一号夜光影像图，包括遮盖陆地灯光和近海岸灯光、遮盖海上油气平台废气燃烧光源以及中值滤波去噪处理过程，图2(c)为船舶检测结果，图2(d)为船舶检测结果分类图；

图3为实例2船舶检测结果；其中3(a)为珞珈一号原始夜光影像图，图3(b)为预处理后的珞珈一号夜光影像图，包括遮盖陆地灯光和近海岸灯光、遮盖海上油气平台废气燃烧光源以及中值滤波去噪处理过程，图3(c)为船舶检测结果，图3(d)为船舶检测结果分类图。

具体实施方式

以下结合附图和具体实例进一步阐述本发明。

本发明采用珞珈一号夜光遥感数据。珞珈一号是全球首颗专业夜光遥感卫星，数据来源精准可靠，该卫星可获取夜间低照度下动态范围高达14bit的高精度夜光亮度辐射信息，并且可在15天内获取全球范围的夜光遥感影像，数据分辨率达到130米。

实例1：如图1所示本发明的实现流程包含以下步骤。

(1) 采集包含海洋区域的珞珈一号夜光遥感数据；

首先采集包含海洋区域的珞珈一号夜光遥感数据，本发明以北部湾局部海域（19.40°N-21.92°N，107.92°E-110.16°E）为例进行说明，所示数据为珞珈一号卫星在2018年9月4日拍摄的影像（见图2(a)），该数据由湖北省高分辨率地球观测系统数据和应用网提供。

(2) 采集海上油气平台燃烧源的位置数据，构建掩膜对夜光遥感数据中的海上油气平台废气燃烧光源进行遮蔽；

由于在夜光影像中不仅记录了船舶灯光，同时也记录了海上油气开采时产生的废气燃烧光源。废气燃烧光源的存在会对船舶检测结果造成影响，并且无法通过灯光辐射值将其与船舶灯光区分开来。而在美国国家海洋和大气管理局（National Oceanic andAtmospheric Administration, NOAA）下的国家环境信息中心（National Centers forEnvironmental Information, NCEI）提供了全球VIIRS夜火检测数据，在2018年9月4日提供的夜间燃烧光源中，实例区域内包含的燃烧光源位置如下：

通过采集的海上油气平台燃烧源的位置数据，使用ArcGIS 10.2软件以燃烧源位置点作为圆心，构建半径为0-2km的圆形掩膜，优选1km，用以遮盖油气平台灯光；掩膜的辐射值设置为0。

(3) 根据陆地掩膜范围，构建缓冲区遮盖陆地灯光以及近海岸灯光；

根据陆地掩膜范围，使用ArcGIS 10.2软件构建宽度为0-4km的缓冲区，优选2km，用以遮盖陆地灯光以及近海岸灯光；掩膜的辐射值设置为0。

(4) 将步骤(3)所得的影像输入宽度为3×3的中值滤波器中进行去噪处理，所用软件为ENVI 5.1，处理结果如图2(b)所示。

(5) 将步骤(4)获得的影像输入双参数恒虚警率检测器中进行船舶检测包含以下子步骤：

(51) 设置双参数恒虚警率检测器参数；

双参数恒虚警率检测器是一种经典的船舶检测工具，其假设背景噪声符合高斯分布模型，采用了三个独立窗口来适应噪声的变化。根据珞珈一号夜光遥感影像的特性，其中目标窗口大小设置为1×1至3×3像素，优选为2×2像素；保护窗口大小设置为8×8至12×12像素，优选为9×9像素；背景窗口大小设置为16×16至24×24像素，优选18×18像素；所述目标窗口、保护窗口和背景窗口均为正方形；具体检测过程如下：

①在目标窗口中，灯光辐射值将与门限值进行比较，进而来判定该辐射值是否为有效目标。门限值的大小可以通过背景噪声的概率密度函数和虚警率确定，计算方法为：

其中

②根据检测门限检测目标：

其中

(52) 对双参数恒虚警率检测的结果进行二值化，将大于检测阀值的目标的DN值赋值为1，小于检测阀值的目标的DN值赋值为0，结果如图2(c)所示；

(53) 由于一个船舶目标产生的灯光区域可能分布在相连通多个栅格中，因此需要判断DN值赋值为1的栅格是否为联通区域，将联通的栅格区域标记为一个船舶目标，之后输出所有船舶目标的坐标信息。

实施例2：采集的数据为珞珈一号卫星在2019年2月15日拍摄的影像（见图3(a)）；

其处理步骤和实例1相同，检测结果如图3(c)和3(d)所示。

图中根据正确检测、遗漏检测和错误检测三种类型对检测结果的分类，从结果可以看出，本发明能够优异地实现大范围夜间海洋船舶的检测，并且具有较高的检测正确率；

检测结果的各类指标统计如下表：

实例1检测结果

实例2检测结果