FME坐标系使用

1、fme 坐标系使用关于beijing54坐标系和xian80坐标系说明的补充前面的博客介绍了在fme中关于beijing54坐标系和xian80的介绍，这里进行一下详细的说明。根据坐标系的定义规则，这两种坐标系都定义了椭球体参数，基准面参数和投影参数。以75度带为中央经度线，高斯投影为例，具体如下：一、beijing54坐标系1. 椭球体参数：ellipsoid parameters：desc_nm: krassovsky - 1940/1948 /椭球体名称：格拉索夫斯基椭球体1940/1948e_rad: 6378245 /椭球体长半径p_rad: 6356863.01877305 /椭球

2、体短半径source: us defense mapping agency, tr-8350.2-b, december 1987 /数据来源这里所用的椭球体参数是标准的参数，和我国选用的一致。2. 基准面参：datum parameters：delta_x: -31.4 /x偏移量delta_y: 144.3 /y偏移量delta_z: 81.2 /z偏移量desc_nm: beijing 1954, /名称ellipsoid: krasov /所用椭球体source: mentor software client /数据来源use: geocentric /使用对于这里的定义，是选择了国际

3、通用的一个三参数法的参数，是因为国家的基准面参数是保密的，而且每个地区的参数不可能是一样的。而对于fme来讲，不可能得到国内的参数，因此是选择了一个统一的参数，对于所有的投影参数，都是使用的这样一个参数，因此利用fme自身带的参数进行空间数据的坐标系变换，会带来一定的误差。只有通过自定义方式，把一个地区的准确的参数，定义一个基准面和坐标系，才有可能完成准确的坐标系变换。3. 投影参数coordinate system parameterscs_name: beijing1954.gk-13 /投影坐标系名称desc_nm: beijing 1954 / gauss-kruger zone 13

4、 /投影坐标系描述dt_name: beijing1954 /选用的基准面group: asia /坐标系分组map_scl: 1 /地图比例parm1: 75 /中央经度线proj: tm /投影类型：tmquad: 1 /scl_red: 1 /source: epsg, v6.6, 21413 /数据来源：epsgunit: meter /坐标单位x_off: 13500000 /x位移量zero_x: 0.0001 /zero_y: 0.0001 /这里就不多介绍这些参数了，基本上都应该了解。二、xian80坐标系1. 椭球体参数：ellipsoid parameters：desc_n

5、m: xian 1980 /椭球体名称：xian 1980e_rad: 6378140 /椭球体长半径p_rad: 6356755.288158 /椭球体短半径source: epsg, v6.3, 7049 epsg /数据来源：epsg这里所用的椭球体参数是标准的参数，和我国选用的一致。2. 基准面参：datum parameters：bwscale: 1.1 /bursa（布尔沙算法）的比例因子delta_x: 24 /x偏移量delta_y: -123 /y偏移量delta_z: -94 /z偏移量desc_nm: xian 1980, /基准面描述ellipsoid: xian80

6、/所引用的椭球体rot_x: -0.02 /x轴旋转rot_y: -0.25 /y轴旋转rot_z: 0.13 /z轴旋转source: mentor software client /数据来源use: 7parameter /基准面的计算方法：7参数对于这里的定义，是选择了国际通用的一个7参数法的参数，是因为国家的基准面参数是保密的，而且每个地区的参数不可能是一样的。而对于fme来讲，不可能得到国内的参数，因此是选择了一个统一的参数，对于所有的投影参数，都是使用的这样一个参数，因此利用fme自身带的参数进行空间数据的坐标系变换，会带来一定的误差。只有通过自定义方式，把一个地区的准确的参数，定

7、义一个基准面和坐标系，才有可能完成准确的坐标系变换。3. 投影参数coordinate system parameterscs_name: xian80.gk-13 /投影坐标系名称desc_nm: xian 1980 / gauss-kruger zone 13 /投影坐标系描述dt_name: xian80 /所引用的基准面名称epsg: 2327 /epsg编号group: asia /坐标系分组map_scl: 1 /地图比例parm1: 75 /中央经度线proj: tm /投影类型：tmquad: 1 /scl_red: 1 /source: epsg, v6.6, 21413 /

8、数据来源：epsgunit: meter /坐标单位x_off: 13500000 /x位移量zero_x: 0.0001 /zero_y: 0.0001 /这里就不多介绍这些参数了，基本上都应该了解。三、对于预定义坐标系的使用对于fme预定义的beijing54坐标系和xian80坐标系，是定义了和wgs84之间的坐标系变换的转换关系，是一个国际上采用的大致的参数，不是准确的参数，因此转换的数据在精度要求不高的情况，是可以使用的。如果，需要准确的转换坐标系，就需要获取准确的转换参数，完成坐标系的自定义，才能进行。关于坐标系变换对于fme进行坐标系变换，涉及两种不同的变换，一个是基准面（dat

9、um）发生变换，另一个基准面不发生变换。本文对第一种变换进行描述。对于基准面发生变换的坐标系变换，在fme workbench中使用reprojector函数，在该函数中要分别选择源数据坐标系和目标数据坐标系，确定后，即可进行坐标系的变换。 或者对于源数据和目标数据，分别设置坐标系，fme在进行数据转换的时候，就可以自动进行坐标系的变换。 但是如何设置一个用户自定义的坐标系，能在fme的坐标系仓库（coordinate system gallery）中选择，然后进行坐标系变换呢。通常，定义一个坐标系是确定这个坐标系的几个参数-椭球体参数，基准面参数以及投影参数。对于fme的坐标系参数定义涉及两

