ChangePoint.js

/*
 * @Author: 叶思铧(292184)
 * @Date: 2022-04-21 11:02:56
 * @Descripttion:坐标转换方法
 * @LastEditors: yesihua
 * @LastEditTime: 2022-08-25 11:35:44
 */
/**
 *  @description  不同类型地图坐标系转换
 *  @param   point   点数
 *  @returns point
 *  */
let ol = window.ol

class ChangePoint {
    constructor(mapObj) {
        this.mapObj = mapObj
        // this.resourceTree = new ResourceTree()
        // // console.log(opts, '测试')
        // // console.log(ol, '测试')
        // this._init(domId, opts)
        // //待定
        // return this
    }
    //默认坐标转换
    transformation(point) {
        return this._changePoint(point)
    }
    //默认的坐标反转
    reversal(point) {
        let mapGlobalData = this.mapObj.mapGlobalData
        if (mapGlobalData.projection == 'EPSG:3857') {
            point = ol.proj.transform(point, 'EPSG:3857', 'EPSG:4326')
        }

        if (mapGlobalData.mapType == 'BAIDU' || mapGlobalData.mapType == 'baidu') {
            point = coordinateConversion.bd09togcj02(point[0], point[1])
            point = coordinateConversion.gcj02towgs84(point[0], point[1])
        } else if (mapGlobalData.mapType == 'amap' || mapGlobalData.mapType == 'aMap') {
            point = coordinateConversion.gcj02towgs84(point[0], point[1])
        }

        //不同厂商之间的坐标系转换
        return point
    }

