2000独立坐标系参数获取方法、装置与计算机可读存储介质制造方法及图纸

本发明专利技术提供了一种2000独立坐标系参数获取方法、装置与计算机可读存储介质，该方法包括根据目标城市坐标系的中央子午线与投影面高程，设置第一投影参数；根据第一投影参数，对预先获取的目标城市的多个控制网中重合控制点的CGCS2000大地坐标进行高斯投影获得第一高斯投影坐标；根据第一高斯投影坐标、预先获取的目标城市重合控制点的现有平面坐标以及预设的参数转换模型，计算第一平面转换参数；根据第一平面转换参数的旋转角参数与尺度参数，确定最佳平移参数与最佳投影参数。通过该方法能够准确获取目标城市2000独立坐标系的投影参数和平移参数，实现目标城市2000独立坐标系与现有坐标系平稳过度，提高目标城市2000独立坐标系与城市现有坐标系之间坐标一致性。

2000 independent coordinate system parameter acquisition method, device and computer readable storage medium

The invention provides a method, device and computer readable storage medium for obtaining parameters of 2000 independent coordinate system. The method includes setting the first projection parameters according to the central meridian and projection plane elevation of the target city coordinate system, and controlling multiple control networks of the pre-acquired target city according to the first projection parameters. The CGCS2000 geodetic coordinates of manufacturing points are projected by Gauss to obtain the first Gauss projection coordinates; the first plane transformation parameters are calculated according to the first Gauss projection coordinates, the existing plane coordinates of the target city coincidence control points obtained in advance and the preset parameter transformation model; and the rotation angle parameters and rulers of the first plane transformation parameters are calculated according to the first plane transformation parameters. The best translation parameters and the best projection parameters are determined. The projection and translation parameters of the 2000 independent coordinate system of the target city can be accurately obtained by this method. The smooth transition between the 2000 independent coordinate system of the target city and the existing coordinate system can be realized, and the coordinate consistency between the 2000 independent coordinate system of the target city and the existing coordinate system of the city can be improved.

