林地保护利用规划中数据的坐标转换精选.pdfVIP

林地保护利用规划数据的坐标转换 在 ArcGIS 中的实现 1 2 蒋之富 ，智长贵 （1 国家林业局昆明勘察设计院，昆明 650216， 2 国家林业局调查规划设计院，北京 100714） 摘要：在此次全国林地保护利用规划过程中部分地区需要对数据进行坐标系的转换， 由于国家测绘局已提供了1∶5 万图幅北京54 坐标系与西安80 坐标系之间的大地坐标 改正量，加上此前坐标系转换大都需要进行复杂的求解过程。基于此，本文利用 ArcGIS 简单实现了数据在北京54 坐标系与西安80 坐标系之间的转换。 关键字：林地保护利用规划 坐标系 ArcGIS 1. 背景 全国林地保护利用规划是保障国土生态安全、建设生态文明的基本要求。它将明 确我国中长期林地保护利用的目标和任务，推动我国林地保护利用管理工作健康有序 开展，为应对全球气候变化和我国经济社会可持续发展奠定坚实的基础。此次全国林 地保护利用规划需要参考部分以前调查成果资料，加之此前林业或其它行业许多调查 成果资料在全国各地所采用坐标系不同，所以在此次全国林地保护利用规划编制实施 中，空间数据坐标系的转换是一个必须解决的问题，也是圆满完成全国林地保护利用 规划中重要的一步。 2. 基本原理 1954 年北京坐标系（下文称54 系）是我国20 世纪50 年代为满足测绘工作的需 [1] 要， 从苏联1942 年普尔科夫坐标系转算过来的 ， 采用克拉索夫斯基椭球体， 其参 数为： 地球椭球长半轴: a= 6 378 246 m 椭球扁率: e= 1/298.3 为清除局部平差和逐级控制产生的不合理影响， 提高大地网的精度， 1982 年我国 完成了全国天文大地网整体平差， 建立了新坐标系即1980 西安坐标系（下文称80 系）， 其参考椭球采用1975 IU GG/IA G 第16 届国际大会推荐的地球参数(简称IA G—75 椭球) [2] [3]。 地球椭球长半轴: a = 6 378 140 m 地球引力场二阶球谐系数: J 2= 1 082.63×10- 6 3 2 地球总质量和引力常数乘积: GM = 3.986 005×1014m /S 地球自转角速度: ω= 7.292 115×10- 5 弧度/s 根据上述常数求得椭球扁率: e= 1/2 98.257 由于此次涉及坐标系转换的地区已取得国家测绘局提供的1∶5 万图幅54 系转换 80 系大地坐标改正量 (d X ， d Y )。基于此，常规所用的坐标系转换方法就不实用 了，因为其需要同名地物点的X、Y 以及Z 值，才能求出三参数或七参数。 而在直接有了坐标改正量之后，对1∶2.5 万、1∶1 万图幅数据可根据1∶5 万坐 标改正量进行双向内插求定54 系与80 系控制点、任意地物点坐标转换存在数学关系: d X = X 80- X 54 d Y = Y 80- Y 54 3. 54 系和80 系在ArcGIS 中的坐标转换方法 3.1 矢量数据 对于矢量数据进行坐标转换时，在跨图幅的情况下不能将其按标准1：5 万图幅裁 切之后按坐标改正量进行平移，那将无形中为后期处理工作增加很大的工作量。 此过程可选用ArcGIS 桌面产品—ArcMap 提供的空间转换工具—Spatial Adjustment 来实现。Spatial Adjustment 空间校正方法中有：空间变换 （Transforma