CN109408656A - 基于MapGIS的地球化学调查矢量数据自动化处理方法及系统 - Google Patents

基于MapGIS的地球化学调查矢量数据自动化处理方法及系统 Download PDF

Info

Publication number
CN109408656A
CN109408656A CN201811133259.8A CN201811133259A CN109408656A CN 109408656 A CN109408656 A CN 109408656A CN 201811133259 A CN201811133259 A CN 201811133259A CN 109408656 A CN109408656 A CN 109408656A
Authority
CN
China
Prior art keywords
point
layer
area
vector
band
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
CN201811133259.8A
Other languages
English (en)
Other versions
CN109408656B (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.)
China University of Geosciences
Original Assignee
China University of Geosciences
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 China University of Geosciences filed Critical China University of Geosciences
Priority to CN201811133259.8A priority Critical patent/CN109408656B/zh
Publication of CN109408656A publication Critical patent/CN109408656A/zh
Application granted granted Critical
Publication of CN109408656B publication Critical patent/CN109408656B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Processing Or Creating Images (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于MapGIS的地球化学调查矢量数据自动化处理方法及系统,首先设置待数据处理区域地球化学调查矢量数据原始点图层、各分带提取框图层、结果矢量区图层,然后通过程序依次使用各分带提取框矢量图层对原始点图层进行裁剪、投影变换和点转区自动化处理,获得最终的待数据处理区域地球化学调查矢量区图层。本发明最明显的不同在于能够实现地球化学调查矢量数据点转区整个过程的自动化,解决了目前地球化学调查矢量数据点转区需要人工处理,工作效率低的问题,本发明可以为地球化学调查成果图件的实时网络共享和发布节约大量的人力成本,提高成果数据发布的时效性。

Description

基于MapGIS的地球化学调查矢量数据自动化处理方法及系统
技术领域
本发明涉及矢量数据自动化处理技术,尤其是一种基于MapGIS的地球化学调查矢量数据自动化处理方法及系统。
背景技术
地球化学调查成果数据制图是成果共享和网络发布的必经过程,而网络发布需要对成果数据进行缩放,而点要素在高级别显示时,无法正确表达地球化学调查每个工作点的真正意义,需要将点要素转换成区要素,并且相邻的两个区还必须是无缝衔接的,中间不能有缝隙。要正确表达地球化学调查的意义,每个区要必须代表2平方公里的空间范围,地球化学调查采样点一般是在平面坐标系中进行的,因此,可以对不同的区域分带进行处理,将每个带内的点要素转换成2公里*2公里的正方形区要素,从而实现整个待数据处理区域的地球化学调查数据的自动化处理。
发明内容
本发明要解决的技术问题在于,针对现有技术中需要大量人工处理工作的缺陷,提供了一种基于MapGIS的地球化学调查矢量数据自动化处理方法及系统。
根据本发明的其中一方面,本发明解决其技术问题所采用的技术方案是:提供一种基于MapGIS的地球化学调查矢量数据自动化处理方法,包括如下步骤:
步骤1:设置原始矢量点图层路径、k1带提取框矢量图层路径、k2带提取框矢量图层路径、...、kn带提取框矢量图层路径以及结果矢量区图层保存路径;其中,k1带、k2带、...、kn带为原始矢量点图层所对应的区域所处的所有的6度带,n为大于1的正整数;
步骤2:根据所述原始矢量点图层路径读取原始矢量点图层,并获取所读取的原始矢量点图层的空间参照系和属性结构;根据所述结果矢量区图层保存路径创建结果矢量区图层,并将所创建的结果矢量区图层的空间参照系、属性结构与原始点图层的空间参照系和属性结构设置为一致;
步骤3:k1带处理步骤,具体包括如下子步骤:
S31、使用k1带提取框矢量图层裁剪原始矢量点图层,得到分布在k1带内的原始点数据,将结果数据保存在临时图层TempOriginLayer_k1中;
S32、将TempOriginLayer_k1投影转换到“高斯大地坐标系_西安80_k1带6_北”空间参照系,生成临时图层TempTransLayer_k1;
S33、获取TempTransLayer_k1的地图范围,得到地图范围的左下角的坐标值(xmin,ymin),将其向下取整,得到取整后的新坐标值(XMIN,YMIN);
S34、查询TempTransLayer_k1中的所有点要素,然后遍历每一个点要素,获取每个点要素的所有属性值、图形颜色,根据每个点要素的坐标,计算包含这个点要素的一个正方形区要素的五个点的坐标,第一个点是正方形的起始左下角点(x0,y0),第二个点是右下角点(x1,y1),第三个点是右上角点(x2,y2),第四个点是左上角点(x3,y3),第五个点是尾部左下角点(x4,y4),计算规则是把点要素当作正方形的中心,正方形的宽度和高度都是2000,推算正方形的所述5个点的坐标;
S35、创建临时结果区图层TempRectLayer_k1,将TempRectLayer_k1的空间参照系设置为“高斯大地坐标系_西安80_k1带6_北”,并设置TempRectLayer_k1的属性结构与TempTransLayer_k1相同;
S36、将计算得到的5个坐标点构成正方形区要素,设置正方形区要素的属性、颜色与当前处理的点要素的属性和颜色相同,然后将设置后的正方形区要素添加到临时结果区图层TempRectLayer_k1中;
S37、将临时结果区图层TempRectLayer_k1投影变换到原始点图层空间参考系SRS_ORIGIN下,生成与原始点图层空间参照系一致的临时结果区图层TempRectTransLayer_k1,然后调用图层追加接口,将TempRectTransLayer_k1中的所有区要素追加到步骤2创建的结果矢量区图层中;
步骤4:采用与步骤3相同的方法,分别使用k2、...、Kn带提取框矢量图层裁剪原始矢量点图层并进行处理,将各个带处理后的区要素都添加到结果矢量区图层中作为最终的处理结果。
进一步地,在本发明的基于MapGIS的地球化学调查矢量数据自动化处理方法中,还包括步骤:
步骤5:获取原始矢量点图层和结果矢量区图层中的点要素及区要素个数,判断两个图层中的各要素的个数是否相等,如果相等,则判定数据处理成功,如果不相等,判定数据处理失败。
进一步地,在本发明的基于MapGIS的地球化学调查矢量数据自动化处理方法,步骤3中,还包括步骤S37之后的步骤:
S38、删除临时图层TempOriginLayer_k1、TempTransLayer_k1、TempRectLayer_k1、TempRectTransLayer_k1。
进一步地,在本发明的基于MapGIS的地球化学调查矢量数据自动化处理方法中,步骤3中,正方形的所述5个点的坐标的计算方法如下:
A1、将当前处理的点要素的坐标向下取整后的位置(px,py)进行校正,得到校正后的点位置(px’,py’),校正计算公式如下:
px’=Math.round((px-XMIN)/2000)*2000+XMIN,
py’=Math.round((py-YMIN)/2000)*2000+YMIN,
其中,Math.round()是四舍五入取整函数;
A2、以(px’,py’)为正方形中心点,计算边长为2000的正方形的所述5
个点的坐标,计算公式如下:
x0=px’-1000;
y0=py’-1000;
x1=px’+1000;
y1=py’-1000;
x2=px’+1000;
y2=py’+1000;
x3=px’-1000;
y3=py’+1000;
x4=px’-1000;
y4=py’-1000。
根据本发明的另一方面,本发明为解决其技术问题,还提供了一种基于MapGIS的地球化学调查矢量数据自动化处理系统,包括如下模块:
路径设置模块,用于设置原始矢量点图层路径、k1带提取框矢量图层路径、k2带提取框矢量图层路径、...、kn带提取框矢量图层路径以及结果矢量区图层保存路径;其中,k1带、k2带、...、kn带为原始矢量点图层所对应的区域所处的所有的6度带,n为大于1的正整数;
预处理模块,用于根据所述原始矢量点图层路径读取原始矢量点图层,并获取所读取的原始矢量点图层的空间参照系和属性结构;根据所述结果矢量区图层保存路径创建结果矢量区图层,并将所创建的结果矢量区图层的空间参照系、属性结构与原始点图层的空间参照系和属性结构设置为一致;
k1带处理模块,具体包括如下单元:
原始点数据获取单元,用于使用k1带提取框矢量图层裁剪原始矢量点图层,得到分布在k1带内的原始点数据,将结果数据保存在临时图层TempOriginLayer_k1中;
投影转换单元,用于将TempOriginLayer_k1投影转换到“高斯大地坐标系_西安80_k1带6_北”空间参照系,生成临时图层TempTransLayer_k1;
坐标预处理单元,用于获取TempTransLayer_k1的地图范围,得到地图范围的左下角的坐标值(xmin,ymin),将其向下取整,得到取整后的新坐标值(XMIN,YMIN);
点要素处理单元,用于查询TempTransLayer_k1中的所有点要素,然后遍历每一个点要素,获取每个点要素的所有属性值、图形颜色,根据每个点要素的坐标,计算包含这个点要素的一个正方形区要素的五个点的坐标,第一个点是正方形的起始左下角点(x0,y0),第二个点是右下角点(x1,y1),第三个点是右上角点(x2,y2),第四个点是左上角点(x3,y3),第五个点是尾部左下角点(x4,y4),计算规则是把点要素当作正方形的中心,正方形的宽度和高度都是2000,推算正方形的所述5个点的坐标;
临时结果区图层处理单元,用于创建临时结果区图层TempRectLayer_k1,将TempRectLayer_k1的空间参照系设置为“高斯大地坐标系_西安80_k1带6_北”,并设置TempRectLayer_k1的属性结构与TempTransLayer_k1相同;
区要素添加单元,用于将计算得到的5个坐标点构成正方形区要素,设置正方形区要素的属性、颜色与当前处理的点要素的属性和颜色相同,然后将设置后的正方形区要素添加到临时结果区图层TempRectLayer_k1中;
区要素追加单元,用于将临时结果区图层TempRectLayer_k1投影变换到原始点图层空间参考系SRS_ORIGIN下,生成与原始点图层空间参照系一致的临时结果区图层TempRectTransLayer_k1,然后调用图层追加接口,将TempRectTransLayer_k1中的所有区要素追加到预处理模块创建的结果矢量区图层中;
最终结果获取模块,用于:采用与k1带处理模块相同的方法,分别使用k2、...、Kn带提取框矢量图层裁剪原始矢量点图层并进行处理,将各个带处理后的区要素都添加到结果矢量区图层中作为最终的处理结果。
进一步地,在本发明的基于MapGIS的地球化学调查矢量数据自动化处理系统,还包括:
结果判断模块,用于获取原始矢量点图层和结果矢量区图层中的点要素及区要素个数,判断两个图层中的各要素的个数是否相等,如果相等,则判定数据处理成功,如果不相等,判定数据处理失败。
进一步地,在本发明的基于MapGIS的地球化学调查矢量数据自动化处理系统,k1带处理模块中,还包括区要素追加单元之后被调用的如下单元:
临时图层删除单元,用于删除临时图层TempOriginLayer_k1、TempTransLayer_k1、TempRectLayer_k1、TempRectTransLayer_k1。
进一步地,在本发明的基于MapGIS的地球化学调查矢量数据自动化处理系统,k1带处理模块中,采用如下单元进行正方形的所述5个点的坐标的计算:
坐标校正单元,用于将当前处理的点要素的坐标向下取整后的位置(px,py)
进行校正,得到校正后的点位置(px’,py’),校正计算公式如下:
px’=Math.round((px-XMIN)/2000)*2000+XMIN,
py’=Math.round((py-YMIN)/2000)*2000+YMIN,
其中,Math.round()是四舍五入取整函数;
坐标计算单元,用于以(px’,py’)为正方形中心点,计算边长为2000的正方形的所述5个点的坐标,计算公式如下:
x0=px’-1000;
y0=py’-1000;
x1=px’+1000;
y1=py’-1000;
x2=px’+1000;
y2=py’+1000;
x3=px’-1000;
y3=py’+1000;
x4=px’-1000;
y4=py’-1000。
实施本发明的基于MapGIS的地球化学调查矢量数据自动化处理方法及系统,能够对待数据处理区域的地球化学调查成果点图层进行分带处理,自动裁剪、投影、转换每个带内的原始点数据,获得最终的地球化学调查成果区图层。通过此方法能够实现点数据,获得最终的地球化学调查成果区图层。通过此方法能够实现地球化学调查成果点图层的全自动化处理,不再需要人工使用GIS软件进行繁琐的处理,极大地提高了地球化学调查成果数据的处理效率,尤其是对大批量的地球化学调查成果数据进行处理时,本发明的效率优势更为明显。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明的基于MapGIS的地球化学调查矢量数据自动化处理方法一实施例的流程图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
下述将以西北五省为例进行说明,其中西北五省具体是指:陕西省、甘肃省、青海省、宁夏回族自治区、新疆维吾尔自治区,西北五省分布在13至19度带。如图1所示,图1是本发明提供的一种西北五省地球化学调查矢量数据自动化处理方法流程图,运行于MapGIS软件为例,所述方法包括以下步骤:
步骤1:设置原始矢量点图层路径、13带提取框矢量图层路径、14带提取框矢量图层路径、15带提取框矢量图层路径、16带提取框矢量图层路径、17带提取框矢量图层路径、18带提取框矢量图层路径、19带提取框矢量图层路径、结果矢量区图层保存路径。
步骤2:根据所述原始矢量点图层路径,通过SFeatureCls.SrID属性和SFeatureCls.Fields,获取原始矢量点图层的空间参照系SRS_ORIGIN和属性结构;根据所述结果矢量区图层保存路径,调用SFeatureCls.Create函数创建结果矢量区图层,设置所创建得结果矢量区图层空间参照系、属性结构与原始点图层空间参照系和属性结构一致。
步骤3:调用SFeatureCls.Create函数创建临时图层TempOriginLayer_13,调用SpatialAnalysis.Clip函数,使用13带提取框矢量图层裁剪原始矢量点图层,得到分布在13带内的原始点数据,将结果数据保存在临时图层TempOriginLayer_13中。
调用SFeatureCls.Create函数创建临时图层TempTransLayer_13,然后调用SFeatureCls.ProjTrans函数将TempOriginLayer_13投影转换到“高斯大地坐标系_西安80_13带6_北”空间参照系,结果保存到临时图层TempTransLayer_13。其中,高斯大地坐标系_西安80_13带6_北为MapGIS软件里面的空间参照系的名字。
通过SFeatureCls.Range属性,获取TempTransLayer_13的地图范围,得到地图范围的左下角的坐标值(xmin,ymin),将其向下取整,得到取整后的新坐标值(XMIN,YMIN)。
调用SFeatureCls.Select函数,查询TempTransLayer_13中的所有点要素,然后遍历每一个点,调用SFeatureCls.GetGeometry函数获取每个点的坐标信息,调用SFeatureCls.GetAtt函数获取点的属性值,调用SFeatureCls.GetInfo函数获取点的图形颜色信息,根据每个点的坐标,计算包含这个点的一个正方形区要素的五个点的坐标,第一个点是正方形的左下角(x0,y0),第二个点是右下角(x1,y1),第三个点是右上角(x2,y2),第四个点是左上角(x3,y3),第五个点是左下角(x4,y4)。计算规则是把点要素当作正方形的中心,正方形的宽度和高度都是2000,推算正方形的各个角的坐标。这里的2000是地球化学调查工作的标准规范,它的单位就是上面提到的高斯平面坐标系下的坐标单位。
以左下角为例,设当前处理的点要素的坐标向下取整后为(px,py),首先要将其位置进行校正,得到校正后的点位置(px’,py’),校正计算公式如下(Math.round是四舍五入取整函数):
px’=Math.round((px-XMIN)/2000)*2000+XMIN
py’=Math.round((py-YMIN)/2000)*2000+YMIN
再以(px’,py’)为正方形中心点,计算边长为2000的正方形的左下角坐标,计算公式如下:
x0=px’-1000
y0=py’-1000
其它几个点的计算方法如下:
x1=px’+1000
y1=py’-1000
x2=px’+1000
y2=py’+1000
x3=px’-1000
y3=py’+1000
x4=px’-1000
y4=py’-1000
调用SFeatureCls.Create函数,创建临时结果区图层TempRectLayer_13,将其空间参照系设置为“高斯大地坐标系_西安80_13带6_北”,并设置其属性结构与TempTransLayer_13相同。
将计算得到的5个坐标点构成正方形区要素,设置其属性、颜色与当前处理的点要素的属性和颜色相同,然后调用SFeatureCls.Append函数将其添加到临时结果区图层TempRectLayer_13中。
调用SFeatureCls.ProjTrans函数,将临时结果区图层TempRectLayer_13投影变换到原始点图层空间参考系SRS_ORIGIN下,生成与原始点图层空间参照系一致的临时结果区图层TempRectTransLayer_13,然后调用图层追加接口,将TempRectTransLayer_13中的所有区要素追加到步骤2创建的结果矢量区图层中。
删除临时图层TempOriginLayer_13、TempTransLayer_13、TempRectLayer_13、TempRectTransLayer_13。
步骤4:与步骤3类似,依次使用14、15、16、17、18、19带提取框矢量图层裁剪原始矢量点图层并进行处理,将各个带处理后的区要素都追加到结果矢量区图层中。在处理过程中涉及到使用高斯大地坐标系空间参考系,每个带需要分别使用相对应的空间参考系,例如,当处理17带时,使用的空间参照系是“高斯大地坐标系_西安80_17带6_北”。
步骤5:通过SFeatureCls.Count属性获取原始矢量点图层和结果矢量区图层中的要素个数,判断两个图层中的要素个数是否相等,如果相等提示“处理成功”,如果不相等,提示“处理失败”,结束。
本发明最明显的不同在于能够将西北五省地球化学调查矢量点图层进行分带处理,针对每个带中的点要素,能够自动生成与之对应的规则的2公里*2公里的网格区要素,通过程序自动化完成整个西北五省所有地球化学调查矢量点图层的处理工作,节省了大量的人工数据处理工作量,极大地提高了数据处理效率,为西北五省所有地球化学调查成果图件的实时网格共享和发布提供了强有力的技术支撑。尤其是对大批量的地球化学调查矢量点图层进行处理时,本发明的优势更为明显。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。

Claims (8)

1.一种基于MapGIS的地球化学调查矢量数据自动化处理方法,其特征在于,包括如下步骤:
步骤1:设置原始矢量点图层路径、k1带提取框矢量图层路径、k2带提取框矢量图层路径、...、kn带提取框矢量图层路径以及结果矢量区图层保存路径;其中,k1带、k2带、...、kn带为原始矢量点图层所对应的区域所处的所有的6度带,n为大于1的正整数;
步骤2:根据所述原始矢量点图层路径读取原始矢量点图层,并获取所读取的原始矢量点图层的空间参照系和属性结构;根据所述结果矢量区图层保存路径创建结果矢量区图层,并将所创建的结果矢量区图层的空间参照系、属性结构与原始点图层的空间参照系和属性结构设置为一致;
步骤3:k1带处理步骤,具体包括如下子步骤:
S31、使用k1带提取框矢量图层裁剪原始矢量点图层,得到分布在k1带内的原始点数据,将结果数据保存在临时图层TempOriginLayer_k1中;
S32、将TempOriginLayer_k1投影转换到“高斯大地坐标系_西安80_k1带6_北”空间参照系,生成临时图层TempTransLayer_k1;
S33、获取TempTransLayer_k1的地图范围,得到地图范围的左下角的坐标值(xmin,ymin),将其向下取整,得到取整后的新坐标值(XMIN,YMIN);
S34、查询TempTransLayer_k1中的所有点要素,然后遍历每一个点要素,获取每个点要素的所有属性值、图形颜色,根据每个点要素的坐标,计算包含这个点要素的一个正方形区要素的五个点的坐标,第一个点是正方形的起始左下角点(x0,y0),第二个点是右下角点(x1,y1),第三个点是右上角点(x2,y2),第四个点是左上角点(x3,y3),第五个点是尾部左下角点(x4,y4),计算规则是把点要素当作正方形的中心,正方形的宽度和高度都是2000,推算正方形的所述5个点的坐标;
S35、创建临时结果区图层TempRectLayer_k1,将TempRectLayer_k1的空间参照系设置为“高斯大地坐标系_西安80_k1带6_北”,并设置TempRectLayer_k1的属性结构与TempTransLayer_k1相同;
S36、将计算得到的5个坐标点构成正方形区要素,设置正方形区要素的属性、颜色与当前处理的点要素的属性和颜色相同,然后将设置后的正方形区要素添加到临时结果区图层TempRectLayer_k1中;
S37、将临时结果区图层TempRectLayer_k1投影变换到原始点图层空间参考系SRS_ORIGIN下,生成与原始点图层空间参照系一致的临时结果区图层TempRectTransLayer_k1,然后调用图层追加接口,将TempRectTransLayer_k1中的所有区要素追加到步骤2创建的结果矢量区图层中;
步骤4:采用与步骤3相同的方法,分别使用k2、...、Kn带提取框矢量图层裁剪原始矢量点图层并进行处理,将各个带处理后的区要素都添加到结果矢量区图层中作为最终的处理结果。
2.根据权利要求1所述的基于MapGIS的地球化学调查矢量数据自动化处理方法,其特征在于,还包括步骤:
步骤5:获取原始矢量点图层和结果矢量区图层中的点要素及区要素个数,判断两个图层中的各要素的个数是否相等,如果相等,则判定数据处理成功,如果不相等,判定数据处理失败。
3.根据权利要求1所述的基于MapGIS的地球化学调查矢量数据自动化处理方法,其特征在于,步骤3中,还包括步骤S37之后的步骤:
S38、删除临时图层TempOriginLayer_k1、TempTransLayer_k1、TempRectLayer_k1、TempRectTransLayer_k1。
4.根据权利要求1所述的基于MapGIS的地球化学调查矢量数据自动化处理方法,其特征在于,步骤3中,正方形的所述5个点的坐标的计算方法如下:
A1、将当前处理的点要素的坐标向下取整后的位置(px,py)进行校正,得到校正后的点位置(px’,py’),校正计算公式如下:
px’=Math.round((px-XMIN)/2000)*2000+XMIN,
py’=Math.round((py-YMIN)/2000)*2000+YMIN,
其中,Math.round()是四舍五入取整函数;
A2、以(px’,py’)为正方形中心点,计算边长为2000的正方形的所述5
个点的坐标,计算公式如下:
x0=px’-1000;
y0=py’-1000;
x1=px’+1000;
y1=py’-1000;
x2=px’+1000;
y2=py’+1000;
x3=px’-1000;
y3=py’+1000;
x4=px’-1000;
y4=py’-1000。
5.一种基于MapGIS的地球化学调查矢量数据自动化处理系统,其特征在于,包括如下模块:
路径设置模块,用于设置原始矢量点图层路径、k1带提取框矢量图层路径、k2带提取框矢量图层路径、...、kn带提取框矢量图层路径以及结果矢量区图层保存路径;其中,k1带、k2带、...、kn带为原始矢量点图层所对应的区域所处的所有的6度带,n为大于1的正整数;
预处理模块,用于根据所述原始矢量点图层路径读取原始矢量点图层,并获取所读取的原始矢量点图层的空间参照系和属性结构;根据所述结果矢量区图层保存路径创建结果矢量区图层,并将所创建的结果矢量区图层的空间参照系、属性结构与原始点图层的空间参照系和属性结构设置为一致;
k1带处理模块,具体包括如下单元:
原始点数据获取单元,用于使用k1带提取框矢量图层裁剪原始矢量点图层,得到分布在k1带内的原始点数据,将结果数据保存在临时图层TempOriginLayer_k1中;
投影转换单元,用于将TempOriginLayer_k1投影转换到“高斯大地坐标系_西安80_k1带6_北”空间参照系,生成临时图层TempTransLayer_k1;
坐标预处理单元,用于获取TempTransLayer_k1的地图范围,得到地图范围的左下角的坐标值(xmin,ymin),将其向下取整,得到取整后的新坐标值(XMIN,YMIN);
点要素处理单元,用于查询TempTransLayer_k1中的所有点要素,然后遍历每一个点要素,获取每个点要素的所有属性值、图形颜色,根据每个点要素的坐标,计算包含这个点要素的一个正方形区要素的五个点的坐标,第一个点是正方形的起始左下角点(x0,y0),第二个点是右下角点(x1,y1),第三个点是右上角点(x2,y2),第四个点是左上角点(x3,y3),第五个点是尾部左下角点(x4,y4),计算规则是把点要素当作正方形的中心,正方形的宽度和高度都是2000,推算正方形的所述5个点的坐标;
临时结果区图层处理单元,用于创建临时结果区图层TempRectLayer_k1,将TempRectLayer_k1的空间参照系设置为“高斯大地坐标系_西安80_k1带6_北”,并设置TempRectLayer_k1的属性结构与TempTransLayer_k1相同;
区要素添加单元,用于将计算得到的5个坐标点构成正方形区要素,设置正方形区要素的属性、颜色与当前处理的点要素的属性和颜色相同,然后将设置后的正方形区要素添加到临时结果区图层TempRectLayer_k1中;
区要素追加单元,用于将临时结果区图层TempRectLayer_k1投影变换到原始点图层空间参考系SRS_ORIGIN下,生成与原始点图层空间参照系一致的临时结果区图层TempRectTransLayer_k1,然后调用图层追加接口,将TempRectTransLayer_k1中的所有区要素追加到预处理模块创建的结果矢量区图层中;
最终结果获取模块,用于:采用与k1带处理模块相同的方法,分别使用k2、...、Kn带提取框矢量图层裁剪原始矢量点图层并进行处理,将各个带处理后的区要素都添加到结果矢量区图层中作为最终的处理结果。
6.根据权利要求5所述的基于MapGIS的地球化学调查矢量数据自动化处理系统,其特征在于,还包括:
结果判断模块,用于获取原始矢量点图层和结果矢量区图层中的点要素及区要素个数,判断两个图层中的各要素的个数是否相等,如果相等,则判定数据处理成功,如果不相等,判定数据处理失败。
7.根据权利要求5所述的基于MapGIS的地球化学调查矢量数据自动化处理系统,其特征在于,k1带处理模块中,还包括区要素追加单元之后被调用的如下单元:
临时图层删除单元,用于删除临时图层TempOriginLayer_k1、TempTransLayer_k1、TempRectLayer_k1、TempRectTransLayer_k1。
8.根据权利要求5所述的基于MapGIS的地球化学调查矢量数据自动化处理系统,其特征在于,k1带处理模块中,采用如下单元进行正方形的所述5个点的坐标的计算:
坐标校正单元,用于将当前处理的点要素的坐标向下取整后的位置(px,py)进行校正,得到校正后的点位置(px’,py’),校正计算公式如下:
px’=Math.round((px-XMIN)/2000)*2000+XMIN,
py’=Math.round((py-YMIN)/2000)*2000+YMIN,
其中,Math.round()是四舍五入取整函数;
坐标计算单元,用于以(px’,py’)为正方形中心点,计算边长为2000的正方形的所述5个点的坐标,计算公式如下:
x0=px’-1000;
y0=py’-1000;
x1=px’+1000;
y1=py’-1000;
x2=px’+1000;
y2=py’+1000;
x3=px’-1000;
y3=py’+1000;
x4=px’-1000;
y4=py’-1000。
CN201811133259.8A 2018-09-27 2018-09-27 基于MapGIS的地球化学调查矢量数据自动化处理方法及系统 Active CN109408656B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811133259.8A CN109408656B (zh) 2018-09-27 2018-09-27 基于MapGIS的地球化学调查矢量数据自动化处理方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811133259.8A CN109408656B (zh) 2018-09-27 2018-09-27 基于MapGIS的地球化学调查矢量数据自动化处理方法及系统

Publications (2)

Publication Number Publication Date
CN109408656A true CN109408656A (zh) 2019-03-01
CN109408656B CN109408656B (zh) 2020-08-07

Family

ID=65465296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811133259.8A Active CN109408656B (zh) 2018-09-27 2018-09-27 基于MapGIS的地球化学调查矢量数据自动化处理方法及系统

Country Status (1)

Country Link
CN (1) CN109408656B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110570385A (zh) * 2019-08-14 2019-12-13 苏州中科天启遥感科技有限公司 基于遥感解译的边界套合方法、系统、存储介质及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103970905A (zh) * 2014-05-27 2014-08-06 重庆大学 一种多源矢量地理信息数据的自动构图整合方法
CN107169080A (zh) * 2017-05-10 2017-09-15 合肥慧图软件有限公司 一种基于gis和空间数据库相结合的地理空间分析系统
CN107391753A (zh) * 2017-08-17 2017-11-24 烟台市公路管理局 一种基于gis的路产矢量化数据自动生成系统及方法
CN108009214A (zh) * 2017-11-15 2018-05-08 中国地质大学(武汉) 基于ArcGIS的矢量数据并行处理方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103970905A (zh) * 2014-05-27 2014-08-06 重庆大学 一种多源矢量地理信息数据的自动构图整合方法
CN107169080A (zh) * 2017-05-10 2017-09-15 合肥慧图软件有限公司 一种基于gis和空间数据库相结合的地理空间分析系统
CN107391753A (zh) * 2017-08-17 2017-11-24 烟台市公路管理局 一种基于gis的路产矢量化数据自动生成系统及方法
CN108009214A (zh) * 2017-11-15 2018-05-08 中国地质大学(武汉) 基于ArcGIS的矢量数据并行处理方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孟海东,张肖云,杨永杰: "基于MapGIS的地球化学数据处理方法", 《现代矿业》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110570385A (zh) * 2019-08-14 2019-12-13 苏州中科天启遥感科技有限公司 基于遥感解译的边界套合方法、系统、存储介质及设备
CN110570385B (zh) * 2019-08-14 2023-05-12 苏州中科天启遥感科技有限公司 基于遥感解译的边界套合方法、系统、存储介质及设备

Also Published As

Publication number Publication date
CN109408656B (zh) 2020-08-07

Similar Documents

Publication Publication Date Title
US11222465B2 (en) Embedded urban design scene emulation method and system
EP2545532B1 (en) Integrated gis system with interactive 3d interface
CN103761397B (zh) 用于面曝光增材成型的3d模型切片及投影面生成方法
CN111737790B (zh) 仿真城市模型的构建方法和设备
CN102855810A (zh) 一种基于利用卫星影像图进行数字线划地图的方法
CN109154503A (zh) 无人机作业航线的规划方法及地面端设备
CN110570491B (zh) 热力图生成方法及装置
CN108648271A (zh) 一种基于gis数据生成复杂地形网格模型的插值方法
CN112686877A (zh) 基于双目相机的三维房屋损伤模型构建测量方法及系统
CN108957507A (zh) 基于增强现实技术的燃气管道泄漏处置方法
JP2014186570A (ja) 3次元地図表示システム
CN106683183B (zh) 一种建设工程模型在线三维可视化方法
CN115270267A (zh) 基于混合hbim模型的智能建筑初级数字孪生系统及构建方法
CN109408656A (zh) 基于MapGIS的地球化学调查矢量数据自动化处理方法及系统
CN113269883B (zh) 基于二维剖面和catia的bim三维地质建模方法
KR100732915B1 (ko) 디지털사진 측량기술 및 인공위성영상을 이용한 기본설계용도로노선의 3차원적 결정 방법
Heipke et al. Semiautomatic extraction of roads from aerial images
CN104574489A (zh) 基于叠层四叉树图集的地形与动态矢量集成方法
CN109657023B (zh) 景区动态电子地图绘制方法、服务器及系统
CN115131511A (zh) 一种基于倾斜摄影技术+Dynamo创建地形的方法
CN112800514B (zh) 一种采用激光点云加bim建模技术运用到换流站可视化管控平台的方法
CN103544675B (zh) 地震波传播展示方法
Zhang et al. A visualization progress management approach of bridge construction based on mixed panoramic and oblique photogrammetry
CN109191556B (zh) 一种从lod分页表面纹理模型提取栅格化数字高程模型的方法
CN117689832B (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