10、个文件，localcoordsysdefs.fme和mycoordsysdefs.fme。这两个文件都位于fme安装目录下的子目录reproject下。在localcoordsysdefs.fme文件中定义基本的参数-椭球体参数和基准面参数。在文件mycoordsysdefs.fme中定义投影参数。椭球体的定义：ellipsoid_def desc_nm source e_rad p_rad 参数描述如下：名称类型描述是否为空 字符串椭球体名称no 字符串椭球体的描述yes 字符串椭球体来源（个人或者组织）yes 浮点数椭球体长半径no 浮点数椭球体短半径no如果一个椭球体已经在fme里面预先

11、定义了，在以后所有的用到的地方都可以引用，不必重新定义，比如krasov，或者xian80等。基准面的定义：datum_def desc_nm source ellipsoid use delta_x delta_y delta_z bwscale rot_x rot_y rot_z 参数描述如下：名称类型描述是否为空字符串基准面的名称no字符串基准面的描述yes字符串椭球体来源（个人或者组织）yes有效的椭球体名称引用已定义的椭球体名称no3parameter4parameter6parameter7parameteragd66agd84ats77bursacsrsdhdngda94hpgn

12、jgd2kmolodenskymulregnad27nad83nzgd2knzgd49wgs72wgs84基准面变换中使用的变换类型（变换公式） no浮点数空间直角坐标系x轴与wgs-84基准面空间直角坐标系x轴的位移量。（单位：米）no浮点数空间直角坐标系y轴与wgs-84基准面空间直角坐标系y轴的位移量。（单位：米）no浮点数空间直角坐标系z轴与wgs-84基准面空间直角坐标系z轴的位移量。（单位：米）no浮点数7参数转换中的比率（bursa类型）。（单位：ppm 百万分之）不是所有类型都需要浮点数7参数转换中的x方向的旋转偏移（bursa类型）。（单位：弧度）不是所有类型都需要浮点数7参

13、数转换中的y方向的旋转偏移（bursa类型）。（单位：弧度）不是所有类型都需要浮点数7参数转换中的z方向的旋转偏移（bursa类型）。（单位：弧度）不是所有类型都需要当已知一个坐标系的基准面和wgs84基准面之间的变换关系，就可以利用上面形式，把一个坐标系定义出来。下面是一个以7参数（布尔莎算法）为例的定义：datum_def dhdn desc_nm “deutsches hauptdreicknetz (dhdn)” source “german government” ellipsoid bessel use bursa delta_x 582.00000000000 delta_y 1

14、05.00000000000 delta_z 414.00000000000 bwscale 8.3000000000000 rot_x -1.0400000000000rot_y -0.35000000000000rot_z 3.0800000000000如果一个基准面已经在fme里面预先定义了，在以后所有的用到的地方都可以引用，不必重新定义，比如wgs84、beijing、xian80等（注意，fme预定义的beijing、xian80定义不是准确的基准面定义，详见前面的博客）。定义中要注意，在localcoordsysdefs.fme文件中最后一行是：include mycoordsys

15、defs.fme所有椭球体和基准面的定义，都要在这一行之前。定义完以上两个参数，就可以在mycoordsysdefs.fme中定义具体的坐标系参数（投影）了。坐标系定义：coordinate_system_def proj unit dt_name |el_name + desc_nm group quad source zero_x zero_y 参数描述如下：名称类型描述是否为空字符串投影坐标系名称no 字符串投影类型（如tm，ll等。具体见fme帮助）no 字符串坐标系来源（个人或者组织）yes字符串坐标单位no字符串分组名称，可以在gallery中进行分类yes字符串基准面名称，引用已

16、经定义的基准面名称yes字符串椭球体名称，引用已定义的名称。椭球体名称和基准面名称必须填写一个。yes 字符串+数字参数+数字（1-24），表示具体参数定义，如中央度线等，取决于投影类型（具体参见fme帮助）。no整型象限区代码yes 非负数设定一个x最小值的非零值yes非负数设定一个y最小值的非零值yes其中注意的是，投影类型的参数要根据需要定义。下面是个例子：coordinate_system_def utm12n83 proj tm unit meter dt_name nad83group custom group parm1 -111.0 scl_red 0.9996 org_lat

17、 0.0 x_off 500000.0 y_off 0.0 map_scl 1.0 zero_x 0.001 zero_y 0.001这个定义所选择的基准面为nad83，投影方式采用的是tm（transverse mercator），参数包括：parm1、org_lat、scl_red、map_scl、x_off、y_off（具体投影类型和参数见fme随机帮助）。简单总结一下，如果已知一个坐标系和wgs84坐标系的变换关系（比如7参数），那么就可以把这个坐标系定义到fme中，和其他已经定义的坐标系进行变换。但是，经常性的一个问题是，我们已知2个非wgs84的坐标系的7参数关系，如何在fme定义，并实现转换呢。例如，