投影参数确定方法、装置、计算机设备、存储介质

技术领域

本申请涉及大地测量技术领域，特别是涉及一种投影参数确定方法、装置、计算机设备、存储介质和计算机程序产品。

背景技术

随着城市化步伐的加速，一方面城区面积不断扩大，由原来的几百平方公里的城市迅速地扩大到几千平方公里甚至上万平方公里；另一方面，大量的人群涌入城市，导致城市基础建设亟待加强，工程领域对城市坐标系也提出了新的要求。按照国家6°带或者3°带的分带方法并以全国参考椭球面为归算面进行投影，会出现在东西方向上随着偏离中央子午线的距离加大投影变形增大的问题，在面积较大的城市，变形量已经不能满足工程建设的需要了，因此有必要建立一个城市的独立坐标系，为生产生活提供方便。建立地方独立坐标系的作用在城市或工程建设地区(如矿山、水库)布设测量控制网时，其成果不仅要满足1：500比列尺测图需要，而且还应该满足一般工程放样的需要，需要减小高程归化与投影变形产生的影响，将它们控制在一个微小的范围，使计算出来的长度在实际利用时(如工程放样)不需要作任何改算。

目前城市发展迅速，城区面积不断扩大，城市工程建设繁多，原有的坐标系及控制网已经不能满足需求，各个城市都相继建立了自己的独立坐标系，为生产生活提供方便。再加上CGCS2000地心坐标系的使用，城市独立坐标系必须与CGCS2000坐标系取得联系，城市独立坐标系的建设也赋予了新的内容。

目前的城市独立坐标系建立过程中大多都采用高斯-克吕格投影即等角横切椭圆柱投影，虽然可以将中央子午线处的投影变形量为0，但中央子午线两边的变形量较大，导致变形量得不到均匀的控制，导致构建出的城市独立坐标系的精度不高。

发明内容

基于此，有必要针对上述技术问题，提供一种能够提高城市独立坐标系构建精度的投影参数确定方法、装置、计算机设备、计算机可读存储介质和计算机程序产品。

第一方面，本申请提供了一种投影参数确定方法。所述方法包括：

获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；

获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；

根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；

将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

在其中一个实施例中，将中央子午线配置于目标区域的中心位置，包括：

对坐标数据进行3度带高斯-克吕格投影，得到目标区域的高斯投影数据；

获取高斯投影数据中各经线的变形量，若存在经线的变形量大于预设变形量，确定目标区域的中心位置；

移动中央子午线，使中央子午线经过目标区域的中心位置。

在其中一个实施例中，投影参数包括q和K，获取多组投影参数，包括：

分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；

根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数。

在其中一个实施例中，方法还包括：

根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据；

根据目标投影数据，建立目标区域的平面坐标系；

根据坐标数据，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中；

对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

在其中一个实施例中，根据目标投影数据，建立目标区域的平面坐标系，包括：

在目标区域中确定一个目标位置，并将目标位置作为平面坐标系的原点；

根据标准坐标系对平面坐标系进行定向，确定平面坐标系的坐标轴方向；

确定目标位置和目标投影数据的对应关系，并根据对应关系、原点和坐标轴方向，建立平面坐标系。

在其中一个实施例中，根据坐标数据，建立目标区域的区域性椭球，包括：

基于目标区域的面积大小确定基准点数量；

根据基准点数量，从坐标数据中获取一个或多个基准点；

根据一个或多个基准点和目标区域的平均高程面，建立目标区域的区域性椭球。

在其中一个实施例中，对区域性椭球进行坐标转换，得到目标区域的独立坐标系，包括：

确定区域性椭球对应的参考坐标系；

采用七参数法计算参考坐标系和目标区域的独立坐标系之间的转换参数；

根据转换参数对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

第二方面，本申请还提供了一种投影参数确定装置。所述装置包括：

获取模块，用于获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；

投影模块，用于获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；

计算模块，用于根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；

比对模块，用于将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

第三方面，本申请还提供了一种计算机设备。所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：

获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；

获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；

根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；

将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

第四方面，本申请还提供了一种计算机可读存储介质。所述计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：

获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；

获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；

根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；

将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

第五方面，本申请还提供了一种计算机程序产品。所述计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以下步骤：

获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；

获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；

根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；

将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

上述投影参数确定方法、装置、计算机设备、存储介质和计算机程序产品，获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。基于目标投影参数获取投影数据，能够提高投影数据的准确性，根据高精确度的投影数据构建城市独立坐标系，能够提高城市独立坐标系构建精度。

