CN109191579B - 坐标转换的方法、装置、计算机设备和存储介质 - Google Patents

坐标转换的方法、装置、计算机设备和存储介质 Download PDF

Info

Publication number
CN109191579B
CN109191579B CN201810706859.2A CN201810706859A CN109191579B CN 109191579 B CN109191579 B CN 109191579B CN 201810706859 A CN201810706859 A CN 201810706859A CN 109191579 B CN109191579 B CN 109191579B
Authority
CN
China
Prior art keywords
coordinate
grid
reduction
point
coordinate system
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
CN201810706859.2A
Other languages
English (en)
Other versions
CN109191579A (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.)
Guangdong Provincial Institute Of Geological Surveying And Mapping
Original Assignee
Guangdong Provincial Institute Of Geological Surveying And Mapping
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 Guangdong Provincial Institute Of Geological Surveying And Mapping filed Critical Guangdong Provincial Institute Of Geological Surveying And Mapping
Priority to CN201810706859.2A priority Critical patent/CN109191579B/zh
Publication of CN109191579A publication Critical patent/CN109191579A/zh
Application granted granted Critical
Publication of CN109191579B publication Critical patent/CN109191579B/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
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明涉及坐标转换的方法、装置、计算机设备和存储介质，属于地理信息技术领域。所述方法包括：根据已知点确定将第一坐标系转换为第二坐标系的特征参数；根据所述特征参数确定所述已知点的转换残差以及格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。上述技术方案，解决了坐标转换过程忽略点位的空间相关性使得转换坐标准确性不高的问题。在考虑已知点改正量的基础上来求解格网改正量，顾及了空间相关性，实现了高精度的坐标转换。

Description

