一种CGCS2000坐标转换方法技术

本发明专利技术提供了一种CGCS2000坐标转换方法，本发明专利技术的目的是提供基础地理信息数据坐标转换过程和方式，批量数据文件坐标转换系统。具体包括以下步骤：步骤1：获取源数据坐标系及坐标类型；步骤2：根据数据类型，选择坐标转换模块；步骤3：设置坐标转换参数(选择数据格式、源坐标系、目标坐标系，设置源坐标系信息、目标坐标信息)；步骤4：选择待转换数据；步骤5：遍历待转换数据和元素集；步骤6：调用高精度高分辨率格网改正量逐点转换；步骤7：更新各数据元素坐标；步骤8：判断是否更新到最后数据的最后一个元素，不是最后一个继续遍历；步骤9：更新到最后数据的最后一个元素，即得到目标坐标系；该坐标转换方法利用全国控制点，计算高精度高分辨率格网改正量，采用逐点转换法，实现基础地理信息成果在全国各个行业的整体统一、高精度、连续无缝转换；可支持批量数据坐标转换，转换效率高。换效率高。换效率高。

【技术实现步骤摘要】
一种CGCS2000坐标转换方法


[0001]本专利技术属于国家基础测绘参心坐标系向地心坐标系转换
，具体涉及一种CGCS2000坐标转换方法。

技术介绍

[0002]为适应国家经济发展的需要，便于资源和成果的共建共享，避免基础设施建设中由于各部门基准不统一造成重大的建设失误与经济损失，经国务院批准，我国自2008年7月1日起启用2000国家大地坐标系—英文名称为China Geodetic Coordinate System 2000，英文缩写为CGCS2000。根据《国家测绘局启用2000国家大地坐标系公告》(2008年第2号)和《关于印发启用2000国家大地坐标系实施方案的通知》(国测国字(2008)24号)，规定2000国家大地坐标系与现行国家大地坐标系转换、衔接的过渡期为8年至10年；2008年7月1日后新生产的各类测绘成果应采用2000国家大地坐标系。自然资源部也于2018年底发布关于停止提供1954年北京坐标系和1980西安坐标系基础测绘成果的公告，致使大量的参心坐标系成果需要转换至2000国家大地坐标系下，以确保国家基础成果在2000国家大地坐标系的正常使用。
[0003]存量参心坐标系成果主要包括DLG数据、栅格数据、DEM数据和文本数据。现有的技术条件下坐标转换针对使用广泛的ArcGIS平台下的gdb数据、栅格数据(img、tiff)和文本数据(txt、xls)。坐标转换参数利用四参数或七参数模型计算。
[0004]现有技术缺陷主要存在以下缺陷：
[0005]1、转换模型采用四参数或七参数，坐标转换适用范围较小；
[0006]2、无法实现全国统一无缝转换，不同区域数据转换后，需要重新接边；
[0007]3、地理信息平台众多，矢量数据格式多样，现有技术适用格式具有局限性；
[0008]4、现有技术不支持DEM数据坐标转换；
[0009]5、针对具有国家标准投影的矢量数据(shp、mdb、gdb、mpj)和栅格数据(img、tiff)，无法实现转后数据自动添加投影。

技术实现思路