【技术实现步骤摘要】
2000独立坐标系参数获取方法、装置与计算机可读存储介质
本专利技术涉及坐标系参数获取
，具体涉及一种2000独立坐标系参数获取方法、装置与计算机可读存储介质。
技术介绍
目前，在国家测绘局的建立城市2000独立坐标系的要求下，国内一部分城市建立了2000独立坐标系，例如：武汉2000城市独立坐标系、宁波2000城市独立坐标系以及南昌2000城市独立坐标系；但是在建立2000城市独立坐标系的过程中存在很多问题，最典型的是无法保证新旧坐标的一致性，例如：武汉2000城市独立坐标系由于城市东、西部边缘地区以及大地高超过145m的地区投影变形超限，且没有考虑新旧资料衔接问题，新旧坐标系坐标有一定的差异；宁波2000城市独立坐标系由于城市东、西两侧投影变形不满足小于1/4万(2.5cm/km)的要求，新老坐标之间只需要平移，但是平移参数为100km的大数A加上一个小数B，但不同地区B不一致，不便于使用；南昌2000城市独立坐标系由于东、西部边缘地区投影变形不满足要求，新旧坐标系之间坐标差异较大，规划区1:500地形图资料需进行坐标转换才可继续使用，需投入大量人力物力财力开展图件转换工作。因此,基于CGCS2000参考椭球新建城市2000独立坐标系的过程中,如何实现既满足投影长度变形不超过1/4万的规范要求，又满足与历史资料无缝衔接，使得已有坐标系下成果可延续使用，避免两坐标系的不一致性导致巨大的转换工作量及费用，成为新建基于CGCS2000参考椭球的城市2000独立坐标系亟待解决的问题。
技术实现思路
本专利技术的目的是提供一种2000独立坐标系参数获取方法、装置与计算机可读存储介质，能够准确获取城市基于2000椭球体的目标城市2000坐标系的投影参数与平移参数，实现目标城市2000独立坐标系与现有坐标系平稳过度，提高目标城市2000独立坐标系与城市现有坐标系之间坐标一致性，以使得新旧资料无缝衔接，同时已有坐标系下成果可延续使用，降低了转换工作量及费用。本专利技术实施例提供了一种2000独立坐标系参数获取方法，包括：获取目标城市的多个控制网中重合控制点的CGCS2000大地坐标以及所述重合控制点对应的目标城市坐标系的平面坐标；根据所述目标城市坐标系的中央子午线以及投影面高程，设置第一投影参数；根据所述第一投影参数，对所述重合控制点的CGCS2000大地坐标进行高斯投影，计算所述重合控制点的第一高斯投影坐标；根据所述重合控制点的第一高斯投影坐标、所述平面坐标以及预设的参数转换模型，计算所述重合控制点的第一高斯投影坐标与所述平面坐标之间的第一平面转换参数；根据所述第一平面转换参数的旋转角参数以及尺度参数，确定最佳平移参数以及最佳投影参数。优选地，所述根据所述第一平面转换参数的旋转角参数以及尺度参数，确定最佳平移参数以及最佳投影参数，具体包括；根据所述第一平面转换参数，设置第二投影参数；根据所述第二投影参数，对所述重合控制点的CGCS2000大地坐标进行高斯投影，计算所述重合控制点的第二高斯投影坐标；根据所述重合控制点的第二高斯投影坐标、所述平面坐标以及所述预设的参数转换模型，计算所述重合控制点的第二高斯投影坐标与所述重合控制点的平面坐标之间的第二平面转换参数；根据所述第二平面转换参数的旋转角参数以及尺度参数，确定所述最佳平移参数以及所述最佳投影参数。优选地，所述根据所述目标城市坐标系的中央子午线以及投影面高程，设置第一投影参数，具体包括：所述第一投影参数包括第一中央子午线以及第一投影面高程；根据所述目标城市坐标系的中央子午线，确定所述第一中央子午线的数值区间，并按照设定的第一间隔设置所述第一中央子午线；根据所述目标城市坐标系的投影高程面，确定所述第一投影高程面的数值区间，并按照设定的第二间隔设置所述第一投影高程面。优选地，所述根据所述第一平面转换参数，设置第二投影参数，具体包括：所述第二投影参数包括第二中央子午线以及第二投影面高程；根据所述第一平面转换参数，确定所述第二中央子午线的数值区间以及所述第二投影面高程的数值区间；按照设定的第三间隔设置所述第二中央子午线以及按照设定的第四间隔设置所述第二投影高程面。优选地，所述按照设定的第三间隔设置所述第二中央子午线以及按照设定的第四间隔设置所述第二投影高程面，具体包括；当所述第二中央子午线为L时，按照所述设定的第四间隔设置所述第一投影高程面，获得N组所述第二中央子午线与所述第二投影高程面的参数组合；其中，L为属于所述第二中央子午线的数值区间内且按照设定的第三间隔设置的中央子午线；当所述第二投影高程面为H时，按照设定的第三间隔设置所述第二中央子午线，获得M组所述第二中央子午线与所述第二投影高程面的参数组合；H为属于所述第二投影高程面的数值区间内且按照设定的第四间隔设置的投影高程面。优选地，所述根据所述第二平面转换参数的旋转角参数以及尺度参数，确定最佳平移参数以及最佳投影参数，具体包括：从M+N个所述第二平面转换参数中获取符合所述第二平面转换参数的旋转角参数α≈0且所述第二平面转换参数的尺度参数m≈1的第i个第二平面转换参数；其中，根据N组所述第二中央子午线与所述第二投影高程面的参数组合以及M组所述第二中央子午线与所述第二投影高程面的参数组合，可获得M+N个所述第二平面转换参数；获取所述第i个第二平面转换参数的平移参数；其中，所述第i个第二平面转换参数的平移参数为所述最佳平移参数；所述最佳投影参数包括最佳中央子午线以及最佳投影高程面；根据所述M+N个所述第二平面转换参数的旋转角参数α与所述第二中央子午线L，拟合出所述旋转角参数与所述第二中央子午线的第一拟合函数；根据所述第一拟合函数，计算当旋转角参数α＝0时对应的最佳中央子午线；根据所述M+N个所述第二平面转换参数的尺度参数m与所述第二投影高程面H，拟合出所述尺度参数与所述第二投影高程面的第二拟合函数；根据所述第二拟合函数，计算当尺度参数m＝1时对应的最佳投影高程面。优选地，所述2000独立坐标系参数获取方法还包括：根据所述重合控制点的第一高斯投影坐标、所述平面坐标以及预设的参数转换模型，计算所述重合控制点的内部精度；获取所述多个控制网的外部检核点；其中，所述外部检核点是所述多个控制网中与所述重合控制点不重合的控制点；根据所述重合控制点以及所述外部检核点，计算所述重合控制点的外部精度；根据所述重合控制点的内部精度以及外部精度，确定所述重合控制点中的异常控制点；剔除所述重合控制点中的异常控制点，并重新计算剔除所述异常控制点后的重合控制点对应的第一平面转换参数。本专利技术实施例还提供了一种2000独立坐标系参数获取装置，包括：坐标获取模块，用于获取目标城市的多个控制网中重合控制点的CGCS2000大地坐标以及所述重合控制点对应的目标城市坐标系的平面坐标；第一投影参数设置模块，用于根据所述目标城市坐标系的中央子午线以及投影面高程，设置第一投影参数；第一高斯投影坐标计算模块，用于根据所述第一投影参数，对所述重合控制点的CGCS2000大地坐标进行高斯投影，计算所述重合控制点的第一高斯投影坐标；第一平面转换参数计算模块，用于根据所述重合控制点的第一高斯投影坐标、所述平面坐标以及预设的参数转换模型，计算所述重合控制点的第一高斯投影坐标与本文档来自技高网...

