public class Gps {
private double wgLat;
private double wgLon;
public Gps(double wgLat, double wgLon) {
setWgLat(wgLat);
setWgLon(wgLon);
}
public double getWgLat() {
return wgLat;
}
public void setWgLat(double wgLat) {
this.wgLat = wgLat;
}
public double getWgLon() {
return wgLon;
}
public void setWgLon(double wgLon) {
this.wgLon = wgLon;
}
@Override
public String toString() {
return wgLat + "," + wgLon;
}
}
  1 /**
  2  * 各地图API坐标系统比较与转换;
  3  * WGS84坐标系:即地球坐标系,国际上通用的坐标系。设备一般包含GPS芯片或者北斗芯片获取的经纬度为WGS84地理坐标系,
  4  * 谷歌地图采用的是WGS84地理坐标系(中国范围除外);
  5  * GCJ02坐标系:即火星坐标系,是由中国国家测绘局制订的地理信息系统的坐标系统。由WGS84坐标系经加密后的坐标系。
  6  * 谷歌中国地图和搜搜中国地图采用的是GCJ02地理坐标系; BD09坐标系:即百度坐标系,GCJ02坐标系经加密后的坐标系;
  7  * 搜狗坐标系、图吧坐标系等,估计也是在GCJ02基础上加密而成的。 chenhua
  8  */
  9 public class PositionUtil {
 10     
 11     public static final String BAIDU_LBS_TYPE = "bd09ll";
 12     
 13     public static double pi = 3.1415926535897932384626;
 14     public static double a = 6378245.0;
 15     public static double ee = 0.00669342162296594323;
 16 
 17     /**
 18      * 84 to 火星坐标系 (GCJ-02) World Geodetic System ==> Mars Geodetic System
 19      * 
 20      * @param lat
 21      * @param lon
 22      * @return
 23      */
 24     public static Gps gps84_To_Gcj02(double lat, double lon) {
 25         if (outOfChina(lat, lon)) {
 26             return null;
 27         }
 28         double dLat = transformLat(lon - 105.0, lat - 35.0);
 29         double dLon = transformLon(lon - 105.0, lat - 35.0);
 30         double radLat = lat / 180.0 * pi;
 31         double magic = Math.sin(radLat);
 32         magic =  - ee * magic * magic;
 33         double sqrtMagic = Math.sqrt(magic);
 34         dLat = (dLat * 180.0) / ((a * ( - ee)) / (magic * sqrtMagic) * pi);
 35         dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
 36         double mgLat = lat + dLat;
 37         double mgLon = lon + dLon;
 38         return new Gps(mgLat, mgLon);
 39     }
 40 
 41     /**
 42      * * 火星坐标系 (GCJ-02) to 84 * * @param lon * @param lat * @return
 43      * */
 44     public static Gps gcj_To_Gps84(double lat, double lon) {
 45         Gps gps = transform(lat, lon);
 46         double lontitude = lon *  - gps.getWgLon();
 47         double latitude = lat *  - gps.getWgLat();
 48         return new Gps(latitude, lontitude);
 49     }
 50 
 51     /**
 52      * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 将 GCJ-02 坐标转换成 BD-09 坐标
 53      * 
 54      * @param gg_lat
 55      * @param gg_lon
 56      */
 57     public static Gps gcj02_To_Bd09(double gg_lat, double gg_lon) {
 58         double x = gg_lon, y = gg_lat;
 59         double z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * pi);
 60         double theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * pi);
 61         double bd_lon = z * Math.cos(theta) + 0.0065;
 62         double bd_lat = z * Math.sin(theta) + 0.006;
 63         return new Gps(bd_lat, bd_lon);
 64     }
 65 
 66     /**
 67      * * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 * * 将 BD-09 坐标转换成GCJ-02 坐标 * * @param
 68      * bd_lat * @param bd_lon * @return
 69      */
 70     public static Gps bd09_To_Gcj02(double bd_lat, double bd_lon) {
 71         double x = bd_lon - 0.0065, y = bd_lat - 0.006;
 72         double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * pi);
 73         double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * pi);
 74         double gg_lon = z * Math.cos(theta);
 75         double gg_lat = z * Math.sin(theta);
 76         return new Gps(gg_lat, gg_lon);
 77     }
 78 
 79     /**
 80      * (BD-09)-->84
 81      * @param bd_lat
 82      * @param bd_lon
 83      * @return
 84      */
 85     public static Gps bd09_To_Gps84(double bd_lat, double bd_lon) {
 86 
 87         Gps gcj02 = PositionUtil.bd09_To_Gcj02(bd_lat, bd_lon);
 88         Gps map84 = PositionUtil.gcj_To_Gps84(gcj02.getWgLat(),
 89                 gcj02.getWgLon());
 90         return map84;
 91 
 92     }
 93 
 94     public static boolean outOfChina(double lat, double lon) {
 95         if (lon < 72.004 || lon > 137.8347)
 96             return true;
 97         if (lat < 0.8293 || lat > 55.8271)
 98             return true;
 99         return false;
     }
 
     public static Gps transform(double lat, double lon) {
         if (outOfChina(lat, lon)) {
             return new Gps(lat, lon);
         }
         double dLat = transformLat(lon - 105.0, lat - 35.0);
         double dLon = transformLon(lon - 105.0, lat - 35.0);
         double radLat = lat / 180.0 * pi;
         double magic = Math.sin(radLat);
         magic =  - ee * magic * magic;
         double sqrtMagic = Math.sqrt(magic);
         dLat = (dLat * 180.0) / ((a * ( - ee)) / (magic * sqrtMagic) * pi);
         dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
         double mgLat = lat + dLat;
         double mgLon = lon + dLon;
         return new Gps(mgLat, mgLon);
     }
 
     public static double transformLat(double x, double y) {
         double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y
                 + 0.2 * Math.sqrt(Math.abs(x));
         ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
         ret += (20.0 * Math.sin(y * pi) + 40.0 * Math.sin(y / 3.0 * pi)) * 2.0 / 3.0;
         ret += (160.0 * Math.sin(y / 12.0 * pi) +  * Math.sin(y * pi / 30.0)) * 2.0 / 3.0;
         return ret;
     }
 
     public static double transformLon(double x, double y) {
         double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1
                 * Math.sqrt(Math.abs(x));
         ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;
         ret += (20.0 * Math.sin(x * pi) + 40.0 * Math.sin(x / 3.0 * pi)) * 2.0 / 3.0;
         ret += (150.0 * Math.sin(x / 12.0 * pi) + 300.0 * Math.sin(x / 30.0
                 * pi)) * 2.0 / 3.0;
         return ret;
     }
 
     public static void main(String[] args) {
 
         // 北斗芯片获取的经纬度为WGS84地理坐标 31.426896,119.496145
         Gps gps = new Gps(31.426896, 119.496145);
         System.out.println("gps :" + gps);
         Gps gcj = gps84_To_Gcj02(gps.getWgLat(), gps.getWgLon());
         System.out.println("gcj :" + gcj);
         Gps star = gcj_To_Gps84(gcj.getWgLat(), gcj.getWgLon());
         System.out.println("star:" + star);
         Gps bd = gcj02_To_Bd09(gcj.getWgLat(), gcj.getWgLon());
         System.out.println("bd  :" + bd);
         Gps gcj2 = bd09_To_Gcj02(bd.getWgLat(), bd.getWgLon());
         System.out.println("gcj :" + gcj2);
     }

