CN116664455A - 一种处理二维矢量图形之间的边界不契合方法 - Google Patents

一种处理二维矢量图形之间的边界不契合方法 Download PDF

Info

Publication number
CN116664455A
CN116664455A CN202310956478.0A CN202310956478A CN116664455A CN 116664455 A CN116664455 A CN 116664455A CN 202310956478 A CN202310956478 A CN 202310956478A CN 116664455 A CN116664455 A CN 116664455A
Authority
CN
China
Prior art keywords
point
folding
building property
dimensional vector
points
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
Application number
CN202310956478.0A
Other languages
English (en)
Other versions
CN116664455B (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.)
Wuhan Tianheng Information Technology Co ltd
Original Assignee
Wuhan Tianheng Information Technology Co ltd
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 Wuhan Tianheng Information Technology Co ltd filed Critical Wuhan Tianheng Information Technology Co ltd
Priority to CN202310956478.0A priority Critical patent/CN116664455B/zh
Publication of CN116664455A publication Critical patent/CN116664455A/zh
Application granted granted Critical
Publication of CN116664455B publication Critical patent/CN116664455B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/77—Retouching; Inpainting; Scratch removal
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/60—Editing figures and text; Combining figures or text
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/80—Geometric correction
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Processing Or Creating Images (AREA)
  • Image Generation (AREA)

Abstract

本发明涉及一种处理二维矢量图形之间的边界不契合方法，包括以下步骤：对分层分户图进行叠加得到建筑产权体模型底图；对建筑产权体模型底图进行属性挂接；对建筑产权体模型底图采用快速排斥法进行筛选，得到二维矢量数据；通过获取单幢建筑产权体模型底图中的所有二维矢量数据的折点数组来获取构成单个二维矢量数据的线型和角度来对筛选后的建筑产权体模型底图二维矢量数据进行数据分类；判断当前单个建筑产权体模型底图二维矢量数据与其它单个建筑产权体模型底图二维矢量数据涉及的边界未契合问题的类型；根据不同类别边界的折点数组和边界未契合问题的类型，完成边界契合修复。本发明有益效果是：大大减少人工量和处理时间。

Description