坐标转换的方法、装置、计算机设备和存储介质
技术领域
本发明涉及地理信息技术领域，特别是涉及坐标转换的方法、装置、计算机设备和存储介质。
背景技术
目前存在多种坐标***，由于不同坐标***之间采用的椭球、框架、历元等不同，这就需要进行坐标的转换，例如将西安1980坐标系向CGCS2000(2000国家大地坐标系)转换。传统的坐标转换往往采用四参数、七参数模型等方式。但是四参数只适合于经度差、纬度差小于2度的区域坐标转换，七参数只适合于省市级的区域坐标转换。近些年，也有研究人员提出采用格网改正模型来实现坐标的转换。在实现本发明过程中，发明人发现现有技术中至少存在如下问题：常规的格网改正模型仅仅根据距离的反加权来实现格网点上改正量的计算，这样忽略了点位的空间相关性，得到的转换坐标准确性不高。
发明内容
基于此，本发明提供了坐标转换的方法、装置、计算机设备和存储介质，能得到准确的转换坐标。
本发明实施例的内容如下：
一种坐标转换的方法，包括以下步骤：根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
在其中一个实施例中，所述特征参数包括四参数；所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤，包括：确定已知点在第一坐标系中的第一坐标以及在第二坐标系中的第二坐标，计算将所述第一坐标转换为第二坐标的四参数。
在其中一个实施例中，根据所述特征参数确定预先建立的格网中的格网点的坐标改正量的步骤，包括：确定所述格网点在第一坐标系中的第一格网坐标，根据所述四参数确定所述格网点在第二坐标系中的第二格网坐标；根据第一格网坐标和第二格网坐标确定所述格网点的坐标改正量。
在其中一个实施例中，所述根据所述已知点的转换残差确定所述格网点的残差改正量的步骤，包括：确定所述已知点与所述格网点的距离；根据所述已知点的转换残差和所述距离，通过距离加权法确定所述格网点的残差改正量。
在其中一个实施例中，所述根据所述坐标改正量和所述残差改正量确定格网改正量的步骤，包括：根据所述坐标改正量和所述残差改正量的和得到所述格网改正量。
在其中一个实施例中，所述转换模型中包含有多个格网点的格网改正量；所述根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型的步骤之后，还包括：从所述转换模型中确定与所述待转换点对应的目标格网点的格网改正量；根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量；根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
在其中一个实施例中，所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤之前，还包括：确定所述格网点预设半径范围内的待分析区域，从所述待分析区域中确定具有所述第一坐标系和第二坐标系的坐标的已知点。
相应的，本发明实施例提供一种坐标转换的装置，包括：参数确定模块，用于根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；改正量确定模块，用于根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；以及，模型建立模块，用于根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
上述坐标转换的方法及装置，根据已知点在第一坐标系和第二坐标系的坐标确定特征参数；根据该特征参数确定已知点的转换残差，并根据该转换残差确定格网点的残差改正量；根据所求的特征参数确定格网点的坐标改正量；根据残差改正量和坐标改正量确定格网改正量并建立转换模型，通过转换模型来确定将第一坐标系中的待转换点转换到第二坐标系后的坐标。在考虑已知点改正量的基础上来求解格网改正量，顾及了空间相关性，实现了高精度的坐标转换。
一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现以下步骤：根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
上述计算机设备，在考虑已知点改正量的基础上来求解格网改正量，顾及了空间相关性，实现了高精度的坐标转换。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
上述计算机可读存储介质，在考虑已知点改正量的基础上来求解格网改正量，顾及了空间相关性，实现了高精度的坐标转换。
附图说明
图1为一个实施例中坐标转换的方法的应用环境图；
图2为一个实施例中坐标转换的方法的流程示意图；
图3为一个实施例中待转换区域的格网的示意图；
图4为一个实施例中转换模型的三维结构示意图；
图5为一个实施例中待转换点投影到格网的示意图；
图6为一个实施例中选取已知点的示意图；
图7为另一个实施例中坐标转换的方法的流程示意图；
图8为一个实施例中坐标转换的装置的结构框图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本申请提供的坐标转换的方法可以应用于如图1所示的计算机设备中。该计算机设备可以是服务器，其内部结构图可以如图1所示。该计算机设备包括通过***总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***、计算机程序和数据库。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的数据库用于存储坐标、坐标改正量、残差改正量、格网改正量等数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时实现一种坐标转换的方法。
本领域技术人员可以理解，图1中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
本发明实施例提供一种坐标转换的方法、装置、计算机设备和存储介质。以下分别进行详细说明。
在一个实施例中，如图2所示，提供了一种坐标转换的方法，能实现对待转换区域内的待转换点的坐标转换，该方法包括以下步骤：
S201、根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标。
在本步骤中，已知点具有第一坐标系和第二坐标系的坐标，根据常规的计算方法计算表征第一坐标系和第二坐标系之间的转换关系的特征参数。
其中，第一坐标系和第二坐标系可以指西安1980坐标系、CGCS2000坐标系、北京54坐标系、WGS-84坐标系等坐标系，本发明实施例对第一坐标系和第二坐标系的类别不做限制，能保证两者不一样且可以进行坐标转换即可。
特征参数指的是表征第一坐标系和第二坐标系之间的转换关系的参数，可以是四参数、七参数等参数。
另外，该已知点可以指待转换区域内的所有具有第一坐标系和第二坐标系的坐标的点，也可以是其中的一部分(例如：预设半径范围内的点)。同时，已知点可以是待转换区域内的随机点(只要它具有第一坐标系和第二坐标系的坐标就选择出来)，也可以是按照一定的规则从待转换区域的多个已知点中选择出来的某些点。
S202、根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量。
本步骤实现的是格网点的坐标改正量以及残差改正量的确定。其中，坐标改正量是根据前一步骤得到的特征参数来确定；而残差改正量是在确定已知点的转换残差后，根据已知点的转换残差以及已知点与格网点的关系来确定。
其中，预先建立的格网指的是在第一坐标系上按照一定的格网间隔绘制的格网图像。其中格网间隔是指：所建立格网的纵横相邻格网的经纬度差异。格网间隔主要根据已知点的平均间隔来确定，一般格网间隔要小于已知点的平均间隔；还可以根据转换图形或待转换区域的大小来确定。其中转换图形指的是需要转换的、依据一定标准进行分幅的地理信息图形。
在某一待转换区域内建立的格网可以如图3所示，图3中的格网由纵横相邻的格网交错构成，每一个格网点都有其在第一坐标系上的坐标。
S203、根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
其中，待转换点指的是具有第一坐标系的坐标但不知其在第二坐标系中的坐标的点。格网改正量指的是某一个格网点的改正量，格网点分布在待转换区域的各个地方，各个待转换点都能根据附近的格网点来确定自己的改正量。
在本步骤中，根据坐标改正量和残差改正量确定格网点的格网改正量，并根据格网改正量建立转换模型，根据该转换模型以及待转换点在第一坐标系中的坐标就能确定待转换点在第二坐标系中的坐标。
转换模型可以根据格网中的一个或多个格网点的格网改正量来确定，例如可以是根据待转换区域中的所有格网点的格网改正量确定。
进一步地，建立转换模型的方法可以是将所确定的格网改正量按照一定的格式存放在一个文件中，该文件就可以作为转换模型。通过读取该文件就可以从中找出所需的目标格网点的格网改正量，进而根据目标格网点的格网改正量确定待转换点在第二坐标系中的坐标。其中，目标格网点可以为多个。
本实施例，格网改正量的获取过程采用周边区域的多个已知点的改正量，即这些已知点都参与了格网改正量的求解，顾及了空间相关性，这样得到的转换模型是连续平滑的，使得大区域的坐标转换具有连续性，从而实现了高精度无缝的大区域坐标转换。
在一个实施例中，所述特征参数包括四参数；所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤，包括：确定已知点在第一坐标系中的第一坐标以及在第二坐标系中的第二坐标，计算将所述第一坐标转换为第二坐标的四参数。
若已知点为多个，可以从中选择一个已知点，根据所选择的已知点的第一坐标和第二坐标来确定四参数。还可以根据全部或者一部分已知点的坐标来确定四参数，若得到有多组四参数，则可以对这些四参数通过一定的运算进行整合。
本实施例限定了特征参数为四参数，四参数的确定过程采用常规的方法实现即可。四参数考虑了小范围区域的坐标转换的偏差，能有效提高所确定的格网改正量的精度。同时四参数的确定过程简单直接，为后续改正量的计算提供数据保障。
在一个实施例中，根据所述特征参数确定预先建立的格网中的格网点的坐标改正量的步骤，包括：确定所述格网点在第一坐标系中的第一格网坐标，根据所述四参数确定所述格网点在第二坐标系中的第二格网坐标；根据第一格网坐标和第二格网坐标确定所述格网点的坐标改正量。
因为格网是在第一坐标系的基础上建立的，因此，各个格网点在第一坐标系中的第一格网坐标能直接获取到。
由于四参数表征了第一坐标系和第二坐标系之间的转换关系，因此，根据四参数就能确定出格网在第二坐标系中的第二格网坐标，将第一格网坐标和第二格网坐标进行比较(例如：求差值)就能确定出格网点的坐标改正量。
本实施例在一般常规格网改正法的基础上，提出采用局部四参数法求解格网点的坐标改正量，根据四参数确定格网的第二格网坐标，并根据第一格网坐标和第二格网坐标来确定格网的坐标改正量。该坐标改正量既表征了格网在第一坐标系的坐标，又表征了格网在第二坐标系中的坐标，因此能用于确定待转换点在第二坐标系中的坐标。
在一个实施例中，所述根据所述已知点的转换残差确定所述格网点的残差改正量的步骤，包括：确定所述已知点与所述格网点的距离；根据所述已知点的转换残差和所述距离，通过距离加权法确定所述格网点的残差改正量。
其中，通过距离加权法确定格网点的残差改正量的过程可以为：通过反距离加权法确定各个已知点的转换残差对于残差改正量的权值(距离越大权值越小)，之后结合转换残差及其对应的权值计算格网点的残差改正量。
本实施例在求解残差改正量的过程中，根据各个已知点的转换残差进行了加权处理，这样使得格网点的解算结果更加符合真实值，提高转换精度。
在一个实施例中，所述根据所述坐标改正量和所述残差改正量确定格网改正量的步骤，包括：根据所述坐标改正量和所述残差改正量的和得到所述格网改正量。
本实施例，根据前述实施例中求得的残差改正量和坐标改正量来确定格网改正量。不仅考虑了转换过程的坐标偏差，还考虑了转换残差，使得偏差的考虑更为全面，能够得到更准确的格网改正量，进而实现更为精确的坐标转换。
在一个实施例中，所述转换模型中包含有多个格网点的格网改正量；所述根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型的步骤之后，还包括：从所述转换模型中确定与所述待转换点对应的目标格网点的格网改正量；根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量；根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
转换模型中包含有多个格网点的格网改正量可以指；转换模型包含待转换区域中所有格网点的格网改正量。根据转换模型绘制的三维结构图如图4所示，在图4中不同坐标的格网改正量可能不同。图4中的格网改正量越大说明偏差越大，图4充分体现了格网改正量的必要性，因为如果不考虑格网改正量，则转换以后的坐标就会存在较大的偏差。
其中，根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标的实现过程可以为：将待转换点在第一坐标系中的坐标加上待转换点的改正量即可得到待转换点在所述第二坐标系中的坐标。
本实施例中，若要对某一待转换点进行坐标转换，首先确定该待转换点所在的格网，通过双线性内插法根据格网点的格网改正量确定待转换点的改正量(在此，待转换点称为内插点)。能准确地确定出待转换点的改正量，进而确定待转换点转换以后的坐标。
在一个实施例中，通过双线性内插法确定内插点的改正量的过程可以为：
如图5所示(在图5中五角星表示内插点，空心圆点表示格网点)，内插点对应的格网内有四个格网点。将插值点向所在格网的x轴、y轴方向分别投影，分别计算投影距离，根据x轴、y轴方向的投影距离确定对应的权值，根据这些权值以及对应四个格网点的格网改正量确定内插点的改正量。
在一个实施例中，内插点对应的格网点还可以不止四个，例如，可以为8个、16个等等。
在一个实施例中，所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤之前，还包括：确定所述格网点预设半径范围内的待分析区域，从所述待分析区域中确定具有所述第一坐标系和第二坐标系的坐标的已知点。
其中，预设半径可以是1m、10m、10km等，该半径可大可小，本发明实施例对预设半径的大小不做限制。
本实施例以格网点为中心按照预设半径确定待分析区域，从待分析区域中确定所有已知点，根据这些已知点来确定格网点的格网改正量，所有已知点均参与了格网改正量的求解，充分考虑了空间相关性，能得到准确的格网改正量。
其中，在待分析区域中选取已知点的示意图见图6，将圆心对应的格网点作为中心，将这些圆的中心以及所有的已知点都选择出来(暂时不考虑其他的格网点)。在图6的右图中，中心点为格网点，其他点为已知点。
在一个实施例中，如图7所示，提供一种坐标转换的方法，包括以下步骤：
S701、确定格网点预设半径范围内的待分析区域，从所述待分析区域中确定具有第一坐标系和第二坐标系的坐标的已知点。
S702、确定已知点在第一坐标系中的第一坐标以及在第二坐标系中的第二坐标，计算将所述第一坐标转换为第二坐标的四参数。
S703、根据所述四参数确定所述已知点的转换残差。
S704、确定所述已知点与所述格网点的距离；根据所述已知点的转换残差和所述距离，通过距离加权法确定所述格网点的残差改正量。
S705、确定所述格网点在第一坐标系中的第一格网坐标，根据所述四参数确定所述格网点在第二坐标系中的第二格网坐标；根据第一格网坐标和第二格网坐标确定所述格网点的坐标改正量。
S706、根据所述坐标改正量和所述残差改正量的和得到格网改正量，根据所述格网改正量建立转换模型。
S707、从所述转换模型中确定与待转换点对应的目标格网点的格网改正量；根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量。
S708、根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
本实施例，在考虑已知点改正量的基础上来求解格网改正量，顾及了空间相关性，实现了高精度的坐标转换。
为了更好地理解上述方法，以下阐述一个本发明坐标转换的方法的应用实例。该应用实例实现的是西安1980坐标和CGCS2000坐标之间的坐标转换，其中待转换区域内有已知点(每个已知点具有西安1980坐标和CGCS2000坐标)。详细说明如下：
一、构建转换模型：
1、已知一张区域地图(待转换区域)中已知点的西安1980坐标系坐标和CGCS2000坐标系坐标。
2、在这张区域地图上制作格网(格网坐标为西安1980坐标系坐标)。
3、以格网中的其中一个格网点A为例，确定格网点A改正量(DX802000，DY802000)，其中，X表示X轴的坐标，Y表示Y轴的坐标，下标为2000表示的是CGCS2000坐标系的X轴坐标，下标为80表示的是西安1980坐标系的X轴坐标，802000表示的是西安1980坐标系向CGCS2000坐标系转换，
改正量的确定过程如下：
S31、在格网点A的适当区域范围内搜索所有已知点；
S32、根据搜索到的各个已知点的坐标(西安1980坐标系坐标和CGCS2000坐标系坐标)计算转换四参数，并基于该四参数计算各个已知点的转换残差；
S33、确定各个已知点与格网点A的距离，根据距离加权法以及已知点的转换残差确定格网点A的残差改正量((dX，dY)；
S34、根据求得的四参数计算格网点A的CGCS2000坐标系坐标，进而计算格网点A的坐标改正量(dX802000，dY802000)；
S35、残差改正量和坐标改正量求和即得到格网点A处的格网改正量：
(DX802000，DY802000)＝(dX802000，dY802000)+(dX，dY)；
S36、确定其他各个格网点的格网改正量，将各个格网改正量按照一定格式放在一个文件里，即为转换模型。
二、根据转换模型进行坐标转换：
4、根据待转换点(B点)的西安1980坐标系坐标确定其位于哪个格网之内。
5、采用双线性内插法以及格网中对应四个格网点的格网改正量确定待转换点B的改正量。
6、确定西安1980坐标系中所有待转换点的改正量，西安1980坐标系坐标加上改正量就得到各个待转换点在CGCS2000中的坐标，即完成西安1980坐标系向CGCS2000坐标系的坐标转换过程。
本实施例，在考虑已知点改正量的基础上来求解格网改正量，顾及了空间相关性，实现了高精度的坐标转换。
需要说明的是，对于前述的各方法实施例，为了简便描述，将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本发明并不受所描述的动作顺序的限制，因为依据本发明，某些步骤可以采用其它顺序或者同时进行。
基于与上述实施例中的坐标转换的方法相同的思想，本发明还提供坐标转换的装置，该装置可用于执行上述坐标转换的方法。为了便于说明，坐标转换的装置实施例的结构示意图中，仅仅示出了与本发明实施例相关的部分，本领域技术人员可以理解，图示结构并不构成对装置的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
如图8所述，坐标转换的装置包括参数确定模块801、改正量确定模块802和模型建立模块803，详细说明如下：
参数确定模块801，用于根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标。
改正量确定模块802，用于根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量。
以及，模型建立模块803，用于根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
本实施例，在考虑已知点改正量的基础上来求解格网改正量，顾及了空间相关性，实现了高精度的坐标转换。
在一个实施例中，参数确定模块801，用于确定已知点在第一坐标系中的第一坐标以及在第二坐标系中的第二坐标，计算将所述第一坐标转换为第二坐标的四参数。
在一个实施例中，改正量确定模块802，包括：坐标确定子模块，用于确定所述格网点在第一坐标系中的第一格网坐标，根据所述四参数确定所述格网点在第二坐标系中的第二格网坐标；改正量确定子模块，用于根据第一格网坐标和第二格网坐标确定所述格网点的坐标改正量。
在一个实施例中，改正量确定模块802，还用于确定所述已知点与所述格网点的距离；根据所述已知点的转换残差和所述距离，通过距离加权法确定所述格网点的残差改正量。
在一个实施例中，模型建立模块803，还用于根据所述坐标改正量和所述残差改正量的和得到所述格网改正量。
在一个实施例中，还包括：目标点改正量确定模块，用于从所述转换模型中确定与所述待转换点对应的目标格网点的格网改正量；待转点改正量确定模块，用于根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量；坐标转换模块，用于根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
在一个实施例中，还包括：已知点确定模块，用于确定所述格网点预设半径范围内的待分析区域，从所述待分析区域中确定具有所述第一坐标系和第二坐标系的坐标的已知点。
需要说明的是，本发明的坐标转换的装置与本发明的坐标转换的方法一一对应，在上述坐标转换的方法的实施例阐述的技术特征及其有益效果均适用于坐标转换的装置的实施例中，具体内容可参见本发明方法实施例中的叙述，此处不再赘述，特此声明。
此外，上述示例的坐标转换的装置的实施方式中，各程序模块的逻辑划分仅是举例说明，实际应用中可以根据需要，例如出于相应硬件的配置要求或者软件的实现的便利考虑，将上述功能分配由不同的程序模块完成，即将所述坐标转换的装置的内部结构划分成不同的程序模块，以完成以上描述的全部或者部分功能。
在一个实施例中，提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现以下步骤：根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：所述特征参数包括四参数；所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤，包括：确定已知点在第一坐标系中的第一坐标以及在第二坐标系中的第二坐标，计算将所述第一坐标转换为第二坐标的四参数。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：根据所述特征参数确定预先建立的格网中的格网点的坐标改正量的步骤，包括：确定所述格网点在第一坐标系中的第一格网坐标，根据所述四参数确定所述格网点在第二坐标系中的第二格网坐标；根据第一格网坐标和第二格网坐标确定所述格网点的坐标改正量。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：所述根据所述已知点的转换残差确定所述格网点的残差改正量的步骤，包括：确定所述已知点与所述格网点的距离；根据所述已知点的转换残差和所述距离，通过距离加权法确定所述格网点的残差改正量。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：所述根据所述坐标改正量和所述残差改正量确定格网改正量的步骤，包括：根据所述坐标改正量和所述残差改正量的和得到所述格网改正量。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：所述转换模型中包含有多个格网点的格网改正量；所述根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型的步骤之后，还包括：从所述转换模型中确定与所述待转换点对应的目标格网点的格网改正量；根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量；根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤之前，还包括：确定所述格网点预设半径范围内的待分析区域，从所述待分析区域中确定具有所述第一坐标系和第二坐标系的坐标的已知点。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：所述特征参数包括四参数；所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤，包括：确定已知点在第一坐标系中的第一坐标以及在第二坐标系中的第二坐标，计算将所述第一坐标转换为第二坐标的四参数。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据所述特征参数确定预先建立的格网中的格网点的坐标改正量的步骤，包括：确定所述格网点在第一坐标系中的第一格网坐标，根据所述四参数确定所述格网点在第二坐标系中的第二格网坐标；根据第一格网坐标和第二格网坐标确定所述格网点的坐标改正量。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：所述根据所述已知点的转换残差确定所述格网点的残差改正量的步骤，包括：确定所述已知点与所述格网点的距离；根据所述已知点的转换残差和所述距离，通过距离加权法确定所述格网点的残差改正量。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：所述根据所述坐标改正量和所述残差改正量确定格网改正量的步骤，包括：根据所述坐标改正量和所述残差改正量的和得到所述格网改正量。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：所述转换模型中包含有多个格网点的格网改正量；所述根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型的步骤之后，还包括：从所述转换模型中确定与所述待转换点对应的目标格网点的格网改正量；根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量；根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤之前，还包括：确定所述格网点预设半径范围内的待分析区域，从所述待分析区域中确定具有所述第一坐标系和第二坐标系的坐标的已知点。
本领域普通技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，作为独立的产品销售或使用。计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(RAM)，只读存储器(ROM)，可擦除可编辑只读存储器(EPROM或闪速存储器)，光纤装置，以及便携式光盘只读存储器(CDROM)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行***执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(PGA)，现场可编程门阵列(FPGA)等。
本发明实施例的术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或(模块)单元的过程、方法、***、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，不能理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (9)

1.一种坐标转换的方法，其特征在于，包括以下步骤：
根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；
根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；
根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标；
所述转换模型中包含有多个格网点的格网改正量；
所述根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量，根据所述格网改正量建立转换模型的步骤之后，还包括：
从所述转换模型中确定与所述待转换点对应的目标格网点的格网改正量；
根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量；
根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
2.根据权利要求1所述的坐标转换的方法，其特征在于，所述特征参数包括四参数；
所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤，包括：
确定已知点在第一坐标系中的第一坐标以及在第二坐标系中的第二坐标，计算将所述第一坐标转换为第二坐标的四参数。
3.根据权利要求2所述的坐标转换的方法，其特征在于，根据所述特征参数确定预先建立的格网中的格网点的坐标改正量的步骤，包括：
确定所述格网点在第一坐标系中的第一格网坐标，根据所述四参数确定所述格网点在第二坐标系中的第二格网坐标；
根据第一格网坐标和第二格网坐标确定所述格网点的坐标改正量。
4.根据权利要求1所述的坐标转换的方法，其特征在于，所述根据所述已知点的转换残差确定所述格网点的残差改正量的步骤，包括：
确定所述已知点与所述格网点的距离；根据所述已知点的转换残差和所述距离，通过距离加权法确定所述格网点的残差改正量。
5.根据权利要求1所述的坐标转换的方法，其特征在于，所述根据所述坐标改正量和所述残差改正量确定格网改正量的步骤，包括：
根据所述坐标改正量和所述残差改正量的和得到所述格网改正量。
6.根据权利要求1至5任一项所述的坐标转换的方法，其特征在于，所述根据已知点确定将第一坐标系转换为第二坐标系的特征参数的步骤之前，还包括：
确定所述格网点预设半径范围内的待分析区域，从所述待分析区域中确定具有所述第一坐标系和第二坐标系的坐标的已知点。
7.一种坐标转换的装置，其特征在于，包括：
参数确定模块，用于根据已知点确定将第一坐标系转换为第二坐标系的特征参数；所述已知点具有所述第一坐标系的坐标和第二坐标系的坐标；
改正量确定模块，用于根据所述特征参数确定所述已知点的转换残差以及预先建立的格网中的格网点的坐标改正量；根据所述已知点的转换残差确定所述格网点的残差改正量；
以及，模型建立模块，用于根据所述坐标改正量和所述残差改正量确定所述格网点的格网改正量建立转换模型；所述转换模型用于确定第一坐标系中的待转换点转换到所述第二坐标系中的坐标；
所述转换模型中包含有多个格网点的格网改正量；
还包括：目标点改正量确定模块，用于从所述转换模型中确定与所述待转换点对应的目标格网点的格网改正量；待转点改正量确定模块，用于根据所述目标格网点的格网改正量，通过双线性内插法确定所述待转换点的改正量；坐标转换模块，用于根据所述待转换点在第一坐标系中的坐标以及所述待转换点的改正量确定所述待转换点在所述第二坐标系中的坐标。
8.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6任一项所述的方法的步骤。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至6任一项所述的方法的步骤。
CN201810706859.2A 2018-07-02 2018-07-02 坐标转换的方法、装置、计算机设备和存储介质 Active CN109191579B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810706859.2A CN109191579B (zh) 2018-07-02 2018-07-02 坐标转换的方法、装置、计算机设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810706859.2A CN109191579B (zh) 2018-07-02 2018-07-02 坐标转换的方法、装置、计算机设备和存储介质

Publications (2)

Publication Number Publication Date
CN109191579A CN109191579A (zh) 2019-01-11
CN109191579B true CN109191579B (zh) 2019-06-28

Family

ID=64948775

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810706859.2A Active CN109191579B (zh) 2018-07-02 2018-07-02 坐标转换的方法、装置、计算机设备和存储介质

Country Status (1)

Country Link
CN (1) CN109191579B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114297324A (zh) * 2021-11-08 2022-04-08 广东省国土资源测绘院 一种地理数据的坐标转换方法、设备、介质及产品
CN116859421B (zh) * 2023-06-27 2024-03-15 国汽大有时空科技(安庆)有限公司 一种多参考框架的定位服务方法及装置
CN117708960B (zh) * 2024-02-04 2024-05-03 武汉大学 平面坐标和正常高的实时转换方法、装置、设备及介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103017653A (zh) * 2012-11-27 2013-04-03 武汉海达数云技术有限公司 一种球面全景影像与三维激光扫描点云的配准及测量方法
CN104537606A (zh) * 2014-12-09 2015-04-22 中国农业大学 地理坐标投影变化方法
CN105300373A (zh) * 2015-09-16 2016-02-03 广州市城市规划勘测设计研究院 一种三维坐标转换方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102339478B (zh) * 2011-09-14 2013-11-13 北京地拓科技发展有限公司 一种从等高线图生成数字高程模型的方法和装置
CN102779231B (zh) * 2012-06-18 2015-11-25 郑州大学 基于邻近的坐标转换参数计算方法
CN102853793B (zh) * 2012-09-27 2015-03-25 中国科学院高能物理研究所 坐标变换数据处理方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103017653A (zh) * 2012-11-27 2013-04-03 武汉海达数云技术有限公司 一种球面全景影像与三维激光扫描点云的配准及测量方法
CN104537606A (zh) * 2014-12-09 2015-04-22 中国农业大学 地理坐标投影变化方法
CN105300373A (zh) * 2015-09-16 2016-02-03 广州市城市规划勘测设计研究院 一种三维坐标转换方法及装置

Also Published As

Publication number Publication date
CN109191579A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN109191579B (zh) 坐标转换的方法、装置、计算机设备和存储介质
CN109189868A (zh) 坐标转换的方法、装置、计算机设备和存储介质
CN110880717B (zh) 基于设备和桥架定位数字化的电缆自动敷设方法
CN107093206A (zh) 利用3d激光扫描技术快速bim建模的方法
CN107391578A (zh) 一种基于格网法的网络地图服务动态坐标转换方法
CN104835202A (zh) 一种三维虚拟场景快速构建方法
CN106780262A (zh) 一种考虑城市道路网络约束的同位模式发现方法及装置
CN109167348A (zh) 基于多参数规划理论的联络线功率可行域确定方法
CN110246095A (zh) 一种图像修复方法、装置、存储介质及终端设备
CN108335354A (zh) 一种基于移动最小二乘法的贪婪投影三角化算法单木重建方法
CN109887073A (zh) 岩心的三维数字模型构建方法及装置
CN106023317B (zh) 一种用于大数据测试的加权Voronoi图生成方法
CN109934866A (zh) 一种图斑划分模式的提取方法、拷贝方法及装置
CN113051637A (zh) 一种关于公共医疗建筑的bim建模方法
CN105809717B (zh) 一种深度估计方法、***及电子设备
CN108921448A (zh) 电能交易方法、装置、设备及存储介质
CN108769996A (zh) 用于定位的信标节点布局方法及定位方法
CN117172134A (zh) 基于融合地形特征的月表多尺度dem建模方法及***
CN109472862A (zh) 一种变电站三维建模***
CN110489510B (zh) 道路数据的处理方法、装置、可读存储介质和计算机设备
CN114329860A (zh) 道路规划方法、装置、存储介质及处理器
CN102324096B (zh) 特征点校正分区拼接方法及装置
CN115293983A (zh) 一种融合多层级互补特征的自适应图像超分辨率复原方法
CN106326492A (zh) 空间矢量数据生成方法及装置
CN102117495B (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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Fu Yan

Inventor after: Wang Jianhui

Inventor after: Li Genger

Inventor after: Han Fei

Inventor after: Wu Lvchuan

Inventor before: Wang Jianhui

Inventor before: Fu Yan

Inventor before: Li Genger

Inventor before: Han Fei

Inventor before: Wu Lvchuan