【技术保护点】
1.一种2000独立坐标系参数获取方法，其特征在于，包括：获取目标城市的多个控制网中重合控制点的CGCS2000大地坐标以及所述重合控制点对应的目标城市坐标系的平面坐标；根据所述目标城市坐标系的中央子午线以及投影面高程，设置第一投影参数；根据所述第一投影参数，对所述重合控制点的CGCS2000大地坐标进行高斯投影，计算所述重合控制点的第一高斯投影坐标；根据所述重合控制点的第一高斯投影坐标、所述平面坐标以及预设的参数转换模型，计算所述重合控制点的第一高斯投影坐标与所述平面坐标之间的第一平面转换参数；根据所述第一平面转换参数的旋转角参数以及尺度参数，确定最佳平移参数以及最佳投影参数。

【技术特征摘要】
1.一种2000独立坐标系参数获取方法，其特征在于，包括：获取目标城市的多个控制网中重合控制点的CGCS2000大地坐标以及所述重合控制点对应的目标城市坐标系的平面坐标；根据所述目标城市坐标系的中央子午线以及投影面高程，设置第一投影参数；根据所述第一投影参数，对所述重合控制点的CGCS2000大地坐标进行高斯投影，计算所述重合控制点的第一高斯投影坐标；根据所述重合控制点的第一高斯投影坐标、所述平面坐标以及预设的参数转换模型，计算所述重合控制点的第一高斯投影坐标与所述平面坐标之间的第一平面转换参数；根据所述第一平面转换参数的旋转角参数以及尺度参数，确定最佳平移参数以及最佳投影参数。2.如权利要求1所述的2000独立坐标系参数获取方法，其特征在于，所述根据所述第一平面转换参数的旋转角参数以及尺度参数，确定最佳平移参数以及最佳投影参数，具体包括；根据所述第一平面转换参数，设置第二投影参数；根据所述第二投影参数，对所述重合控制点的CGCS2000大地坐标进行高斯投影，计算所述重合控制点的第二高斯投影坐标；根据所述重合控制点的第二高斯投影坐标、所述平面坐标以及所述预设的参数转换模型，计算所述重合控制点的第二高斯投影坐标与所述重合控制点的平面坐标之间的第二平面转换参数；根据所述第二平面转换参数的旋转角参数以及尺度参数，确定所述最佳平移参数以及所述最佳投影参数。3.如权利要求1所述的2000独立坐标系参数获取方法，其特征在于，所述根据所述目标城市坐标系的中央子午线以及投影面高程，设置第一投影参数，具体包括：所述第一投影参数包括第一中央子午线以及第一投影面高程；根据所述目标城市坐标系的中央子午线，确定所述第一中央子午线的数值区间，并按照设定的第一间隔设置所述第一中央子午线；根据所述目标城市坐标系的投影高程面，确定所述第一投影高程面的数值区间，并按照设定的第二间隔设置所述第一投影高程面。4.如权利要求2所述的2000独立坐标系参数获取方法，其特征在于，所述根据所述第一平面转换参数，设置第二投影参数，具体包括：所述第二投影参数包括第二中央子午线以及第二投影面高程；根据所述第一平面转换参数，确定所述第二中央子午线的数值区间以及所述第二投影面高程的数值区间；按照设定的第三间隔设置所述第二中央子午线以及按照设定的第四间隔设置所述第二投影高程面。5.如权利要求4所述的2000独立坐标系参数获取方法，其特征在于，所述按照设定的第三间隔设置所述第二中央子午线以及按照设定的第四间隔设置所述第二投影高程面，具体包括；当所述第二中央子午线为L时，按照所述设定的第四间隔设置所述第一投影高程面，获得N组所述第二中央子午线与所述第二投影高程面的参数组合；其中，L为属于所述第二中央子午线的数值区间内且按照设定的第三间隔设置的中央子午线；当所述第二投影高程面为H时，按照设定的第三间隔设置所述第二中央子午线，获得M组所述第二中央子午线与所述第二投影高程面的参数组合；H为属于所述第二投影高程面的数值区间内且按照设定的第四间隔设置的投影高程面。6.如权利要求5所述的2000独立坐标系参数获取方法，其特征在于，所述根据所述第二平面转换参数的旋转角参数...

【专利技术属性】
技术研发人员：林鸿，欧海平，秦亮军，喻永平，刘洋，余锐，王天应，张荣，叶玲洁，王鹏，娄俊萍，
申请(专利权)人：广州市城市规划勘测设计研究院，
类型：发明
国别省市：广东,44