CN108427824B - 一种利用船位数据提取刺网网长与角度的方法 - Google Patents
一种利用船位数据提取刺网网长与角度的方法 Download PDFInfo
- Publication number
- CN108427824B CN108427824B CN201810083391.6A CN201810083391A CN108427824B CN 108427824 B CN108427824 B CN 108427824B CN 201810083391 A CN201810083391 A CN 201810083391A CN 108427824 B CN108427824 B CN 108427824B
- Authority
- CN
- China
- Prior art keywords
- net
- ship
- point
- frequency
- lat
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Abstract
本发明涉及一种利用船位数据提取刺网网长与角度的方法,包括以下步骤:统计各个点回报的即时航速的频数,回报的经度和纬度,回报时间和即时航向,并根据经度和纬度使用半正矢公式计算各点之间的空间距离,使用回报时间和反正切公式计算各点之间的时间差和航向差,根据各点的航速、空间距离、时间差和航向差,判断刺网捕捞的收网和放网状态,根据开始收网和结束收网的点的距离,提取单次刺网网次的长度。本发明能够适用于从北斗船位数据中提取刺网渔船的捕捞状态。
Description
技术领域
本发明涉及船舶船位数据应用技术领域,特别是涉及一种利用船位数据提取刺网网长与角度的方法。
背景技术
渔船监测系统(VMS,Vessel Monitoring System)能够提供渔船时间、位置和动态信息,目前我国渔船安装船舶自动识别设备(AIS)、北斗卫星船位监控系统,CDMA公众移动通讯设备等,初步实现对海洋渔船船位的实时联络及跟踪监控。国际上研究现有渔船监测位置分析渔船捕捞状态,其方法可以概括为三种:(1)通过分析船速的变化判断捕捞状态;(2)通过分析轨迹变化判断捕捞状态;(3)通过分析船速、航向等特征数据组成向量判断捕捞状态。由于北斗卫星传送的经纬度船位数据的时间分辨率为3分钟记录一次,空间分辨率约为10米,时空精度高,实时性强。已有的方法不适用于高时间分辨率的北斗卫星的船位数据分析,目前没有专门针对北斗卫星的渔船捕捞状态判断方法。需要找到针对特定网具的计算方法,用来判断基于高时间分辨率的北斗船位数据的渔船捕捞状态。
发明内容
本发明所要解决的技术问题是提供一种利用船位数据提取刺网网长与角度的方法,能够适用于从北斗船位数据中提取刺网渔船的捕捞状态。
本发明解决其技术问题所采用的技术方案是:提供一种利用船位数据提取刺网网长与角度的方法,包括以下步骤:统计各个点回报的即时航速的频数,回报的经度和纬度,回报时间和即时航向,并根据经度和纬度使用半正矢公式计算各点之间的空间距离,使用回报时间和反正切公式计算各点之间的时间差和航向差,根据各点的航速、空间距离、时间差和航向差,判断刺网捕捞的收网和放网状态,根据开始收网和结束收网的点的距离,提取单次刺网网次的长度。
根据刺网渔船航速统计船位点在各航速出现的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的航速阈值的最小值与最大值;利用航速阈值筛选出可能处于收网状态的船位点,然后计算相邻船位点的空间距离,统计船位点在各空间距离间隔长度出现的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的空间距离间隔阈值的最小值与最大值;计算相邻船位点的时间间隔,统计在各时间间隔出现的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的时间间隔阈值的最小值与最大值;计算相邻船位点的角度变化,统计在各角度变化值的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的角度变化值的最小值与最大值;根据刺网空间间隔、时间间隔、角度变化间隔的阈值提取刺网起始点和终止点,然后利用获取起始点和终止点计算刺网长度和角度。
在统计各个点回报的即时航速的频数时,当出现连续的航速为0m/s时,将这些数据舍去。
所述各点之间的航向差通过θ=atan2(sin(Blon-Alon)×cosBlat,cosAlat×sinBlat-sinAlat×cosBlat×cos(Blon-Alon))计算得到,其中,Alon和Alat表示船位点A的经度和纬度,Blon和Blat表示船位点B的经度和纬度,
有益效果
由于采用了上述的技术方案,本发明与现有技术相比,具有以下的优点和积极效果:本发明利用船位数据提取刺网网长与角度的方法,精度较高,实时性强。
附图说明
图1是刺网出海作业轨迹点示意图;
图2是刺网船位点在各航速出现的频数统计图;
图3是航速阈值筛选可能的刺网收网点(粗的黑色点)示意图;
图4是刺网各网次间存在距离间隔示意图;
图5是相邻船位点各距离间隔出现的频数统计图;
图6是刺网相邻网次存在距离较小的现象示意图;
图7是相邻点各时间间隔出现的频数统计图;
图8是符合距离间隔和时间间隔但不是网次的点示意图;
图9是相邻点的各角度变化值频数统计图;
图10是刺网各网次提取示意图;
图11是刺网各网次放网长度分布统计图;
图12是刺网各网次放网角度分布统计图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明的实施方式涉及一种利用船位数据提取刺网网长与角度的方法,包括以下步骤:统计各个点回报的即时航速的频数,回报的经度和纬度,回报时间和即时航向,并根据经度和纬度使用半正矢公式计算各点之间的空间距离,使用回报时间和反正切公式计算各点之间的时间差和航向差,根据各点的航速、空间距离、时间差和航向差,判断刺网捕捞的收网和放网状态,根据开始收网和结束收网的点的距离,提取单次刺网网次的长度。当得到单次刺网网次的长度后可以通过求和计算单艘渔船刺网作业的总长度和所有刺网渔船作业的总网长。
下面以具体的实施例来进一步说明本发明,一艘刺网渔船在2016年全年船位点的空间分布如图1,按照时间排序组成空间轨迹点。刺网作业过程包括开始放网、结束放网、等待渔获进入刺网、开始收网、结束收网。其中放网状态航速较快,等待时渔船处于抛锚或漂流,因此航速基本为0,收网渔船速度较慢。一个完整的刺网作业过程放网速度快,放网状态点少,收网速度慢,放网状态点多,因此刺网网次的起止点,用收网起止点来表示。
1.航速频数统计
根据刺网渔船航速统计船位点在各航速出现的频数,在渔船进行收网作业时,处于一个相对接近的较低的航速,因此,找到航速频数的峰值,获取峰值两侧的谷值可以作为刺网收网的航速阈值的最小值与最大值;可以认为,只有航速处于航速阈值内的船位点,才是正在进行收网作业的船位点。图2是渔船的航速频数曲线,第一个峰值为渔船收网状态,第二峰值为渔船放网状态,第一个峰值两侧获取收网状态的阈值,即最小值为0m/s,最大值为1.5m/s,通过航速阈值筛选可能的刺网收网点(图3中粗的黑色点);
2.空间距离频数统计
根据北斗导航系统回报的数据,可以获取每个船位点的实时经纬度,根据前后两个点之间的经纬度,能计算出这两个点的空间距离。
由于地球是一个两级部位略扁的不规则球体,在进行大地坐标系的计算时,一般可将其作为一个球体计算。因此,地球上已知经纬度的两点之间,可以使用球面距离公式计算其距离。
假设两个船位点A,B,其经纬度分别为Alon,Alat,Blon,Blat。
那么A、B两点之间的距离D可以用半正矢公式(Haversine formula)求得。
其公式为:
其中R为地球平均半径6378145m
利用航速阈值筛选出可能处于收网状态的点之后,计算相邻船位点的空间距离,统计船位点在各距离间隔长度出现的频数。根据对于数据的统计和刺网作业的流程,刺网作业收网时的船位点存在空间距离相对较近的特点,而每个网次之间存在一定的空间间隔。因此,找到频数峰值,获取峰值两侧的谷值作为刺网收网的空间距离间隔阈值的最小值与最大值;只有和其他点处在空间距离间隔阈值范围内的船位点才是处于收网状态的船位点。而刺网的各网次间存在较大的距离(图4),在实例中(图5)各刺网的空间距离间隔阈值的最小值为178m,最大值为193m。
3.时间间隔频数统计
北斗船位数据回报的数据集中包括每个船位点的时间,因此可以计算出相邻船位点的时间间隔,统计在各时间间隔出现的频数,和空间间隔相似,根据刺网作业的特点,收网作业时,存在船位点以较小的时间间隔密集存在。找到频数峰值,获取峰值两侧的谷值,即时间间隔阈值的最小值与最大值作为刺网各网次的分界值,刺网的各网次有时间距离较近,利用空间距离有时难以区分(图6),但各网次间的都有一定的时间差值,在实例中(图7)刺网各网次的时间间隔阈值的最小值为30分钟,最大值为46分钟。
4.刺网角度频数统计
在计算航速阈值、空间间隔和时间间隔后,已经可以排除绝大多数的非收网状态船位点,但依然存在少数符合距离间隔和时间间隔但不是网次的船位点(图8),为排除这类点,根据刺网作业流程,可以根据船只的航向变化进行判断。
虽然北斗船位数据回报的数据中包括各船位点的即时航向数据,但是由于计算方式的原因,在航速为0的时候所统计的航向值一般为0,无法反映此时的渔船方向。而由于刺网作业的操作规程的原因,在收网过程中,有多个船位点航速可能为0,如果将其全部删除,可能影响对于单个网次长度的统计。因此,在此采用大圆航向(Great-circle navigation)的方法计算船只的数据,即计算前一船位点和至后一船位点二点连线的距离作为该点的实时航向,而非使用船位数据系统回报的数据。
假设两个船位点A,B,其经纬度分别为Alon,Alat,Blon,Blat。
那么从A至B的航向方位角θ可以如下公式求得。
θ=atan2(sin(Blon-Alon)×cosBlat,cosAlat×sinBlat-sinAlat×cosBlat×cos(Blon-Alon))
其中,atan2是一个专门的函数,用以计算两个值之间的反正切,将直角坐标系中的坐标转换为极坐标。表示坐标原点为起点,指向(x,y)的射线在坐标平面上与x轴正方向之间的角的角度。其结果可以用以下等式说明:
对于两个坐标值x,y
该函数的值域为(-π,π],可以通过对负数结果加2π的方法,将函数的结果映射到[0,2π)范围内。
其余的算式则是将处于球面坐标系中的Alon,Alat,Blon,Blat转换为以A为原点的平面直角坐标系上的坐标。
得到每个船位点的航向角后,计算相邻船位点的角度变化,统计在各角度变化值的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的角度变化值的最小值与最大值;根据操作规范,在收网时应该大致保证航向不变。频数峰值也正如同规程所述,峰值应当出现在变化值为0的点处附近。在实例中(图9)角度变化值阈值的最小值为-8°,最大值为10°。
5.刺网网次长度提取
根据刺网航速阈值先提取可能的刺网收网状态点,再用空间间隔、时间间隔的阈值区分各网次,然后用角度变化间隔去除非刺网的网次,得出每次作业中开始收网和结束收网的船位点(图10)。这两个点之间的空间距离作为渔船刺网作业中的一个网次的长度(图11),这两个点之间的空间航向角可以认为是该次刺网作业中收网的方向(图12)。
根据收网的船位点,可以逆推出放网的过程。根据刺网的操作流程,放网的开始点应该邻近收网的结束点,放网的结束点应该邻近收网的开始点。因此,在此可以取大于收网航速阈值(本例中为1.5m/s)的点中最靠近收网起始/结束点的船位点作为放网结束、开始的船位点。
单次刺网作业的总时间,可以简单使用放网开始点的时间减去收网结束点的时间进行计算。
表1刺网网次提取表
6.渔船刺网网次长度提取
在计算出单艘渔船每次作业的网次长度后,将该航次中所有刺网网次的长度求和并提取,得到单艘渔船刺网作业的总作业时间和总网次长度。
Claims (4)
1.一种利用船位数据提取刺网网长与角度的方法,其特征在于,包括以下步骤:统计各个点回报的即时航速的频数,回报的经度和纬度,回报时间和即时航向,并根据经度和纬度使用半正矢公式计算各点之间的空间距离,使用回报时间和反正切公式计算各点之间的时间差和航向差,根据各点的航速、空间距离、时间差和航向差,判断刺网捕捞的收网和放网状态,根据开始收网和结束收网的点的距离,提取单次刺网网次的长度;具体为:根据刺网渔船航速统计船位点在各航速出现的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的航速阈值的最小值与最大值;利用航速阈值筛选出可能处于收网状态的船位点,然后计算相邻船位点的空间距离,统计船位点在各空间距离间隔长度出现的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的空间距离间隔阈值的最小值与最大值;计算相邻船位点的时间间隔,统计在各时间间隔出现的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的时间间隔阈值的最小值与最大值;计算相邻船位点的角度变化,统计在各角度变化值的频数,找到频数峰值,获取峰值两侧的谷值作为刺网收网的角度变化值的最小值与最大值;根据刺网空间间隔、时间间隔、角度变化间隔的阈值提取刺网起始点和终止点,然后利用获取起始点和终止点计算刺网长度和角度。
2.根据权利要求1所述的利用船位数据提取刺网网长与角度的方法,其特征在于,在统计各个点回报的即时航速的频数时,当出现连续的航速为0m/s时,将这些数据舍去。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810083391.6A CN108427824B (zh) | 2018-01-29 | 2018-01-29 | 一种利用船位数据提取刺网网长与角度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810083391.6A CN108427824B (zh) | 2018-01-29 | 2018-01-29 | 一种利用船位数据提取刺网网长与角度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108427824A CN108427824A (zh) | 2018-08-21 |
CN108427824B true CN108427824B (zh) | 2021-08-31 |
Family
ID=63156086
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810083391.6A Active CN108427824B (zh) | 2018-01-29 | 2018-01-29 | 一种利用船位数据提取刺网网长与角度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108427824B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109389186B (zh) * | 2018-11-23 | 2022-03-18 | 中国水产科学研究院东海水产研究所 | 基于dbscan的vms定置刺网渔船网次提取方法 |
CN109918608B (zh) * | 2019-01-10 | 2023-07-14 | 中国水产科学研究院东海水产研究所 | 一种提取渔船作业航次中离岸和离港距离的方法 |
CN110135559A (zh) * | 2019-04-28 | 2019-08-16 | 中国水产科学研究院东海水产研究所 | 一种基于船位数据获取张网网位和捕捞努力量的方法 |
CN110047331A (zh) * | 2019-05-22 | 2019-07-23 | 镇江市地方海事局 | 一种基于ais的船舶与目标物之间距离测算方法 |
CN112434465B (zh) * | 2020-11-19 | 2021-11-09 | 江苏省海洋水产研究所 | 基于船位数据提取毛虾网有效布网网长的方法 |
CN112699531B (zh) * | 2020-12-08 | 2023-06-09 | 浩亚信息科技有限公司 | 一种通用航空低空飞行目视参考点建立方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104115795A (zh) * | 2014-04-10 | 2014-10-29 | 湛江胜浪海洋捕捞研究所 | 一种大洋性流刺网渔船 |
CN104463433A (zh) * | 2014-11-19 | 2015-03-25 | 中国水产科学研究院东海水产研究所 | 一种基于渔船船位数据提取刺网捕捞状态的方法 |
CN204659993U (zh) * | 2015-05-22 | 2015-09-23 | 中国水产科学研究院东海水产研究所 | 一种船用打捞弯钩 |
CN106021675A (zh) * | 2016-05-13 | 2016-10-12 | 国家海洋局第三海洋研究所 | 一种基于多源船舶自动识别系统数据的数据融合方法 |
CN106372094A (zh) * | 2015-07-24 | 2017-02-01 | 通用电气公司 | 一种路径数据库的构建方法及基于该数据库的车辆运行控制系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100242851A1 (en) * | 2009-03-26 | 2010-09-30 | Lee Carstensen | Systems and Methods for Aquatic Electrified Barriers |
CN103927858B (zh) * | 2014-04-09 | 2016-09-14 | 广东海洋大学 | 基于北斗定位通信单元的船载手持渔获信息采集装置 |
-
2018
- 2018-01-29 CN CN201810083391.6A patent/CN108427824B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104115795A (zh) * | 2014-04-10 | 2014-10-29 | 湛江胜浪海洋捕捞研究所 | 一种大洋性流刺网渔船 |
CN104463433A (zh) * | 2014-11-19 | 2015-03-25 | 中国水产科学研究院东海水产研究所 | 一种基于渔船船位数据提取刺网捕捞状态的方法 |
CN204659993U (zh) * | 2015-05-22 | 2015-09-23 | 中国水产科学研究院东海水产研究所 | 一种船用打捞弯钩 |
CN106372094A (zh) * | 2015-07-24 | 2017-02-01 | 通用电气公司 | 一种路径数据库的构建方法及基于该数据库的车辆运行控制系统 |
CN106021675A (zh) * | 2016-05-13 | 2016-10-12 | 国家海洋局第三海洋研究所 | 一种基于多源船舶自动识别系统数据的数据融合方法 |
Non-Patent Citations (3)
Title |
---|
Interaction of Harsh Weather Operation and Collision Avoidance in Autonomous Navigation;H.‐C.Burmeister等;《the International Journal on Marine Navigation and Safety of Sea Transportation》;20150331;第9卷(第1期);31-40 * |
一种大椭圆航法导航参数计算方法研究;房慎冲等;《测绘科学》;20150520;第40卷(第5期);21-25 * |
基于神经网络和VMS的渔船捕捞类型辨别;郑巧玲等;《南方水产科学》;20160405;第12卷(第2期);81-87 * |
Also Published As
Publication number | Publication date |
---|---|
CN108427824A (zh) | 2018-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108427824B (zh) | 一种利用船位数据提取刺网网长与角度的方法 | |
Podesta et al. | Exploring the association between swordfish catch rates and thermal fronts on US longline grounds in the western North Atlantic | |
Hu et al. | Monitoring coastline variations in the Pearl River Estuary from 1978 to 2018 by integrating Canny edge detection and Otsu methods using long time series Landsat dataset | |
CN109389186B (zh) | 基于dbscan的vms定置刺网渔船网次提取方法 | |
CN102243298B (zh) | 一种基于dem的机载气象雷达地杂波剔除方法 | |
CN103292792B (zh) | 一种适用海底探测与假地形处理的实测svp重构方法 | |
CN107977635B (zh) | 一种格状水系识别方法 | |
CN110147716B (zh) | 一种基于频域与空域相结合的sar图像内波区域检测方法 | |
CN103971370A (zh) | 一种针对遥感大图像的海洋溢油智能检测方法 | |
Zimmermann | Calculation of untrawlable areas within the boundaries of a bottom trawl survey | |
CN109188430A (zh) | 一种基于地面监视雷达系统的目标提取方法 | |
CN113093119A (zh) | 一种时频恒虚警高频雷达目标检测方法及系统 | |
CN116645789A (zh) | 一种基于无人机测绘的土地自然灾害预警管控系统 | |
CN112784180A (zh) | 一种金枪鱼围网渔船捕捞强度空间信息提取方法 | |
Wang et al. | Characterizing the capability of mesoscale eddies to carry drifters in the northwest Pacific | |
CN104297753B (zh) | 一种基于自适应缩减算子的导航雷达图像反演海面风向方法 | |
Shao et al. | Automated searching of ground points from airborne lidar data using a climbing and sliding method | |
CN114063616A (zh) | 基于三维激光扫描探测规划林区路径的方法及装置 | |
CN110988888B (zh) | 海底信息的获取方法和装置 | |
CN113269078B (zh) | 一种基于无人机的江河鱼类资源量统计方法及系统 | |
CN111080536A (zh) | 一种机载激光雷达点云的自适应滤波方法 | |
CN104615643A (zh) | 一种基于渔船船位数据提取耙刺捕捞状态的方法 | |
CN112541487B (zh) | 基于亚米级高分卫星遥感影像的牧区羊群的自动监测方法 | |
CN109918968A (zh) | 一种舰船目标检测方法 | |
CN104063614A (zh) | 一种基于odp的海底沙波特征自动识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |