CN115731284B - 一种界址点图斑绘制及面积计算算法 - Google Patents

一种界址点图斑绘制及面积计算算法 Download PDF

Info

Publication number
CN115731284B
CN115731284B CN202211477418.2A CN202211477418A CN115731284B CN 115731284 B CN115731284 B CN 115731284B CN 202211477418 A CN202211477418 A CN 202211477418A CN 115731284 B CN115731284 B CN 115731284B
Authority
CN
China
Prior art keywords
coordinate
area
coordinates
boundary
longitude
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
CN202211477418.2A
Other languages
English (en)
Other versions
CN115731284A (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.)
Zhongke Xingtu Smart Technology Anhui Co ltd
Original Assignee
Zhongke Xingtu Smart Technology Anhui 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 Zhongke Xingtu Smart Technology Anhui Co ltd filed Critical Zhongke Xingtu Smart Technology Anhui Co ltd
Priority to CN202211477418.2A priority Critical patent/CN115731284B/zh
Publication of CN115731284A publication Critical patent/CN115731284A/zh
Application granted granted Critical
Publication of CN115731284B publication Critical patent/CN115731284B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Processing Or Creating Images (AREA)

Abstract

本发明公开了一种界址点图斑绘制及面积计算算法,涉及技术领域,基于非GIS的报盘导入坐标数据,通过分割方式识别坐标类型,进行分组,通过坐标经纬度进行图斑绘制,再经过坐标圈号分析图斑,结合坐标数据源,计算图斑面积;图斑绘制是通过界址点坐标数据信息进行分组,坐标经纬度进行60进制转10进制转换,根据传入的分割坐标方式对坐标数据源进行分组,通过缩放、偏移处理分组后的坐标数据,根据坐标经纬度构造图形路径,遍历坐标数据源设置拐点号,根据坐标画图,并且通过偏移处理不至于使图形过于紧密;解决了计算复杂度以及误差大的问题。

Description

一种界址点图斑绘制及面积计算算法
技术领域
本发明属于界址点图斑领域,涉及界址点图斑绘制技术,具体是一种界址点图斑绘制及面积计算算法。
背景技术
随着时代的发展和科技的进步,各行各业都面临着行业相关数据爆炸增长的态势,有效、快速的整理复杂的数据就成为了行业急需解决的痛点。而用地图形绘制以及面积计算,一直都是从业人员面临的难题之一。
传统方法虽然可以借助GIS软件,绘制界址点图斑图形与计算面积,但是仍然费时费力,需用到矢量数据重新分类、分析、统计数据,才能绘制出图斑图形及计算面积,出错率高,且仅能绘制图形,再通过复杂方法计算面积,并且不包含坐标经纬度,不能实现数据与图斑的直接联系与应用,分析方法手段落后;
图斑图形绘制与面积计算一直以来都是很重要的工作,图斑绘制与面积计算为城镇建设用地统计、耕地保护面积统计以及各类其他用地的统计分析等应用非常重要,能够有效的帮助相关从业人员以及政企事业单位对于地形图形绘制和面积计算统计。目前,界址点图斑图形绘制以及面积计算的主要方式有GIS建模与图解分析、GIS类地图斑净面积计算、GIS矢量数据界址点图斑计算以及GIS几何计算与Shp转GDB图斑边界及面积计算,这些方式均基于GIS手段才能实现,需要一定专业能力与软件操作熟练度能力,均需要GIS数据才可实现;
具体来说,基于GIS手段解析界址点图斑的方法需要大量矢量数据以及兼容GIS相关软件,并且无法做到分析界址点图斑用地类型,无法及时检测、计算需要的用地类型图斑面积。并且以上方法及手段并无相关经纬度数据,需要手动换算成需要的坐标类型,其相关计算过程与手段方法比较复杂;
为此,提出一种界址点图斑绘制及面积计算算法。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。为此,本发明提出一种界址点图斑绘制及面积计算算法,该一种界址点图斑绘制及面积计算算法基于非GIS的报盘导入坐标数据,通过分割方式识别坐标类型,进行分组,通过坐标经纬度进行图斑绘制,再经过坐标圈号分析图斑,结合坐标数据源,计算图斑面积;图斑绘制是通过界址点坐标数据信息进行分组,坐标经纬度进行60进制转10进制转换,根据传入的分割坐标方式对坐标数据源进行分组,通过缩放、偏移处理分组后的坐标数据,根据坐标经纬度构造图形路径,遍历坐标数据源设置拐点号,根据坐标画图,并且通过偏移处理不至于使图形过于紧密;解决了计算复杂度以及误差大的问题。
为实现上述目的,根据本发明的第一方面的实施例提出一种界址点图斑绘制及面积计算算法,包括以下步骤:
步骤一:报盘导入坐标信息,根据传入的坐标类型参数判断待生成图形的地块坐标类型;
步骤二:根据传入的坐标类型判断要生成的图形的地块坐标类型,并根据坐标类型进行对应的坐标转换操作;
步骤三:根据坐标信息中的分割坐标方式对坐标数据源进行分组,并对分组处理后的坐标数据根据实际经验进行缩放以及偏移处理;
步骤四:根据坐标信息中的坐标经纬度构造图形路径,遍历坐标数据源设置拐点号,并根据坐标画图,通过原点将多边形拆分成三角形,再通二维向量叉乘法求三角形面积;
步骤五:获得坐标信息中最小的经度值和纬度值对坐标进行首次偏移,保证坐标从图片的左上角开始,以及保证各个点之间的距离不变,待放大之后各点间的距离变大解决各个坐标在图形上显示的过于紧密的问题;获得坐标信息能放大的最大倍数并放大,将放大后的坐标偏移到图片的中央;获得界址点图斑会制图;
步骤六:使用界址点图斑面积算法求解界址点图斑面积;
其中,所述坐标信息包括但不限于项目ID、坐标圈号、生成图形的宽度、生成图形的宽度、坐标信息的外键、拐点号显示标记、坐标分割方式、坐标类型以及增列矿种标记;其中,坐标信息的外键用于查询导入的地块的坐标信息;坐标分割方式包括0与1两个值;其中,0代表*号分割方式,1代表圈分割方式;其中,坐标类型包括0、1、2、3、4五类,其中,0代表油气非油气探矿与油气采矿,1代表建筑用地,2代表非油气采矿,3代表增列矿种坐标专用,4代表城市边界规模调整shape文件图形预览;
其中,所述根据坐标类型进行对应的坐标转换操作包括以下方式:
如果坐标类型是油气非油气探矿与油气采矿,进行第一坐标转换操作;
如果坐标类型为建设用地,则进行第二坐标转换操作;
如果坐标类型是非油气采矿,则进行第三坐标转换操作;
如果坐标类型是增列矿种坐标专用,则进行第一坐标转换操作;
如果坐标类型是城市边界规模调整shape文件图形预览:则进行第二坐标转换操作;
其中,所述第一坐标转换操作包括以下步骤:
步骤S1:根据坐标信息的外键查询出坐标数据源;
步骤S2:把坐标数据源纵坐标和横坐标进行60进制转换为10进制坐标点;
其中,所述第二坐标转换操作为根据坐标信息的外键查询出坐标数据源;
其中,所述第三坐标转换操作包括以下步骤:
步骤P1:根据坐标信息的外键查询出坐标数据源;
步骤P2:判断坐标类型是直角坐标还是经纬坐标;若坐标类型为直角坐标,将坐标数据的源经度数据和维度数据进行调换;
其中,所述使用界址点图斑面积算法求解界址点图斑面积包括以下步骤:
步骤Z1:根据报盘导入坐标信息,根据传入的所属项目ID和坐标圈号查询出坐标数据源信息,以及坐标信息和分割方式对坐标数据源进行分组处理坐标信息;
步骤Z2:循环所有地块下面的坐标后根据多边形各个点坐标,计算出当前地块对应坐标的多边形的面积;
步骤Z3:根据圈号判断是主区域,还是挖空地块,圈号1为主区域、其他都为挖空地块;
步骤Z4:如果是主区域则多个坐标面积相加,其他都为挖空地块则多个坐标面积相减;
步骤Z5:将计算出的面积单位换算平方米转换成公顷;
其中,所述计算出当前地块对应坐标的多边形的面积包括以下步骤:
步骤X1:将面积标记为S,计算面积的初始值,计算公式为S=w1*(jN-j 1);其中,w1为坐标集合的第一个坐标的维度,N为坐标集合中坐标的总量;j为坐标集合中每个坐标的经度;jN即为第N个坐标的经度;j 1为第1个坐标的经度;
步骤X2:将面积S进行迭代更新;迭代更新的方式为:设置变量n=1;对于每个n=1,2,3,…N,将面积S依次更新为S+wn*(j(n-1)-j((n+1)%N));
步骤X3:将面积S更新为|S/2|,获得当前地块对应坐标的多边形的面积。
与现有技术相比,本发明的有益效果是:
本发明基于非GIS的报盘导入坐标数据,通过分割方式识别坐标类型,进行分组,通过坐标经纬度进行图斑绘制,再经过坐标圈号分析图斑,结合坐标数据源,计算图斑面积;图斑绘制是通过界址点坐标数据信息进行分组,坐标经纬度进行60进制转10进制转换,根据传入的分割坐标方式对坐标数据源进行分组,通过缩放、偏移处理分组后的坐标数据,根据坐标经纬度构造图形路径,遍历坐标数据源设置拐点号,根据坐标画图,并且通过偏移处理不至于使图形过于紧密;解决了计算复杂度以及误差大的问题。
附图说明
图1为本发明的流程图。
具体实施方式
下面将结合实施例对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
如图1所示,一种界址点图斑绘制及面积计算算法,包括以下步骤:
步骤一:报盘导入坐标信息,根据传入的坐标类型参数判断待生成图形的地块坐标类型;
步骤二:根据传入的坐标类型判断要生成的图形的地块坐标类型,并根据坐标类型进行对应的坐标转换操作;
步骤三:根据坐标信息中的分割坐标方式对坐标数据源进行分组,并对分组处理后的坐标数据根据实际经验进行缩放以及偏移处理;
步骤四:根据坐标信息中的坐标经纬度构造图形路径,遍历坐标数据源设置拐点号,并根据坐标画图,通过原点将多边形拆分成三角形,再通二维向量叉乘法求三角形面积;
步骤五:获得坐标信息中最小的经度值和纬度值对坐标进行首次偏移,保证坐标从图片的左上角开始,以及保证各个点之间的距离不变,待放大之后各点间的距离变大解决各个坐标在图形上显示的过于紧密的问题;获得坐标信息能放大的最大倍数并放大,将放大后的坐标偏移到图片的中央;获得界址点图斑会制图;
步骤六:使用界址点图斑面积算法求解界址点图斑面积;
可以理解的是,基于GIS手段解析界址点图斑的方法需要大量矢量数据以及兼容GIS相关软件,并且无法做到分析界址点图斑用地类型,无法及时检测、计算需要的用地类型图斑面积。并且以上方法及手段并无相关经纬度数据,需要手动换算成需要的坐标类型,其相关计算过程与手段方法比较复杂;
在一个优选的实施例中,所述坐标信息包括但不限于项目ID、坐标圈号、生成图形的宽度、生成图形的宽度、坐标信息的外键、拐点号显示标记、坐标分割方式、坐标类型以及增列矿种标记;其中,坐标信息的外键用于查询导入的地块的坐标信息;坐标分割方式包括0与1两个值;其中,0代表*号分割方式,1代表圈分割方式;其中,坐标类型包括0、1、2、3、4五类,其中,0代表油气非油气探矿与油气采矿,1代表建筑用地,2代表非油气采矿,3代表增列矿种坐标专用,4代表城市边界规模调整shape文件图形预览;
在一个优选的实施例中,所述根据坐标类型进行对应的坐标转换操作包括以下方式:
如果坐标类型是油气非油气探矿与油气采矿,进行第一坐标转换操作;
如果坐标类型为建设用地,则进行第二坐标转换操作;
如果坐标类型是非油气采矿,则进行第三坐标转换操作;
如果坐标类型是增列矿种坐标专用,则进行第一坐标转换操作;
如果坐标类型是城市边界规模调整shape文件图形预览:则进行第二坐标转换操作;
在一个优选的实施例中,所述第一坐标转换操作包括以下步骤:
步骤S1:根据坐标信息的外键查询出坐标数据源;
步骤S2:把坐标数据源纵坐标和横坐标进行60进制转换为10进制坐标点;
在一个优选的实施例中,所述第二坐标转换操作为根据坐标信息的外键查询出坐标数据源;
在一个优选的实施例中,所述第三坐标转换操作包括以下步骤:
步骤P1:根据坐标信息的外键查询出坐标数据源;
步骤P2:判断坐标类型是直角坐标还是经纬坐标;若坐标类型为直角坐标,将坐标数据的源经度数据和维度数据进行调换;
在一个优选的实施例中,所述使用界址点图斑面积算法求解界址点图斑面积包括以下步骤:
步骤Z1:根据报盘导入坐标信息,根据传入的所属项目ID和坐标圈号查询出坐标数据源信息,以及坐标信息和分割方式对坐标数据源进行分组处理坐标信息;
步骤Z2:循环所有地块下面的坐标后根据多边形各个点坐标,计算出当前地块对应坐标的多边形的面积;
步骤Z3:根据圈号判断是主区域,还是挖空地块,圈号1为主区域、其他都为挖空地块;
步骤Z4:如果是主区域则多个坐标面积相加,其他都为挖空地块则多个坐标面积相减;
步骤Z5:将计算出的面积单位换算平方米转换成公顷;
在一个优选的实施例中,所述计算出当前地块对应坐标的多边形的面积包括以下步骤:
步骤X1:将面积标记为S,计算面积的初始值,计算公式为S=w1*(jN-j1);其中,w1为坐标集合的第一个坐标的维度,N为坐标集合中坐标的总量;j为坐标集合中每个坐标的经度;jN即为第N个坐标的经度;j1为第1个坐标的经度;
步骤X2:将面积S进行迭代更新;迭代更新的方式为:设置变量n=1;对于每个n=1,2,3,…N,将面积S依次更新为S+wn*(j(n-1)-j((n+1)%N));
步骤X3:将面积S更新为|S/2|,获得当前地块对应坐标的多边形的面积。
以上实施例仅用以说明本发明的技术方法而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方法进行修改或等同替换,而不脱离本发明技术方法的精神和范围。