152 }

WGS84,GCJ02, BD09坐标转换的更多相关文章

  1. GIS坐标系测绘原理:大地水准面/基准面/参考椭球体/EPSG/SRI/WKT

    预热文章系列:<GIS历史概述与WebGis应用开发技术浅解>.<GIS坐标系:WGS84,GCJ02,BD09,火星坐标,大地坐标等解析说与转换>.<OGC标准WMTS ...

  2. BD09坐标(百度坐标) WGS84(GPS坐标) GCJ02(国测局坐标) 的相互转换

    BD09坐标(百度坐标) WGS84(GPS坐标) GCJ02(国测局坐标) 的相互转换 http://www.cnphp6.com/archives/24822 by root ⋅ Leave a ...

  3. 坐标转换——GCJ-02

    WGS84(World Geodetic System 1984),是为GPS 全球定位系统 使用而建立的坐标系统GCJ-02,我国在WGS84的基础上加密得到BD-09,百度坐标在GCJ-02基础上 ...

  4. wgs84坐标系与gcj02坐标系转换误差分布图 | Mapping the Error in Transformation between WGS84 and GCJ02 Coordinations

    国际上通用的是wgs84坐标系,而我国对于境内的坐标进行了加密,采用了gcj02坐标系,或者称为火星坐标系.亢孟军老师带的一门课<多媒体电子地图设计>要求我们从wgs84坐标系转换为gcj ...

  5. java中WGS84坐标(ios)转换BD-09坐标(百度坐标)

    iPhone的GPS定位(CLLocationManager)获得的经纬坐标是基于WGS-84坐标系(世界标准),Google地图使用的是GCJ-02坐标系(中国特色的火星坐标系),百度的经纬坐标在G ...

  6. 简谈百度坐标反转至WGS84的三种思路

    文章版权由作者李晓晖和博客园共有,若转载请于明显处标明出处:http://www.cnblogs.com/naaoveGIS/ 1.背景 基于百度地图进行数据展示是目前项目中常见场景,但是因为百度地图 ...

  7. 各类坐标系相互之间的转换(84互转GC02,GC02互转BD09)

    在遥感行业我们经常会用到各类的坐标系相互之间的转换,常见的度分秒转化为度很简单,直接上代码: //经纬度 ////118度48分54.152秒=118+(48/60)+(54.152/3600)=11 ...

  8. 【PHP版】火星坐标系 (GCJ-02) 与百度坐标系 (BD-09ll)转换算法

    首先感谢java版作者@宋宋宋伟,java版我是看http://blog.csdn.net/coolypf/article/details/8569813 然后根据java代码修改成了php代码. & ...

  9. 使用ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标的转换

    [摘 要] 本文针对从事测绘工作者普遍遇到的坐标转换问题,简要介绍ArcGIS实现WGS84经纬度坐标到北京54高斯投影坐标转换原理和步骤. [关键词] ArcGIS 坐标转换 投影变换 1 坐标转换 ...