一种处理二维矢量图形之间的边界不契合方法
技术领域
本发明涉及二维图像边界处理领域，尤其涉及一种处理二维矢量图形之间的边界不契合方法。
背景技术
在三维地籍管理系统的建设中数据是重中之重，其中产权体模型数据生产的基本流程为：
生成房产分层分户图；将《房产测绘报告》（CAD文件或PDF文件）中的每一层每一户的图形数据转换为二维矢量面图形数据；叠加分层分户图；即将同一幢的分层分户图叠加在一起生成一幢的分层分户图；空间数据校正：参照地籍图或地形图确定叠加后的分层分户图的空间位置和建筑朝向；属性挂接；将每一户对应的标识码，地籍区编号，地籍区名称，不动产单元号，基底高程，层高，备注，户号，所在层，坐落，幢坐落等属性信息挂载在每一户上；5.产权体模型生产；通过叠加后的分层分户图以及挂接属性中的基底高程，层高属性转换为三维白模。
目前数据生产中生成房产分层分户图的基础数据通常为《房产测绘报告》。
《房产测绘报告》有两个难以避免的问题：一是《房产测绘报告》多为PDF格式需要进行矢量化后由人工进行二次处理后才能继续使用；二是《房产测绘报告》是人工测绘数据，会有一些细微的误差。在叠加分层分户图和空间数据校正的步骤中由于同一幢的分层分户图有可能来自于多个《房产测绘报告》，在叠加合并后就会因为人工测量的一些小误差造成不同层的分层分户图矢量二维图形的边界之间出现边界未契合的情况，这些边界未契合的错误往往误差在0.05米之内，但却非常多需要人工进行一些修复，而在三维地籍的研发和运用上，数据生产范围动则一个区甚至一个市，因为细微数据误差问题造成的边界未契合问题数量会急剧增加，处理这些问题会占用大量的人力和时间成本。
发明内容
为解决目前二维矢量数据存在因为数据源细微数据误差问题造成的边界未契合问题，本发明提出了一种处理二维矢量图形之间的边界不契合方法，该方法基于同一幢的分层分户图二维矢量数据的每一个折点与其他分层分户图二维矢量数据折点及折线的距离为基本依据来判断哪些折点需要处理并以此来批量修复二维矢量数据边界未契合的问题。
本发明所述的一种处理二维矢量图形之间的边界不契合方法，包括以下步骤：
S1、对分层分户图进行叠加得到建筑产权体模型底图；
S2、对建筑产权体模型底图进行属性挂接；
S3、对建筑产权体模型底图采用快速排斥法进行筛选，得到筛选后的单幢建筑产权体模型底图二维矢量数据；
S4、通过获取单幢建筑产权体模型底图中的所有二维矢量数据的折点数组来获取构成单个二维矢量数据的线型和角度来对筛选后的建筑产权体模型底图二维矢量数据进行数据分类；
S5、判断当前单个建筑产权体模型底图二维矢量数据与其它单个建筑产权体模型底图二维矢量数据涉及的边界未契合问题的类型；
S6、根据不同类别边界的折点数组和边界未契合问题的类型，完成边界契合修复。
本发明提供的有益效果是：可以最大限度的保证原图形的完整性下批量处理因为人工测量的一些小误差造成不同层的分层分户图矢量二维图形的边界之间出现边界未契合的情况，可以大大减少人工量和处理时间。
附图说明
图1是本发明方法流程示意图；
图2是产权体模型底图折线类型示意图；
图3是产权体模型底图边界未契合现象示意图；
图4是产权体模型底图边界未契合类型示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。
请参考图1，图1是本发明方法流程示意图；本发明提供的一种处理二维矢量图形之间的边界不契合方法，具体包括以下步骤：
S1、对分层分户图进行叠加得到建筑产权体模型底图；
S2、对建筑产权体模型底图进行属性挂接；
S3、对建筑产权体模型底图采用快速排斥法进行筛选，得到筛选后的单幢建筑产权体模型底图二维矢量数据；
需要说明的是，步骤S3中，采用快速排斥法筛选的过程如下：
S31、通过已有的楼栋幢坐落属性，筛选需要处理的建筑产权体模型底图；
需要说明的是，幢是指房屋的计量单位，通常是指一座独立的、包括不同结构和不同层次的房屋。
多个不同功能、不同结构、不同层次的房屋相互组合与交织组成一个建筑综合体(如同一裙房上矗立多个主楼或同一地下室上矗立多个单体建筑等) 。
在三维地籍的研发和运用上，数据生产范围动则一个区甚至一个市，数据非常之繁杂，一次性处理这些数据效率很低。
因此本发明在处理数据的时候以幢为条件设置了一个范围，幢是房产分层分户图之间出现边界未契合问题最大的一个范围，因此把数据按幢坐落去分区处理是比较准确的处理范围，同时也可以减少一次性处理数据的数量，提高处理效率。
在一些其它实施例中，可能会没有幢属性来区分，那么可以根据实际情况跳过此步骤。
S32、设置步长阈值r，所述步长阈值用以约束多个建筑产权体模型底图的折点或边界的距离，当距离超过所述步长阈值时，进入快速排除法；
作为一种实施例而言，在叠加分层分户图和空间数据校正的步骤中由于同一幢的分层分户图有可能来自于多个《房产测绘报告》，在叠加合并后就会因为人工测量的一些小误差造成不同层的分层分户图矢量二维图形的边界之间出现边界未契合的情况，这些边界未契合的错误往往误差在0.05米之内。
而《房产测绘报告》中是不包含建筑外立面细节部位以及柱子之类的建筑内部细节的，它所记录的边界折线长度除了弧形结构不会小于0.3米。因此为了批量处理分层分户图矢量二维图形之间的边界未契合错误，本发明需要设置一个阈值来保证包含图形可能出现的边界未契合误差的最大范围（>=0.05米），同时也需要保证不会影响到建筑原本结构(<=0.3米)。需要说明的是，弧形情况会在后续步骤中分类处理。
本发明在处理数据的流程中设置阈值r为0.1米来确认数据的处理范围。当然，在一些其它实施例中，此值根据具体数据情况可以进行细微修改。
S33、获取使用幢坐落字段筛选的建筑产权体模型底图二维矢量数据的矩形包围框，将矩形包围框扩大一圈，扩大的值为所述步长r/2；
需要说明的是，在获得了以幢为条件筛选出来的数据后，需要遍历这些分层分户图数据（也可以称作建筑产权体模型底图，方便阐述，后文部分处以分层分户图进行表述），将这些分层分户图数据和可能出现边界未契合现象的分层分户图数据进行匹配筛选。
可能会出现边界未契合现象的分层分户图可以理解为两个分户图的边界（折线，折点）之间的距离小于设置的阈值，简单来说就是二维矢量图形的边界扩展阈值的一半后两个二维矢量图形边界存在相切或相交，这两个分层分户图数据可能会存在边界未契合现象。
在实际情况中分成分户图的数据形状会非常的多，直接去计算两个二维矢量图形边界存在相切或交集非常麻烦，因此本发明使用图形的矩形包围框配合快速排斥试验快速的筛选匹配数据。
矩形包围框点集合获取方法：1.收集需要获取矩形包围框点集合的分成分户图的所有坐标，排序筛选出所有坐标中的X值和Y值的最大值（X max,Y max），X值Y值的最小值（X min,Y min）。2.包围框的四个点的坐标为（X min,Y max），（X max,Y max），（X max,Y min），（X min,Y min）。
如果直接使用矩形包围框进行快速排斥试验那么计算的是两个图形本身是否相切或相交，为判断两个分层分户图数据是否会存在边界未契合现象，就需要加上阈值，由于是两个图形的边界之间的距离为阈值，所以单个图形的矩形包围框的四个点坐标需要加上阈值的一半，矩形包围框的四个点坐标更新为（X min-r/2,Y max+r/2），（X max+r/2,Y max+r/2），（X max+r/2,Y min），（X min-r/2,Y min-r/2）；
S34、遍历使用幢坐落数据筛选的建筑产权体模型底图，使用当前建筑产权体模型底图的包围框和其它建筑产权体模型底图的包围框进行一次快速排斥试验，初步筛选出与当前建筑产权体模型底图边界可能需要契合的其他建筑产权体模型底图，作为筛选后的建筑产权体模型底图。
步骤S34具体如下：
S341、获取当前建筑产权体模型底图的包围框和其它建筑产权体模型底图的包围框的坐标集合；
所述当前建筑产权体模型底图的包围框坐标为：（X 1min-r/2,Y 1max+r/2），（X 1max+r/2,Y 1max+r/2），（X 1max+r/2,Y 1min），（X 1min-r/2,Y 1min-r/2）；
所述其它建筑产权体模型底图的包围框的坐标为:（X 2min-r/2,Y 2max+r/2），（X 2max+r/2,Y 2max+r/2），（X 2max+r/2,Y 2min），（X 2min-r/2,Y 2min-r/2）；
S342、根据坐标集合判断两个包围框是否相切或相交，若是则判定其它建筑产权体模型底图的边界需要与当前建筑产权体模型底图的边界相契合；否则表示不需要契合。
作为一种实施例而言，当获取了两个图形加了一半阈值的坐标集合后，就可以使用快速排斥试验判断两个图形分别加上阈值距离一半以后是否有相切或者相交，从而得出两个图形的边界是否可能出现未契合现象。
快速排斥试验步骤：
1.首先获取两个图形矩形包围框的坐标集合，图形1（X 1min-r/2,Y 1max+r/2），（X 1max+r/2,Y 1max+r/2），（X 1max+r/2,Y 1min），（X 1min-r/2,Y 1min-r/2），图形2（X 2min-r/2,Y 2max+r/2），（X 2max+r/2,Y 2max+r/2），（X 2max+r/2,Y 2min），（X 2min-r/2,Y 2min-r/2）；
2.判断两个矩形包围框是否相切或相交需要同时符合这四个条件，条件一：X 1max+r/2与X 1min-r/2中小的那个值小于等于X 2max+r/2与X 2min-r/2中大的那个值；
条件二：X 2max+r/2与X 2min-r/2中小的那个值小于等于X 1max+r/2与X 1min-r/2中大的那个值；
条件三：Y 1max+r/2与Y 1min-r/2中小的那个值小于等于Y 2max+r/2与Y 2min-r/2中大的那个值；
条技四：Y 2max+r/2与Y 2min-r/2中小的那个值小于等于Y 1max+r/2与Y 1min-r/2大的那个值。
本发明使用快速排斥试验筛选匹配分层分户图和可能与分层分户图有边界未契合图形的集合用来继续下一阶段的处理。
S4、通过获取单幢建筑产权体模型底图中的所有二维矢量数据的折点数组来获取构成单个二维矢量数据的线型和角度来对筛选后的建筑产权体模型底图二维矢量数据进行数据分类；
需要说明的是，在三维地籍的研发和运用上，分层分户图的图形并都是不是一个个规则的矩形或规整的多边形，通常是由多种线形（弧线，直线），多种角度（直角，钝角，锐角）构成的复杂图形，因此本申请中不可能直接用一种方法处理所有类型的图形，故而需要对图形进行分类，执行不同的策略。
由于常规二维矢量图形数据格式比如SHP是不支持弧线这种类型的，通常弧线类型是由大量钝角角度的直线来描述的。因此本发明使用角度来进行图形内的折点折线分类执行不同的修复策略，为了方便本发明会把角度和折点挂钩。
具体参考如图2所示的图形一；图形一中，P代表的是点，L代表的是线，折线L 1是由点P 1和点P 2构成的，折线L 2是由点P 2和点P 3构成的，点P2的角度就是点P 2涉及到的折线L 1与L 2之间的角度，而折线L 1 L 2涉及到的折点为P 1 P 2 P 3，因此计算折点P 2的角度只需要获取这个折点和这个折点涉及到的折线的两个折点P 1 P 3的坐标就可以求出这个折点的角度。
步骤S4的具体过程如下：
步骤S4中的数据分类过程具体如下：
S41、收集需要契合的建筑产权体模型底图的所有折点的坐标集合(P 1，P 2，P 3，……)；需要说明的是，点的坐标集合的顺序是按照构成这个图形的折点的顺序；
S42、获取折点P i涉及到的折线L i和L i+1；其中i为折点编号；
S43、折线L i由折点P i和P i+1构成；折线L i+1由折点P i+1和P i+2构成；
S44、获取折线L i和L i+1的向量v i和v i+1；
作为一种实施例，比如点P 2，折点P 2涉及到的折线为L 1和L 2，我们可以通过构成折线L 1的折点P 1 P 2，构成折线L 2的折点P 2 P 3来得出折线L 1 L 2的向量；
设向量L 1为（b 1，b 2），向量L 2为(d 1,d 2)，则：
b 1=P 1坐标的X值减去P 2坐标的X值；
b 2=P 1坐标的Y值减去P 2坐标的Y值；
d 1=P 2坐标的X值减去P 3坐标的X值；
d 2=P 2坐标的Y值减去P 3坐标的Y值；
S45、计算点P i+1涉及到的向量v i和v i+1之间的夹角α；
仍以点P2为例，则P2的α=arccos((b 1*d 1+b 2*d 2) / (sqrt(b 1^2 +b 2^2) *sqrt(d 1^2 +d 2^2)))；这里arccos表示反余弦函数，sqrt表示平方根；^表示幂；
S46、当α处于大于第一预设角度或小于第二预设角度且和这个连续相关的两个折点夹角角度不是全部都大于等于第三预设角度时，认为点P i+1涉及到的夹角α属于第一类夹角；
当α小于等于第一预设角度或大于等于第二预设角度，则认为点P i+1涉及到的夹角α属于第二类夹角；
当α大于等于第三预设角度且连续相关的两个折点夹角角度全都大于第三预设角度或前两个折点全都大于第三预设角度或后两个折点全都大于第三预设角度，则认为夹角α属于点P i+1涉及到的第三类夹角；
其中，第二预设角度小于第一预设角度小于第三预设角度。
作为一种实施例，通过步骤S45求出的角度是弧度制的，需要以度数表示，将其结果乘以180/π进行转换；
由于两个线段夹角可能会大于180度，为了方便后续计算，需要把超过180度的夹角转换为小于180度的夹角，方法就是加一个判断条件，如果计算出来的角度大于180度则此角度改为360度减去此角度。
当获得了分层分户图后需要按照点集合对图形所有夹角按照角度区间进行分类。
本发明将这些夹角按照角度为三类去进行修复：
1.第一类夹角（普通夹角）：即折点夹角角度符合大于91或小于89度且和这个连续相关的两个折点夹角角度不是全部都大于等于110度。比如说折点P2，P2的夹角小于89度和折点P2连续相关的两个折点P1和P3的夹角都不大于等于110度，所以它是普通夹角。
2.第二类夹角（直角夹角）：即折点夹角角度小于等于91或大于等于89度，比如折点P 3。
3.弧线夹角：即折点角度大于等于110度且连续相关的两个折点夹角角度全都大于110度或前两个折点全都大于110度或后两个折点全都大于110度。比如折点P 8，它的夹角角度大于110度，它的后面两个折点P 9 P 10的角度也全都大于110度，所以它是弧线夹角。
S5、判断当前单个建筑产权体模型底图二维矢量数据与其它单个建筑产权体模型底图二维矢量数据涉及的边界未契合问题的类型；
步骤S5具体为：当建筑产权体模型底图的折点到另一个建筑产权体模型底图的折线或者折点的距离如果小于设置好的阈值且该折点到另一建筑产权体模型底图的折线或者折点的距离大于零，则该折点涉及到的边界存在边界未契合问题；
所述边界未契合问题分为：折点到折线的未契合问题和折点到折点的未契合。
需要说明的是，本发明将判断分层分户图的哪些边界与其他分层分户图未契合的条件由分层分户图边界分解到了分层分户图的折点和折线，即分户图的折点到另一个分户图的折线或者折点的距离如果小于设置好的阈值且该折点到另一个分户图的折线或者折点的距离大于零则该折点涉及到的边界有边界未契合问题。
以此为核心算法配合其他方法来完善整个流程。
举例说明，请参考图3，图3中的建筑一层和建筑二层之间的边界契合错误，建筑一层的折点D 2很明显应该处于建筑折二层折点P 2上就属于是折点到折点的问题；
而建筑二层的折点P 0很明显应该在建筑一层折线LL 0上，可以理解为是折点到折线的问题。
需要说明的是，所述折点到折线的未契合问题的计算，具体如下:
获取折点P i(x 1,y 1)与折点P i+1(x 2,y 2)所构成的直线L；
获取线段LL j；线段LL j两个端点的坐标为(x 3,y 3) 和 (x 4,y 4)；
直线L的斜率设为m 1，m 1= (y 2-y 1) / (x 2-x 1)
线段LL j的斜率设为m 2，m2 = (y 4-y 3) / (x 4-x 3)；
求直线L与线段LL j的斜距，其中直线L的斜距为b 1，b 1=y 1-m 1*x 1；线段LL j的斜距为b 2，b 2=y 3-m 2*x 3；
计算直线L与线段LL j延长后的直线的交点，设交点的坐标为（x，y），x= (b 2-b 1) /(m 1-m 2)，y=m 1*x+b 1；
判断交点是否在线段LL j上，若是，则计算交点到折点P i的距离，若距离小于阈值，就将折点的坐标修改为交点的坐标；
其中，交点在线段LL j上需要满足下面条件：min(x 3,x 4)<=x<=max(x 3,x 4)并且min(y 3,y 4)<=y<=max(y 3,y 4)；其中min是指括号中的两个值当中最小的那个值，max是指括号中的两个值当中最大的那个值。
作为一种实施例而言，比如建筑二层折点P 0到建筑一层折线LL 0的边界未契合错误修复后应该在的位置为折点P 0 P 1构成的直线与线段LL 0的交点；
同理可知P 0 P 1构成的直线与线段LL 0的交点与P 0之间的距离如果小于阈值，则折点P0为需要修复的边界未契合问题的折点。
以下为直线与线段是否相交，如果相交求交点的具体算法：
首先可以将直线和线段当成两个直线去求他们的斜率，设直线1的两个折点的坐标为(x 1,y 1) 和 (x 2,y 2)；线段1的两个折点的坐标为(x 3,y 3) 和 (x 4,y 4)；直线1的斜率设为m 1，m 1= (y 2-y 1) / (x 2-x 1)；线段1的斜率设为m 2，m 2= (y 4-y 3) / (x 4-x 3)。如果直线或线段的斜率都为0或者相等则直线和线段平行，那么就不需要使用直线到线段检测直接使用点到点的检测即可；求斜距，设直线1的斜距为b 1，b 1=y 1-m 1*x 1；设线段1的斜距为b 2，b 2=y 3-m 2*x 3；
计算两个直线的交点，设交点的坐标为（x，y），x= (b 2-b 1) / (m 1-m 2)，y=m 1*x+b 1；
判断交点是否在线段1上，交点在线段1上需要满足下面条件：min(x 3,x 4)<=x<=max(x 3,x 4)并且min(y 3,y 4)<=y<=max(y 3,y 4)；其中min是指括号中的两个值当中最小的那个值，max是指括号中的两个值当中最大的那个值。
需要说明的是，所述折点到折点的未契合问题的计算，具体如下：遍历求取建筑产权体模型底图所有折点与另外一个建筑产权体模型底图的所有折点之间的距离，如果这个建筑产权体模型底图的一个折点与另外一个建筑产权体模型底图的折点之间的距离小于阈值则将小于阈值另外一个建筑产权体模型底图的折点坐标按照距离进行排序只保留距离最近的那个，并将这个折点的坐标更新为保留的那个折点。
S6、根据不同类别边界的折点数组和边界未契合问题的类型，完成边界契合修复。
步骤S6中，当折点为第一类夹角的修复过程如下：先进行折点到折点的未契合问题的计算，再进行折点到折点的未契合问题的计算。
步骤S6中，折点为第二类夹角的修复过程如下：先进行折点到折点的未契合问题的计算，根据计算结果，再判定是否折点到折点的未契合问题的计算。
折点为第三类夹角的修复过程如下：进行折点到折点的未契合问题的计算。
作为一种实施例而言，为了满足各种形状的二维矢量图形的边界未契合错误的处理，本发明按照折点涉及到的折线的夹角角度来对数据进行分类修复。
本发明在修复前的处理流程是先对折点进行分类；
1.如果折线为弧线夹角则直接进行点到点的数据计算，如果低于阈值则把折点的坐标修改为最近的那个低于阈值的其他图形的折点坐标。
2.如果折线为普通夹角,则先进行点到线的计算阶段，如果判定为交点在线段上则计算交点到折点的距离，如果距离小于阈值，就将折点的坐标修改为交点的坐标；然后再进行点到点的数据计算，如果低于阈值则把折点的坐标修改为最近的那个低于阈值的其他图形的折点坐标。
3.如果折点为直角夹角，那么先进行点到线的计算阶段，如果判定为交点在线段上则计算交点到折点的距离，如果距离小于阈值，就将折点的坐标修改为交点的坐标，如果折点涉及到的续相关的两个折点中不为进行数据计算的那个折点也是直角夹角的话，那个折点偏移原折点坐标到新折点坐标的偏移量；然后再进行点到点的数据计算，如果低于阈值则把折点的坐标修改为最近的那个低于阈值的其他图形的折点坐标，如果折点涉及到的续相关的两个折点中不为进行数据计算的那个折点也是直角夹角的话，那个折点偏移原折点坐标到新折点坐标的偏移量。之所以直角夹角需要相邻点一起偏移是因为折点和相邻点都是直角，一般建筑结构直角的结构大部分为准确建筑结构，如果只偏移其中一个点就会破坏它的整体结构如图4中的P 5折点往D 5 D 6中间的折线偏移的时候，如果只偏移P 5，那么整体结构就变了所以在P 6和P 5都是直角夹角的时候，需要一起往P 5偏移的量偏移。
综合来看，本发明的有益效果是：可以最大限度的保证原图形的完整性下批量处理因为人工测量的一些小误差造成不同层的分层分户图矢量二维图形的边界之间出现边界未契合的情况，可以大大减少人工量和处理时间。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种处理二维矢量图形之间的边界不契合方法，其特征在于：包括以下步骤：
S1、对分层分户图进行叠加得到建筑产权体模型底图；
S2、对建筑产权体模型底图进行属性挂接；
S3、对建筑产权体模型底图采用快速排斥法进行筛选，得到筛选后的单幢建筑产权体模型底图二维矢量数据；
S4、通过获取单幢建筑产权体模型底图中的所有二维矢量数据的折点数组来获取构成单个二维矢量数据的线型和角度来对筛选后的建筑产权体模型底图二维矢量数据进行数据分类；
S5、判断当前单个建筑产权体模型底图二维矢量数据与其它单个建筑产权体模型底图二维矢量数据涉及的边界未契合问题的类型；
S6、根据不同类别边界的折点数组和边界未契合问题的类型，完成边界契合修复；
步骤S3中，采用快速排斥法筛选的过程如下：
S31、通过已有的楼栋幢坐落属性，筛选需要处理的建筑产权体模型底图；
S32、设置步长阈值r，所述步长阈值用以约束多个建筑产权体模型底图的折点或边界的距离，当距离超过所述步长阈值时，进入快速排除法；
S33、获取使用幢坐落字段筛选的建筑产权体模型底图二维矢量数据的矩形包围框，将矩形包围框扩大一圈，扩大的值为所述步长r/2；
S34、遍历使用幢坐落数据筛选的建筑产权体模型底图，使用当前建筑产权体模型底图的包围框和其它建筑产权体模型底图的包围框进行一次快速排斥试验，初步筛选出与当前建筑产权体模型底图边界可能需要契合的其他建筑产权体模型底图，作为筛选后的建筑产权体模型底图。
2.如权利要求1所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：步骤S34中，所述快速排斥试验如下：
S341、获取当前建筑产权体模型底图的包围框和其它建筑产权体模型底图的包围框的坐标集合；
所述当前建筑产权体模型底图的包围框坐标为：（X 1min-r/2,Y 1max+r/2），（X 1max+r/2,Y 1max+r/2），（X 1max+r/2,Y 1min），（X 1min-r/2,Y 1min-r/2）；
所述其它建筑产权体模型底图的包围框的坐标为:（X 2min-r/2,Y 2max+r/2），（X 2max+r/2,Y 2max+r/2），（X 2max+r/2,Y 2min），（X 2min-r/2,Y 2min-r/2）；
S342、根据坐标集合判断两个包围框是否相切或相交，若是则判定其它建筑产权体模型底图的边界需要与当前建筑产权体模型底图的边界相契合；否则表示不需要契合。
3.如权利要求1所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：步骤S4中的数据分类过程具体如下：
S41、收集需要契合的建筑产权体模型底图的所有折点的坐标集合(P 1，P 2，P 3，……)；
S42、获取折点P i涉及到的折线L i和L i+1；其中i为折点编号；
S43、折线L i由折点P i和P i+1构成；折线L i+1由折点P i+1和P i+2构成；
S44、获取折线L i和L i+1的向量v i和v i+1；
S45、计算点P i+1涉及到的向量v i和v i+1之间的夹角α；
S46、当α处于大于第一预设角度或小于第二预设角度且和这个连续相关的两个折点夹角角度不是全部都大于等于第三预设角度时，认为点P i+1涉及到的夹角α属于第一类夹角；
当α小于等于第一预设角度或大于等于第二预设角度，则认为点P i+1涉及到的夹角α属于第二类夹角；
当α大于等于第三预设角度且连续相关的两个折点夹角角度全都大于第三预设角度或前两个折点全都大于第三预设角度或后两个折点全都大于第三预设角度，则认为夹角α属于点P i+1涉及到的第三类夹角；
其中，第二预设角度小于第一预设角度小于第三预设角度。
4.如权利要求3所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：步骤S5具体为：当建筑产权体模型底图的折点到另一个建筑产权体模型底图的折线或者折点的距离如果小于设置好的阈值且该折点到另一建筑产权体模型底图的折线或者折点的距离大于零，则该折点涉及到的边界存在边界未契合问题；
所述边界未契合问题分为：折点到折线的未契合问题和折点到折点的未契合问题。
5.如权利要求4所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：步骤S6中，当折点为第一类夹角的修复过程如下：先进行折点到折点的未契合问题的计算，再进行折点到折点的未契合问题的计算。
6.如权利要求4所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：步骤S6中，折点为第二类夹角的修复过程如下：先进行折点到折点的未契合问题的计算，根据计算结果，再判定是否折点到折点的未契合问题的计算。
7.如权利要求4所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：折点为第三类夹角的修复过程如下：进行折点到折点的未契合问题的计算。
8.如权利要求4所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：所述折点到折线的未契合问题的计算，具体如下:
获取折点P i(x 1, y 1)与折点P i+1(x 2, y 2)所构成的直线L；
获取线段LL j；线段LL j两个端点的坐标为(x 3, y 3) 和 (x 4, y 4)；
直线L的斜率设为m 1，m 1 = (y 2 - y 1) / (x 2 - x 1)
线段LL j的斜率设为m 2，m2 = (y 4 - y 3) / (x 4 - x 3)；
求直线L与线段LL j的斜距，其中直线L的斜距为b 1，b 1 = y 1 - m 1 * x 1；线段LL j的斜距为b 2，b 2 = y 3 - m 2 * x 3；
计算直线L与线段LL j延长后的直线的交点，设交点的坐标为（x，y），x = (b 2 - b 1) /(m 1 - m 2)，y = m 1 * x + b 1；
判断交点是否在线段LL j上，若是，则计算交点到折点P i的距离，若距离小于阈值，就将折点的坐标修改为交点的坐标；
其中，交点在线段LL j上需要满足下面条件：min(x 3, x 4) <= x <= max(x 3, x 4)并且min(y 3, y 4) <= y <= max(y 3, y 4)；其中min是指括号中的两个值当中最小的那个值，max是指括号中的两个值当中最大的那个值。
9.如权利要求4所述的一种处理二维矢量图形之间的边界不契合方法，其特征在于：所述折点到折点的未契合问题的计算，具体如下：遍历求取建筑产权体模型底图所有折点与另外一个建筑产权体模型底图的所有折点之间的距离，如果这个建筑产权体模型底图的一个折点与另外一个建筑产权体模型底图的折点之间的距离小于阈值则将小于阈值另外一个建筑产权体模型底图的折点坐标按照距离进行排序只保留距离最近的那个，并将这个折点的坐标更新为保留的那个折点。
CN202310956478.0A 2023-08-01 2023-08-01 一种处理二维矢量图形之间的边界不契合方法 Active CN116664455B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310956478.0A CN116664455B (zh) 2023-08-01 2023-08-01 一种处理二维矢量图形之间的边界不契合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310956478.0A CN116664455B (zh) 2023-08-01 2023-08-01 一种处理二维矢量图形之间的边界不契合方法

