GIS地图工作记录：1、坐标系

坐标系

1、地理坐标系和投影坐标系

地理坐标系：使用三维球面来定义地球表面位置，以实现通过经纬度对地球表面点位引用的坐标系。

1. 一个地理坐标系包括角度测量单位、本初子午线和参考椭球体三部分。
2. 地理坐标系依据其所选用的本初子午线、参考椭球的不同而略有区别。
3. 参考平面地是椭球面，地理坐标系的坐标单位为经纬度

投影坐标系 ：使用基于X,Y值的坐标系统来描述地球上某个点所处的位置。这个坐标系是从地球的近似椭球型投影得到的，它对应于某个地理坐标系。

1. 投影坐标系的生成是以地理坐标系为基准的，所以每个投影坐标系前面都会挂有地理坐标系。
2. 参考平面地是水平面，平面坐标系统地图单位通常为米

两者区别：

1. 地理坐标系统是一种球面坐标，而投影坐标系统是平面坐标

2. 投影坐标系统在二维平面上有恒定的长度、角度和面积

3. 投影坐标系统是由地理坐标投影到平面坐标上形成的

上述内容可忽略，很多内容还在学习整理中

---

2.常用坐标系

WGS84：一种国际上采用的地心坐标系。坐标原点为地球质心，其地心空间直角坐标系的Z轴指向BIH （国际时间服务机构）1984.O定义的协议地球极（CTP)方向，X轴指向BIH 1984.0的零子午面和CTP赤道的交点，Y轴与Z轴、X轴垂直构成右手坐标系，称为1984年世界大地坐标系统。

GCJ02：由中国国家测绘局（G表示Guojia国家，C表示Cehui测绘，J表示Ju局）制订的地理信息系统的坐标系统。它其实就是对真实坐标系统进行人为的加偏处理，按照特殊的算法，将真实的坐标加密成虚假的坐标，而这个加偏并不是线性的加偏，所以各地的偏移情况都会有所不同。而加密后的坐标也常被大家称为“火星坐标系统”。

BD09：经纬度投影属于百度坐标系，它是在标准经纬度的基础上进行GCJ-02加偏之后，再加上百度自身的加偏算法，也就是在标准经纬度的基础之上进行了两次加偏。

CGCS2000：2000中国大地坐标系(China Geodetic Coordinate System 2000，CGCS2000)，又称之为2000国家大地坐标系，是中国新一代大地坐标系，21世纪初已在中国正式实施。

1. GCJ02坐标系是在WGS84坐标系进行非线性加密计算形成的
2. BD09坐标系是在GCJ02坐标系进行非线性加密计算形成的

 

地图	坐标系
高德地图、天地图	GCJ02
百度地图	BD09
谷歌地图	BD09（国内）、WGS84

 

 

 

 

 

坐标系转换的一些代码

public class MapCoordinateChangeUtil {

    //定义一些常量
    public static double x_PI = 3.14159265358979324 * 3000.0 / 180.0;
    public static double  PI = 3.1415926535897932384626;
    public static double  a = 6378245.0;
    public static double  ee = 0.00669342162296594323;
    /**
     * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
     * 即 百度 转 谷歌、高德
     * @param bd_lon
     * @param bd_lat
     * @returns {*[]}
     */
    public static  String bd09togcj02(float bd_lon, float bd_lat) {
        double x_pi = 3.14159265358979324 * 3000.0 / 180.0;
        double x = bd_lon - 0.0065;
        double y = bd_lat - 0.006;
        double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
        double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
        double gg_lng = z * Math.cos(theta);
        double gg_lat = z * Math.sin(theta);
        return gg_lng+","+gg_lat;
    }
    /**
     * GCJ02 转换为 WGS84
     * @param lng
     * @param lat
     * @returns {*[]}
     */
    public static String  gcj02towgs84(double lng, double lat) {
        if (out_of_china(lng, lat)) {
            return lng+","+lat;
        }
        else {
            double dlat = transformlat(lng - 105.0, lat - 35.0);
            double dlng = transformlng(lng - 105.0, lat - 35.0);
            double radlat = lat / 180.0 * PI;
            double magic = Math.sin(radlat);
            magic = 1 - ee * magic * magic;
            double sqrtmagic = Math.sqrt(magic);
            dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI);
            dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
            double mglat  = lat + dlat;
            double mglng  = lng + dlng;
            return (lng * 2 - mglng)+ ","+ (lat * 2 - mglat);
        }
    }
    public static double transformlat(double lng, double lat) {
        double ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng));
        ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(lat * PI) + 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0;
        ret += (160.0 * Math.sin(lat / 12.0 * PI) + 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0;
        return ret;
    }
    public static double  transformlng(double lng,double lat) {
        double ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));
        ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(lng * PI) + 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(lng / 12.0 * PI) + 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0;
        return ret;
    }
    /**
     * 判断是否在国内，不在国内则不做偏移
     * @param lng
     * @param lat
     * @returns {boolean}
     */
    public static boolean out_of_china(double lng, double lat) {
        return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false);
    }
}