本次投影变换坐标的源数据采用的是采用1980西安的地理坐标系统,1985国家高程基准的1:50000的DLG数据。
二、投影变换基础知识准备
北京54坐标系和西安80坐标系之间的转换其实是两种不同的椭球参数之间的转换。
在ArcGIS中定义了两套坐标系:地理坐标系(Geographic coordinate system)和投影坐标系(Projected coordinate system)。
1、地理坐标系,是以经纬度为地图的存储单位的,是球面坐标系统。地球是一个不规则的椭球,为了将数据信息以科学的方法放到椭球上,这就需要有一个可以量化计算的椭球体。具有长半轴,短半轴,偏心率。一下几行是GCS_Xian_1980椭球及其相应的参数。
Geographic Coordinate System: GCS_Xian_1980
Datum: D_Xian_1980
Prime Meridian: Greenwich
Angular Unit: Degree
每个椭球体都需要一个大地基准面将这个椭球定位,因此可以看到在坐标系统中有Datum: D_Xian_1980的描述,表示,大地基准面是D_Xian_1980。
2、有了椭球体和基准面这两个基本条件,地理坐标系便可以定义投影坐标系统了。以下是已定义Beijing_1954坐标的投影坐标系统的参数:
Projected Coordinate System: Beijing_1954_GK_Zone_19
Projection: Gauss_Kruger
False_Easting: 19500000.00000000
False_Northing: 0.00000000
Central_Meridian: 111.00000000
Scale_Factor: 1.00000000
Latitude_Of_Origin: 0.00000000
Linear Unit: Meter
Geographic Coordinate System: GCS_Beijing_1954
Datum: D_Beijing_1954
Prime Meridian: Greenwich
Angular Unit: Degree
投影坐标系统,实质上是平面坐标系统,其地图单位是米。将球面坐标转化为平面坐标的过程便称为投影,即投影的条件一是有球面坐标,二是要有转化的算法。因此,从参数中可以看出,每一个投影坐标系统都必定会有Geographic Coordinate System。
3、关于坐标偏移量的问题
(1)偏移量的由来
不同国家由于采用的参考椭球及定位方法不同,因此同一地面点在不同坐标系中大地坐标值也不相同。北京1954坐标系的原点在原苏联西部的普尔科夫,采用的是克拉索夫斯基椭球体;西安1980坐标系选用的是1975年国际大地测量协会推荐的参考椭球,其坐标原点设在我国中部的西安市附近的泾阳县境内。
因此,通常情况下,直接转换过来的数据会有一定的误差存在,所以为了保证数据的精度,在转换的过程中通过设置横坐标和纵坐标的偏移量来修正转换后的坐标值。
由西安1980坐标系转换成北京1954坐标系,那么它们的偏移量就是北京1954坐标系相对于WGS84椭球体的偏移量减去西安1980坐标系相对于WGS84偏移量。
(2)偏移量的计算方法
在测区附近选择一国家已知点(X1,Y1),在该已知点上用GPS测定WGS84坐标经纬度,将此坐标视为有误的西安80坐标系,并将其转换为西安80的平面直角坐标X,Y,然后与已知坐标相比较则课计算出偏移量。
即△X1=X- X1
△Y1= Y- Y1
同理可求得北京54坐标系相对于WGS84坐标的偏移量△X2,△Y2,所以由西安80坐标转换成北京54坐标的偏移量即是:△X=△X2-△X1,△Y=△Y2-△Y1
三、“西安80坐标系”转“北京54坐标系”的操作步骤
1、启动ArcMAP,载入coverage数据层,加载arctoolbox工具箱,选择Data Management Tools—>projections and transformations—>feature—>project,打开project对话框,a、在Input Dataset or Feature Class中选择需要进行转换的数据,b、在Output Dataset or Feature Class中选择输出路径和输出的文件名,c、在Output Coordinate System中输入需要定义的地理坐标类型Xian 1980.prj。
其中地理坐标系统在Geographic Coordinate Systems中定义,投影坐标系在Projected Coordinate Systems中选择。
2、上述的coverage数据在定义了西安1980的地理坐标后就转换成了.shp格式的文件,如同上述操作打开project对话框,选择此.shp格式的文件进行投影,投影类型选择Xian 1980 GK Zone 19.prj。
为了区分不同带间的点位,在每个点位的横坐标前加上所在的带号,如Xian 1980 GK Zone 19.prj,即是表示六度分带法的西安80坐标系,分带号为19,横坐标前加带号。
3、为了将数据能够正确的转换为北京54坐标系,需要对以定义西安80坐标系的数据进行平移纠偏。
载入第二步的结果数据,加载Editor工具条,打开Editor—>start editing,让此数据处于编辑状态。加载Spatial Adjustment工具条,选择new displacement link图标,在图像上选择四个Link点,然后打开Link Tabel表修正坐标,纠正公式是X.destination=X.source+66,Y.destination=Y.source+53,66和53是已经计算出来的坐标偏移量,若计算正确,则计算完后可以看到residual error全部变为零,关闭Link Tabel表,
4、已进行了平移纠偏的数据就可以直接转换成北京54坐标了,同样是在arctoolbox工具箱中选择Data Management Tools—>projections and transformations—>Define Projection,当数据量比较大时,可以用批处理操作,速度会快很多,在Samples工具中选择Data Management—>Projections—>Batch Define Coordinate System。