BD09坐标系转换到CGCS2000坐标转换函数【JS】

   // BD09到CGCS2000坐标转换函数
        function bd09ToCgcs2000(bd_lng, bd_lat) {
            // BD09 -> WGS84 (通过GCJ02中转)
            var x_pi = 3.14159265358979324 * 3000.0 / 180.0;
            var x = bd_lng - 0.0065;
            var y = bd_lat - 0.006;
            var z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
            var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
            var gcj_lng = z * Math.cos(theta);
            var gcj_lat = z * Math.sin(theta);
            
            // GCJ02 -> WGS84 (火星坐标系转WGS84)
            var a = 6378245.0;
            var ee = 0.00669342162296594323;
            var dLat = transformLat(gcj_lng - 105.0, gcj_lat - 35.0);
            var dLng = transformLng(gcj_lng - 105.0, gcj_lat - 35.0);
            var radLat = gcj_lat / 180.0 * Math.PI;
            var magic = Math.sin(radLat);
            magic = 1 - ee * magic * magic;
            var sqrtMagic = Math.sqrt(magic);
            dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * Math.PI);
            dLng = (dLng * 180.0) / (a / sqrtMagic * Math.cos(radLat) * Math.PI);
            
            var wgs_lng = gcj_lng - dLng;
            var wgs_lat = gcj_lat - dLat;
            
            // WGS84 -> CGCS2000 (在中国境内差异很小，可以直接使用)
            return {lng: wgs_lng, lat: wgs_lat};
        }
        
        function transformLat(lng, lat) {
            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 * Math.PI) + 20.0 * Math.sin(2.0 * lng * Math.PI)) * 2.0 / 3.0;
            ret += (20.0 * Math.sin(lat * Math.PI) + 40.0 * Math.sin(lat / 3.0 * Math.PI)) * 2.0 / 3.0;
            ret += (160.0 * Math.sin(lat / 12.0 * Math.PI) + 320 * Math.sin(lat * Math.PI / 30.0)) * 2.0 / 3.0;
            return ret;
        }
        
        function transformLng(lng, lat) {
            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 * Math.PI) + 20.0 * Math.sin(2.0 * lng * Math.PI)) * 2.0 / 3.0;
            ret += (20.0 * Math.sin(lng * Math.PI) + 40.0 * Math.sin(lng / 3.0 * Math.PI)) * 2.0 / 3.0;
            ret += (150.0 * Math.sin(lng / 12.0 * Math.PI) + 300.0 * Math.sin(lng / 30.0 * Math.PI)) * 2.0 / 3.0;
            return ret;
        }

                           // 百度拾取的经纬度（BD09坐标系）
                            var bd_lng = 118.378381;  
                            var bd_lat = 35.101901;   
                            
                            // 转换为天地图坐标系（CGCS2000）
                            var converted = bd09ToCgcs2000(bd_lng, bd_lat);
                            var lng = converted.lng;
                            var lat = converted.lat;