Claims (8)

1.一种界址点图斑绘制及面积计算算法,其特征在于,包括以下步骤:
步骤一:报盘导入坐标信息,根据传入的坐标类型参数判断待生成图形的地块坐标类型;
步骤二:根据传入的坐标类型判断要生成的图形的地块坐标类型,并根据坐标类型进行对应的坐标转换操作;
步骤三:根据坐标信息中的分割坐标方式对坐标数据源进行分组,并对分组处理后的坐标数据根据实际经验进行缩放以及偏移处理;
步骤四:根据坐标信息中的坐标经纬度构造图形路径,遍历坐标数据源设置拐点号,并根据坐标画图,通过原点将多边形拆分成三角形,再通二维向量叉乘法求三角形面积;
步骤五:获得坐标信息中最小的经度值和纬度值对坐标进行首次偏移,保证坐标从图片的左上角开始,以及保证各个点之间的距离不变,待放大之后各点间的距离变大解决各个坐标在图形上显示的过于紧密的问题;获得坐标信息能放大的最大倍数并放大,将放大后的坐标偏移到图片的中央;获得界址点图斑会制图;
步骤六:使用界址点图斑面积算法求解界址点图斑面积。
2.根据权利要求1所述的一种界址点图斑绘制及面积计算算法,其特征在于,所述坐标信息包括项目ID、坐标圈号、生成图形的宽度、生成图形的宽度、坐标信息的外键、拐点号显示标记、坐标分割方式、坐标类型以及增列矿种标记;其中,坐标信息的外键用于查询导入的地块的坐标信息;坐标分割方式包括0与1两个值;其中,0代表*号分割方式,1代表圈分割方式;其中,坐标类型包括0、1、2、3、4五类,其中,0代表油气非油气探矿与油气采矿,1代表建筑用地,2代表非油气采矿,3代表增列矿种坐标专用,4代表城市边界规模调整shape文件图形预览。
3.根据权利要求1所述的一种界址点图斑绘制及面积计算算法,其特征在于,所述根据坐标类型进行对应的坐标转换操作包括以下方式:
如果坐标类型是油气非油气探矿与油气采矿,进行第一坐标转换操作;
如果坐标类型为建设用地,则进行第二坐标转换操作;
如果坐标类型是非油气采矿,则进行第三坐标转换操作;
如果坐标类型是增列矿种坐标专用,则进行第一坐标转换操作;
如果坐标类型是城市边界规模调整shape文件图形预览:则进行第二坐标转换操作。
4.根据权利要求3所述的一种界址点图斑绘制及面积计算算法,其特征在于,所述第一坐标转换操作包括以下步骤:
步骤S1:根据坐标信息的外键查询出坐标数据源;
步骤S2:把坐标数据源纵坐标和横坐标进行60进制转换为10进制坐标点。
5.根据权利要求3所述的一种界址点图斑绘制及面积计算算法,其特征在于,所述第二坐标转换操作为根据坐标信息的外键查询出坐标数据源。
6.根据权利要求3所述的一种界址点图斑绘制及面积计算算法,其特征在于,所述第三坐标转换操作包括以下步骤:
步骤P1:根据坐标信息的外键查询出坐标数据源;
步骤P2:判断坐标类型是直角坐标还是经纬坐标;若坐标类型为直角坐标,将坐标数据的源经度数据和维度数据进行调换。
7.根据权利要求1所述的一种界址点图斑绘制及面积计算算法,其特征在于,使用界址点图斑面积算法求解界址点图斑面积包括以下步骤:
步骤Z1:根据报盘导入坐标信息,根据传入的所属项目ID和坐标圈号查询出坐标数据源信息,以及坐标信息和分割方式对坐标数据源进行分组处理坐标信息;
步骤Z2:循环所有地块下面的坐标后根据多边形各个点坐标,计算出当前地块对应坐标的多边形的面积;
步骤Z3:根据圈号判断是主区域,还是挖空地块,圈号1为主区域、其他都为挖空地块;
步骤Z4:如果是主区域则多个坐标面积相加,其他都为挖空地块则多个坐标面积相减;
步骤Z5:将计算出的面积单位换算平方米转换成公顷。
8.根据权利要求1所述的一种界址点图斑绘制及面积计算算法,其特征在于,计算出当前地块对应坐标的多边形的面积包括以下步骤:
步骤X1:将面积标记为S,计算面积的初始值,计算公式为S=w1*(jN-j 1);其中,w1为坐标集合的第一个坐标的维度,N为坐标集合中坐标的总量;j为坐标集合中每个坐标的经度;jN即为第N个坐标的经度;j 1为第1个坐标的经度;
步骤X2:将面积S进行迭代更新;迭代更新的方式为:设置变量n=1;对于每个n=1,2,3,…N,将面积S依次更新为S+wn*(j(n-1)-j((n+1)%N));
步骤X3:将面积S更新为|S/2|,获得当前地块对应坐标的多边形的面积。
CN202211477418.2A 2022-11-23 2022-11-23 一种界址点图斑绘制及面积计算算法 Active CN115731284B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211477418.2A CN115731284B (zh) 2022-11-23 2022-11-23 一种界址点图斑绘制及面积计算算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211477418.2A CN115731284B (zh) 2022-11-23 2022-11-23 一种界址点图斑绘制及面积计算算法