附图说明

图1为一个实施例中投影参数确定方法的流程示意图；

图2为一个实施例中七参数法的坐标转换示意图；

图3为一个实施例中投影参数确定装置的结构框图；

图4为一个实施例中计算机设备的内部结构图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

在一个实施例中，如图1所示，提供了一种投影参数确定方法，本实施例以该方法应用于计算机设备进行举例说明，可以理解的是，该计算机设备具体可以是终端或服务器。其中，终端可以但不限于是各种个人计算机、笔记本电脑、智能手机、平板电脑、物联网设备、便携式可穿戴设备，物联网设备可为智能音箱、智能电视、智能空调、智能医用设备等。便携式可穿戴设备可为智能手表、智能手环、头戴设备等。服务器可以用独立的服务器或者是多个服务器组成的服务器集群来实现。本实施例中，该方法包括以下步骤：

步骤102，获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置。

其中，标准坐标系是指CGCS2000坐标系(国家大地坐标系，英文名称为ChinaGeodetic Coordinate System 2000，英文缩写为CGCS2000)或WGS-84坐标系(WorldGeodetic System—1984Coordinate System，是一种国际上采用的地心坐标系)。坐标数据是指基于一种标准坐标系获取的地理信息数据。目标区域可以但不限于是城市区域。

可选的，在CGCS2000坐标系下，获取城市区域的相关坐标数据，然后对坐标数据进行3度带高斯-克吕格投影，若城市区域内存在经线的变形量大于预设变形量，将中央子午线移动至城市区域的中心位置。中心位置可以是对城市区域测量确定的中心点，则中央子午线需要穿过中心点；中心位置也可以是对城市区域测量确定的中心区域，则中央子午线需要穿过中心区域；中心位置也可以是对城市区域测量确定的中心区域，则中央子午线需要穿过中心区域；中心位置还可以是对城市区域测量确定的南北方向的中心线，则可以移动中央子午线与中心线重合。

步骤104，获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据。

可选的，投影参数包括q和K，分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数，按照每一组q和K对城市区域的坐标数据进行投影，得到城市区域的多组样本投影数据。

步骤106，根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差。

其中，均方误差(mean-square error,MSE)是反映估计量与被估计量之间差异程度的一种度量。设t是根据子样确定的总体参数θ的一个估计量，(θ-t)2的数学期望，称为估计量t的均方误差。它等于σ2+b2，其中σ2与b分别是t的方差与偏倚。

步骤108，将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

可选的，将经线变形量的MSE最小的一组样本投影数据对应的一组q和K，作为目标投影参数。

上述投影参数确定方法中，获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。基于目标投影参数获取投影数据，能够提高投影数据的准确性，根据高精确度的投影数据构建城市独立坐标系，能够提高城市独立坐标系构建精度。

在一个实施例中，将中央子午线配置于目标区域的中心位置，包括：对坐标数据进行3度带高斯-克吕格投影，得到目标区域的高斯投影数据；获取高斯投影数据中各经线的变形量，若存在经线的变形量大于预设变形量，确定目标区域的中心位置；移动中央子午线，使中央子午线经过目标区域的中心位置。

进一步的，投影参数包括q和K，获取多组投影参数，包括：分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数。根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据。根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差。将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

可选的，在国家标准的3°带高斯投影下，城市区域内长度变形值大于2.5cm/km时，将中央子午线确定在城区中心，以q＝0，K＝0为初值，以q＝0.0005和K＝2作为上限值，分别以0.0001和0.5为步距对q和K进行变换，并进行不同组合，分别进行投影，以MSE最小为基准，选择能够有效控制变形量的最优投影参数组合。

本实施例中，获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置，获取多组投影参数。适用于城市独立坐标系的建立方法，能够有效防止城区边缘或坐标系边缘的投影变形量超限。

在一个实施例中，方法还包括：根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据。在目标区域中确定一个目标位置，并将目标位置作为平面坐标系的原点；根据标准坐标系对平面坐标系进行定向，确定平面坐标系的坐标轴方向；确定目标位置和目标投影数据的对应关系，并根据对应关系、原点和坐标轴方向，建立平面坐标系。基于目标区域的面积大小确定基准点数量；根据基准点数量，从坐标数据中获取一个或多个基准点；根据一个或多个基准点和目标区域的平均高程面，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中。确定区域性椭球对应的参考坐标系；采用七参数法计算参考坐标系和目标区域的独立坐标系之间的转换参数；根据转换参数对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

