CN111008598B - 一种基于缓冲区生长模型的岛屿自动选取方法 - Google Patents

一种基于缓冲区生长模型的岛屿自动选取方法 Download PDF

Info

Publication number
CN111008598B
CN111008598B CN201911240476.1A CN201911240476A CN111008598B CN 111008598 B CN111008598 B CN 111008598B CN 201911240476 A CN201911240476 A CN 201911240476A CN 111008598 B CN111008598 B CN 111008598B
Authority
CN
China
Prior art keywords
island
buffer
row
width
islands
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
Application number
CN201911240476.1A
Other languages
English (en)
Other versions
CN111008598A (zh
Inventor
唐露露
张立华
贾帅东
戴泽源
杨一曼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
PLA Dalian Naval Academy
Original Assignee
PLA Dalian Naval Academy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by PLA Dalian Naval Academy filed Critical PLA Dalian Naval Academy
Priority to CN201911240476.1A priority Critical patent/CN111008598B/zh
Publication of CN111008598A publication Critical patent/CN111008598A/zh
Application granted granted Critical
Publication of CN111008598B publication Critical patent/CN111008598B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V30/00—Character recognition; Recognising digital ink; Document-oriented image-based pattern recognition
    • G06V30/40—Document-oriented image-based pattern recognition
    • G06V30/42—Document-oriented image-based pattern recognition based on the type of document
    • G06V30/422—Technical drawings; Geographical maps
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume

Abstract

本发明公开了一种基于缓冲区生长模型的岛屿自动选取方法，包括以下步骤：读取待综合区域的全部海图数据并提取其中的岛屿数据；计算各岛屿之间的最小距离，构建岛屿距离矩阵，记录各岛屿属性；为各岛屿设置是否可继续生长的缓冲区标识和生长因子，并根据上述构建的岛屿距离矩阵及属性要素生成各岛屿缓冲区的宽度，根据岛屿缓冲区宽度，构建岛屿缓冲区；按照面积删除其中首个岛屿；根据剩余的待选取岛屿及其宽度，对模型进行更新；重复上述步骤直至剩余待选取岛屿数目达到阈值；输出选取出的岛屿编号。本方法通过构建缓冲区生长模型，识别并删除岛屿群中密度最大的岛屿，并设计了缓冲区生长模型的动态更新方法，可明显提高海图岛屿选取质量。

Description

一种基于缓冲区生长模型的岛屿自动选取方法
技术领域
本发明涉及海洋测绘技术领域，尤其涉及一种岛屿自动选取方法。
背景技术
岛屿作为海图上一类重要的地理要素，通常以点状要素或面状要素的形式在海图上大量存在，其综合的优劣程度直接影响到最终海图的质量。长期以来，岛屿综合通常由制图作业人员根据规范要求、结合自身经验，采用手工方式来完成的。
近年来，随着计算机制图技术的不断发展，有学者开始探索岛屿的自动综合算法，但由于岛屿的自动选取具有极大的复杂性和难度，一直未在方法上取得突破。陈子澎探讨了Voronoi图用于岛屿进行分区域自动选取的可能性；陆毅进一步在用Voronoi图对岛屿密集程度进行评估的基础上，优先删除密度最大的内部岛屿，并对边界岛屿借鉴曲线化简的原理进行取舍，实现了岛屿的自动选取。黄亚峰基于Gestalt原则，对岛屿群的骨架岛屿以及边界岛屿进行识别的基础上，对剩余岛屿采用岛屿群空间等剖分结构(类似于多边形的Voronoi图)对“密度”较大的岛屿进行识别并删除。
当前这些算法大多是在Voronoi图对岛屿的密集程度进行量化评估的基础上，优先对分布较为密集区域的岛屿进行综合，以达到选取的目的。但经过研究发现，使用Voronoi图在某些特殊区域并不能对要素的密集程度进行有效描述，从而导致选取结果难以满足生产实际的需求，制约了海图自动综合的研究。
发明内容
为了克服当前岛屿的自动综合方法不能满足生产实际需要的不足，本发明提供了一种基于缓冲区动态生长的岛屿自动选取方法。
本发明为实现上述目的所采用的技术方案是：一种基于缓冲区生长模型的岛屿自动选取方法，包括以下步骤：
a、读取海图上待选取区域的岛屿平面坐标数据；
b、利用读取的岛屿平面坐标数据，计算各岛屿之间的最小距离，并构建岛屿距离矩阵；
c、为各个待选取岛屿设置缓冲区标识以及增长因子，并依据距离矩阵生成缓冲区宽度；
d、提取各缓冲区宽度，为各个岛屿构建缓冲区，计算个岛屿缓冲区的面积；
e、根据上述步骤所计算的缓冲区面积，将待选取岛屿按照构建的缓冲区面积由大到小排序；
f、判断是否存在面积大小相等的缓冲区，若是则进入步骤g，否则进入步骤h；
g、提取其中缓冲区面积大小相等的岛屿，获取其相邻缓冲区数目，并判断缓冲区数目是否相等，若相等，则计算岛屿与其周围岛屿间平均间距，将其平均间距较大者置于前部，否则将其中相邻缓冲区数目较多的岛屿放于前侧，随后转步骤h；
h、删除其中第一个岛屿；
i、判断剩余岛屿是否达到选取指标，若是，则进入步骤j，否则，将剩余岛屿设置为待选取岛屿，并进入步骤c；
j、提取选取的岛屿，更新海图数据，输出选取结果。
所述步骤b中，构建的岛屿距离矩阵概念为：求取各岛屿与其他所有岛屿之间的最小距离，并将其放在距离矩阵S中，
Figure BDA0002306073570000021
以元素为例，sij表示岛屿i与岛屿j之间的最小距离。
所述步骤c中，各个待选取岛屿缓冲区标识以及生长因子为：
(1)岛屿缓冲区标识
岛屿缓冲区标识是对岛屿的缓冲区宽度以及宽度是否可以更改进行标识的式子，其表达形式为：
Figure BDA0002306073570000022
其中sym_changei用于标识第i个岛屿的缓冲区宽度是否可更改，其中di表示第i个岛屿的缓冲区宽度，在初始状态时，将宽度是否可更改标识全部设置为可更改，将全部岛屿的缓冲区宽度设置为0，sym_elementi表示第i个岛屿的属性；
(2)岛屿缓冲区宽度矩阵
为了之后方便生成缓冲区，设置两个宽度矩阵，分别为行主导矩阵，列主导矩阵：
Figure BDA0002306073570000023
Figure BDA0002306073570000024
其中行主导矩阵Drow的任意一行的所有元素相同，以第i行为例，这一行的值均为di，di为SYM中第i个岛屿的缓冲区宽度，同理，对于列主导矩阵而言，其任意一列的所有元素均相同，其值为第i个岛屿的缓冲区宽度；
(3)缓冲区宽度增长因子
由于需要使岛屿的影响域宽度不断地增大，因此设置一个增长因子w，使w以一定的速率逐渐增大，w的初始值设置为0。
所述步骤c中，基于岛屿距离矩阵生成缓冲区的方法为：
第一步：为岛屿缓冲区标识SYM、岛屿缓冲区宽度矩阵Drow，Dcol赋值。其过程为：使增长因子w增大，获取全部SYM中sym_change设置为可更改的岛屿，将这些岛屿的缓冲区宽度d设置为w；将宽度矩阵中的行主导矩阵Drow每一行的值都设置为SYM中的对应岛屿缓冲区宽度值，例如，取行主导矩阵Drow的第i行的全部元素，记为Drow(i,),令Drow(i,)＝di，其中di为SYM中第i个岛屿的缓冲区宽度，即将第i行的全部元素大小均设置为di；对于列主导矩阵Dcol赋值过程与行主导矩阵类似；
第二步：判断当前缓冲区宽度下，是否有岛屿的缓冲区相交，利用第一步设置好的Drow，Dcol，进行矩阵运算S-Drow-Dcol，其结果矩阵记为Sres，在此结果矩阵中，任意一元素Sres(i,j)表达了在当前SYM中所述的岛屿缓冲区宽度下，岛屿i与岛屿j缓冲区之间的最小距离，遍历结果矩阵Sres，获取其中最小的元素Sres(row,col)，判断Sres(row,col)是否大于0，若是，则表明在当前缓冲区宽度下，各岛屿的缓冲区没有相切或相交，返回第一步；若小于0，则表明在此缓冲区宽度下，岛屿row与岛屿col的缓冲区相交，此时将row与col进行记录，进入第三步；若等于0，此时无法判断是否是已经完成宽度求取的两岛屿缓冲区，因此在这里转入第一步；
第三步：为出现相交的岛屿缓冲区设置最终的宽度值，并对增长因子重新赋值，利用第二步获取的row与col，获取这两个岛屿的相关信息，即第row个岛屿与第col个岛屿之间的最小距离srowcol，岛屿缓冲区标识SYM中的sym_changerow，sym_changecol，sym_elementrow，sym_elementcol；针对不同的sym_changerow，sym_changecol，分以下3中情况对增长因子w，SYM中的岛屿缓冲区宽度d进行如下赋值：
①当sym_changerow，sym_changecol均标识为可更改时
增长因子：
Figure BDA0002306073570000041
岛屿缓冲区宽度：
Figure BDA0002306073570000042
并将sym_changerow，sym_changecol设置为不可更改；
②当sym_changerow标识为可更改，sym_changecol标识为不可更改时
增长因子：
w＝srowcol-dcol (8)
岛屿缓冲区宽度：
drow＝w (9)
将sym_changerow设置为不可更改；
③当sym_changerow标识为不可更改，sym_changecol标识为可更改时
该情形与②中所述类似；
在赋值结束后，在拓扑表中将岛屿row与岛屿col的缓冲区相切的信息进行记录，之后转到第四步；
第四步：判断最终的岛屿缓冲区宽度是否设置完毕。遍历岛屿缓冲区标识SYM，判断是否其中的sym_change是否均为不可更改，若是,则设置完毕，进入第五步；否则转到第一步；
第五步：构建缓冲区生长模型中各岛屿的缓冲区，利用前述步骤生成的SYM中的各岛屿缓冲区宽度，求取岛屿的缓冲区，输出结果，模型构建完成。
所述步骤g中，提取相邻缓冲区数目的方法为：提取缓冲区距离矩阵中该岛屿编号所在的行S(col,)以及缓冲区标识矩阵的第二行SYM(2,)，做矩阵运算S(col,)-SYM(2,)，统计该结果矩阵中0的数目即为所求相邻缓冲区数目。
本发明的基于缓冲区动态生长的岛屿自动选取方法，克服了以往岛屿自动选取方法采用Voronoi图表达岛屿密集程度时的缺陷问题，使用本方法能通过计算机作业自动选取出海图上的岛屿要素，在满足海图编绘规范要求的同时，可明显提高岛屿自动选取质量。
附图说明
图1是本发明基于缓冲区动态生长的岛屿自动选取方法主流程框图。
图2是本发明缓冲区生长模型的构建流程图。
图3是本发明缓冲区生长模型构造过程中判断岛屿缓冲区是否相交的流程图。
图4是本发明提取相邻缓冲区数目的流程图。
具体实施方式
本发明的实现过程是采用计算机实现海图制图过程中岛屿的自动选取。假设海图制图综合中有某区域A(海图待综合区域)需进行岛屿选取，利用综合区域之内资料图岛屿坐标数据，采用本发明方法进行岛屿自动选取，包括以下步骤，如图1所示：
a、读取海图上待选取区域的岛屿平面坐标数据；
b、利用读取的岛屿平面坐标数据，计算各岛屿之间的最小距离，并构建岛屿距离矩阵，构建的岛屿距离矩阵概念为：求取各岛屿与其他所有岛屿之间的最小距离，并将其放在距离矩阵S中。
Figure BDA0002306073570000051
以元素为例，sij表示岛屿i与岛屿j之间的最小距离。在这里为了之后的计算方便，将距离矩阵S对角线上的元素值设置为DBL_MAX(浮点型数值的最大值)；
c、为各个待选取岛屿设置缓冲区标识以及生长因子，并依据距离矩阵生成缓冲区宽度，各个待选取岛屿缓冲区标识以及生长因子的概念为：
(1)岛屿缓冲区标识。岛屿缓冲区标识是对岛屿的缓冲区宽度以及宽度是否可以更改进行标识的式子，其表达形式为：
Figure BDA0002306073570000052
其中sym_changei用于标识第i个岛屿的缓冲区宽度是否可更改，其中di表示第i个岛屿的缓冲区宽度。在初始状态时，将宽度是否可更改标识全部设置为可更改，将全部岛屿的缓冲区宽度设置为0，sym_elementi表示第i个岛屿的属性，例如对面状岛屿、点状岛屿分别设置不同的值。
(2)岛屿缓冲区宽度矩阵。为了之后方便生成缓冲区，设置两个宽度矩阵，分别为行主导矩阵，列主导矩阵：
Figure BDA0002306073570000061
Figure BDA0002306073570000062
其中行主导矩阵Drow的任意一行的所有元素相同，以第i行为例，这一行的值均为di，di为SYM中第i个岛屿的缓冲区宽度。同理，对于列主导矩阵而言，其任意一列的所有元素均相同，其值为第i个岛屿的缓冲区宽度。
(3)缓冲区宽度增长因子。由于在之前的设计中，需要使岛屿的影响域宽度不断地增大，因此在这里设置一个增长因子w，使w以一定的速率逐渐增大，w的初始值设置为0。
基于岛屿距离矩阵生成缓冲区的方法为：
第一步：为岛屿缓冲区标识SYM、岛屿缓冲区宽度矩阵Drow，Dcol赋值。其过程为，使增长因子w增大，获取全部SYM中sym_change设置为可更改的岛屿，将这些岛屿的缓冲区宽度d设置为w；将宽度矩阵中的行主导矩阵Drow每一行的值都设置为SYM中的对应岛屿缓冲区宽度值，例如，取行主导矩阵Drow的第i行的全部元素，记为Drow(i,),令Drow(i,)＝di(其中di为SYM中第i个岛屿的缓冲区宽度)，即将第i行的全部元素大小均设置为di；对于列主导矩阵Dcol赋值过程与行主导矩阵类似。
第二步：判断当前缓冲区宽度下，是否有岛屿的缓冲区相交，如图3所示，利用第一步设置好的Drow，Dcol，进行矩阵运算S-Drow-Dcol，其结果矩阵记为Sres(在此结果矩阵中，任意一元素Sres(i,j)表达了在当前SYM中所述的岛屿缓冲区宽度下，岛屿i与岛屿j缓冲区之间的最小距离)，遍历结果矩阵Sres，获取其中最小的元素Sres(row,col)，判断Sres(row,col)是否大于0，若是，则表明在当前缓冲区宽度下，各岛屿的缓冲区没有相切或相交，返回第一步；若小于0，则表明在此缓冲区宽度下，岛屿row与岛屿col的缓冲区相交，此时将row与col进行记录，进入第三步；若等于0，此时无法判断是否是已经完成宽度求取的两岛屿缓冲区，因此在这里转入第一步。
第三步：为出现相交的岛屿缓冲区设置最终的宽度值，并对增长因子重新赋值。利用第二步获取的row与col，获取这两个岛屿的相关信息，即第row个岛屿与第col个岛屿之间的最小距离srowcol，岛屿缓冲区标识SYM中的sym_changerow，sym_changecol，sym_elementrow，sym_elementcol；针对不同的sym_changerow，sym_changecol，分以下3中情况对增长因子w，SYM中的岛屿缓冲区宽度d进行如下赋值：
①当sym_changerow，sym_changecol均标识为可更改时
增长因子：
Figure BDA0002306073570000071
岛屿缓冲区宽度：
Figure BDA0002306073570000072
并将sym_changerow，sym_changecol设置为不可更改。
②当sym_changerow标识为可更改，sym_changecol标识为不可更改时增长因子：
w＝srowcol-dcol (8)
岛屿缓冲区宽度：
drow＝w (9)
将sym_changerow设置为不可更改。
③当sym_changerow标识为不可更改，sym_changecol标识为可更改时该情形与②中所述类似，在这里不再赘述。
在赋值结束后，在拓扑表中将岛屿row与岛屿col的缓冲区相切的信息进行记录，之后转到第四步。
第四步：判断最终的岛屿缓冲区宽度是否设置完毕。遍历岛屿缓冲区标识SYM，判断是否其中的sym_change是否均为不可更改，若是,则设置完毕，进入第五步；否则转到第一步。第五步：构建缓冲区生长模型中各岛屿的缓冲区。利用前述步骤生成的SYM中的各岛屿缓冲区宽度，求取岛屿的缓冲区，输出结果，模型构建完成，如图2所示；
d、提取各缓冲区宽度，为各个岛屿构建缓冲区，计算个岛屿缓冲区的面积，如图4所示；
e、根据上述步骤所计算的缓冲区面积，将待选取岛屿按照构建的缓冲区面积由大到小排序；
f、判断是否存在面积大小相等的缓冲区，若是则进入步骤g，否则进入步骤h；
g、提取其中缓冲区面积大小相等的岛屿，获取其相邻缓冲区数目，并判断缓冲区数目是否相等，若相等，则计算岛屿与其周围岛屿间平均间距，将其平均间距较大者置于前部，否则将其中相邻缓冲区数目较多的岛屿放于前侧，随后转步骤h，提取相邻缓冲区数目的方法为：遍历缓冲区距离矩阵中该岛屿编号所在的行，以其中任意两岛屿为例，判断两岛屿之间的距离是否等于步骤(c)中生成的岛屿缓冲区长度之和，若是则认为该两岛屿缓冲区相邻；统计该距离矩阵中与该岛屿缓冲区相邻的岛屿数目；
h、删除其中第一个岛屿；
i、判断剩余岛屿是否达到选取指标，若是，则进入步骤j，否则，将剩余岛屿设置为待选取岛屿，并进入步骤c；
j、提取选取的岛屿，更新海图数据，输出选取结果。
具体实施例：
步骤a，在计算机上打开编绘海图，输入需要进行岛屿选取区域A的边界点平面坐标(XAi,YAi)(i＝1,2,Λ,nA)，提取该区域中的全体待选取岛屿集合A(Island)；
步骤b，采用多边形之间最小距离计算方法计算上述岛屿集合中各岛屿之间的最小距离，并将其放入距离矩阵
Figure BDA0002306073570000081
中，同时构建初始缓冲区标识矩阵
Figure BDA0002306073570000082
并将其中的宽度是否可更改标识全部设置为可更改，将全部岛屿的缓冲区宽度设置为0，依据岛屿属性为sym_elementi表设置初值；
步骤c，令生长因子w＝w+0.01，依据此将SYM中sym_change为可更改的列中的d设置为w，随后采用图2所述方法生成岛屿缓冲区生长模型；
步骤d，提取已生成岛屿缓冲区生长模型的岛屿缓冲区标识矩阵中的第二行，该行即为所求各岛屿的缓冲区宽度集合D(di,i＝0,1,L,n)，并依据此为各岛屿设置缓冲区，可得岛屿缓冲区集合Buffers(Island)；
步骤e，计算上述构建的岛屿缓冲区集合Buffers(Island)中各缓冲区的面积集合Area(Buffersi)(i＝0,1,L,n)。
步骤f，将上述各岛屿按照缓冲区面积集合Area(Buffersi)中的面积大小进行由小到大进行排序，得岛屿序列T；
步骤g，遍历面积集合Area(Buffersi)，判断是否存在面积大小相等的缓冲区，若存在，则采用图3所示方法求取这些缓冲区面积大小相等的岛屿的相邻缓冲区数目，将其中数目较多的岛屿放于前部，更新岛屿序列T，若相邻缓冲区数目也相同，则计算该岛屿与周围一定范围内的岛屿间平均间距，将其中间距较大者放于前侧，并更新岛屿序列T。反之，若不存在面积大小相同的面积大小相同的缓冲区，则进入步骤h；
步骤h，删除经步骤g更新后的岛屿序列T中的第一个岛屿；
步骤i，判断剩余的岛屿数量是否达到指标，若达到指标，则将剩余的岛屿进行输出，并更新海图上选取区域A中的岛屿数据，反之，将剩余的岛屿集合A′(Island)设置为待选取岛屿，重复步骤b到步骤i，直到剩余岛屿的数目达到指标。
本发明是通过实施例进行描述的，本领域技术人员知悉，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本申请的权利要求范围内的实施例都属于本发明的保护范围。

Claims (5)

1.一种基于缓冲区生长模型的岛屿自动选取方法，其特征在于：包括以下步骤：
a、读取海图上待选取区域的岛屿平面坐标数据；
b、利用读取的岛屿平面坐标数据，计算各岛屿之间的最小距离，并构建岛屿距离矩阵；
c、为各个待选取岛屿设置缓冲区标识以及增长因子，并依据距离矩阵生成缓冲区宽度；
d、提取各缓冲区宽度，为各个岛屿构建缓冲区，计算个岛屿缓冲区的面积；
e、根据上述步骤所计算的缓冲区面积，将待选取岛屿按照构建的缓冲区面积由大到小排序；
f、判断是否存在面积大小相等的缓冲区，若是则进入步骤g，否则进入步骤h；
g、提取其中缓冲区面积大小相等的岛屿，获取其相邻缓冲区数目，并判断缓冲区数目是否相等，若相等，则计算岛屿与其周围岛屿间平均间距，将其平均间距较大者置于前部，否则将其中相邻缓冲区数目较多的岛屿放于前侧，随后转步骤h；
h、删除其中第一个岛屿；
i、判断剩余岛屿是否达到选取指标，若是，则进入步骤j，否则，将剩余岛屿设置为待选取岛屿，并进入步骤c；
j、提取选取的岛屿，更新海图数据，输出选取结果。
2.根据权利要求1所述的一种基于缓冲区生长模型的岛屿自动选取方法，其特征在于：所述步骤b中，构建的岛屿距离矩阵概念为：求取各岛屿与其他所有岛屿之间的最小距离，并将其放在距离矩阵S中，
Figure FDA0002306073560000011
以元素为例，sij表示岛屿i与岛屿j之间的最小距离。
3.根据权利要求1所述的一种基于缓冲区生长模型的岛屿自动选取方法，其特征在于：所述步骤c中，各个待选取岛屿缓冲区标识以及生长因子为：
(1)岛屿缓冲区标识
岛屿缓冲区标识是对岛屿的缓冲区宽度以及宽度是否可以更改进行标识的式子，其表达形式为：
Figure FDA0002306073560000012
其中sym_changei用于标识第i个岛屿的缓冲区宽度是否可更改，其中di表示第i个岛屿的缓冲区宽度，在初始状态时，将宽度是否可更改标识全部设置为可更改，将全部岛屿的缓冲区宽度设置为0，sym_elementi表示第i个岛屿的属性；
(2)岛屿缓冲区宽度矩阵
为了之后方便生成缓冲区，设置两个宽度矩阵，分别为行主导矩阵，列主导矩阵：
Figure FDA0002306073560000021
Figure FDA0002306073560000022
其中行主导矩阵Drow的任意一行的所有元素相同，以第i行为例，这一行的值均为di，di为SYM中第i个岛屿的缓冲区宽度，同理，对于列主导矩阵而言，其任意一列的所有元素均相同，其值为第i个岛屿的缓冲区宽度；
(3)缓冲区宽度增长因子
由于需要使岛屿的影响域宽度不断地增大，因此设置一个增长因子w，使w以一定的速率逐渐增大，w的初始值设置为0。
4.根据权利要求1所述的一种基于缓冲区生长模型的岛屿自动选取方法，其特征在于：所述步骤c中，基于岛屿距离矩阵生成缓冲区的方法为：
第一步：为岛屿缓冲区标识SYM、岛屿缓冲区宽度矩阵Drow，Dcol赋值， 其过程为：使增长因子w增大，获取全部SYM中sym_change设置为可更改的岛屿，将这些岛屿的缓冲区宽度d设置为w；将宽度矩阵中的行主导矩阵Drow每一行的值都设置为SYM中的对应岛屿缓冲区宽度值，例如，取行主导矩阵Drow的第i行的全部元素，记为Drow(i,),令Drow(i,)＝di，其中di为SYM中第i个岛屿的缓冲区宽度，即将第i行的全部元素大小均设置为di；对于列主导矩阵Dcol赋值过程与行主导矩阵类似；
第二步：判断当前缓冲区宽度下，是否有岛屿的缓冲区相交，利用第一步设置好的Drow，Dcol，进行矩阵运算S-Drow-Dcol，其结果矩阵记为Sres，在此结果矩阵中，任意一元素Sres(i,j)表达了在当前SYM中所述的岛屿缓冲区宽度下，岛屿i与岛屿j缓冲区之间的最小距离，遍历结果矩阵Sres，获取其中最小的元素Sres(row,col)，判断Sres(row,col)是否大于0，若是，则表明在当前缓冲区宽度下，各岛屿的缓冲区没有相切或相交，返回第一步；若小于0，则表明在此缓冲区宽度下，岛屿row与岛屿col的缓冲区相交，此时将row与col进行记录，进入第三步；若等于0，此时无法判断是否是已经完成宽度求取的两岛屿缓冲区，因此在这里转入第一步；
第三步：为出现相交的岛屿缓冲区设置最终的宽度值，并对增长因子重新赋值，利用第二步获取的row与col，获取这两个岛屿的相关信息，即第row个岛屿与第col个岛屿之间的最小距离srowcol，岛屿缓冲区标识SYM中的sym_changerow，sym_changecol，sym_elementrow，sym_elementcol；针对不同的sym_changerow，sym_changecol，分以下3中情况对增长因子w，SYM中的岛屿缓冲区宽度d进行如下赋值：
①当sym_changerow，sym_changecol均标识为可更改时
增长因子：
Figure FDA0002306073560000031
岛屿缓冲区宽度：
Figure FDA0002306073560000032
并将sym_changerow，sym_changecol设置为不可更改；
②当sym_changerow标识为可更改，sym_changecol标识为不可更改时
增长因子：
w＝srowcol-dcol (8)
岛屿缓冲区宽度：
drow＝w (9)
将sym_changerow设置为不可更改；
③当sym_changerow标识为不可更改，sym_changecol标识为可更改时
该情形与②中所述类似；
在赋值结束后，在拓扑表中将岛屿row与岛屿col的缓冲区相切的信息进行记录，之后转到第四步；
第四步：判断最终的岛屿缓冲区宽度是否设置完毕， 遍历岛屿缓冲区标识SYM，判断是否其中的sym_change是否均为不可更改，若是,则设置完毕，进入第五步；否则转到第一步；
第五步：构建缓冲区生长模型中各岛屿的缓冲区，利用前述步骤生成的SYM中的各岛屿缓冲区宽度，求取岛屿的缓冲区，输出结果，模型构建完成。
5.根据权利要求1所述的基于缓冲区生长模型的岛屿自动选取方法，其特征在于：所述步骤g中，提取相邻缓冲区数目的方法为：提取缓冲区距离矩阵中该岛屿编号所在的行S(col,)以及缓冲区标识矩阵的第二行SYM(2,)，做矩阵运算S(col,)-SYM(2,)，统计该结果矩阵中0的数目即为所求相邻缓冲区数目。
CN201911240476.1A 2019-12-06 2019-12-06 一种基于缓冲区生长模型的岛屿自动选取方法 Active CN111008598B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911240476.1A CN111008598B (zh) 2019-12-06 2019-12-06 一种基于缓冲区生长模型的岛屿自动选取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911240476.1A CN111008598B (zh) 2019-12-06 2019-12-06 一种基于缓冲区生长模型的岛屿自动选取方法

Publications (2)

Publication Number Publication Date
CN111008598A CN111008598A (zh) 2020-04-14
CN111008598B true CN111008598B (zh) 2022-10-28

Family

ID=70114928

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911240476.1A Active CN111008598B (zh) 2019-12-06 2019-12-06 一种基于缓冲区生长模型的岛屿自动选取方法

Country Status (1)

Country Link
CN (1) CN111008598B (zh)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107255801B (zh) * 2017-06-13 2019-08-02 中国电子科技集团公司第二十八研究所 一种岸基对海监视雷达信号的处理方法
CN108428220B (zh) * 2018-03-05 2020-12-01 武汉大学 静止轨道卫星序列遥感影像海岛礁区域自动几何校正方法
CN108959348B (zh) * 2018-04-20 2021-07-09 中国人民解放军海军大连舰艇学院 一种融合多幅海图数据的碍航区自动生成方法

Also Published As

Publication number Publication date
CN111008598A (zh) 2020-04-14

Similar Documents

Publication Publication Date Title
WO2020233152A1 (zh) 基于城市建筑空间数据的建成区边界识别方法及设备
CN102289991B (zh) 一种基于视觉变量的地图注记自动分类配置方法
CN112070769B (zh) 一种基于dbscan的分层点云分割方法
CN101901503B (zh) 使用光谱分析的伸展驱动的网格参数化
JP3796658B2 (ja) 解析モデル作成装置
CN108537274A (zh) 一种基于网格的空间多尺度快速聚类方法
CN110866922B (zh) 基于强化学习和迁移学习的图像语义分割模型及建模方法
CN109102538B (zh) 利用等值线数据提取槽脊特征点与槽脊特征线的方法、天气图分析方法
CN115757604B (zh) 一种基于夜光影像数据的gdp时空演变分析方法
CN102165487A (zh) 前景区提取程序、前景区提取装置以及前景区提取方法
CN111858810A (zh) 一种面向道路dem构建的建模高程点筛选方法
CN110992481B (zh) 基于最邻近连线的建筑物白模一致性合并方法
KR20140024590A (ko) 연속지적도의 다축척 모델 생성 방법
CN103793549A (zh) 采用模糊聚类与随机游走的计算机辅助乱针绣制作方法
CN111008598B (zh) 一种基于缓冲区生长模型的岛屿自动选取方法
CN112346409B (zh) 一种快速判断点云包含问题的方法
KR20160059452A (ko) 이미지의 파티셔닝
CN117413988A (zh) 一种非常规的大翻领服装制版配领方法
CN109190255B (zh) 一种面向城市三维产权空间立体重构方法
CN115187127B (zh) 基于空间分析的详细规划分级管理智能检测方法
CN114758104B (zh) 一种适用于三维设计钢筋图的抗干涉散布筋自动标注方法
CN110110685A (zh) 数据建档方法、装置及服务器
CN106599044A (zh) 一种道路网目标信息的识别和处理方法
CN113887737B (zh) 一种基于机器学习的样本集自动生成方法
CN115310237B (zh) 一种基于图形文件的型钢孔型自动参数化方法

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