    _changePoint = function (point) {
        let mapGlobalData = this.mapObj.mapGlobalData
        if (mapGlobalData.mapType == 'BAIDU' || mapGlobalData.mapType == 'baidu' || mapGlobalData.mapType == 'amap' || mapGlobalData.mapType == 'aMap' || mapGlobalData.projection == 'EPSG:3857') {
            if (mapGlobalData.mapType == 'aMap' || mapGlobalData.mapType == 'amap') {
                point = coordinateConversion.wgs84togcj02(point[0], point[1])
                point = ol.proj.transform(point, 'EPSG:4326', 'EPSG:3857')
            } else if (mapGlobalData.mapType == 'BAIDU' || mapGlobalData.mapType == 'baidu') {
                point = coordinateConversion.wgs84togcj02(point[0], point[1])
                point = coordinateConversion.gcj02tobd09(point[0], point[1])
                point = ol.proj.transform(point, 'EPSG:4326', 'EPSG:3857')
            } else {
                point = ol.proj.transform(point, 'EPSG:4326', 'EPSG:3857')
            }
        }
        // console.log(point, '测试')
        return point
    }
    wgs84togcj02(lng, lat) {
        return coordinateConversion.wgs84togcj02(lng, lat)
    }
    gcj02tobd09(lng, lat) {
        return coordinateConversion.gcj02tobd09(lng, lat)
    }
    bd09togcj02(lng, lat) {
        return coordinateConversion.bd09togcj02(lng, lat)
    }
    gcj02towgs84(lng, lat) {
        return coordinateConversion.gcj02towgs84(lng, lat)
    }
}
//不同坐标系坐标转换算法
const coordinateConversion = {
    x_PI: (3.14159265358979324 * 3000.0) / 180.0,
    PI: 3.1415926535897932384626,
    a: 6378245.0,
    ee: 0.00669342162296594323,

    /**
     * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
     * 即 百度 转 谷歌、高德
     * @param bd_lon
     * @param bd_lat
     * @returns {*[]}
     */
    bd09togcj02: function (bd_lon, bd_lat) {
        bd_lon = +bd_lon
        bd_lat = +bd_lat
        var x = bd_lon - 0.0065
        var y = bd_lat - 0.006
        var z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * coordinateConversion.x_PI)
        var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * coordinateConversion.x_PI)
        var gg_lng = z * Math.cos(theta)
        var gg_lat = z * Math.sin(theta)
        return [gg_lng, gg_lat]
    },
    /**
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
     * 即谷歌、高德 转 百度
     * @param lng
     * @param lat
     * @returns {*[]}
     */
    gcj02tobd09: function (lng, lat) {
        lat = +lat
        lng = +lng
        var z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * coordinateConversion.x_PI)
        var theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * coordinateConversion.x_PI)
        var bd_lng = z * Math.cos(theta) + 0.0065
        var bd_lat = z * Math.sin(theta) + 0.006
        return [bd_lng, bd_lat]
    },
    /**
     * WGS84转GCj02
     * @param lng
     * @param lat
     * @returns {*[]}
     */
    wgs84togcj02: function (lng, lat) {
        lat = +lat
        lng = +lng
        if (coordinateConversion.outOfChina(lng, lat)) {
            return [lng, lat]
        } else {
            var dlat = coordinateConversion.transformlat(lng - 105.0, lat - 35.0)
            var dlng = coordinateConversion.transformlng(lng - 105.0, lat - 35.0)
            var radlat = (lat / 180.0) * coordinateConversion.PI
            var magic = Math.sin(radlat)
            magic = 1 - coordinateConversion.ee * magic * magic
            var sqrtmagic = Math.sqrt(magic)
            dlat = (dlat * 180.0) / (((coordinateConversion.a * (1 - coordinateConversion.ee)) / (magic * sqrtmagic)) * coordinateConversion.PI)
            dlng = (dlng * 180.0) / ((coordinateConversion.a / sqrtmagic) * Math.cos(radlat) * coordinateConversion.PI)
            var mglat = lat + dlat
            var mglng = lng + dlng
            return [mglng, mglat]
        }
    },
    /**
     * GCJ02 转换为 WGS84
     * @param lng
     * @param lat
     * @returns {*[]}
     */
    gcj02towgs84: function (lng, lat) {
        lat = +lat
        lng = +lng
        if (coordinateConversion.outOfChina(lng, lat)) {
            return [lng, lat]
        } else {
            var dlat = coordinateConversion.transformlat(lng - 105.0, lat - 35.0)
            var dlng = coordinateConversion.transformlng(lng - 105.0, lat - 35.0)
            var radlat = (lat / 180.0) * coordinateConversion.PI
            var magic = Math.sin(radlat)
            magic = 1 - coordinateConversion.ee * magic * magic
            var sqrtmagic = Math.sqrt(magic)
            dlat = (dlat * 180.0) / (((coordinateConversion.a * (1 - coordinateConversion.ee)) / (magic * sqrtmagic)) * coordinateConversion.PI)
            dlng = (dlng * 180.0) / ((coordinateConversion.a / sqrtmagic) * Math.cos(radlat) * coordinateConversion.PI)
            var mglat = lat + dlat
            var mglng = lng + dlng
            return [lng * 2 - mglng, lat * 2 - mglat]
        }
    },
    transformlat: function (lng, lat) {
        lat = +lat
        lng = +lng
        var 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 * coordinateConversion.PI) + 20.0 * Math.sin(2.0 * lng * coordinateConversion.PI)) * 2.0) / 3.0
        ret += ((20.0 * Math.sin(lat * coordinateConversion.PI) + 40.0 * Math.sin((lat / 3.0) * coordinateConversion.PI)) * 2.0) / 3.0
        ret += ((160.0 * Math.sin((lat / 12.0) * coordinateConversion.PI) + 320 * Math.sin((lat * coordinateConversion.PI) / 30.0)) * 2.0) / 3.0
        return ret
    },
    transformlng: function (lng, lat) {
        lat = +lat
        lng = +lng
        var 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 * coordinateConversion.PI) + 20.0 * Math.sin(2.0 * lng * coordinateConversion.PI)) * 2.0) / 3.0
        ret += ((20.0 * Math.sin(lng * coordinateConversion.PI) + 40.0 * Math.sin((lng / 3.0) * coordinateConversion.PI)) * 2.0) / 3.0
        ret += ((150.0 * Math.sin((lng / 12.0) * coordinateConversion.PI) + 300.0 * Math.sin((lng / 30.0) * coordinateConversion.PI)) * 2.0) / 3.0
        return ret
    },
    /**
     * 判断是否在国内，不在国内则不做偏移
     * @param lng
     * @param lat
     * @returns {boolean}
     */
    outOfChina: function (lng, lat) {
        lat = +lat
        lng = +lng
        // 纬度3.86~53.55,经度73.66~135.05
        return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55)
    },
}
window.coordinateConversion = coordinateConversion
export default ChangePoint