Publications (2)

Publication Number Publication Date
CN116664455A true CN116664455A (zh) 2023-08-29
CN116664455B CN116664455B (zh) 2023-11-14

Family

ID=87714029

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310956478.0A Active CN116664455B (zh) 2023-08-01 2023-08-01 一种处理二维矢量图形之间的边界不契合方法

Country Status (1)

Country Link
CN (1) CN116664455B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109544687A (zh) * 2018-11-21 2019-03-29 南京师范大学 一种定向分层分户图扩展方法
US20190251209A1 (en) * 2018-02-09 2019-08-15 Guangzhou University Method for Automatic Modeling of an Architecture Based on an Architectural Drawing
CN110517353A (zh) * 2019-08-30 2019-11-29 西南交通大学 基于二维矢量图及少量高程点的快速建筑物三维重建方法
CN112489212A (zh) * 2020-12-07 2021-03-12 武汉大学 一种基于多源遥感数据的建筑物智能化三维测图方法
CN115222873A (zh) * 2022-05-19 2022-10-21 南京信息工程大学 一种基于标注和空间语义的三维产权体图属批量关联方法
CN115525951A (zh) * 2022-09-28 2022-12-27 土豆数据科技集团有限公司 建筑物轮廓优化及特征边重构的规则化方法及装置
CN116310251A (zh) * 2023-02-03 2023-06-23 北斗伏羲中科数码合肥有限公司 房产ifc模型与建筑矢量边界自动坐标配准方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190251209A1 (en) * 2018-02-09 2019-08-15 Guangzhou University Method for Automatic Modeling of an Architecture Based on an Architectural Drawing
CN109544687A (zh) * 2018-11-21 2019-03-29 南京师范大学 一种定向分层分户图扩展方法
CN110517353A (zh) * 2019-08-30 2019-11-29 西南交通大学 基于二维矢量图及少量高程点的快速建筑物三维重建方法
CN112489212A (zh) * 2020-12-07 2021-03-12 武汉大学 一种基于多源遥感数据的建筑物智能化三维测图方法
CN115222873A (zh) * 2022-05-19 2022-10-21 南京信息工程大学 一种基于标注和空间语义的三维产权体图属批量关联方法
CN115525951A (zh) * 2022-09-28 2022-12-27 土豆数据科技集团有限公司 建筑物轮廓优化及特征边重构的规则化方法及装置
CN116310251A (zh) * 2023-02-03 2023-06-23 北斗伏羲中科数码合肥有限公司 房产ifc模型与建筑矢量边界自动坐标配准方法及装置

