CN107748793B - 一种基础地理数据的共点检查的方法及系统 - Google Patents
一种基础地理数据的共点检查的方法及系统 Download PDFInfo
- Publication number
- CN107748793B CN107748793B CN201711061203.1A CN201711061203A CN107748793B CN 107748793 B CN107748793 B CN 107748793B CN 201711061203 A CN201711061203 A CN 201711061203A CN 107748793 B CN107748793 B CN 107748793B
- Authority
- CN
- China
- Prior art keywords
- node
- coordinate
- equation
- representing
- formula
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Instructional Devices (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种共点检查的方法及系统,属于地理信息技术领域,具体是涉及一种基础地理数据的共点检查的方法及系统。本发明采用四叉树对数据进行分块,并采用多线程的处理机制,能高效的检查基础地理数据的共点错误,准确度高速度快,大大降低了基础地理数据共点错误检查的时间成本和人力成本。
Description
技术领域
本发明涉及一种共点检查的方法及系统,属于地理信息技术领域,具体是涉及一种基础地理数据的共点检查的方法及系统。
背景技术
常规基础地理数据的采集通常有实测法、图解法、综合法三种方法。实际作业生产中,这三种方法都可能因为机器测量误差和人工误差产生共点错误:
i.实测法。由于测量的机器,如GPS、全站仪,本身在测量过程中会存在一定误差,所以采集的基础地理数据也会对应着存在共点错误;
ii.图解法。由人工在数字正射影像上采集基础地理数据。由于人工勾绘本身存在误差,所以采集的基础地理数据也会存在共点错误;
iii.综合法是结合前两个方法的一种综合的方法,所以也会产生共点错误;
共点错误指的是因为一定原因,同一个点位的地方,不同地物的点位实际上有偏差。常见共点错误如图1、图2所示。
目前,检查共点错误的最主要的技术方案是人工检查的方法。其具体方法步骤如下:
①用户加载需要检查共点错误的数据;
②用户打开软件的节点查看功能;
③用户选择一个基础地理图形,查看图形的所有节点,然后定位到一个节点;
④用户缩放节点范围,以大比例尺进行查看;
⑤用户测量该节点和周边图形最近节点或者边的距离L,如图3所示。
⑥设置计算容差A(比如0.0001米),和测绘精度要求B(比如GPS一般的精度0.05米);
⑦进行比对:
i.如果距离L大于计算容差A,且小于测绘精度要求B,则该节点存在共点错误;
ii.如果距离L小于A或者大于B,则不存在共点错误;
⑧重复第三步至第第七步,完成一个图形的共点检查;
传统手工共点错误检查的技术方法,耗时耗力,效率低下,下面针对传统方法和本发明进行对比说明。
以农村承包地块为例来说明人工检查共点错误的工作量。以湖北省的某县下的一个乡镇作为数据来源,这个乡镇有9437户,共80048块地,面积为47868亩。一个县往往有十几个乡镇,这里按照10个乡镇计算,那么一个县的地块总数为:80048*10=800480。假定这些地块平均一个地块4个节点,那么,总节点数为4*800480=3201920个。经过实际测试,一个内业人员,假定平均检查一个节点的时间为1分钟,则总时间需要1*3201920=3201920分钟,即53365.33小时。成本上,按5000元一月的内业作业人员成本计算,每小时成本5000/22/8=28.41元每小时,总费用就是28.41*53365.33=1516109.03元。即,一个县的共点错误检查时间成本为53365.33小时,人力成本为1516109.03元。
从时间成本上,使用本发明的算法,基本上一个地块的共点检查只需要0.01秒的时间就可以完成,一个县的共点检查,只需要800480*0.01/60/60=2.22小时即可;从人力成本上,一个县的共点检查,人力成本2.22*28.41=63.07元,相对于1516109.03元的传统人工检查成本,可以忽略不计。
综合以上,本发明可以降低共点错误检查的时间成本和人力成本,提高共点检查效率,更好的满足实际作业需要。
发明内容
本发明主要是解决现有技术所存在的人工检查共点错误需要大量的时间成本、人力成本,不能满足实际基础地理数据生产的需求的问题,提供了一种基础地理数据的共点检查的方法。该方法使用计算机遍历所有基础地理数据的图形,依次进行共点检查,利用算法的高效性和计算机的计算性能,从而实现高效快速的共点检查。
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种基础地理数据的共点检查的方法,包括:
包络分块步骤,对包络范围内的数据进行分块,每分块内包括的要素个数不超过预设阈值;
对象收集步骤,查找与各分块单元相交的所有几何对象,收集几何对象的节点信息和/或折线信息;
共点查找步骤,遍历所述节点信息和/或折线信息以查找共点数据。
优选的,上述的一种基础地理数据的共点检查的方法,所述包络分块步骤中采用四叉树算法对包络范围进行数据分块。
优选的,上述的一种基础地理数据的共点检查的方法,数据的分块基于以下公式:
式中:e1,e2,e3,e4分别为由矩形框四分得到的四个子矩形框;(x1,y1,x2,y2)为新矩形框的包络范围;xmin、xmax、ymin、ymax分别为原包络范围的坐标值。
优选的,上述的一种基础地理数据的共点检查的方法,所述对象收集步骤中,几何对象的查找基于下式:
G={(Xmn,Ymn)|Xmn∈[Xmin,Xmax],Ymn∈[Ymin,Ymax],m∈N,n∈N} 公式5
式中:G为分块单元内所有点、线、面几何对象集合;Xmn为第m个要素中第n个节点的X坐标值;Ymn为为第m个要素中第n个节点的Y坐标值;Xmin、Xmax、Ymin、Ymax分别为分块单元的范围值;N为自然数。
优选的,上述的一种基础地理数据的共点检查的方法,节点信息的收集基于以下公式:
N={(x,y,i,j)|x=X(i,j),y=Y(i,j),i>=0,j>=0} 公式6
式中:(x,y,i,j)为第i个几何对象的第j个节点的节点信息;x为x坐标值;y为y坐标值;i为几何对象编号;j为节点编号;X(i,j)为第i个几何对象的第j个节点的x坐标值;Y(i,j)为第i个几何对象的第j个节点的y坐标值。
和/或
基于下式收集几何对象的折线信息。
L={(n1,n2)|n2=f(n1),n1∈N,n2∈N} 公式6
式中:(n1,n2)为收集的线段;L为收集的线段集合,n2=f(n1)表示节点顺序的映射关系。
优选的,上述的一种基础地理数据的共点检查的方法,所述共点查找步骤中,遍历节点信息,查找本几何对象内部节点重叠的情况,其中,基于下式查找重复节点:
N={(n1,n2)||X(n1)-X(n2)|<e,|Y(n1)-Y(n2)|<e,e>0} 公式7
式中,N为重复节点对集合;(n1,n2)为重复节点对;X(n1)为n1的X坐标;X(n2)为n2的X坐标;Y(n1)为n1的Y坐标;Y(n2)为n2的Y坐标;e为容差值。
优选的,上述的一种基础地理数据的共点检查的方法,所述共点查找步骤中,遍历节点信息,查找不同几何对象之间节点距离过近的情况,其中,基于下式查找过近节点:
N={(n1,n2)|0<|X(n1)-X(n2)|<e,0<|Y(n1)-Y(n2)|<e,O(n1)≠O(n2),e>0} 公式8
式中,N表示检查结果的节点对集合,(n1,n2)表示检查结果的节点对,X(n2)为n2的X坐标;Y(n1)为n1的Y坐标;Y(n2)为n2的Y坐标;e为容差值;O(n1)表示n1节点所在几何对象的编号;O(n2)表示n2节点所在几何对象的编号。
优选的,上述的一种基础地理数据的共点检查的方法,所述共点查找步骤中,遍历节点信息,检查节点与线段的悬挂问题,其中,节点悬挂基于下式获得:
a=Y(n2)-Y(n1) 公式9
b=X(n1)-X(n2) 公式10
c=-aX(n1)-bY(n1) 公式11
式中,X(n1)表示线段起节点的X坐标;X(n2)表示线段末节点的X坐标;Y(n1)表示线段起节点的Y坐标;Y(n2)表示线段末节点的Y坐标;X(n3)表示当前节点的X坐标;Y(n3)表示当前节点的Y坐标;L表示当前节点到线段的最短距离。
因此,本发明具有如下优点:采用四叉树对数据进行分块,并采用多线程的处理机制,能高效的检查基础地理数据的共点错误,准确度高速度快,大大降低了基础地理数据共点错误检查的时间成本和人力成本。
附图说明
附图1是共点错误情况1;
附图2是共点错误情况2;
附图3是测量距离L;
附图4是技术流程图;
附图5是数据分块图;
附图6是6内部节点检查图;
附图7是外部节点距离小于容差图;
附图8是悬挂节点检查图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:
本方案用于检查基础地理数据的共点问题,工作流程如图4所示,包括:
(1)加载需要进行共点检查的数据。
(2)采用四叉树的算法对包络范围内的数据进行分块,要素个数超出阈值时,继续分块。
式中:em1,em2,em3,em4分别为由矩形框四分得到的四个子矩形框;(xmin,xmax,ymin,ymax)为新矩形框的包络范围;xmin、xmax、ymin、ymax分别为原包络范围的坐标值。
(3)根据每个分块单元,查找与本分块单元相交的所有几何对象。
式中:G为分块单元内所有点、线、面几何对象集合;Xmn为第m个要素中第n个节点的X坐标值;Ymn为为第m个要素中第n个节点的Y坐标值;Xmin、Xmax、Ymin、Ymax分别为分块单元的范围值。
(4)收集所有几何对象的节点信息。
N={(x,y,i,j)|x=X(i,j),y=(i,j),i>=0,j>=0} 公式6
式中:(x,y,i,j)为第i个几何对象的第j个节点的节点信息;x为x坐标值;y为y坐标值;i为几何对象编号;j为节点编号。
(5)收集几何对象的折线信息。
L={(n1,n2)|n2=f(n1),n1∈N,n2∈N} 公式6
式中:(n1,n2)为收集的线段;L为收集的线段集合,n2=f(n1)表示节点顺序的映射关系。
(6)遍历节点信息,查找本几何对象内部节点重叠的情况,如图所示,节点距离小于容差e,则节点重复。
几何内部节点重复的公式如下:
N={(n1,n2)||X(n1)-X(n2)|<e,|Y(n1)-Y(n2)|<e,e>0} 公式7
式中,N为重复节点对集合;(n1,n2)为重复节点对;X(n1)为n1的X坐标;X(n2)为n2的X坐标;Y(n1)为n1的Y坐标;Y(n2)为n2的Y坐标;e为容差值。
(6)遍历节点信息,查找不同几何对象之间节点距离过近的情况,如图所示,节点之间有距离小于e,则该处存在共点错误。
节点距离过近的公式如下:
N={(n1,n2)|0<|X(n1)-X(n2)|<e,0<|Y(n1)-Y(n2)|<e,O(n1)≠O(n2),e>0} 公式8
式中,N表示检查结果的节点对集合,(n1,n2)表示检查结果的节点对,X(n2)为n2的X坐标;Y(n1)为n1的Y坐标;Y(n2)为n2的Y坐标;e为容差值;O(n1)表示n1的OID;O(n2)表示n2的OID。
(7)遍历节点信息,检查节点与线段的悬挂问题,如图所示,节点在线段上投影的距离小于e,而找不到对应的节点,则存在共点错误。
节点悬挂的计算公式如下:
a=Y(n2)-Y(n1) 公式9
b=X(n1)-X(n2) 公式10
c=-aX(n1)-bY(n1) 公式11
式中,X(n1)表示线段起节点的X坐标;X(n2)表示线段末节点的X坐标;Y(n1)表示线段起节点的Y坐标;Y(n2)表示线段末节点的Y坐标;X(n3)表示当前节点的X坐标;Y(n3)表示当前节点的Y坐标;L表示当前节点到线段的最短距离。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (4)
1.一种基础地理数据的共点检查的方法,其特征在于,包括:
包络分块步骤,采用四叉树算法对包络范围内的数据进行分块,每分块内包括的要素个数不超过预设阈值;
对象收集步骤,查找与各分块单元相交的所有几何对象,收集几何对象的节点信息和/或折线信息;其中,基于以下公式查找几何对象:
G={(Xmn,Ymn)|Xmn∈[Xmin,Xmax],Ymn∈[Ymin,Ymax],m∈N,n∈N} 公式5
式中:G为分块单元内所有点、线、面几何对象集合;Xmn为第m个要素中第n个节点的X坐标值;Ymn为为第m个要素中第n个节点的Y坐标值;Xmin、Xmax、Ymin、Ymax分别为分块单元的范围值;N为自然数;
共点查找步骤,遍历所述节点信息和/或折线信息以查找共点数据;
数据的分块基于以下公式:
式中:e1,e2,e3,e4分别为由矩形框四分得到的四个子矩形框;(x1,y1,x2,y2)为新矩形框的包络范围;xmin、xmax、ymin、ymax分别为原包络范围的坐标值;
所述共点查找步骤中,遍历节点信息,检查节点与线段的悬挂问题,其中,节点悬挂基于下式获得:
a=Y(n2)-Y(n1) 公式9
b=X(n1)-X(n2) 公式10
c=-aX(n1)-bY(n1) 公式11
式中,X(n1)表示线段起节点的X坐标;X(n2)表示线段末节点的X坐标;Y(n1)表示线段起节点的Y坐标;Y(n2)表示线段末节点的Y坐标;X(n3)表示当前节点的X坐标;Y(n3)表示当前节点的Y坐标;L表示当前节点到线段的最短距离。
2.根据权利要求1所述的一种基础地理数据的共点检查的方法,其特征在于,节点信息的收集基于以下公式:
N={(x,y,i,j)|x=X(i,j),y=Y(i,j),i>=0,j>=0} 公式6
式中:(x,y,i,j)为第i个几何对象的第j个节点的节点信息;x为x坐标值;y为y坐标值;i为几何对象编号;j为节点编号;X(i,j)为第i个几何对象的第j个节点的x坐标值;Y(i,j)为第i个几何对象的第j个节点的y坐标值;
和/或
基于下式收集几何对象的折线信息:
L={(n1,n2)|n2=f(n1),n1∈N,n2∈N} 公式6
式中:(n1,n2)为收集的线段;L为收集的线段集合,n2=f(n1)表示节点顺序的映射关系。
3.根据权利要求1所述的一种基础地理数据的共点检查的方法,其特征在于,所述共点查找步骤中,遍历节点信息,查找本几何对象内部节点重叠的情况,其中,基于下式查找重复节点:
N={(n1,n2)||X(n1)-X(n2)|<e,|Y(n1)-Y(n2)|<e,e>0} 公式7
式中,N为重复节点对集合;(n1,n2)为重复节点对;X(n1)为n1的X坐标;X(n2)为n2的X坐标;Y(n1)为n1的Y坐标;Y(n2)为n2的Y坐标;e为容差值。
4.根据权利要求1所述的一种基础地理数据的共点检查的方法,其特征在于,所述共点查找步骤中,遍历节点信息,查找不同几何对象之间节点距离过近的情况,其中,基于下式查找过近节点:
N={(n1,n2)|0<|X(n1)-X(n2)|<e,0<|Y(n1)-Y(n2)|<e,O(n1)≠O(n2),e>0} 公式8
式中,N表示检查结果的节点对集合,(n1,n2)表示检查结果的节点对,X(n2)为n2的X坐标;Y(n1)为n1的Y坐标;Y(n2)为n2的Y坐标;e为容差值;O(n1)表示n1节点所在几何对象的编号;O(n2)表示n2节点所在几何对象的编号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711061203.1A CN107748793B (zh) | 2017-11-02 | 2017-11-02 | 一种基础地理数据的共点检查的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711061203.1A CN107748793B (zh) | 2017-11-02 | 2017-11-02 | 一种基础地理数据的共点检查的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107748793A CN107748793A (zh) | 2018-03-02 |
CN107748793B true CN107748793B (zh) | 2020-04-14 |
Family
ID=61254325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711061203.1A Active CN107748793B (zh) | 2017-11-02 | 2017-11-02 | 一种基础地理数据的共点检查的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107748793B (zh) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7949668B2 (en) * | 2001-08-20 | 2011-05-24 | Pardalis, Inc. | Common point authoring system for the complex sharing of hierarchically authored data objects in a distribution chain |
CN103399730B (zh) * | 2013-07-29 | 2016-08-10 | 中国科学院地理科学与资源研究所 | 一种基于cpu+gpu架构的空间几何体线段相交判断并行处理方法 |
MX2016000639A (es) * | 2013-08-16 | 2016-11-10 | Landmark Graphics Corp | Generacion de representaciones de estructuras geologicas reconocibles a partir de una recopilacion de puntos comunes. |
CN106204305A (zh) * | 2016-07-13 | 2016-12-07 | 武大吉奥信息技术有限公司 | 一种地块自动批量切分方法及装置 |
CN106202433A (zh) * | 2016-07-13 | 2016-12-07 | 武大吉奥信息技术有限公司 | 一种基础地理数据接边入库方法 |
-
2017
- 2017-11-02 CN CN201711061203.1A patent/CN107748793B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107748793A (zh) | 2018-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102147260B (zh) | 电子地图匹配方法和装置 | |
CN102800052B (zh) | 非标准地图的半自动数字化方法 | |
CN104703143A (zh) | 一种基于wifi信号强度的室内定位方法 | |
CN103247062B (zh) | 采集农田关键点测绘成图的方法 | |
CN107044852B (zh) | 不平整状态下全站仪测量方法 | |
CN107529644B (zh) | 一种电力系统静态电压稳定域边界线性近似方法 | |
CN109947881B (zh) | 一种poi判重方法、装置、移动终端及计算机可读存储介质 | |
CN109215338B (zh) | 一种基于网格化的北斗gps车辆轨迹管理系统及其方法 | |
CN105203120B (zh) | 导航路线评测方法及装置 | |
CN104809724A (zh) | 多波段遥感影像的自动精配准方法 | |
CN107885790A (zh) | 一种道路空间网络多因子自动更新方法 | |
CN105674926A (zh) | 点云拼接方法及系统 | |
CN107133236A (zh) | 地理位置轨迹优化的方法和装置 | |
CN104408055A (zh) | 一种激光雷达点云数据的存储方法及装置 | |
JP2014132711A (ja) | 複数の測位方式に基づく移動端末の位置情報から対象エリアを特定するエリア管理サーバ、プログラム及び方法 | |
CN111079773A (zh) | 一种基于Mask R-CNN网络的砾石参数获取方法、装置、设备及存储介质 | |
CN102609501A (zh) | 一种基于实时历史数据库的数据清洗方法 | |
CN108776660A (zh) | 一种基于ArcGIS的批量匹配道路属性的方法 | |
CN104636401A (zh) | 一种scada系统数据回滚的方法及装置 | |
CN103760586B (zh) | 一种在gps姿态测量中快速探测与修复周跳的方法 | |
CN107748793B (zh) | 一种基础地理数据的共点检查的方法及系统 | |
CN102506812A (zh) | 一种变形监测基准点稳定性判断的vt检验方法 | |
CN113516406A (zh) | 一种基于实时观测分析的高铁沿线降水测点布置方法 | |
CN110348483B (zh) | 基于空间信息与统计学习的空间信息变化检测与分类方法 | |
CN112381941B (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 | ||
CP03 | Change of name, title or address |
Address after: 430000 Wuda science and Technology Park, Jiangxia Avenue, Miaoshan District, Donghu Development Zone, Wuhan City, Hubei Province Patentee after: Geospace Information Technology Co., Ltd. Address before: 430223 Wuda science and Technology Park, Jiangxia Avenue, Miaoshan community, Donghu Development Zone, Wuhan City, Hubei Province Patentee before: WUDA GEOINFORMATICS Co.,Ltd. |
|
CP03 | Change of name, title or address |