CN117971959A - 内河船舶大气污染物排放清单编制及可视化方法 - Google Patents
内河船舶大气污染物排放清单编制及可视化方法 Download PDFInfo
- Publication number
- CN117971959A CN117971959A CN202410356993.XA CN202410356993A CN117971959A CN 117971959 A CN117971959 A CN 117971959A CN 202410356993 A CN202410356993 A CN 202410356993A CN 117971959 A CN117971959 A CN 117971959A
- Authority
- CN
- China
- Prior art keywords
- ship
- point
- track
- folding
- inland
- 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.)
- Granted
Links
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 53
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000010586 diagram Methods 0.000 claims abstract description 11
- 238000004326 stimulated echo acquisition mode for imaging Methods 0.000 claims abstract description 7
- 230000000007 visual effect Effects 0.000 claims abstract description 5
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 description 7
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 6
- 238000012937 correction Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 235000010269 sulphur dioxide Nutrition 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000004291 sulphur dioxide Substances 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Landscapes
- Navigation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种内河船舶大气污染物排放清单编制及可视化方法，包括S1：结合AIS数据和STEAM模型，计算出船舶大气污染物的排放量；S2：船舶轨迹点空间位置转换，计算船舶轨迹点对应的里程坐标；S3：建立可视化模型，包括S3.1：根据里程坐标判断船舶的航向；S3.2：沿长度方向将内河航道对应的折线分为若干线段，判断船舶轨各个迹点对应里程坐标具体属于哪一线段；S3.3：按照时间间隔，根据步骤S1计算得到的船舶轨迹段污染物排放量，计算步骤S3.2中各个线段上对应的污染物排放量，绘制热力图展现出污染物排放量的时空分布情况。该方法将船舶的空间轨迹用一维的里程坐标表示，结合时间维度，能够在一张图上展现内河船舶污染物排放的时空分布特征。
Description
技术领域
本发明属于船舶污染物排放清单编制技术领域，具体涉及一种内河船舶大气污染物排放清单编制及可视化方法。
背景技术
船舶自动识别系统（Automatic Identification System，AIS），由岸基设施和船载设备共同组成，通过对全球船舶进行实时动态数据采集，其中广播的每一帧AIS报文通常是由动态信息（例如船载台向岸基台自动播发的本船的经纬度坐标、航速、航向等）及静态信息（船名、船的大小、目的港等）构成的。
公布号CN105590024A的中国专利申请，公开了一种基于活动的港区大气污染物排放清单编制方法，该方法中结合AIS和船舶档案数据库信息构建了本土化STEAM模型，计算得到南京龙潭港靠港船舶排放清单。
通过采用STEAM模型，可以计算相邻轨迹点之间的大气污染物排放量，这种方法计算出来的排放清单是一种复杂的空间数据，具有数据量大、维度高、层次多的特点，但是难以直观地进行可视化。
目前的排放清单可视化技术大多通过网格划分方法将大量矢量数据点栅格化，然后使用克里金插值或核密度图方法绘制船舶排放的空间分布图。然而，这种方法存在一些明显的问题。首先，由于数据点的不均匀性，这种方法得出的结果往往不够准确。其次，这种方法不能很好地分析船舶排放在时间上的变化情况，需要设定统一的标准，分别绘制多个时间点的分布图进行对比，因此效率较低。此外，这种方法一般只适用于海洋上，内河船舶分布受河道宽度限制，绘制出来的效果不佳。
发明内容
本发明的内河船舶大气污染物排放清单编制及可视化方法，采用AIS数据和STEAM模型，准确地计算出内河船舶的排放清单，并采用时空可视化模型将排放清单可视化。
为了实现上述目的，本发明的内河船舶大气污染物排放清单编制及可视化方法，包括以下步骤：
S1：结合AIS数据和STEAM模型，计算出船舶大气污染物的排放量；
S2：船舶轨迹点空间位置转换，并计算船舶轨迹点对应的里程坐标，具体包括如下步骤：
S2.1：将内河航道的遥感影像和船舶AIS数据中船舶轨迹点转换为统一的平面坐标；
S2.2：矢量化内河航道，记录内河航道的折点坐标；
S2.3：根据折点位置将内河航道划分为不同的区域，计算船舶轨迹点所属区域；
S2.4：根据所属区域，计算船舶轨迹点对应的里程坐标；
S3：建立可视化模型，包括如下步骤：
S3.1：根据里程坐标判断船舶的航向；
S3.2：沿长度方向将内河航道对应的折线分为若干个等间距的线段，并沿航向依次编号，并判断各个船舶轨迹点对应里程坐标具体属于哪一个线段；
S3.3：按照时间间隔，根据步骤S1计算得到的船舶轨迹段排放量，计算步骤S3.2中各个线段上所有船舶对应的污染物排放量大小，绘制热力图展现出污染物排放量的时空分布情况。
进一步地，步骤S2.1中，将内河航道的遥感影像导入到软件ArcGIS中，使用投影关系和变换工具，将遥感影像的坐标系转换为高斯平面直角坐标系；利用高斯正算，将船舶AIS数据中船舶轨迹点坐标，转化为高斯平面直角坐标系中的坐标。
进一步地，步骤S2.3中，将内河航道对应的折线上各个折点从左至右依次进行编号，记为1、2、3、……，M；M表示折线折点的数量；根据各个折点与其两侧相邻两个折点的连线所形成夹角的角平分线，将内河航道划分为不同的区域，角平分线作为区域的边界；对所有区域进行编号，从左至右依次记为1、2、3、……，M+1。
第i个折点处的角平分线方程为：
y-k i x-b i =0;
；
b i= y i -k i x i ;；
其中，（x i ，y i ）为第i个折点的平面直角坐标；（x i-1，y i-1）是第i-1个折点的平面直角坐标；（x i+1，y i+1）是第i+1个折点的平面直角坐标；k i 为第i个折点处的角平分线的斜率，即第i-1个折点、第i个折点和第i+1个折点这三点的连线所形成夹角的角平分线斜率，b i 为第i个折点处的角平分线截距。
进一步地，步骤S2.3中，判断各个船舶轨迹点具体属于哪个区域的方法如下：
S2.3.1: 设船舶轨迹上某一点的平面直角坐标为(x 船，y 船)；
S2.3.1.1:设第1折点处的角平分线方程斜率为k 1，第1折点处的角平分线方程为：y 1-k 1 x-b 1 =0；
如果满足条件一：k 1<0，且同时，y 船≤k 1 x 船 + b 1；则，船舶轨迹点(x 船，y 船)归属于第1区域；
或者，如果满足条件二：k 1>0，且同时，y 船＞k 1 x 船 +b 1；则，船舶轨迹点(x 船，y 船)归属于第1区域；
S2.3.1.2：如果既不满足条件一也不满足条件二，则根据第2折点处的角平分线方程来判断；即，
如果满足条件三：k 2<0，且同时，y 船≤k 2 x 船 +b 2；则，船舶轨迹点(x 船，y 船)归属于第2区域；
或者，如果满足条件四：k 2>0，且同时，y 船＞k 2 x 船 +b 2；则，船舶轨迹点(x 船，y 船)归属于第2区域；
……
以此，遍历其余所有折点对应的角平分线方程；
如果不满足前面的所有条件，则船舶轨迹上的该点 (x 船，y 船)归属于第M+1区域；
S2.3.2:重复步骤S2.3.1，对整个船舶轨迹上的所有点进行区域划分。
进一步地，步骤S2.4中，根据船舶轨迹点所属区域计算该点对应的里程坐标，具体包括如下步骤：
S2.4.1：将内河航道对应折点的平面直角坐标转化为一维的里程坐标，设内河航道对应折线起点的里程坐标为0，并依次计算所有折点对应的里程坐标；
其中，第1个折点的里程坐标S 1：
；
其余折点的里程坐标按照如下方式计算，第i个折点的里程坐标S i ：
；
S2.4.2：当船舶轨迹点(x 船，y 船)投影位于第i个折点和第i+1个折点之间，即船舶轨迹点(x 船，y 船)落入上述计算出的对应区域内时，该点(x 船，y 船)的里程坐标S 船的计算公式如下：
；
进一步地，步骤S3.2具体包括如下步骤：
S3.2.1：沿长度方向将内河航道对应的折线分为N个线段，并计算出各个线段终点处的里程坐标S n ，n∈[1，N]；
；
其中，S M 表示内河航道折线最末端的折点的里程坐标，l为划分线段的间隔值，为向上取整运算；
S3.2.2：根据各个船舶轨迹点投影在折线上的位置，判断其具体属于哪个线段。具体为：根据对应船舶轨迹点的里程坐标与最近的线段终点里程坐标进行比较，通过比较大小很容易得知船舶轨迹点该线段终点的左边或右边。
进一步地，步骤S3.3，计算船舶轨迹在各个线段上对应的污染物排放量，具体为:
设定船舶轨迹的起点位置对应的里程坐标S 船起，船舶轨迹的终点位置对应的里程坐标S 船终，该船舶整个轨迹点共落入K个线段中；
设定船舶轨迹的起点位于第n线段中时，对应第n线段上的污染物排放量E n 为：
；
则船舶轨迹的终点位于第n+K线段中，对应第n+K线段上的污染物排放量E n+K 为：
；
其余中间K-2个线段上的污染物排放量均为E n+1 ；
；
其中，E为在船舶整个轨迹段的污染物排放量；
重复该步骤，计算航道上所有船舶在对应线段上的污染物排放量。
与现有技术相比，本发明具有如下有益效果：
1、本发明中将船舶的空间轨迹转换为一维的里程坐标，结合时间维度，能够在一张图上展现内河船舶排放的时空分布特征，并且可以区分不同航向的排放。这使得排放清单的时间和空间分布情况更加直观、清晰。
2、本发明中先计算出轨迹点间的排放量，再进行轨迹点空间上的转换，并且将一个轨迹段的排放量按比例分配到不同线段上，保障了时空分布图的准确性。
3、本发明所述方法能够提供有关内河航运大气污染的全面数据，以帮助改善污染物监管和治理政策，并改善沿线城市空气环境质量，为内河航运的绿色可持续发展提供支持和方向。
附图说明
图1是本发明所述内河船舶大气污染物排放清单编制及可视化方法的流程图；
图2是南京长江段主航道的折线图；
图3是南京长江段主航道对应折线划分区域后的示意图；
图4是将内河航道线段划分后，将相邻轨迹点间排放量分配到各个航线段上的示意图；
图5是从上游到下游的船舶二氧化硫排放量的示意图；
图6是从下游到上游的船舶二氧化硫排放量的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1，内河船舶大气污染物排放清单编制及可视化方法，包括以下步骤：
S1：计算船舶大气污染物的排放量，具体包括如下步骤：
S1.1：采用AIS数据和STEAM模型结合，计算大气污染物的排放量，公式为
；
其中，E i 为第g种大气污染物的排放量；j为工作引擎，包括主机和辅机；k为航行工况，包括停泊、操纵和巡航； Load j 为船舶主机和辅机的负荷功率；MCR j 为船舶主机和辅机的额定功率；LF j 为船舶主机和辅机负载系数；EF j 为污染物在不同工作引擎工作条件下的排放因子；CF为主机低负荷修正系数；hr jk 为不同工作引擎在不同航行工况下的运行时间。
S1.2：将上式中指标本地化处理：
以长江南京段为例，船舶特点分为：客船、工作船和其他船三大类；船舶的工作引擎一般包括主机和辅机；船舶的航行工况为停泊、操纵和巡航。船舶航行状态的划分标准表1所示：
表1：船舶航行状态的划分标准
表1中， LF 1为主机负载系数。
式中，AS为船舶的实际运行速度，用轨迹段的平均速度代替；MS为船舶的最大设计航速。
参考《川江及三峡库区运输船舶标准船型主尺度系列》(2010年修订版），得到符合本地的船舶最大设计航速，如表2所示：
表2：船舶最大设计航速及辅机额定功率与主机额定功率的比值
选用生态环境部发布的《移动源（船舶）排放系数手册》推荐数据，代表我国船舶辅助发动机的实际情况，如表3所示：
表3：船舶不同行驶状态下辅机负载系数
船舶主机在航行中的排放因子不是固定的。当主机负荷功率降低时，气缸内温下降，燃油燃烧质量变差，污染物排放浓度逐渐升高。低负荷状态被定义为主机负荷功率低于20%，并引入主机低负荷修正系数CF，具体如表4所示。
不同污染物对应的船舶主机低负荷修正系数：
;
式中，CF为主机低负荷修正系数，X为主机负荷系数。
表4：内河船舶发动机排放因子
本实施例中，长江南京段船舶大气污染物排放量测算实例的部分数据如表5所示。
表5：长江南京段船舶大气污染物排放量测算的部分数据
S2：轨迹点空间位置转换。
根据步骤S1计算出同一船舶相邻轨迹点间的排放量，但是这样的排放量只是在时间上统计区域整体的排放量情况，不能分析污染物排放量的空间分布特征。
为此，传统的方法是将研究区域网格划分，然后将网格内所有轨迹段的排放汇总到网格中心点，再使用空间插值的方法进行可视化。
该方式因为存在背景技术中所提到的问题，所以在此不使用传统的区域网格划分方法。本申请提出另一种方法。
由于内河船舶分布受河道限制，只需研究船舶排放沿航道方向上的空间分布。本申请中步骤S2具体包括以下步骤：
S2.1：将内河航道的遥感影像和船舶AIS数据中轨迹点转换为统一的平面坐标。
本实施例中，内河航道的遥感影像是从地理空间数据云（地理空间数据云(gscloud.cn)）下载的Landsat 8 南京地区的遥感影像。将遥感影像导入软件ArcGIS中，利用投影和变换工具，将遥感影像的坐标系转换为高斯平面直角坐标系。
AIS轨迹点的坐标采用的是大地坐标系，将其转为高斯平面直角坐标系，该转换称为高斯投影的正算。
S2.2：矢量化内河航道。
首先在遥感影像中裁剪出包含内河航道的研究区域，然后根据内河航道的弯曲程度，人为选取关键节点，绘制一条折线表示内河航道，记录各个折点的坐标。
S2.3：根据折点位置将内河航道划分为不同的区域，计算船舶轨迹点所属区域，具体包括以下步骤：
S2.3.1：根据折点位置将内河航道划分为不同的区域。
将内河航道对应的折线上各个折点从左至右依次进行编号，记为1、2、3、……，M；M表示折线折点的数量；根据各个折点与其两侧相邻两个折点的连线所形成夹角的角平分线，将内河航道划分为不同的区域，角平分线作为区域的边界；对所有区域进行编号，从左至右依次记为1、2、3、……，M+1。
本实例中，南京长江段主航道提取后的矢量折线如图2所示，共有5个折点，并编号1、2、3、4、5。南京长江段主航道被分为6个区域，并编号1、2、3、4、5、6，如图3所示。各折点的高斯平面直角坐标如表6所示。
第i个折点处的角平分线方程为：
y-k i x-b i =0;
；
b i= y i -k i x i ;
其中，（x i ，y i ）为第i个折点的平面直角坐标；（x i-1，y i-1）是第i-1个折点的平面直角坐标；（x i+1，y i+1）是第i+1个折点的平面直角坐标；k i 为第i个折点处的角平分线的斜率，即第i-1个折点、第i个折点和第i+1个折点这三点的连线所形成夹角的角平分线斜率，b i 为第i个折点处的角平分线截距。
S2.3.2：判断各个船舶轨迹点具体属于哪个区域。
S2.3.2.1: 设船舶轨迹上某一点的平面直角坐标为(x 船，y 船)；
S2.3.2.1.1:设第1折点处的角平分线方程斜率为k 1，第1折点处的角平分线方程为：y-k 1 x-b 1 =0；
如果满足条件一：k 1<0，且同时，y 船≤k 1 x 船 + b 1；则，船舶轨迹点(x 船，y 船)归属于第1区域；
或者，如果满足条件二：k 1>0，且同时，y 船＞k 1 x 船 +b 1；则，船舶轨迹点(x 船，y 船)归属于第1区域；
S2.3.2.1.2：如果既不满足条件一也不满足条件二，则根据第2折点处的角平分线方程来判断；即，
如果满足条件三：k 2<0，且同时，y 船≤k 2 x 船 +b 2；则，船舶轨迹点(x 船，y 船)归属于第2区域；
或者，如果满足条件四：k 2>0，且同时，y 船＞k 2 x 船 +b 2；则，船舶轨迹点(x 船，y 船)归属于第2区域；
……
以此，遍历其余所有折点对应的角平分线方程；
如果不满足前面的所有条件，则船舶轨迹上的该点 (x 船，y 船)归属于第M+1区域。
S2.3.2.2:重复步骤S2.3.2.1，对整个船舶轨迹上所有点进行区域划分。
S2.4：根据所属区域，计算船舶轨迹点对应的里程坐标；
S2.4.1：将内河航道对应折点的平面直角坐标转化为一维的里程坐标，设内河航道对应折线起点的里程坐标为0，并依次计算所有折点对应的里程坐标；
其中，第1个折点的里程坐标S 1：
；
其余折点的里程坐标按照如下方式计算，第i个折点的里程坐标S i ：
；
本实施例中，南京长江段主航道对应折线各个折点的高斯平面坐标，以及里程坐标如表6所示。
表6：折点的坐标
S2.4.2：当船舶轨迹点(x 船，y 船)投影位于第i个折点和第i+1个折点之间，即船舶轨迹点(x 船，y 船)落入上述计算出的对应区域内时，该点(x 船，y 船)的里程坐标S 船的计算公式如下：
；
本实施例中，南京长江段主航道上某只船舶位于折点4与折点5之间的部分轨迹点的数据，如表7所示，包括船舶轨迹点的AIS数据、高斯平面坐标系和里程坐标。由于整个轨迹点的数据量太大，限于篇幅，此处仅提供一部分。
表7：AIS轨迹数据坐标转换
S3：建立可视化模型，具体包括如下步骤：
S3.1：根据里程坐标判断船舶的航向；
内河中船舶只有两个航向，上游至下游或下游至上游。本实施例中船舶轨迹点都换算成了里程，设里程越大，距离上游越远。通过对比轨迹段起点和终点的里程，可以判断出船舶的航向。
S3.2：沿长度方向将内河航道对应的折线分为若干个线段，并沿航向依次编号，并判断船舶轨迹点对应里程坐标具体属于哪一个线段；
S3.2.1：沿长度方向将内河航道对应的折线分为N个线段，并计算出各个线段终点处的里程坐标S n ，n∈[1，N];
；
其中，S M 表示内河航道折线最末端的折点的里程坐标，l为划分线段的间隔值，为向上取整运算；
本实施例中，以l=3(km)划分内河航道。长江南京段主航道共96.485794km，故将其对应的折线划分为33个线段，其中前32个线段为等间隔，第33个线段为其末段的0.485794km。
S3.2.2：根据船舶轨迹点投影在折线上的位置，判断其具体属于哪个线段。
S3.3：按照时间间隔，根据步骤S1计算得到的船舶轨迹段排放量，计算步骤S3.2中各个线段上所有船只对应的排放量大小，绘制热力图展现出排放量的时空分布情况。具体包括:
S3.3.1：设定船舶轨迹的起点位置对应的里程坐标S 船起，船舶轨迹的终点位置对应的里程坐标S 船终，以及该船舶整个轨迹点共落入K个线段中。如图4所示，图中线型较粗的线条表示其中某一只船舶轨迹。
设定船舶轨迹的起点位于第n线段中时，对应第n线段上的污染物排放量E n 为：
；
则船舶轨迹的终点位于第n+K-1线段中，对应第n+K-1线段上的污染物排放量E n+K-1为：
；
其余中间K-2个线段上的排放量均为E n+1 ；
；
其中，E为在船舶整个轨迹段的污染物排放量。
重复该步骤，计算航道上所有船舶在对应线段上的污染物排放量。
S3.3.2：按照一定的时间间隔，汇总S3.2中各个线段上所有船舶排放量大小，绘制热力图即可展现出排放量的时空分布情况。本实施例中，以二氧化硫为例，分别展示了从上游到下游，以及从下游到上游的排放量分布图，如图5和图6所示，图中右侧条形图表示排放量对应的颜色，单位为105g，在一张图上展现内河船舶排放的时空分布特征，并且可以区分不同航向的排放，非常清晰、直观。
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。
Claims (7)
1.内河船舶大气污染物排放清单编制及可视化方法，其特征在于，包括以下步骤，
S1：结合AIS数据和STEAM模型，计算出船舶大气污染物的排放量；
S2：船舶轨迹点空间位置转换，并计算船舶轨迹点对应的里程坐标，具体包括如下步骤：
S2.1：将内河航道的遥感影像和船舶AIS数据中船舶轨迹点转换为统一的平面坐标；
S2.2：矢量化内河航道，记录内河航道的折点坐标；
S2.3：根据折点位置将内河航道划分为不同的区域，计算船舶轨迹点所属区域；
S2.4：根据所属区域，计算船舶轨迹点对应的里程坐标；
S3：建立可视化模型，包括如下步骤：
S3.1：根据里程坐标判断船舶的航向；
S3.2：沿长度方向将内河航道对应的折线分为若干个线段，并沿航向依次编号，判断各个船舶轨迹点对应里程坐标具体属于哪一个线段；
S3.3：按照时间间隔，根据步骤S1计算得到的船舶轨迹段排放量，计算步骤S3.2中各个线段上所有船舶对应的排放量，绘制热力图展现出污染物排放量的时空分布情况。
2.根据权利要求1所述的内河船舶大气污染物排放清单编制及可视化方法，其特征在于，步骤S2.1中，将内河航道的遥感影像导入到软件ArcGIS中，使用投影关系和变换工具，将遥感影像的坐标系转换为高斯平面直角坐标系；利用高斯正算，将船舶AIS数据中船舶轨迹点坐标，转化为高斯平面直角坐标系中的坐标。
3.根据权利要求2所述的内河船舶大气污染物排放清单编制及可视化方法，其特征在于，步骤S2.3中，将内河航道对应的折线上各个折点从左至右依次进行编号，记为1、2、3、……，M；M表示折线折点的数量；根据各个折点与其两侧相邻两个折点的连线所形成夹角的角平分线，将内河航道划分为不同的区域，角平分线作为区域的边界；对所有区域进行编号，从左至右依次记为1、2、3、……，M+1；
第i个折点处的角平分线方程为：
y-k i x-b i =0; （1）
；（2）
b i= y i -k i x i ; （3）
其中，（x i ，y i ）为第i个折点的平面直角坐标；（x i-1，y i-1）是第i-1个折点的平面直角坐标；（x i+1，y i+1）是第i+1个折点的平面直角坐标；k i 为第i个折点处的角平分线的斜率，即第i-1个折点、第i个折点和第i+1个折点这三点的连线所形成夹角的角平分线斜率，b i 为第i个折点处的角平分线截距。
4.根据权利要求3所述的内河船舶大气污染物排放清单编制及可视化方法，其特征在于，步骤S2.3中，判断各个船舶轨迹点具体属于哪个区域的方法如下：
S2.3.1: 设船舶轨迹上某一点的平面直角坐标为(x 船，y 船)；
S2.3.1.1:设第1折点处的角平分线方程斜率为k 1，第1折点处的角平分线方程为：y 1-k 1 x-b 1 =0；
如果满足条件一：k 1<0，且同时，y 船≤k 1 x 船 + b 1；则，船舶轨迹点(x 船，y 船)归属于第1区域；
或者，如果满足条件二：k 1>0，且同时，y 船＞k 1 x 船 +b 1；则，船舶轨迹点(x 船，y 船)归属于第1区域；
S2.3.1.2：如果既不满足条件一也不满足条件二，则根据第2折点处的角平分线方程来判断；即，
如果满足条件三：k 2<0，且同时，y 船≤k 2 x 船 +b 2；则，船舶轨迹点(x 船，y 船)归属于第2区域；
或者，如果满足条件四：k 2>0，且同时，y 船＞k 2 x 船 +b 2；则，船舶轨迹点(x 船，y 船)归属于第2区域；
……
以此，遍历其余所有折点对应的角平分线方程；
如果不满足前面的所有条件，则船舶轨迹上的该点 (x 船，y 船)归属于第M+1区域；
S2.3.2:重复步骤S2.3.1，对整个船舶轨迹上的所有点进行区域划分。
5.根据权利要求4所述的内河船舶大气污染物排放清单编制及可视化方法，其特征在于，步骤S2.4中，根据舶轨迹点所属区域计算该点对应的里程坐标，具体包括如下步骤：S2.4.1：将内河航道对应折点的平面直角坐标转化为一维的里程坐标，设内河航道对应折线起点的里程坐标为0，并依次计算所有折点对应的里程坐标；
其中，第1个折点的里程坐标S 1：
；（4）
其余折点的里程坐标按照如下方式计算，第i个折点的里程坐标S i ：
；（5）
S2.4.2：当船舶轨迹点(x 船，y 船)投影位于第i个折点和第i+1个折点之间，即船舶轨迹点(x 船，y 船)落入上述计算出的对应区域内时，该点(x 船，y 船)的里程坐标S 船的计算公式如下：
；（6）
（7）。
6.根据权利要求5所述的内河船舶大气污染物排放清单编制及可视化方法，其特征在于，步骤S3.2具体包括如下步骤：
S3.2.1：沿长度方向将内河航道对应的折线分为N个线段，并计算出各个线段终点处的里程坐标S n ，n∈[1，N];
；（8）
其中，S M 表示内河航道折线最末端的折点的里程坐标，l为划分线段的间隔值，为向上取整运算；
S3.2.2：根据各个船舶轨迹点投影在折线上的位置，判断其具体属于哪个线段 。
7.根据权利要求6所述的内河船舶大气污染物排放清单编制及可视化方法，其特征在于，步骤S3.3，计算船舶轨迹在各个线段上对应的排放量，具体为:
设定船舶轨迹的起点位置对应的里程坐标为S 船起、船舶轨迹的终点位置对应的里程坐标为S 船终，该船舶整个轨迹点共落入K个线段中；
设定船舶轨迹的起点位于第n线段中时，对应第n线段上的污染物排放量E n 为：
； （9）
则船舶轨迹的终点位于第n+K线段中，对应第n+K线段上的污染物排放量E n+K 为：
；（10）
其余中间K-2个线段上的污染物排放量均为E n+1 ；
；（11）
其中，E为在船舶整个轨迹段的污染物排放量；
重复该步骤，计算航道上所有船舶在对应线段上的污染物排放量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410356993.XA CN117971959B (zh) | 2024-03-27 | 2024-03-27 | 内河船舶大气污染物排放清单编制及可视化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410356993.XA CN117971959B (zh) | 2024-03-27 | 2024-03-27 | 内河船舶大气污染物排放清单编制及可视化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117971959A true CN117971959A (zh) | 2024-05-03 |
| CN117971959B CN117971959B (zh) | 2024-06-07 |
Family
ID=90866012
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410356993.XA Active CN117971959B (zh) | 2024-03-27 | 2024-03-27 | 内河船舶大气污染物排放清单编制及可视化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117971959B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104680850A (zh) * | 2015-03-19 | 2015-06-03 | 南京大学 | 基于志愿观测船数据的舰船轨迹线生成方法 |
| CN105590024A (zh) * | 2015-12-16 | 2016-05-18 | 河海大学 | 一种基于活动的港区大气污染物排放清单编制方法 |
| WO2022252398A1 (zh) * | 2021-05-31 | 2022-12-08 | 武汉理工大学 | 基于船舶轨迹特征点提取的时空dp方法 |
-
2024
- 2024-03-27 CN CN202410356993.XA patent/CN117971959B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104680850A (zh) * | 2015-03-19 | 2015-06-03 | 南京大学 | 基于志愿观测船数据的舰船轨迹线生成方法 |
| CN105590024A (zh) * | 2015-12-16 | 2016-05-18 | 河海大学 | 一种基于活动的港区大气污染物排放清单编制方法 |
| WO2022252398A1 (zh) * | 2021-05-31 | 2022-12-08 | 武汉理工大学 | 基于船舶轨迹特征点提取的时空dp方法 |
Non-Patent Citations (1)
| Title |
|---|
| 苑帅;封学军;朱逸凡;: "内河航道船舶尾气快速排放清单研究――以长江江苏段为例", 交通运输研究, no. 02, 15 April 2020 (2020-04-15), pages 94 - 103 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117971959B (zh) | 2024-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111060109B (zh) | 一种基于改进a星算法的无人艇全局路径规划方法 | |
| Ma et al. | Method for simultaneously optimizing ship route and speed with emission control areas | |
| CN111814367A (zh) | 港口船舶排放监测监管云服务系统 | |
| CN105675448B (zh) | 一种面向水环境质量监测的动态模拟污染物扩散的方法 | |
| CN106546244B (zh) | 一种针对凸形障航物确定船舶航行局部最优路径的方法 | |
| Wu et al. | Long-voyage route planning method based on multi-scale visibility graph for autonomous ships | |
| Huang et al. | Estimation and spatio-temporal analysis of ship exhaust emission in a port area | |
| CN105539797A (zh) | 一种基于ecdis的风力助航船舶的航行方法及系统 | |
| CN109655063B (zh) | 大型水陆两栖飞机海上搜索航路规划方法 | |
| CN112214721B (zh) | 一种基于ais数据的动态化船舶排放清单建立方法 | |
| CN117971959B (zh) | 内河船舶大气污染物排放清单编制及可视化方法 | |
| CN114239426A (zh) | 基于水流资料同化的长江干线船舶排放清单生成方法 | |
| Zaccone et al. | An optimization tool for ship route planning in real weather scenarios | |
| He et al. | A calculation algorithm for ship pollutant gas emissions and diffusions based on real-time meteorological conditions and its application | |
| Kao et al. | Utilizing the fuzzy IoT to reduce Green Harbor emissions | |
| CN115328903A (zh) | 船舶停留状态识别方法和装置、电子设备和存储介质 | |
| CN115719550B (zh) | 船舶骨干航路网的构建方法、装置和可读储存介质 | |
| CN115077523A (zh) | 一种基于改进Dijkstra的变速航迹规划方法 | |
| Ma et al. | Ship emission accounting based on AIS data mining | |
| CN115905770A (zh) | 一种基于ais数据的船舶污染排放量轨迹测算方法 | |
| CN118035606A (zh) | 一种拖网渔船温室气体排放时空分布计算方法 | |
| CN115146836A (zh) | 基于a星算法的气象航线动态优化方法 | |
| Ma et al. | A two-step optimization method for ship navigation decision-making with emission control area | |
| CN116244557B (zh) | 航道船舶流量监测方法、系统及计算机可读存储介质 | |
| Kamal | Inland electronic navigation charts for enhancing the efficiency of navigational waterway in the Nile River Case Study-El-Wasta Reach |
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 |