proj4js 坐标转换

时间：2021-07-02 15:29:04 收藏：0 阅读：0

一、概述

Proj4js 是一个开源的 JavaScript 库，用于将点坐标从一个坐标系转换到另一个坐标系，包括基准转换。

git代码库地址：https://github.com/proj4js/proj4js

另一个坐标系在线查询和坐标转换地址：https://epsg.io/

在PostGIS中有一个表 spatial_ref_sys ，可以查询所有的坐标系信息。

在Geoserver中也有所有坐标系的信息。在ArcGIS中也有。

以下介绍在node中使用proj4，以及在Cesium三维开发中使用。

技术图片

二、node中使用Proj4

1、引入库，npm install proj4

2、引用库，const proj4 = require(‘proj4‘)

3、转换方法，proj4(fromProjection, toProjection, [x, y])

4、坐标系举例说明

EPSG:4549 CGCS2000 / 3-degree Gauss-Kruger CM 120E

2000大地坐标系，3度分带，高斯克吕格投影，中央经线120。

如果坐标是加了带号的，就用4528。具体看需要转换的数据坐标。加了带号的X坐标前面两位就是带号。

东西向X坐标不加带号6位数（坐标向东偏移，了500km，因此中央子午线的X坐标就是500000，加了带号就是365000000），加了带号8位数。南北方向7位数

EPSG:4528 CGCS2000 / 3-degree Gauss-Kruger zone 40

3度带中央经线=3乘以带号

技术图片

5、代码，在下面的代码中实现将点坐标从CGCS2000坐标系EPSG:4549转换为WGS84坐标系

函数convert111_和convert111结果相同。因为proj4函数可以传入坐标系的wkt，也可以传入proj.4。

/*
 * @Author: 苹果园dog
 * @Date: 2021-07-01 11:48:07
 * @LastEditTime: 2021-07-01 17:44:56
 * @LastEditors: Please set LastEditors
 * @Description: In User Settings Edit
 * @FilePath: 
 */
const proj4 = require(‘proj4‘);
function convert111_(x, y) {
  const fromProjection = `PROJCS["CGCS2000 / 3-degree Gauss-Kruger CM 120E",GEOGCS["China Geodetic Coordinate System 2000",DATUM["China_2000",SPHEROID["CGCS2000",6378137,298.257222101,AUTHORITY["EPSG","1024"]],AUTHORITY["EPSG","1043"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4490"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",120],PARAMETER["scale_factor",1],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","4549"]]
  `;
  const toProjection = `GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG", "6326"]], PRIMEM["Greenwich", 0, AUTHORITY["EPSG", "8901"]],
        UNIT["degree", 0.0174532925199433, AUTHORITY["EPSG", "9122"]], AUTHORITY["EPSG", "4326"]]`;
  const newCoordinates = proj4(fromProjection, toProjection, [x, y]);
  return newCoordinates;
}

function convert111(x, y) {
  const fromProjection = `+proj=tmerc +lat_0=0 +lon_0=120 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs`;
  const toProjection = `+proj=longlat +datum=WGS84 +no_defs`;
  const newCoordinates = proj4(fromProjection, toProjection, [x, y]);
  return newCoordinates;
}

module.exports = {
  convert111: convert111
}

6、坐标转换测试

三、Cesium中使用Proj4

由于Cesium三维球默认是WGS84的，要想加载一些CGCS2000坐标系的一些矢量文件。

比如geojson，kml等，需要在加载数据的时候动态转换坐标。

以下介绍加载2000坐标系的Geojson方法。

在下面的代码中实现将点坐标从CGCS2000坐标系EPSG:4524转换为WGS84坐标系。

/*
 * @Author: 苹果园dog
 * @Date: 2021-04-12 14:42:01
 * @LastEditTime: 2021-04-12 14:46:41
 * @LastEditors: Please set LastEditors
 * @Description: In User Settings Edit
 * @FilePath: 
 */
import { default as proj4 } from "proj4";

Cesium.GeoJsonDataSource.crsNames[
  "urn:ogc:def:crs:EPSG::4524"
] = Cesium.GeoJsonDataSource.crsNames["EPSG:4524"] = function(coordinates) {
  const fromProj = `PROJCS["CGCS2000 / 3-degree Gauss-Kruger zone 36", 
        GEOGCS["China Geodetic Coordinate System 2000", 
          DATUM["China 2000", 
            SPHEROID["CGCS2000", 6378137.0, 298.257222101, AUTHORITY["EPSG","1024"]], 
            AUTHORITY["EPSG","1043"]], 
          PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], 
          UNIT["degree", 0.017453292519943295], 
          AXIS["Geodetic longitude", EAST], 
          AXIS["Geodetic latitude", NORTH], 
          AUTHORITY["EPSG","4490"]], 
        PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]], 
        PARAMETER["central_meridian", 108.0], 
        PARAMETER["latitude_of_origin", 0.0], 
        PARAMETER["scale_factor", 1.0], 
        PARAMETER["false_easting", 36500000.0], 
        PARAMETER["false_northing", 0.0], 
        UNIT["m", 1.0], 
        AXIS["Easting", EAST], 
        AXIS["Northing", NORTH], 
        AUTHORITY["EPSG","4524"]]`;
  const toProj = `GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG", "6326"]], PRIMEM["Greenwich", 0, AUTHORITY["EPSG", "8901"]],
        UNIT["degree", 0.0174532925199433, AUTHORITY["EPSG", "9122"]], AUTHORITY["EPSG", "4326"]]`;

  const x = coordinates[0];
  const y = coordinates[1];

  const newCoordinates = proj4(fromProj, toProj, [x, y]);
  return Cesium.Cartesian3.fromDegrees(newCoordinates[0], newCoordinates[1], 0);
};