CN109801552A - 一种人工海岸线化简方法 - Google Patents
一种人工海岸线化简方法 Download PDFInfo
- Publication number
- CN109801552A CN109801552A CN201910068952.XA CN201910068952A CN109801552A CN 109801552 A CN109801552 A CN 109801552A CN 201910068952 A CN201910068952 A CN 201910068952A CN 109801552 A CN109801552 A CN 109801552A
- Authority
- CN
- China
- Prior art keywords
- artificial
- coastline
- line
- point
- directed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000008602 contraction Effects 0.000 claims description 27
- 238000005452 bending Methods 0.000 claims description 24
- 238000012545 processing Methods 0.000 claims description 11
- 238000012217 deletion Methods 0.000 claims description 7
- 230000037430 deletion Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims 1
- 238000011946 reduction process Methods 0.000 claims 1
- 230000015572 biosynthetic process Effects 0.000 abstract description 4
- 238000003786 synthesis reaction Methods 0.000 abstract description 4
- 238000003909 pattern recognition Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 9
- 239000013598 vector Substances 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Processing Or Creating Images (AREA)
Abstract
本发明涉及一种人工海岸线化简方法，属于空间数据中模式识别及线要素自动综合技术领域。本发明通过识别构成突出式海岸人工建筑的平行弧段组，从中提取有向中心线，对提取出的有向中心线进行合并、自由点和收缩点移动，得到突出式海岸线人工建筑中心线，实现对突出式海岸人工建筑进行化简；从除突出式海岸建筑外的剩余人工海岸线中提取海侧弯曲单元，并进行化简；对化简结果与突出式海岸线人工建筑中心线存在拓扑冲突进行处理，最终实现人工海岸线的自动化简。整个化简过程，不仅确保了化简结果符合目标尺度下制图要求，还保持了丰富、准确的人工海岸特征，解决了现有人工海岸线化简方法导致的形态失真问题。
Description
技术领域
本发明涉及一种人工海岸线化简方法，属于空间数据中模式识别及线要素自动综合技术领域。
背景技术
线要素化简是制图综合中的重要研究内容和经典研究问题，当前线要素自动化简方法有很多，主要包括节点压缩式、弯曲取舍式和智能算法式等。但通用的线化简方法忽略了不同线要素地理特征差异性，难以保证化简后地理特征完整性以及化简前后地理特征一致性，必须研制专门的海岸线化简方法。顾及海图保障航行的安全性需求，需要研制遵循“扩陆缩海”的原则的海图海岸线自动化简方法。如发表于《武汉大学学报(信息科学版)》2011年第12期上一篇名称为《以弯曲骨架线为化简指标的海岸线综合方法》的论文通过改进单调链法识别弯曲，以弯曲骨架线为综合指标取舍海侧弯曲，专门用于海图海岸线化简。不同类型海岸线间形态差异明显、地理特征各异，针对不同类型海岸线也应研制专门的化简方法，如发表于《Marine Geodesy》2014年第2期中一篇名为《A Simplification of RiaCoastline with Geomorphologic Characteristics Preserved》的论文则针对海图中树枝状溺谷海岸线研制专门的化简方法。
上述都是针对自然海岸线的化简，人工海岸线是一种区别于自然海岸线、受人类活动影响显著的常见海岸线。通过人工构筑海岸工程形成的人造海岸线都属于人工海岸线范畴。几何层次看，人工海岸线相对规则、平直；地理层次看，人工海岸包含港口、码头、防波堤等典型海岸人工建筑。依据中国海图图示(GB12319-1998)，常见海岸人工建筑在海图中依比例尺表示和半依比例尺表示形式(其中，不具有半依比例尺表达形式的海岸人工建筑称为“顺岸式”海岸人工建筑，顺岸式海岸人工建筑不能降维表达；具有半依比例尺表达形式的海岸人工建筑称为“突出式”海岸人工建筑，突出式海岸人工建筑能够降维表达，如图1-a和图1-b所示。依据降维表达形态差异，包含引桥的突出式海岸人工建筑称为“T”型，不包含引桥的突出式海岸人工建筑称为折线型)。而目前针对人工海岸线化简方法的研究十分匮乏，仅武汉大学2011年博士论文《海岛(礁)要素地图综合数据模型与方法研究》中建议借用多边形化简方法中的最小二乘法化简人工海岸线，此方法过于强调直角特征而容易导致化简后形态失真。
发明内容
本发明的目的是提供一种人工海岸线的化简方法，以解决目前采用多边形化简方法化简人工海岸线导致的化简后形态失真的问题。
本发明为解决上述技术问题而提供一种人工海岸线化简方法，该化简方法包括以下步骤：
1)人工海岸线描述：
以顶点顺序递增方向为人工海岸线方向，确保人工海岸线方向一侧为海洋，另一侧为陆地，将人工海岸线描述成相邻顶点构成的有向弧段集合；
2)提取有向中心线：
从有向弧段集合中筛选出构成突出式海岸人工建筑的平行弧段组，提取出每个平行弧段组内的两弧段间的中心线段，并以两弧段中顶点顺序相对靠前弧段的方向作为中心线段的方向，得到各平行弧段组的有向中心线；
3)有向中心线初次合并：
对有向中心线按照设定的约束条件进行初次合并，所述的设定的约束条件包括：
若两连续的有向中心线有且只有一个公共顶点，将该公共顶点作为连接点将两者合并；
若两连续的有向中心线没有公共顶点，但具有相交、重叠或者近似重叠的关系，将这两连续的有向中心线进行全连接，并选择全连接中最长的线段作为初次合并的结果；
4)对初次合并后的有向中心线进行端点移动：
所述初次合并后的有向中心线还包含有未参加合并的有向中心线，对初次合并的中心线的自由点和收缩点进行移动，所述自由点指的是有向中心线的终止顶点，所述收缩点指的是有向中心线的起始顶点；
5)利用端点移动后的有向中心线化简突出式海岸人工建筑：
若构成突出式海岸人工建筑的平行弧段组间最小距离小于制图要求的宽度阈值，则将该平行弧段组对应的经过步骤3)和4)处理后的有向中心线作为对应突出式海岸人工建筑的降维表达，实现对突出式海岸人工建筑的化简；
6)识别弯曲单元：
利用拐点法从除突出式海岸建筑外的剩余人工海岸线中提取海侧弯曲单元，删除面积阈值或基线长度小于对应阈值的海侧弯曲单元，得到剩余人工海岸线的化简结果；
7)拓扑冲突处理：
将剩余人工海岸线化简结果与对应原海岸线部分叠加，以删除部分构成的面状多边形集合记为Bm，若存在有向中心线包含于Bm，则说明存在拓扑冲突，若存在拓扑冲突的有向中心线于Bm外的部分不满足制图视觉要求时，将该有向中心线删除，若存在拓扑冲突的有向中心线于Bm外的部分满足制图视觉要求时，将该有向中心线移动至Bm外。
本发明通过识别构成突出式海岸人工建筑的平行弧段组，从中提取有向中心线，对提取出的有向中心线进行合并、自由点和收缩点移动，得到突出式海岸线人工建筑中心线，实现对突出式海岸人工建筑进行化简；从除突出式海岸建筑外的剩余人工海岸线中提取海侧弯曲单元，并进行化简；对化简结果与突出式海岸线人工建筑中心线存在拓扑冲突进行处理，最终实现人工海岸线的自动化简。整个过程，不仅确保了海岸线化简结果符合目标尺度下制图要求，还保持了丰富、准确的人工海岸特征，解决了现有人工海岸线化简方法导致的形态失真问题。
进一步地，为了更好的描述突出式海岸人工建筑的几何特征，本发明还给出了平行弧段组的筛选方法，提高了对突出式海岸人工建筑刻画的精准性，所述步骤2)中筛选出的构成突出式海岸人工建筑的平行弧段组满足以下条件：
a.平行弧段组内的两弧段平行或近似平行；
b.平行弧段组内只包含陆地，不包含海洋；
c.平行弧段组内两弧段间的间距小于两弧段的平均长度。
进一步地，为了保证有向中心线更加完整，还需对有向中心线中的自由点进行移动，所述步骤4)中自由点的移动方式为：
自由点向邻近的初次合并后的中心线的收缩点移动，并以该移动轨迹作为连接线段，与相邻的中心线合并；
自由点沿初次合并的中心线内移动至与人工海岸线交点处停止；
自由点沿初次合并的中心线延长线方向移动至海岸线交点处停止；
自由点沿初次合并的中心线延长线方向移动至另一中心线终止，并以此移动轨迹作为连接线段，合并两中心线。
进一步地，为了保证有向中心线更加完整，还需对有向中心线中的收缩点进行移动，所述步骤4)中收缩点的移动方式为：
当收缩点不满足设定条件时，将收缩点沿有向中心线方向移动至lk1或lk2与有向中心线的交点位置；其中包含收缩点的有向中心线记为lm_s，其对应的平行弧段组的两弧段分别为pl(lm_s)1和pl(lm_s)2，人工海岸线中包含pl(lm_s)1起始顶点且不是pl(lm_s)1的弧段记为lk1，包含pl(lm_s)2终止顶点且不是pl(lm_s)2的弧段记为lk2，lk1-2为pl(lm_s)1起始顶点与pl(lm_s)2终止顶点连结的线段；
所述的设定条件为：(1)lk1、lk2及lk1-2近似共线；(2)lk1、lk2及lk1-2不近似共线，沿lk1方向延长交于pl(lm_s)2的延长线且沿lk2反向延长交于pl(lm_s)1的延长线，或沿lk1方向延长交于pl(lm_s)2内且沿lk2反向延长交于pl(lm_s)1内。
进一步地，为了避免转折点、冗余点的干扰，提高有向弧段对海岸线刻画的准确性，所述步骤1)还包括对所述有向弧段集合进行预处理，预处理过程为：将相邻且近似共线的有向弧段进行合并，所述近似共线指的是两个弧段之间的夹角小于等于设定角度阈值。
进一步地，为了保证化简后的突出式海岸人工建筑的长度满足制图要求，所述步骤5)中对突出式海岸人工建筑的化简还包括若降维处理后的突出式海岸人工建筑小于制图要求的长度阈值，将该突出式海岸人工建筑删除的步骤。
进一步地，为了保证有向中心线的初次合并，本发明还给出了相交、重叠和近似重叠的判断依据，所述步骤3)中的相交指的是两连续的有向中心线有且只有一个交点且交点不是公共顶点；重叠指的是两连续的有向中心线间包含重叠弧段；近似重叠指的是两连续的有向中心线的全连接均包含于海岸线陆地一侧。
附图说明
图1-a是现有常见海岸人工建筑的依比例尺表达；
图1-b是图1-a中常见海岸人工建筑的半依比例尺表达；
图2是本发明实施例中人工海岸线构成的示意图；
图3是本发明实施例中构成突出式海岸人工建筑的平行组的示意图；
图4-a是本发明实施例中具有相接关系的有向中心线初次合并的示意图；
图4-b是本发明实施例中具有相交关系的有向中心线初次合并的示意图；
图4-c是本发明实施例中具有重叠关系的有向中心线初次合并的示意图；
图4-d是本发明实施例中具有近似重叠关系的有向中心线初次合并的示意图；
图5是本发明实施例中两有向中心线间全连接及利用全连接合并有向中心线的示意图；
图6-a是本发明实施例中第一种自由点移动的示意图；
图6-b是本发明实施例中第二种自由点移动的示意图；
图6-c是本发明实施例中第三种自由点移动的示意图；
图6-d是本发明实施例中第四种自由点移动的示意图；
图7是本发明实施例中收缩点移动得到准确降维表达的示意图；
图8-a是本发明实施例中收缩点无需移动情况的示意图；
图8-b是本发明实施例中收缩点移动的示意图；
图9-a是本发明实施例中删除降维表达协调拓扑冲突的示意图；
图9-b是本发明实施例中移动降维表达协调拓扑冲突的示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步地说明。
本发明通过识别构成突出式海岸人工建筑的平行弧段组，从中提取有向中心线，对提取出的有向中心线进行合并、自由点和收缩点移动，得到突出式海岸线人工建筑中心线，实现对突出式海岸人工建筑进行化简(降维)；从除突出式海岸建筑外的剩余人工海岸线中提取海侧弯曲单元，并进行化简；对化简结果与突出式海岸线人工建筑中心线存在拓扑冲突进行处理，最终实现人工海岸线的自动化简。整个化简过程，不仅确保了化简结果符合目标尺度下制图要求，还保持了丰富、准确的人工海岸特征，解决了现有人工海岸线化简方法导致的形态失真问题。
本发明人工海岸线化简方法的实施例
为方便表述，作如下定义：
定义l的起始顶点记为l_f，终止顶点记为l_t。l为l对应的向量，|l|＝l，l的方向是从l的起始顶点指向终止顶点的方向，顶点va向顶点vb的连线记为dl(va,vb)，线要素l的长度记为Len(l)；
某条人工海岸线的有序顶点依次为v1、…、vn，顶点序号递增方向定义为该条海岸线的方向，L内相邻顶点vi-1、vi构成的向量记为li，其中2≤i≤n，L可离散成{li}；同理L'可离散成{li'}，其中2≤i'≤n'；
定义函数Ang(li,lj)用于计算li到lj的角度，设α为近似共线的角度阈值，β为近似平行的角度阈值。
有向中心线lm对应的平行组弧段记为pl(lm)1和pl(lm)2，规定对应有序集合{lx}中，pl(lm)1先于pl(lm)2；
lm、lm+1端点间的全连接记为CC(lm,lm+1)＝{dl(lm_f,lm+1_f),dl(lm_f,lm+1_t),dl(lm_t,lm+1_f),dl(lm_t,lm+1_t)}，CC(lm,lm+1)中距离最大的连接记为Max{CC(lm,lm+1)}。
化简前后比例尺分别记为S1、S2，图上最小可视距离记为svo，线宽记为b。thl表示突出式人工建筑半依比例尺表达形式的长度阈值，thw表示图上线划间宽度阈值，thb表示图上弯曲基线阈值，ths表示图上弯曲面积阈值。通常情况下，thl＝1mm，thw＝2·b+svo，thb＝2·(b+svo)，ths＝thb2，svo＝0.2mm，b＝0.1mm。
下面以图2中的人工海岸线为例对本发明的化简方法进行详细说明，该方法的具体实现过程如下。
1.人工海岸线描述。
以顶点递增方向为人工海岸线方向，确保人工海岸线方向一侧为海洋，另一侧为陆地，将人工海岸线描述成相邻顶点构成的有向弧段集合。
1.1定义顶点递增方向为岸线方向，如图2所示，多引桥的人工建筑的存在使得人工海岸线由岸线主体(L)和闭合环(L')构成，更改与L方向相同的L'，使其方向与L相反。确保沿岸线方向一侧为海洋，另一侧为陆地，其中，L'也可为空。本实施例而言，如图2所示，该人工海岸线的有序顶点依次为v1、…、vn，顶点序号递增方为人工海岸线的方向为从左至右将L与L'离散成相邻顶点构成的有向弧段(相邻顶点vi-1、vi构成的向量记为li)集合。
1.2为避免各有向弧段存在冗余点、假转折点的干扰，还需将有向弧段集合中的多段邻接且近似共线的li合并成lk，其中L、L'中邻接且近似共线向量合并后的有向弧段集合分别记为{lk}、{lk’}。其中，近似共线指的是两个有向弧段的夹角小于α，α为近似共线的角度阈值，一般情况下0°≤α≤10；当li到li+1的角度不大于α时，则二者邻接且近似共线，需要将其合并。
步骤1通过计算机实现的过程如下：
Step1：若存在与L同向的L'，则更新L'内各个顶点下标，令i＝n-i+1，重复此步骤；否则，对L及所有L'依次执行Step2；
Step2：若i＝n(或i'＝n')，则此L(或L')预处理结束；若i<n(或i'<n'，令i＝i'，n＝n')，则执行Step3；
Step3：若Ang(li,li+1)>α，则令lk＝li，i＝i+1，执行Step2；若Ang(li,li+1)≤α，则li与li+1近似共线，新建链表List，li与li+1加入List，执行Step4；
Step4：若不存在Ang(lj,li+2)≤α(其中lj∈List)，则li+2加入List，令i＝i+1，重复Step4；若存在Ang(lj,li+2)>α或i+2＝n，令lk＝List，则lk即为近似共线的弧段，令i＝i+2，转入Step2。
L预处理后，被依次划分为有序弧段集合{lk}，其中1≤k≤u；所有L'预处理后划分的弧段集合依次存入弧段集合{lk’}，其中1≤k'≤u'；lx∈{lk}或{lk'}，顶点按原线对应顶点序号递增的顺序存储，lx方向与原线方向一致，lx为由lx首顶点向末顶点连结成的向量。
2.提取有向中心线。
2.1识别平行弧段组。
识别预处理后的人工海岸线中构成突出式海岸人工建筑的一组平行或近似平行的弧段，称为平行弧段组，构成突出式海岸人工建筑的平行弧段组应满足以下条件：1)平行组内的两弧段平行或近似平行；2)平行组间只包含陆地，不包括海洋；3)平行组的两弧段间间距小于两弧段的平均长度。
对本实施例而言，如图3所示，l3与l7平行，满足条件1)，但不满足条件2)、3)，不能构成平行弧段组；l14与l16满足条件1)、2)，但不满足条件3)，二者间距较大，对应顺岸式人工建筑，也不是构成突出式海岸人工建筑的平行组。其中，l2与l4、l6与l12、l7与l9、l9与l11都是构成突出式海岸人工建筑的平行弧段组。
该步骤通过计算机实现的过程如下：
根据前述原理，针对构成突出式海岸人工建筑的平行组的三个条件构造三个函数F1()～F3()，用于平行组的识别。
F1(lp,lq)：若β<Ang(lp,lq)≤180°，170°≤β≤180°，则返回True；否则，返回False。
F2(lp,lq)：构建辅助线dl(lp_f,lq_t)和dl(lp_t,lq_f)，若两辅助线均包含于陆地内，则返回True；否则，返回False。
F3(lp,lq)：dl(lp_f,lq_t)中点和dl(lp_t,lq_f)中点连结的线段为lm。对于近似平行的lp和lq：对应的lm为lp和lq的中位线，Len(lm)为lp和lq的平均长度；lp、dl(lp_f,lq_t)、lq、dl(lp_t,lq_f)所围成近似梯形的多边形，面积为S；lp与lq为该梯形的平行边，二者间平均距离dis(lp,lq)可以看作此梯形的高，由dis(lp,lq)＝S/Len(lm)求出；若dis(lp,lq)<Len(lm)，则返回lm，|lm|＝lm，lm方向与lp一致，lm称为有向中心线；否则，返回Null。
2.2提取有向中心线：
平行弧段组内两弧段间的中心线段(平行弧段组内不同弧段的起止顶点连线中点所连结而成的线段)即为有向中心线段，有向中心线lm对应的平行组弧段记为pl(lm)1和pl(lm)2，规定对应有序集合{lx}中，pl(lm)1先于pl(lm)2；有向中心线只有起始顶点(标记为收缩点)和终止顶点(标记为自由点)两个顶点，起始顶点向终止顶点的方向即为该有向中心线的方向，该有向中心线是构成突出式海岸人工建筑中心线的基本单元。
通过计算机手段的提取过程为：
从{lk}中l1开始，按如下步骤探测L内、L与L'间的平行组弧段，并依次得到平行组对应的有向中心线集合{lm}。
Step1：若u>2时，则执行下一步；若u＝1且u'>0，则执行Step5；若u＝1且u'＝0，终止；若u＝0，则终止；
Step2：生成计步器y，令y＝2。顾及地理学第一定律，避免计算量过大，可以规定步数阈值为常量thr；
Step3：若F1(l1,l1+y)＝True且F2(l1,l1+y)＝True且F3(l1,l1+y)≠Null，则l1和l1+y为一组平行组弧段，将F3(l1,l1+y)返回的lm加入集合{lm}，将l1从{lk}中移除后，执行Step1；否则，执行下一步；
Step4：若y<thr，则令y＝y+1，重复Step3；否则，执行下一步；
Step5：遍历{lk'}，若存在lk’使得F1(l1,lk')＝True且F2(l1,lk')＝True且F3(l1,lk')≠Null，则l1和lk'构成一个平行组，将F3(l1,lk')返回的lm加入集合{lm}，将对应lk'移除，l1也从{lk}中移除，执行Step1；否则，将l1从{lk}中移除后执行Step1。
类似的，重复以上步骤再探测剩余{lk'}内的平行组及其对应的有向中心线，初步得到有向中心线集合{lm}。
按如下步骤确定有向中心线方向，得到最终的有向中心线集合{lm}。
Step1：将所有lm方向确定为与pl(lm)1方向一致：若Ang(lj,li+2)≤90°，则lm方向与pl(lm)1方向一致；否则，交换lm的端点，改变方向；执行下一步；
Step2：依次判断所有lm对应的平行组弧段是否与L'相离，若均相离，则终止；否则，执行下一步；
Step3：若存在lm对应的平行组弧段与L'中闭合环不相离，则lm与对应闭合环相关。
3.有向中心线初次合并。
3.1若两连续的有向中心线有且只有一个公共顶点，即二者相接，则将公共顶点作为连接点直接将二者合并，初次合并后的中心线的起止顶点的标记与合并前有向中心线端点的标记一致，示意过程如图4-a所示。
3.2若两连续的有向中心线具有相交、重叠或近似重叠的的关系，有两向中心线进行全连接(即两有向中心间所有顶点间进行连接)，并选择全连接中最长的线段作为初次合并的结果，实现初次合并，初次合并后的中心线的起止顶点的标记与合并前有向中心线端点的标记一致。
相交指的是两连续的有向中心线有且只有一个交点且交点不是公共顶点，则二者相交；重叠指的是两连续的有向中心线间包含重叠弧段，则二者重叠；近似重叠指的是两连续的有向中心线的全连接均包含于海岸线陆地一侧，则二者进行重叠，具有相交、重叠、近似重叠关系两有向中心线初次合并的示意过程分别如图4-b、图4-c、图4-d所示。其原理如图5所示，获取两有向中心线全连接，并从中选取最长的线段作为合并结果。
初次合并后的中心线(包含未参加合并的有向中心线)的端点只存在两种标记情况：端点均被标记为自由点(终止顶点)、一个端点被标记为自由点和一个端点被标记为收缩点(起始顶点)。
通过计算机的实现手段为：
Step1：新建集合List{}用于存储待合并的有向中心线，将lm加入到List{}，执行下一步。
Step2：若lm与lm+1有公共顶点，lm+1加入List{}，令m＝m+1，重复此步骤；若lm与lm+1无公共顶点，执行下一步。
Step3：若lm与lm+1有且只有一个交点，且交点不是公共顶点，则令lm＝Max{CC(lm,lm+1)}且lm+1＝Max{CC(lm,lm+1)}，更新List{}中的lm，令m＝m+1，执行Step2；否则，执行下一步。
Step4：若lm与lm+1存在重叠部分，则令lm＝Max{CC(lm,lm+1)}且lm+1＝Max{CC(lm,lm+1)}，更新List{}中的lm，令m＝m+1，执行Step2；否则，执行下一步。
Step5：若CC(lm,lm+1)均包含于海岸线陆地侧，则令lm＝Max{CC(lm,lm+1)}且lm+1＝Max{CC(lm,lm+1)}，更新List{}中的lm，令m＝m+1，执行Step2；否则，执行下一步。
Step6：若List{}中包含不止一个元素，则以公共顶点作为连接点，得到初次合并的中心线；否则，List{}中的唯一元素即为初次合并后的中心线。令m＝m+1，执行Step1。
初次合并后的中心线集合记为{lm'}，lm'端点的标记与合并前有向中心线起止顶点标记相同。
4.对初次合并后的有向中心线进行端点移动。
4.1自由点移动
被标记为自由点的端点的移动方式如下：
1)自由点向邻近的初次合并的中心线的收缩点移动，并以此移动轨迹作为连结线段，与相邻的中心线合并，如图6-a所示；
2)自由点沿初次合并的中心线内移动至与海岸线交点处停止，如图6-b所示；
3)自由点沿初次合并的中心线延长线方向移动至与海岸线交点处停止，如图6-c所示；
4)自由点沿初次合并的中心线延长线方向移动至另一中心线终止，且以此移动轨迹作为连结线段，合并两中心线，如图6-d所示。
自由点移动后，中心线合并完成，合并后的中心线对应的平行弧段组中下标最小的顶点与最大的顶点间所夹的弧段即为初步识别的突出式海岸人工建筑。自由点移动过程中伴随有有向中心线的再次合并。通过计算机的实现过程为：
Step1：若lm'自由点与l(m±1)'收缩点连线位于陆侧，且是CC(lm',l(m±1)')中包含于陆地一侧连结中距离最小的，则连接lm'自由点与l(m±1)'收缩点，且合并lm'与l(m±1)'；否则，执行下一步；
Step2：若lm'与L相交，则自由点沿lm'向内移动至lm'与L的交点处，终止移动；否则，执行下一步；
Step3：自由点沿lm'延长线向外移动，若与l(m±1)'相交，则终止移动，执行下一步；否则，自由点沿lm'延长线向外移动，直至与L相交，终止；
Step4：若交点是l(m±1)'的端点或l(m±1)'中不含收缩点，则合并lm'与l(m±1)'即可；若交点不是l(m±1)'的端点且l(m±1)'包含收缩点，则l(m±1)'收缩点沿l(m±1)'向内移动至交点时终止，且合并lm'与l(m±1)'。
4.2收缩点移动
自由点移动引发合并后的中心线集合记为{lm”}，lm”对应一个初步识别的、独立的突出式人工建筑的中心线。在此基础上，还要通过移动收缩点，如图7所示，处理合并后的中心线lm”，得到最终的中心线结果。
其中包含收缩点的有向中心线记为lm_s，其对应平行组的弧段分别为pl(lm_s)1和pl(lm_s)2，预处理后的人工海岸线中包含pl(lm_s)1起始顶点且不是pl(lm_s)1的弧段记为lk1，包含pl(lm_s)2终止顶点且不是pl(lm_s)2的弧段记为lk2，lk1-2为pl(lm_s)1起始顶点与pl(lm_s)2终止顶点连结的线段。收缩点在以下两种情况下无需移动，如图8-a所示，图8-a中的图例与图8-b中的一样。
1)lk1、lk2及lk1-2近似共线，收缩点不移动；
2)lk1、lk2及lk1-2不近似共线时，若沿lk1方向延长交于pl(lm_s)2的延长线且沿lk2反向延长交于pl(lm_s)1的延长线，或沿lk1方向延长交于pl(lm_s)2内且沿lk2反向延长交于pl(lm_s)1内，也不用移动收缩点。
当不满足上述两种情况时，收缩点沿有向中心线方向移动至与该有向中心线的交点位置，并向初步识别的突出式海岸人工建筑中插入新定点用于更新对应平行组弧段的起始(或终止)顶点，进而对突出式海岸建筑识别结果进行更新，得到最终识别结果，其过程如图8-b所示。通过计算机实现收缩点移动的具体过程为：
Step1：连结pl(lm_s)1_f和pl(lm_s)2_t得到dl(pl(lm_s)1_f,pl(lm_s)2_t)，若Ang(lk1,dl(pl(lm_s)1_f,pl(lm_s)2_t))≤α且Ang(dl(pl(lm_s)1_f,pl(lm_s)2_t),lk2)≤α，即lk1、dl(pl(lm_s)1_f,pl(lm_s)2_t)、lk2近似共线，则不做处理；否则，执行下一步；
Step2：若沿lk1方向延长交于pl(lm_s)2的延长线且沿lk2反向延长交于pl(lm_s)1的延长线，或沿lk1方向延长交于pl(lm_s)2内且沿lk2反向延长交于pl(lm_s)1内，则不做处理；否则，执行下一步；
Step3：沿lk1方向(或lk2反向)延长交于pl(lm_s)2(或pl(lm_s)1)内，交点作为L中新插入的顶点记为vm_i，且沿lk2反向(或lk1方向)延长交于pl(lm_s)1反向延长线(或pl(lm_s)2延长线)上，将lm_s收缩点移动至lm_s与dl(pl(lm_s)1_f,vm_i)(或dl(pl(lm_s)2_t,vm_i))的交点，则改变原中心线。将pl(lm_s)2_t(或pl(lm_s)1_f)替换为vm_i，对应平行组的弧段也得到处理。
有向中心线合并过程中，对应的平行组加入同一集合；处理后得到最终中心线lm”'及其对应的平行组集合；获取平行组集合中所有线要素端点对应于L中顶点的下标，下标最小值与最大值顶点记为vm_min、vm_max，L中vm_min、vm_max及其间所有顶点构成lm”’对应的突出式人工建筑岸线。
5.突出式海岸人工建筑化简。
突出式海岸人工建筑中所有平行组内最小距离小于制图要求的宽度阈值时，将通过上述过程提取的有向中心线作为该突出式人工建筑的降维表达，实现降维处理；当降维表达后的突出式海岸人工建筑存在小于制图要求长度阈值的部分，则删除此突出式人工建筑。
提取的突出式人工建筑岸线集合记为{LT}，其余部分岸线依次记为{LS}，LT与LS交替出现。一般情况下突出式人工建筑间相互独立，化简处理互不影响，{LT}化简只考虑LT化简即可。
LT的端点分别为vm_min、vm_max，LT对应中心线为lm”'，lm”'由有向中心线集合{lm}构成，lm”'对应平行组集合{{pl(lm)1,pl(lm)2}}，pl(lm)1与pl(lm)2间的最小距离为dis_min(pl(lm)1,pl(lm)2)。通常情况下，thl＝1mm，b＝0.1mm。
dis_min(pl(lm)1,pl(lm)2)
＝Min{Dis(pl(lm)1_f,pl(lm)2),Dis(pl(lm)1_t,pl(lm)2),Dis(pl(lm)2_f,pl(lm)1),Dis(pl(lm)2_t,pl(lm)1}
该步骤通过计算机的实现过程为：
Step1：若Min{dis_min(pl(lm)1,pl(lm)2)}/S1≥thw/S2，则LT不做处理；否则，删除vm_min与vm_max间所有顶点，将lm”'作为LT半依比例尺的表达结果，执行Step2。
Step2：lm””进行拓扑化处理得到{lm””}。若{lm””}数量为1时，对应突出式人工建筑为折线型，执行Step3；当{lm””}数量不为1时，对应突出式人工建筑为T型，利用垂直关系将{lm””}分为引桥集合{lm””_b}和堤岸集合{lm””_d}，执行Step4。
Step3：若Min{Len(dl(lm””_f,lm””_t))，Len(lm””)}/S1<thl/S2，则删除lm”'；否则，不作处理。
Step4：若Min{Min{Len(lm””_b)},∑{lm””_d}}/S1<thl/S2，则删除lm”'；否则，不作处理。
6.识别弯曲单元。
通过上述过程，可实现对人工海岸线中突出式海岸建筑的化简，但是人工海岸线还包括除突出式海岸建筑外的剩余人工海岸线，剩余人工海岸线形态相对规则简单，与海湾、岬角等地理要素一样均可作为一般弯曲单元进行识别、取舍。首先利用拐点法从除突出式海岸建筑外的剩余人工海岸线中识别弯曲单元，顾及“扩陆缩海”约束，提取海侧弯曲单元即可；探后通过弯曲面积、弯曲基线长度等量化衡量所识别的海侧弯曲单元，删除小于制图要求面积阈值或基线长度阈值的海侧弯曲单元；重复此识别、删除过程，直至所有不符合目标尺度下制图要求的细小海侧弯曲单元均被删除为止。
对本实例而言，可采用Springer出版的第773卷CCIS中《A Progressive Methodof Simplifying Polylines with Multi-bends Simplification Templates》中改进的拐点法识别弯曲单元，若识别弯曲单元的基线(起止顶点连结的线段)包含于海洋内时，该单元为海侧单元，识别的海侧单元集合记为{U(i)}。该步骤通过计算机实现时采用的具体手段为：
构造函数Base(U(i))返回海侧单元U(i)的基线长度；函数Area(U(i)))返回海侧单元U(i)及其基线所围成的多边形的面积。若Base(U(i))/S1<thb/S2或Area(U(i))/S1<thb/S2，则删除U(i)起止顶点间的所有顶点，即可实现该单元删除。遍历{U(i)}，删除所有需要删除的海侧弯曲，得到初次的化简结果。对上一次化简后的部分重复海侧单元识别、弯曲删除，直至没有弯曲被删除为止，完成化简。
7.拓扑冲突处理。
步骤5对突出式海岸人工建筑的化简和步骤6对其余人工海岸线的化简可能会产生拓扑冲突。将除突出式海岸人工建筑外的其余人工海岸线的化简结果与对应原海岸线进行叠加，已删除部分构成的面状多边形集合记为{Bm}，若Bm与突出式海岸人工建筑降维表达的结果间存在重叠部分，则化简结果存在拓扑冲突。
突出式海岸人工建筑的降维表达结果较多地包含于Bm内(降维表达结果于Bm外的部分不满足制图视觉要求)时，则将该降维表达删除，示意过程如图9-a所示。降维表达结果较少地包含于Bm内(降维表达结果于Bm外的部分足够满足制图视觉要求)时，则移动该降维表达，示意过程如图9-b所示。
将包含于Bm内的部分记为l1m”'，其余部分为l2m”'。对l1m”'、l2m”'进行拓扑化处理得到{l1m””}、{l2m””}，化简前后比例尺分别记为S1、S2，图上最小可视距离记为svo，通常情况下，svo＝0.2mm，lm”'与Bm的冲突处理步骤如下：
Step1：若l2m”'为空，则lm”'包含于Bm，直接删除lm”'；否则，执行下一步；
Step2：若Min{Len(l2m””)}/S1<svo/S2，则删除lm”'；否则，执行下一步；
Step3：Min{Len(l1m””)}记为l1m””_min，l1m””_min与lm”'有公共端点，lm”'从公共端点开始沿l1m””_min移动至l1m””_min的另一端点为止。移动结束后，若仍存在包含于Bm内的部分，则用弯曲基线截断移动后的lm”'，删除包含于Bm内的部分。
通过上述化简过程可知，本发明能够确保化简结果符合目标尺度下制图要求的同时，保持了丰富、准确的海岸人工特征，顾及航行安全，具有良好的拓扑一致性，是可以应用于海图综合的自动化方法。
Claims (7)
1.一种人工海岸线化简方法，其特征在于，该化简方法包括以下步骤：
1)人工海岸线描述：
以顶点顺序递增方向为人工海岸线方向，确保人工海岸线方向一侧为海洋，另一侧为陆地，将人工海岸线描述成相邻顶点构成的有向弧段集合；
2)提取有向中心线：
从有向弧段集合中筛选出构成突出式海岸人工建筑的平行弧段组，提取出每个平行弧段组内的两弧段间的中心线段，并以两弧段中顶点顺序相对靠前的弧段的方向作为中心线段的方向，得到各平行弧段组的有向中心线；
3)有向中心线初次合并：
对有向中心线按照设定的约束条件进行初次合并，所述的设定的约束条件包括：
若两连续的有向中心线有且只有一个公共顶点，将该公共顶点作为连接点将两者合并；
若两连续的有向中心线没有公共顶点，但具有相交、重叠或者近似重叠的关系，将这两连续的有向中心线进行全连接，并选择全连接中最长的线段作为初次合并的结果；
4)对初次合并后的有向中心线进行端点移动：
所述初次合并后的有向中心线还包含有未参加合并的有向中心线，对初次合并的中心线的自由点和收缩点进行移动，所述自由点指的是有向中心线的终止顶点，所述收缩点指的是有向中心线的起始顶点；
5)利用端点移动后的有向中心线化简突出式海岸人工建筑：
若构成突出式海岸人工建筑的平行弧段组间最小距离小于制图要求的宽度阈值，则将该平行弧段组对应的经过步骤3)和4)处理后的有向中心线作为对应突出式海岸人工建筑的降维表达，实现对突出式海岸人工建筑的化简；
6)识别弯曲单元：
利用拐点法从除突出式海岸建筑外的剩余人工海岸线中提取海侧弯曲单元，删除面积阈值或基线长度小于对应阈值的海侧弯曲单元，得到剩余人工海岸线的化简结果；
7)拓扑冲突处理：
将剩余人工海岸线化简结果与对应原海岸线部分叠加，将删除部分构成的面状多边形集合记为Bm，若存在有向中心线包含于Bm，则说明存在拓扑冲突，若存在拓扑冲突的有向中心线于Bm外的部分不满足制图视觉要求时，将该有向中心线删除，若存在拓扑冲突的有向中心线于Bm外的部分满足制图视觉要求时，将该有向中心线移动至Bm外。
2.根据权利要求1所述的人工海岸线化简方法，其特征在于，所述步骤2)中筛选出的构成突出式海岸人工建筑的平行弧段组满足以下条件：
a.平行弧段组内的两弧段平行或近似平行；
b.平行弧段组内只包含陆地，不包含海洋；
c.平行弧段组内两弧段间的间距小于两弧段的平均长度。
3.根据权利要求1所述的人工海岸线化简方法，其特征在于，所述步骤4)中自由点的移动方式为：
自由点向邻近的初次合并后的中心线的收缩点移动，并以该移动轨迹作为连接线段，与相邻的中心线合并；
自由点沿初次合并的中心线内移动至与人工海岸线交点处停止；
自由点沿初次合并的中心线延长线方向移动至海岸线交点处停止；
自由点沿初次合并的中心线延长线方向移动至另一中心线终止，并以此移动轨迹作为连接线段，合并两中心线。
4.根据权利要求1或3所述的人工海岸线化简方法，其特征在于，所述步骤4)中收缩点的移动方式为：
当收缩点不满足设定条件时，将收缩点沿有向中心线方向移动至lk1或lk2与有向中心线的交点位置；其中包含收缩点的有向中心线记为lm_s，其对应的平行弧段组的两弧段分别为pl(lm_s)1和pl(lm_s)2，人工海岸线中包含pl(lm_s)1起始顶点且不是pl(lm_s)1的弧段记为lk1，包含pl(lm_s)2终止顶点且不是pl(lm_s)2的弧段记为lk2，lk1-2为pl(lm_s)1起始顶点与pl(lm_s)2终止顶点连结的线段；
所述的设定条件为：(1)lk1、lk2及lk1-2近似共线；(2)lk1、lk2及lk1-2不近似共线，沿lk1方向延长交于pl(lm_s)2的延长线且沿lk2反向延长交于pl(lm_s)1的延长线，或沿lk1方向延长交于pl(lm_s)2内且沿lk2反向延长交于pl(lm_s)1内。
5.根据权利要求1所述的人工海岸线化简方法，其特征在于，所述步骤1)还包括对所述有向弧段集合进行预处理，预处理过程为：将相邻且近似共线的有向弧段进行合并，所述近似共线指的是两个弧段之间的夹角小于等于设定角度阈值。
6.根据权利要求4所述的人工海岸线化简方法，其特征在于，所述步骤5)中对突出式海岸人工建筑的化简还包括若降维处理后的突出式海岸人工建筑小于制图要求的长度阈值，将该突出式海岸人工建筑删除的步骤。
7.根据权利要求1所述的人工海岸线化简方法，其特征在于，所述步骤3)中的相交指的是两连续的有向中心线有且只有一个交点且交点不是公共顶点；重叠指的是两连续的有向中心线间包含重叠弧段；近似重叠指的是两连续的有向中心线的全连接均包含于海岸线陆地一侧。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910068952.XA CN109801552B (zh) | 2019-01-24 | 2019-01-24 | 一种人工海岸线化简方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910068952.XA CN109801552B (zh) | 2019-01-24 | 2019-01-24 | 一种人工海岸线化简方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109801552A true CN109801552A (zh) | 2019-05-24 |
| CN109801552B CN109801552B (zh) | 2020-11-27 |
Family
ID=66560313
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910068952.XA Active CN109801552B (zh) | 2019-01-24 | 2019-01-24 | 一种人工海岸线化简方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109801552B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110569532A (zh) * | 2019-07-19 | 2019-12-13 | 中国地质大学（武汉） | 一种室内边界要素矩形的拓扑一致性优化方法 |
| CN112926465A (zh) * | 2021-03-02 | 2021-06-08 | 中国人民解放军战略支援部队信息工程大学 | 基于点云类型的海岸线性质识别方法及装置 |
| CN115100322A (zh) * | 2022-07-23 | 2022-09-23 | 中国人民解放军战略支援部队信息工程大学 | 一种监督学习支持的多化简算法协同的线要素自适应化简方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101751682A (zh) * | 2009-12-31 | 2010-06-23 | 武汉大学 | 一种通用地图的制备方法 |
| CN102609898A (zh) * | 2012-02-09 | 2012-07-25 | 武汉大学 | 一种顾及地理特征的溺谷海岸线化简方法 |
| CN102708405A (zh) * | 2012-04-20 | 2012-10-03 | 武汉大学 | 一种基于克隆选择算法的线要素化简方法 |
| CA3001119A1 (en) * | 2015-11-09 | 2017-05-18 | Landmark Graphics Corporation | Modelling complex geological sequences using geologic rules and paleographic maps |
| CN108537116A (zh) * | 2018-03-05 | 2018-09-14 | 中国地质大学（武汉） | 一种基于多尺度特征的海岸线二级类型提取方法及系统 |
| CN108765440A (zh) * | 2018-05-31 | 2018-11-06 | 大连海事大学 | 一种单极化sar图像的线引导超像素海岸线提取方法 |
-
2019
- 2019-01-24 CN CN201910068952.XA patent/CN109801552B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101751682A (zh) * | 2009-12-31 | 2010-06-23 | 武汉大学 | 一种通用地图的制备方法 |
| CN102609898A (zh) * | 2012-02-09 | 2012-07-25 | 武汉大学 | 一种顾及地理特征的溺谷海岸线化简方法 |
| CN102708405A (zh) * | 2012-04-20 | 2012-10-03 | 武汉大学 | 一种基于克隆选择算法的线要素化简方法 |
| CA3001119A1 (en) * | 2015-11-09 | 2017-05-18 | Landmark Graphics Corporation | Modelling complex geological sequences using geologic rules and paleographic maps |
| CN108537116A (zh) * | 2018-03-05 | 2018-09-14 | 中国地质大学（武汉） | 一种基于多尺度特征的海岸线二级类型提取方法及系统 |
| CN108765440A (zh) * | 2018-05-31 | 2018-11-06 | 大连海事大学 | 一种单极化sar图像的线引导超像素海岸线提取方法 |
Non-Patent Citations (3)
| Title |
|---|
| 刘鹏程: "贝叶斯模型下面向地图表达的海岸线分段模型", 《计算机工程与应用》 * |
| 杜佳威: "一种面向化简的海岸线弯曲识别与结构化表达方法", 《测绘科学技术学报》 * |
| 杜佳威: "采用多元弯曲组划分的线要素化简方法", 《计算机辅助设计与图形学学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110569532A (zh) * | 2019-07-19 | 2019-12-13 | 中国地质大学（武汉） | 一种室内边界要素矩形的拓扑一致性优化方法 |
| CN110569532B (zh) * | 2019-07-19 | 2023-04-21 | 中国地质大学（武汉） | 一种室内边界要素矩形的拓扑一致性优化方法 |
| CN112926465A (zh) * | 2021-03-02 | 2021-06-08 | 中国人民解放军战略支援部队信息工程大学 | 基于点云类型的海岸线性质识别方法及装置 |
| CN115100322A (zh) * | 2022-07-23 | 2022-09-23 | 中国人民解放军战略支援部队信息工程大学 | 一种监督学习支持的多化简算法协同的线要素自适应化简方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109801552B (zh) | 2020-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109801552B (zh) | 一种人工海岸线化简方法 | |
| CN110176018B (zh) | 一种保持结构化地物轮廓特征的图斑合并方法 | |
| CN107845139B (zh) | 一种狭长图斑分裂线的处理方法 | |
| CN110222642B (zh) | 一种基于全局图聚类的平面建筑构件点云轮廓提取方法 | |
| CN110543885B (zh) | 一种交互式提取高分辨率遥感影像道路并生成路网的方法 | |
| CN107330901B (zh) | 一种基于骨架的物体构件分解方法 | |
| Shen et al. | A polygon aggregation method with global feature preservation using superpixel segmentation | |
| CN102682479A (zh) | 一种在不规则三角网上进行三维地形特征点生成的方法 | |
| CN101251929B (zh) | 基于Delaunay三角网技术的区域拓扑关系构建方法 | |
| CN107092978A (zh) | 一种面向虚拟地球的最短路径分层规划方法 | |
| CN111858810A (zh) | 一种面向道路dem构建的建模高程点筛选方法 | |
| CN109034234A (zh) | 一种毗邻区特征识别及自动处理方法 | |
| Kastrisios et al. | A cohesive methodology for the delimitation of maritime zones and boundaries | |
| CN115359149A (zh) | 一种泛化拓扑结构约束下的连续图斑化简方法 | |
| CN111881919A (zh) | 一种基于追踪式格网剖分的线要素智能化简方法与装置 | |
| CN111428811A (zh) | 一种单环多边形自相交模式识别及处理方法 | |
| US7181377B1 (en) | Method of modifying a volume mesh using sheet extraction | |
| CN107341823A (zh) | 一种基于特征融合的最小分支立体匹配方法 | |
| CN106296822B (zh) | 一种网格模型的表面边界恢复方法 | |
| Oh et al. | The geodesic 2-center problem in a simple polygon | |
| CN113763563B (zh) | 一种基于平面识别的三维点云几何网格结构生成方法 | |
| CN112362073B (zh) | 一种基于凸包特征的导航路径建模及最短路径求解方法 | |
| Kenzhebay | Planar roof structure extraction from Very High-Resolution aerial images and Digital Surface Models using deep learning | |
| Dunkars | Multiple representation databases for topographic information | |
| Bildirici et al. | A generic approach to building typification |
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 |