Publications (2)

Publication Number Publication Date
CN115731284A CN115731284A (zh) 2023-03-03
CN115731284B true CN115731284B (zh) 2024-06-25

Family

ID=85297732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211477418.2A Active CN115731284B (zh) 2022-11-23 2022-11-23 一种界址点图斑绘制及面积计算算法

Country Status (1)

Country Link
CN (1) CN115731284B (zh)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111062658B (zh) * 2019-12-04 2022-11-18 南京苏宁加电子商务有限公司 配送范围的生成方法、装置、计算机设备和存储介质
CN112131329A (zh) * 2020-09-11 2020-12-25 北京数字政通科技股份有限公司 基于椭球面梯形剖分的国土空间面积计算方法
CN113706604B (zh) * 2021-08-20 2024-02-09 园测信息科技股份有限公司 一种基于两个凸多边形交集求解算法的地类图斑分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
国家水土保持重点工程项目管理系统拐点坐标输入方法――AUTOCAD基点坐标转换经纬度;王连磊;陈新军;;中国水土保持;20180805(第08期);全文 *
在AutoCAD中实现面积的快速量测;龚有亮, 曹歆宏, 翟翊;测绘通报;20020425(第04期);全文 *

Also Published As

Publication number Publication date
CN115731284A (zh) 2023-03-03

