土地确权登记中数据七参数坐标转换方法.docx

土地确权登记中数据七参数坐标转换方法摘要农村土地承包经营权确权登记颁证工作中涉及数据坐标转换的事项较多对确权、颁证过程中土地信息的准确性以及承包经营权的核准均有重要影响。本文基于上述背景对农村土地承包经营权确权登记颁证工作中数据坐标转换工作进行概述并探讨数据七参数坐标转换方法以期能为农村土地确权工作提供参考。关键词农村土地确权颁证七参数坐标转换ACRGIS一、开展农村土地确权颁证工作的背景在确权颁证过程中要以第二轮土地承包合同以及相应的农村土地承包经营权权属资料为主要依据同时要结合国土二调成果资料、林权发证资料和图件成果、农村集体土地所有权数据、基本农田划定数据、行政区域勘界资料以及相关基础地理信息资料等。由于涉及多部门的多类资料且坐标基准和投影方式各不相同其中包含了1980西安坐标系、1954北京坐标系而本次确权颁证工作明确了坐标系统为2000国家大地坐标系选择高斯-克吕格投影采用标准三度分带。坐标基准不统一将对数据资料的使用造成极大不便还会引起不必要的权属纠纷和部门纠纷。因此对数据坐标进行转换处理在农村土地确权颁证工作中有重要意义。以下将从实际数据出发对数据坐标转换的特点进行简要介绍希望能提高农村土地确权颁证的工作效率。二、坐标转换方法确权颁证工作中大量数据需从80坐标系转换到2000国家大地坐标系由于采用不一样的椭球基准因此转换是不严密的转换一般选用七参数法即X平移Y平移Z平移X旋转Y旋转Z旋转尺度变化K。当然如果区域范围不大最远点间的距离小于30Km（经验值）就可以使用三参数法即X平移Y平移Z平移而将X旋转Y旋转Z旋转尺度变化K视为0三参数只是七参数的一种特例因此本文只对七参数转换法进行介绍。（一）七参数坐标转换原理两个空间直角坐标系它们的原点不一致相应的坐标轴相互不平行两个坐标轴间除了三个平移参数还有三个欧勒角即三个旋转参数又考虑到两个坐标系的尺度不尽相同还需设一个尺度变化参数k总计共有七个参数。用七参数进行空间直角坐标系转换有布尔莎公式、莫洛琴斯基公式和范式公式等。下面给出布尔莎七参数公式：式中的（X0Y0Z0）为三个平移参数（εXεYεZ）为三个旋转角参数m为尺度参数。（二）七参数求取方法由于坐标转换的计算复杂为避免不必要的错误、提高效率通常采用软件进行计算采用布尔莎模型进行求解七参数至少需3个同名点。考虑到确权颁证工作中面对的经常是一个县（区）的资料转换故应多选取几个同名点一般10到15个为宜过少会导致离同名点较远的误差偏大过多则会因为约束条件太多导致整体精度不佳。此外所选同名点分布要尽量均匀并留出2～3个不参与七参数计算用于转换参数校核。以我市红塔区为例确权颁证项目涉及11个乡（镇）行政区域面积1004平方公里村（社区居）委会104个其中：社区94个、村委会10个；村（居）民小组1106个其中：社区居民小组1035个、村民小组71个、自然村436个。人口约40.29万人农户约65670户二调面积耕地为262934.7亩、二轮承包面积120351亩。项目中采用七参数坐标转换方法应用COORDGX参数计算软件步骤如下：第一打开COORDGX软件“工具”菜单中“坐标转换”在“文件”菜单下对“椭球管理”定义好需要的椭球参数。分别输入椭球名称、长半径a和扁率1/e的具体参数。第二将“投影设置”选为高斯三度带投影并设定红塔区域所对应的中央子午线经度。第三同名点输入将源坐标和目标坐标类型均选择为平面坐标椭球分别设置成1980西安和2000国家大地坐标系。选择“计算七参数”按软件格式要求分别输入点的源坐标和目标坐标点击计算。红塔区本次选取12个同名点其中参与计算七参数的同名点9个3个作为校核点。计算过程中未出现误差超限的点当然如果出现超限情况应及时将超限点剔除再次计算。通过计算我们得到（X0、Y0、Z0、εX、εY、εZ、m）7个转换参数。第四将剩余3个同名点带入COORDGX软件用已经求得的七参数进行80坐标到2000大地坐标的转换得出结果与3个点原有的2000大地坐标进行比较差值均在厘米以内满足确权颁证工作的要求。（三）利用七参数对图件进行转换第一打开ARCGIS10.1在“图层”中将图形窗口的地理空间参考参数设置为“国家2000”坐标系。第二在ArcToolbox中选择CreateCustomGeographicTransformation工具在弹出的窗口中输入一个转换的名字如xian80To2000。在定义地理转换方法下面在Method中选择合适的转换方法如COORDIN