可选的，首先，确定城市平面坐标系。城市平面坐标系定向用于确定整个坐标系的起始方位角，与国家坐标系的定向一致，城市坐标原点一般选择在城市区域的中心或易于保存的的、稳定可靠的城市高等级控制点。坐标系原点坐标可加入常数，该常数选择应保证城市平面坐标系所有坐标值为正值。

进一步的，确定区域性椭球，大面积城市采用多点法，小面积城市采用单点法，使区域性椭球与城市的平均高程面充分接近，然后将目标投影数据归算到该椭球面，再将该面的值归算到高斯平面。

最后，进行坐标成果转换。若确定新椭球元素时是以WGS-84椭球为参考椭球，将WGS-84下的坐标，全部转化为新坐标系下的坐标。若确定椭球元素时是使用的CGCS2000坐标，则可以先转换为WGS-84的坐标，然后再转化为城市独立坐标系下的坐标。城市独立坐标转换为CGCS2000或者WGS-84坐标时，可以采用七参数法求解转换参数。

其中，七参数法求解转换参数的方法如图2所示，设坐标系A和坐标系B间有七个转换参数——3个平移参数、3个旋转参数和1个尺度参数。例如，坐标系A为CGCS2000坐标系，坐标系B为城市独立坐标系。

若：(X

(X

(ΔX

(ω

m为坐标系A转换到坐标系B的尺度参数。

则由坐标系A转换到坐标系B的转换关系为：

其中：

一般ω

cosω≈1

sinω≈ω(5)

则有：

在不知转换参数的情况下，可以通过三个或三个以上已知点的两套坐标逆向求解转换参数，方法如下。

式(1)可以写作：

式中，a＝1+m,b＝(1+m)ω

式(7)变形得到：

如果把坐标系A中的坐标当做准确值，坐标系B中的坐标当做观测值，则可以列出误差方程：

当有三个或三个以上的点时，就可以列出九个或九个以上的方程，而参数只有七个，此时就要按照最小二乘原则采用间接平差的方法求解。至于平面坐标的转换，就可以根据同一点的新旧坐标投影成果，采用相似变换的方法求解转换参数。相似变换公式如下：

式中，X

为了便于处理，可以将式(10)简化为：

式中，

本实施例中，根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据；根据目标投影数据，建立目标区域的平面坐标系；根据坐标数据，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中；对区域性椭球进行坐标转换，得到目标区域的独立坐标系。可以获得最优的投影参数组合以均衡不同经度处的投影变形量，减小投影变形，有效防止城区边缘或坐标系边缘的投影变形量超限，提高了城市独立坐标系的构建精度。

在一个实施例中，一种投影参数确定方法，包括：

获取标准坐标系下目标区域的坐标数据，并对坐标数据进行3度带高斯-克吕格投影，得到目标区域的高斯投影数据；获取高斯投影数据中各经线的变形量，若存在经线的变形量大于预设变形量，确定目标区域的中心位置；移动中央子午线，使中央子午线经过目标区域的中心位置。

分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数。根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据。根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差。将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据。

在目标区域中确定一个目标位置，并将目标位置作为平面坐标系的原点；根据标准坐标系对平面坐标系进行定向，确定平面坐标系的坐标轴方向；确定目标位置和目标投影数据的对应关系，并根据对应关系、原点和坐标轴方向，建立平面坐标系。基于目标区域的面积大小确定基准点数量；根据基准点数量，从坐标数据中获取一个或多个基准点；根据一个或多个基准点和目标区域的平均高程面，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中。确定区域性椭球对应的参考坐标系；采用七参数法计算参考坐标系和目标区域的独立坐标系之间的转换参数；根据转换参数对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

例如，选取某市的实测资料作为研究数据，该市的整体区域范围位于东经东经120°55'至122°16'，北纬28°51'至30°33'之间，国家统一的3°带第40带的范围内(中央子午线经度为120°)，所以该市整个城区都位于中央子午线的西侧，经计算得到该市最西侧距离中央子午线的距离为91Km，最东侧距离中央子午线的距离为204km。本次试验所采用的数据是均匀分布在某市的各等级控制点，总计106个，其坐标系统为CGCS2000。该市的平均大地高约为53.65m，平均曲率半径取6371km，将中央子午线移动至城区中心，中心经度为121.53°，选一个基准点gj18(29.8119545，121.4822751，16.9)。将中央子午线取至城区中心，将投影参数q和K初始值分别设置为0和0，分别按照0.0001和0.5的变化步长，将q和K值进行变换，并进行不同组合，以MSE最小为基准，选择最优方案。得到结果发现：若采用国家标准三度带投影，由于城市距离中央子午线的距离过大，106个控制点处的投影变形均无法达到城市测量规范的要求；当只把中央子午线移动至城市中心时，有66个点处的综合投影变形达到了城市测量规范的要求；此外，经过多次组合，发现q值的变化比K值的变化对结果影响更大，而q值在十的负四次方量级上当q值为0.0001，K＝1时，控制形变的效果最好。

将采用高斯投影的计算结果和采用q＝0.0001,K＝1投影参数的结果进行对比，如表1和表2所示，可以看到采用q＝0.0001,K＝1投影参数的双经线割椭圆柱投影，此时双经线在距离中央子午线约±0.8°处，这种投影方式能够更加均匀地限制形变量，防止在城区边缘处(或离中央子午线较远距离处)的投影变形量超限。由此可以说明，把中央子午线移动至城市中央，采用q＝0.0001,K＝1投影参数确定的投影方案的城市独立坐标系的建立方法可以更加有效地抑制投影变形。

表1 q＝0时某市区域的变形量

表1 q＝0时某市区域的变形量

表2 q＝0.0001 K＝1时该市区域的变形量

应该理解的是，虽然如上所述的各实施例所涉及的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，如上所述的各实施例所涉及的流程图中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。

基于同样的发明构思，本申请实施例还提供了一种用于实现上述所涉及的投影参数确定方法的投影参数确定装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现方案相似，故下面所提供的一个或多个投影参数确定装置实施例中的具体限定可以参见上文中对于投影参数确定方法的限定，在此不再赘述。

在一个实施例中，如图3所示，提供了一种投影参数确定装置300，包括：获取模块301、投影模块302、计算模块303和比对模块304，其中：

获取模块301，用于获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；

投影模块302，用于获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；

计算模块303，用于根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；

比对模块304，用于将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

在一个实施例中，获取模块301还用于对坐标数据进行3度带高斯-克吕格投影，得到目标区域的高斯投影数据；获取高斯投影数据中各经线的变形量，若存在经线的变形量大于预设变形量，确定目标区域的中心位置；移动中央子午线，使中央子午线经过目标区域的中心位置。

在一个实施例中，投影模块302还用于分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数。

在一个实施例中，装置还包括：

构建模块305，用于根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据；根据目标投影数据，建立目标区域的平面坐标系；根据坐标数据，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中；对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

在一个实施例中，构建模块305还用于在目标区域中确定一个目标位置，并将目标位置作为平面坐标系的原点；根据标准坐标系对平面坐标系进行定向，确定平面坐标系的坐标轴方向；确定目标位置和目标投影数据的对应关系，并根据对应关系、原点和坐标轴方向，建立平面坐标系。

在一个实施例中，构建模块305还用于基于目标区域的面积大小确定基准点数量；根据基准点数量，从坐标数据中获取一个或多个基准点；根据一个或多个基准点和目标区域的平均高程面，建立目标区域的区域性椭球。

在一个实施例中，构建模块305还用于确定区域性椭球对应的参考坐标系；采用七参数法计算参考坐标系和目标区域的独立坐标系之间的转换参数；根据转换参数对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

上述投影参数确定装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。

在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图4所示。该计算机设备包括处理器、存储器、输入/输出接口、通信接口、显示单元和输入装置。其中，处理器、存储器和输入/输出接口通过系统总线连接，通信接口、显示单元和输入装置通过输入/输出接口连接到系统总线。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质和内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的输入/输出接口用于处理器与外部设备之间交换信息。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信，无线方式可通过WIFI、移动蜂窝网络、NFC(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种投影参数确定方法。该计算机设备的显示单元用于形成视觉可见的画面，可以是显示屏、投影装置或虚拟现实成像装置。显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。

本领域技术人员可以理解，图4中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：对坐标数据进行3度带高斯-克吕格投影，得到目标区域的高斯投影数据；获取高斯投影数据中各经线的变形量，若存在经线的变形量大于预设变形量，确定目标区域的中心位置；移动中央子午线，使中央子午线经过目标区域的中心位置。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据；根据目标投影数据，建立目标区域的平面坐标系；根据坐标数据，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中；对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：在目标区域中确定一个目标位置，并将目标位置作为平面坐标系的原点；根据标准坐标系对平面坐标系进行定向，确定平面坐标系的坐标轴方向；确定目标位置和目标投影数据的对应关系，并根据对应关系、原点和坐标轴方向，建立平面坐标系。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：基于目标区域的面积大小确定基准点数量；根据基准点数量，从坐标数据中获取一个或多个基准点；根据一个或多个基准点和目标区域的平均高程面，建立目标区域的区域性椭球。

在一个实施例中，处理器执行计算机程序时还实现以下步骤：确定区域性椭球对应的参考坐标系；采用七参数法计算参考坐标系和目标区域的独立坐标系之间的转换参数；根据转换参数对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：对坐标数据进行3度带高斯-克吕格投影，得到目标区域的高斯投影数据；获取高斯投影数据中各经线的变形量，若存在经线的变形量大于预设变形量，确定目标区域的中心位置；移动中央子午线，使中央子午线经过目标区域的中心位置。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据；根据目标投影数据，建立目标区域的平面坐标系；根据坐标数据，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中；对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：在目标区域中确定一个目标位置，并将目标位置作为平面坐标系的原点；根据标准坐标系对平面坐标系进行定向，确定平面坐标系的坐标轴方向；确定目标位置和目标投影数据的对应关系，并根据对应关系、原点和坐标轴方向，建立平面坐标系。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：基于目标区域的面积大小确定基准点数量；根据基准点数量，从坐标数据中获取一个或多个基准点；根据一个或多个基准点和目标区域的平均高程面，建立目标区域的区域性椭球。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：确定区域性椭球对应的参考坐标系；采用七参数法计算参考坐标系和目标区域的独立坐标系之间的转换参数；根据转换参数对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

在一个实施例中，提供了一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以下步骤：获取标准坐标系下目标区域的坐标数据，并将中央子午线配置于目标区域的中心位置；获取多组投影参数，并根据每一组投影参数，对坐标数据进行投影，得到目标区域的多组样本投影数据；根据每一组样本投影数据中各经线的变形量，计算每一组样本投影数据的经线变形量的均方误差；将经线变形量的均方误差最小的一组样本投影数据对应的一组投影参数，作为目标投影参数。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：对坐标数据进行3度带高斯-克吕格投影，得到目标区域的高斯投影数据；获取高斯投影数据中各经线的变形量，若存在经线的变形量大于预设变形量，确定目标区域的中心位置；移动中央子午线，使中央子午线经过目标区域的中心位置。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：分别配置q对应的第一初始值和K对应的第二初始值，以及配置q对应的第一步距和K对应的第二步距；根据第一初始值和第一步距对q进行变换，并根据第二初始值和第二步距对K进行变换，获取多个q和K的组合，将每一个q和K的组合作为一组投影参数。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据目标投影参数对坐标数据进行投影，得到目标区域的目标投影数据；根据目标投影数据，建立目标区域的平面坐标系；根据坐标数据，建立目标区域的区域性椭球，并将平面坐标系中的坐标数据映射到区域性椭球中；对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：在目标区域中确定一个目标位置，并将目标位置作为平面坐标系的原点；根据标准坐标系对平面坐标系进行定向，确定平面坐标系的坐标轴方向；确定目标位置和目标投影数据的对应关系，并根据对应关系、原点和坐标轴方向，建立平面坐标系。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：基于目标区域的面积大小确定基准点数量；根据基准点数量，从坐标数据中获取一个或多个基准点；根据一个或多个基准点和目标区域的平均高程面，建立目标区域的区域性椭球。

在一个实施例中，计算机程序被处理器执行时还实现以下步骤：确定区域性椭球对应的参考坐标系；采用七参数法计算参考坐标系和目标区域的独立坐标系之间的转换参数；根据转换参数对区域性椭球进行坐标转换，得到目标区域的独立坐标系。

需要说明的是，本申请所涉及的用户信息(包括但不限于用户设备信息、用户个人信息等)和数据(包括但不限于用于分析的数据、存储的数据、展示的数据等)，均为经用户授权或者经过各方充分授权的信息和数据，且相关数据的收集、使用和处理需要遵守相关国家和地区的相关法律法规和标准。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-OnlyMemory，ROM)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(ReRAM)、磁变存储器(Magnetoresistive Random Access Memory，MRAM)、铁电存储器(Ferroelectric Random Access Memory，FRAM)、相变存储器(Phase Change Memory，PCM)、石墨烯存储器等。易失性存储器可包括随机存取存储器(Random Access Memory，RAM)或外部高速缓冲存储器等。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器(Static Random Access Memory，SRAM)或动态随机存取存储器(Dynamic RandomAccess Memory，DRAM)等。本申请所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等，不限于此。本申请所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等，不限于此。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请的保护范围应以所附权利要求为准。