Also Published As

Publication number Publication date
CN116664455B (zh) 2023-11-14

Similar Documents

Publication Publication Date Title
TWI773063B (zh) 圖像產生系統及圖像產生方法
US20180025541A1 (en) Method for automatic modeling of complex buildings with high accuracy
WO2010103942A1 (ja) 表面加工データの作成方法および装置
CN102812463A (zh) 能够进行三维对象模型的3d印刷的方法和系统
US20140063038A1 (en) Method and device for forming surface treatment data
CN113538373A (zh) 一种基于三维点云的施工进度自动检测方法
CN112700529A (zh) 根据规范文档生成三维模型的方法和系统
CN116051533B (zh) 一种电力杆塔的施工质量检测方法及装置
Li et al. A deep learning-based indoor acceptance system for assessment on flatness and verticality quality of concrete surfaces
CN113868904B (zh) 一种基于三维点云自动化生成既有结构有限元模型的方法
CN116664455B (zh) 一种处理二维矢量图形之间的边界不契合方法
CN102506753A (zh) 基于十四点球面小波变换的不规则零件形状差异检测方法
CN102012964A (zh) 激光扫描的采样数据的处理方法和装置
Gökgöz Generalization of contours using deviation angles and error bands
Sindram et al. Voluminator 2.0–speeding up the approximation of the volume of defective 3D building models
CN114170160A (zh) 一种基于几何模型的钢筋分类扫描规划方法
CN118229670A (zh) 一种智能的轮毂模型尺寸快速测量方法
CN107526859A (zh) 设计复合零件的外表面
CN118229918A (zh) 基于工业ct图像的三维建模与孔洞修复方法及系统
CN111089539A (zh) 晶圆轮廓图的构造方法
Sugihara et al. Roof report from automatically generated 3D building models by straight skeleton computation
CN112017294B (zh) 基于可见性判断的三维cad网格模型简化方法和系统
Steinbrenner et al. Procedural CAD model edge tolerance negotiation for surface meshing
KR20110037015A (ko) 연속지적도 오류 발견 프로그램 및 그 방법
CN112948919B (zh) 一种基于图像细化的bim模型跨楼层路网提取方法

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