gps84转换gcj02公式_WGS84-GCJ-02坐标转化
GCJ-02火星坐标系和WGS-84坐标系转换关系
WGS-84:GPS坐标系
GCJ-02:火星坐标系,国测局02年发布的坐标体系,高德,腾讯等使用。
BD-09:百度坐标系,百度自研,百度地图使用。
GCJ-02 转WGS-84
公式
下面的改进后的代码
package com.google.android.apps.location.gps.gnsslogger;
import java.math.BigDecimal;
import java.math.RoundingMode;
/**
* gps纠偏算法,适用于google,高德体系的地图
*/
public abstract class GpsUtilNew {
private final static double a = 6378245.0;
private final static double pi = 3.1415926535897932384626;
private final static double ee = 0.00669342162296594323;
/**
* 计算地球上任意两点(经纬度)距离
*
* @param lat1
* 第一点纬度
* @param lng1
* 第一点经度
* @param lat2
* 第二点纬度
* @param lng2
* 第二点经度
* @return 返回距离 单位:米
*/
public static double distance(double lat1, double lng1, double lat2, double lng2) {
double a, b, R;
R = 6378137; // 地球半径
lat1 = lat1 * Math.PI / 180.0;
lat2 = lat2 * Math.PI / 180.0;
a = lat1 - lat2;
b = (lng1 - lng2) * Math.PI / 180.0;
double d;
double sa2, sb2;
sa2 = Math.sin(a / 2.0);
sb2 = Math.sin(b / 2.0);
d = 2 * R * Math.asin(Math.sqrt(sa2 * sa2 + Math.cos(lat1) * Math.cos(lat2) * sb2 * sb2));
return d;
}
/**
* Description: WGS-84 to GCJ-02
*
* @author dsn
* @date 2017年10月24日 下午2:09:27
* @param latitude
* 纬度
* @param longitude
* 经度
* @return [纬度,经度]
* @version 1.0
*/
public static double[] toGCJ02Point(double latitude, double longitude) {
double[] dev = calDev(latitude, longitude);
double retLat = latitude + dev[0];
double retLon = longitude + dev[1];
return new double[] { retLat, retLon };
}
/**
* Description: WGS-84 to GCJ-02
*
* @author dsn
* @date 2017年10月24日 下午2:09:27
* @param latitude
* 纬度
* @param longitude
* 经度
* @param scale
* 经纬度保留小数位数
* @return [纬度,经度]
* @version 1.0
*/
public static double[] toGCJ02Point(double latitude, double longitude, int scale) {
double[] dev = calDev(latitude, longitude);
double retLat = latitude + dev[0];
double retLon = longitude + dev[1];
return new double[] { new BigDecimal(retLat).setScale(scale, RoundingMode.DOWN).doubleValue(),
new BigDecimal(retLon).setScale(scale, RoundingMode.DOWN).doubleValue() };
}
/**
* Description:GCJ-02 to WGS-84
*
* @author dsn
* @date 2017年10月24日 下午2:09:54
* @param latitude
* 纬度
* @param longitude
* 经度
* @return [纬度,经度]
* @version 1.0
*/
public static double[] toWGS84Point(double latitude, double longitude) {
double[] dev = calDev(latitude, longitude);
double retLat = latitude - dev[0];
double retLon = longitude - dev[1];
dev = calDev(retLat, retLon);
retLat = latitude - dev[0];
retLon = longitude - dev[1];
return new double[] { retLat, retLon };
}
private static double[] calDev(double wgLat, double wgLon) {
if (isOutOfChina(wgLat, wgLon)) {
return new double[] { 0, 0 };
}
double dLat = calLat(wgLon - 105.0, wgLat - 35.0);
double dLon = calLon(wgLon - 105.0, wgLat - 35.0);
double radLat = wgLat / 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);
dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
return new double[] { dLat, dLon };
}
private static boolean isOutOfChina(double lat, double lon) {
if (lon < 72.004 || lon > 137.8347)
return true;
if (lat < 0.8293 || lat > 55.8271)
return true;
return false;
}
private static double calLat(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) + 320 * Math.sin(y * pi / 30.0)) * 2.0 / 3.0;
return ret;
}
private static double calLon(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;
}
}
gps84转换gcj02公式_WGS84-GCJ-02坐标转化相关推荐
- gps84转换gcj02公式_【测绘基础】关于投影坐标和坐标转换,这几点你不得不知道!...
[测绘基础]关于投影坐标和坐标转换,这几点你不得不知道! 本文要点: 1.地理坐标系是什么?常见的地理坐标系有哪些? 2.投影坐标系是什么?常见的投影坐标系有哪些? 3.关于坐标转换. 4.如何在Ar ...
- gps84转换gcj02公式_百度坐标(BD09)、国测局坐标(火星坐标,GCJ02)、和WGS84坐标系之间的转换...
//定义一些常量 var x_PI = 3.14159265358979324 * 3000.0 / 180.0; var PI = 3.1415926535897932384626; var a = ...
- gps84转换gcj02公式_地球坐标系 (WGS-84) 相互转火星坐标系 (GCJ-02) 的转换算法
import java.util.HashMap; import java.util.Map; /** * 地球坐标系 (WGS-84) 相互转火星坐标系 (GCJ-02) 的转换算法 * * @au ...
- PHP 不同地图坐标系经纬度转换 GCj02 WGS84 BD-09
PHP 不同地图坐标系经纬度转换 GCj02 WGS84 BD-09 在项目中需要根据坐标计算距离,结果发现地图坐标不一致需要转换一下,在网上找到大佬写的一篇文章,亲测有效,原文地址: https:/ ...
- 题目内容: 写一个将华氏温度转换成摄氏温度的程序,转换的公式是: °F = (9/5)*°C + 32 其中C表示摄氏温度,F表示华氏温度。 程序的输入是一个整数,表示华氏温度。输出对
#题目内容: 写一个将华氏温度转换成摄氏温度的程序,转换的公式是: °F = (9/5)*°C + 32 其中C表示摄氏温度,F表示华氏温度. 程序的输入是一个整数,表示华氏温度.输出对应的摄氏温度, ...
- arcmap中有火星坐标码_在ArcMap中使用坐标值转换一个CAD文件中的坐标位置
(2013-04-13 20:51:19) 转 载 ▼ 标签: 股票 分类: GIS 技术研究 您现在的位置 : 首页 -> 知识库 -> ArcGIS Desktop -> 技术文 ...
- 【10位时间戳转换日期公式】
10位时间戳转换日期公式 时间戳TimeStamp是时间序列中的最基本的数据类型,Pandas可以将整型或浮点型表示的时间转换为时间戳.程序如下: import pandas as pd print( ...
- 写一个将华氏温度转换成摄氏温度的程序,转换的公式是:°F = (9/5)*°C + 32 其中C表示摄氏温度F表示华氏温度程序的输入是一个整数表示华氏温度。输出对应的摄氏温度也是一个整数(java)
1温度转换(5分) (此题出自:中国大学MOOC零基础学java,翁恺,第1周编程题) 题目内容: 写一个将华氏温度转换成摄氏温度的程序,转换的公式是: °F = (9/5)*°C + 32 其中C表 ...
- 将Materials Studio导出的pdb文件转换成LAMMPS所用的data坐标文件(含程序)
将Materials Studio导出的pdb文件转换成LAMMPS所用的data坐标文件的程序使用说明 程序在此 1.程序可以识别的文件格式和输出的文件格式示例如下: 输入的PDB文件应该包含模拟边 ...
- 点云处理:Kitti数据集中点云坐标转化到彩色摄像机坐标
一.实现目标 激光雷达探测到的点云目标,将目标的激光雷达下的坐标转化到摄像机的坐标下. 实现结果: 雷达检测到车辆目标. 找到对应的2D坐标点,并框出目标. 二.实现原理 其中,x为点云坐标系坐标,y ...
最新文章
- 财务学习mysql_从零开始学 MySQL - 前言
- CIKM 2020 | 知识库问答复杂问题的分层查询图生成方法
- SpringCloud熔断器介绍
- 为什么不可以使用哈曼顿距离_请对比下欧式距离和曼哈顿距离的差别
- git push时报错fatal: Could not read from remote repository.
- PHP操作MongoDB GridFS 存储文件
- 微型通用程序框架-规划
- java的Teacher类_java类的继承 - osc_tauwfamo的个人空间 - OSCHINA - 中文开源技术交流社区...
- swfobject.js 2.2简单使用方法
- 为LogAnalyzer添加源IP
- python tkinter界面 控件位置随窗体自动调整适应_如何实现python tkinter控件随窗口变化?...
- 9. Browser 对象 - Location 对象
- [Python+Django]Web图书管理系统毕业设计之源码+论文篇
- python 判断素数
- 第九组 通信3班 063 防火墙配置
- 用python做一个木马_Python编程简单的木马程序(转载于乌云中)
- 江苏省高等学校计算机等级考试操作题,江苏省高等学校计算机等级考试09春一级操作题汇总...
- Pid控制算法-专家PID与模糊PID的C++实现
- Honor—Python note(2)——PyCharm创建新项目
- 浏览器页面性能分析指南(chrome)
热门文章
- java 处理汉字枚举类_浅谈java枚举类(附代码)
- spurious wakeup 虚假唤醒
- 【小家java】JUC并发编程之:虚假唤醒(spurious wakeup)以及推荐的解决方案
- 怎么让模糊的数字变清楚_Re: 请问将模糊的图片中数字变的清晰方法是?多谢!!!...
- Face2Face: Real-time Face Capture and Reenactment of RGB Videos 论文翻译
- 腾讯云即时通讯IM_ 小贴士
- mysql中怎样把字段名改为中文_mysql修改字段名-Mysql,修改
- DHSNet: Deep Hierarchical Saliency Network for Salient Object Detection阅读总结
- 计算机word题目设置标题,word文档标题怎么设置?word设置标题样式的方法
- arduinouno的地是相连的吗_垫圈锁紧,靠谱吗?