Similar Documents

Publication Publication Date Title
CN113724279B (zh) 路网自动划分交通小区的系统、方法、设备及存储介质
CN115018249B (zh) 一种基于激光扫描技术的地铁站施工质量评价方法
CN117557681B (zh) 一种基于多源测绘数据的高精度地形图生成方法及装置
CN115731284B (zh) 一种界址点图斑绘制及面积计算算法
CN113486429B (zh) 一种基于插值算法的空间汇交结构自动化建造方法
CN117237567B (zh) 一种数字孪生洪水流场可视化渲染方法
CN112525208B (zh) 城市道路地图快速构建方法和装置及设备
CN112069445A (zh) 一种2d slam算法评估及量化方法
CN111475592A (zh) 一种面向矢量图形文件的流域地图制作方法
CN114239096B (zh) 一种基于ArcGIS的土壤侵蚀类型图自动绘制模型的建立方法
CN115906717A (zh) 一种水土保持的侵蚀计算方法及系统
CN107220615B (zh) 一种融合兴趣点大数据的城市不透水面信息提取方法
CN109509234B (zh) 一种基于城区建筑群高清栅格地图矢量化的方法
JPWO2010092680A1 (ja) 地図情報処理装置、地図情報処理方法、地図情報処理プログラム、及び記録媒体
CN109493419B (zh) 一种倾斜摄影数据获取数字表面模型的方法及装置
CN112487124A (zh) 一种使用VBA将CorelDraw地图中点状要素转换到SuperMap的方法
CN112434081A (zh) 地下管线数据检查方法、装置、电子设备及存储介质
CN111241221A (zh) 一种受损地形坐标数据的自动化匹配及高精度修复方法
CN116012296B (zh) 基于超分辨率和半监督学习的装配式预制构件检测方法
Jin et al. Analysis of the Modeling Method and Application of 3D City Model based on the CityEngine
CN111400891B (zh) 管道巡检点偏差程度获取方法、装置、设备及存储介质
CN118230192B (zh) 基于实例分割模型耕地影像中耕地图斑共线的处理方法
CN115328373B (zh) 一种确定性系数计算方法及系统
CN112464153A (zh) 一种土地性质变化分析方法
Nikolić IS GIS PRESENT IN LOGISTICS 4.0 AND 5.0? EXPLORING THE POTENTIAL OF USING GIS IN SMART LOGISTICS

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