随机推荐

  1. IndentationError 这个错误是缩进的问题

    1.自定义库导入时,提示这个错误,是缩进的问题

  2. object detection[SSD]

    0. 背景 经过了rcnn,spp,fast rcnn, faster rcnn,yolo,这里又到了ssd模型. faster rcnn的贡献是将候选框区域提取的部分也集成到CNN中去,并且与对象的 ...

  3. A2D JS框架 - DES加密解密 与 Cookie的封装(C#与js互相加密解密)

    这次实现了JS端的DES加密与解密,并且C#端也能正确解析DES的密文(反之也实现了) 使用的代码如下,非常方便: <script src="A2D.js" type=&qu ...

  4. 2019 年起如何开始学习 ABP 框架系列文章-开篇有益

    2019 年起如何开始学习 ABP 框架系列文章-开篇有益 [[TOC]] 本系列文章推荐阅读地址为:52ABP 开发文档 https://www.52abp.com/Wiki/52abp/lates ...

  5. 用Flask+Redis维护代理池

    GitHub:https://github.com/LXL-YAN/ProxyPool 视频讲解:https://www.bilibili.com/video/av19057145/?p=17

  6. set和multiset的用法

    set #include<iostream> #include<algorithm> #include<cstdio> #include<string.h&g ...

  7. 1000/problem/A

    传送门: [http://codeforces.com/contest/1000/problem/A] 题意: 一个比赛颁奖,要准备T-Shirt给获奖者,但有的去年获奖过,衣服尺寸可以不改,有的需要 ...

  8. Java工具类——UUIDUtils

    借用一下百度百科的解释,来看一下UUID是什么. UUID含义是通用唯一识别码 (Universally Unique Identifier),这 是一个软件建构的标准,也是被开源软件基金会 (Ope ...

  9. c++局部变量在外可用的方法

    C++的局部变量在作用域结束后,一般都会被回收.如下面这段代码 map<a*, b*> _map; void fun() { a _a; b _b; _map[&_a] = &am ...

  10. #Leetcode# 989. Add to Array-Form of Integer

    https://leetcode.com/problems/add-to-array-form-of-integer/ For a non-negative integer X, the array- ...