[0010](1)要解决的技术问题
[0011]本专利技术为了实现存量基础地理信息数据从参心坐标系向2000国家大地坐标系的转换，本专利技术要解决的技术问题是针对全国各类型、各格式的数据，确定坐标转换模型和坐标转换方法，提供一种快速、高效、无缝的基础地理信息数据坐标转换系统。
[0012](2)技术方案
[0013]本专利技术的目的是提供一种CGCS2000坐标转换方法，包括如下步骤：
[0014]步骤1：获取源数据坐标系及坐标类型；
[0015]步骤2：根据数据类型，选择坐标转换模块；
[0016]步骤3：设置坐标转换参数(选择数据格式、源坐标系、目标坐标系，设置源坐标系
信息、目标坐标信息)；
[0017]步骤4：选择待转换数据；
[0018]步骤5：遍历待转换数据和元素集；
[0019]步骤6：调用全国统一的高精度高分辨率格网改正量逐点转换；
[0020]步骤7：更新各数据元素坐标；
[0021]步骤8：判断是否更新到最后数据的最后一个元素，不是最后一个继续遍历；
[0022]步骤9：更新到最后数据的最后一个元素，即得到目标坐标系。
[0023]进一步的，所述步骤6：调用全国统一的高精度高分辨率格网改正量逐点转换，采用的是移动拟合转换法。所述坐标转换改正量，用反距离加权平均模型计算：
[0024][0025]其中，(B0,L0)为格网点大地坐标，(ΔB
i
,ΔL
i
)为各控制点坐标转换改正量，用下式计算：
[0026]ΔB
i
＝B
T
‑
B
S
[0027]ΔL
i
＝L
T
‑
L
S
[0028]式中，(B
T
,L
T
)分别为目标坐标系大地坐标，(B
S
,L
S
)分别为源坐标系大地坐标。
[0029]进一步的，所述步骤3中坐标转换类型包括1954年北京坐标系、1980西安坐标系向2000国家大地坐标系转换；同一坐标系下坐标正反算、高斯平面坐标换带(3度分带、6度分带、任意带)。
[0030]进一步的，所述步骤4中待转换数据包括文本数据、DLG数据、DEM数据、栅格数据，具体数据格式包括以下内容：
[0031]文本数据：格式包括txt、xlsx和xls。
[0032]DLG数据：格式包括shp、mdb、gdb、dwg、dxf、mpj、sdb。具有正确投影信息的shp、mdb、gdb数据，可自动添加转换后坐标系投影信息。
[0033]DEM数据：格式包括tif、img、asc、grid、dem。此外，具备asc、dem两种格式坐标转换后相邻图形间的接边检查功能。
[0034]栅格数据：格式包括tif和img。具有正确投影信息的tif、img数据，可自动添加转换后坐标系投影信息
[0035]进一步的，所述文本数据转换包括如下步骤：
[0036]步骤501、获取源坐标系文本数据；
[0037]步骤502、按要求整理坐标；
[0038]步骤503、读取源坐标系坐标，利用全国统一的高精度高分辨率格网改正量计算目标坐标系坐标；
[0039]步骤504、将目标坐标系坐标写入文本；
[0040]步骤505、得到目标坐标系文本数据。
[0041]进一步的，所述DLG数据转换包括如下步骤：
[0042]步骤601、获取源坐标系DLG数据；
[0043]步骤602、获取源要素类的投影信息；
[0044]步骤603、获取各要素的节点坐标；
[0045]步骤604、将源坐标系高斯平面坐标转换到大地坐标；
[0046]步骤605、逐点计算其目标坐标系下各要素的节点坐标；
[0047]步骤606、将目标坐标系大地坐标转换为高斯坐标；
[0048]步骤607、目标坐标系下逐要素写入要素类；
[0049]步骤608、目标坐标系数据图廓整饰；
[0050]步骤609、添加目标坐标系投影信息；
[0051]步骤610、得到目标坐标系DLG数据。
[0052]进一步的，所述DEM数据转换包括如下步骤：
[0053]步骤701、读取源坐标系DEM数据；
[0054]步骤702、将DEM数据转换为ASCII文本格式；
[0055]步骤703、获取源坐标系下DEM坐标，计算其目标坐标系下坐标；
[0056]步骤704、将主块与周边相邻块DEM数据进行拼接；
[0057]步骤705、计算起算点格网点坐标，参考像素分辨率，重新生成DEM数据；
[0058]本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种CGCS2000坐标转换方法，其特征在于，包括如下步骤：步骤1：获取源数据坐标系及坐标类型；步骤2：根据数据类型，选择坐标转换模块；步骤3：设置坐标转换参数(选择数据格式、源坐标系、目标坐标系，设置源坐标系信息、目标坐标信息)；步骤4：选择待转换数据；步骤5：遍历待转换数据和元素集；步骤6：调用全国统一的高精度高分辨率格网改正量逐点转换；步骤7：更新各数据元素坐标；步骤8：判断是否更新到最后数据的最后一个元素，不是最后一个继续遍历；步骤9：更新到最后数据的最后一个元素，即得到目标坐标系。2.根据权利要求1所述一种CGCS2000坐标转换方法，其特征在于：所述步骤6，调用全国统一的高精度高分辨率格网改正量逐点转换，采用的是移动拟合转换法。格网改正量用反距离加权平均模型计算：其中，(B0,L0)为格网点大地坐标，(ΔB
i
,ΔL
i
)为各控制点坐标转换改正量，用下式计算：ΔB
i
＝B
T
‑
B
S
ΔL
i
＝L
T
‑
L
S
式中，(B
T
,L
T
)分别为目标坐标系大地坐标，(B
S
,L
S
)分别为源坐标系大地坐标。利用全国统一的高精度高分辨率改正量，保证坐标转换后仍为全国一张图。3.根据权利要求1所述一种CGCS2000坐标转换方法，其特征在于：所述步骤3中坐标转换类型包括1954年北京坐标系、1980西安坐标系向2000国家大地坐标系转换；同一坐标系下坐标正反算、高斯平面坐标换带(3度分带、6度分带、任意带)。4.根据权利要求1所述一种CGCS2000坐标转换方法，其特征在于：所述步骤4中待转换数据具体包括以下内容：文本数据：格式包括txt、xlsx和xls。DLG数据：格式包括shp、mdb、gdb、dwg、dxf、mpj、sdb。具有正确投影信息的shp、mdb、gdb数据，可自动添加转换后坐标系投影信息。DEM数据：格式包括tif、img、asc、grid、dem。此外，具备asc、dem两种格式坐标转换后相邻图形间的接边检查。栅格数据：格式包括tif和img。具有正确投影信息的tif、img数据，可自动添加转换后坐标系投影信息。5.根据权利要求5所述一种CGCS2000坐标转换方法，其特征在于：所述文本数据转换包括如下步骤：
步骤501、获取源坐标系文本数据；步骤502、按要求整理坐标；步骤503、读取源坐标系坐标，利用高精度、高分辨率格网改正量计算目标坐标系坐标；步骤504、将目标坐标系坐标写入文本；步骤505、得到目标坐标系文本数据。6.根据权利要求5所述一...

【专利技术属性】
技术研发人员：范宏涛，张颜秋，樊琦，
申请(专利权)人：自然资源部大地测量数据处理中心，
类型：发明
国别省市：