WGS ddmm.mmmm 与度 互相转化

2024-06-13 11:20
文章标签 互相转化 wgs mmmm ddmm 与度

本文主要是介绍WGS ddmm.mmmm 与度 互相转化，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

1、搜索半天都是WEG转为度的，很少有度转为WGS的这里我完善了下

public class GpsToUtil {private static double pi = 3.14159265358979324;private static double a = 6378245.0;private static double ee = 0.00669342162296594323;private final static double x_pi = 3.14159265358979324 * 3000.0 / 180.0;/*** WGS转2GCJ* @param lat* @param lon* @return* @Author : pjh. create at 2021年1月21日 上午10:42:42*/public static double[] wgs2gcj(double lat, double 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 = 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);double mgLat = lat + dLat;double mgLon = lon + dLon;double[] loc = { mgLat, mgLon };return loc;}/*** WGS先转为dd.dddd 再转为2gcj* @param latStr  DDFF.FFFF* @param lonStr  DDFF.FFFF* @return   latStr   lonStr* @Author : pjh. create at 2021年1月21日 上午10:42:57*/public static double[] wgs2gcj_dddddd(String latStr, String lonStr) {String wmm = latStr.substring(latStr.indexOf(".") - 2);String jmm = lonStr.substring(lonStr.indexOf(".") - 2);String wdd = latStr.substring(0, latStr.indexOf(".") - 2);String jdd = lonStr.substring(0, lonStr.indexOf(".") - 2);double wmm_a = Double.parseDouble(wmm) / 60;double jmm_a = Double.parseDouble(jmm) / 60;// 转成dd.dddd格式double lat = Double.parseDouble(wdd) + wmm_a;double lon = Double.parseDouble(jdd) + jmm_a;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 = 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);double mgLat = lat + dLat;double mgLon = lon + dLon;double[] loc = { mgLat, mgLon };return loc;}private static double transformLat(double lat, double lon) {double ret = -100.0 + 2.0 * lat + 3.0 * lon + 0.2 * lon * lon + 0.1 * lat * lon+ 0.2 * Math.sqrt(Math.abs(lat));ret += (20.0 * Math.sin(6.0 * lat * pi) + 20.0 * Math.sin(2.0 * lat * pi)) * 2.0 / 3.0;ret += (20.0 * Math.sin(lon * pi) + 40.0 * Math.sin(lon / 3.0 * pi)) * 2.0 / 3.0;ret += (160.0 * Math.sin(lon / 12.0 * pi) + 320 * Math.sin(lon * pi / 30.0)) * 2.0 / 3.0;return ret;}private static double transformLon(double lat, double lon) {double ret = 300.0 + lat + 2.0 * lon + 0.1 * lat * lat + 0.1 * lat * lon + 0.1 * Math.sqrt(Math.abs(lat));ret += (20.0 * Math.sin(6.0 * lat * pi) + 20.0 * Math.sin(2.0 * lat * 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 += (150.0 * Math.sin(lat / 12.0 * pi) + 300.0 * Math.sin(lat / 30.0 * pi)) * 2.0 / 3.0;return ret;}/*** 高德转百度* @param gd_lat* @param gd_lon* @return* @Author : pjh. create at 2021年1月21日 上午10:43:30*/public static double[] gaoDeToBaidu(double gd_lat, double gd_lon) {double PI = 3.14159265358979324 * 3000.0 / 180.0;double x = gd_lon, y = gd_lat;double z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * PI);double theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * PI);double wgLat = z * Math.sin(theta) + 0.006;// 纬度double wgLon = z * Math.cos(theta) + 0.0065;// 经度double[] loc = { wgLat, wgLon };return loc;}/*** 度转为   ddmm.mmmm* @param angleInDegrees* @return*/public static String degreesToDMSLon(double angleInDegrees) {// 转换为度、分钟int degrees = (int) angleInDegrees;double temp = (angleInDegrees - degrees) * 60 ;int minutes = (int) temp;double minutesPoint = (temp - minutes);minutesPoint = Math.round(minutesPoint * 10000);// 返回度分格式的字符串return String.format("%3d%02d.%04d", degrees, minutes, (int) Math.round(minutesPoint));}/*** 度转为   ddmm.mmmm* @param angleInDegrees* @return*/public static String degreesToDMSLat(double angleInDegrees) {// 转换为度、分钟int degrees = (int) angleInDegrees;double temp = (angleInDegrees - degrees) * 60 ;int minutes = (int) temp;double minutesPoint = (temp - minutes);minutesPoint = Math.round(minutesPoint * 10000);// 返回度分格式的字符串return String.format("%2d%02d.%04d", degrees, minutes, (int) Math.round(minutesPoint));}/*** WGS 转为 度* @param dm* @return*/public static double df2LatLng(String dm) {if (dm == null) {return 0;}try {dm = dm.replace(" ", "");int d = Integer.valueOf(dm.substring(0, dm.lastIndexOf(".") - 2));// 兼容经纬度的转换double fen = Double.parseDouble(dm.substring(String.valueOf(d).length()));double lat = (fen / 60) + Math.abs(d);if (lat < 0) {lat = -lat;}return Double.parseDouble(String.format("%.7f", lat));} catch (Exception e) {e.printStackTrace();}return 0;}public static void main(String[] args) {// 北斗芯片获取的经纬度为WGS84地理坐标ddmm.mmmm格式String w = "3004.1628";String j = "10403.5574";System.out.println( "WGS84数据："+ df2LatLng(w) + ";" + df2LatLng(j));double[] a = wgs2gcj(df2LatLng(w), df2LatLng(j));System.out.println("WGS84转为高德数据："+ a[0] +";" + a[1]);double[] c = wgs2gcj_dddddd(w, j);System.out.println("WGS84转为高德数据："+ c[0] +";" + c[1]);MapPoint p = PositionUtil.gcj_To_Gps84(c[1],c[0]);System.out.println("高德转为WGS84："+ p.getLat() +";" + p.getLng());//dd.dddd 转为 ddmm.mmmmSystem.out.println("WGS原始数据：" + degreesToDMSLat(p.getLat()) + ";" + degreesToDMSLon(p.getLng()));}
}
public class PositionUtil {public static final String BAIDU_LBS_TYPE = "bd09ll";public static double pi = 3.1415926535897932384626;public static double a = 6378245.0;public static double ee = 0.00669342162296594323;/*** 84 to 火星坐标系 (GCJ-02)** @param lat* @param lon* @return*/public static MapPoint gps84_To_Gcj02(double lon, double lat) {if (outOfChina(lon, lat)) {return null;}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 = 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);double mgLat = lat + dLat;double mgLon = lon + dLon;return new MapPoint(mgLon, mgLat);}/*** * 火星坐标系 (GCJ-02) to 84 * * @param lon * @param lat * @return* */public static MapPoint gcj_To_Gps84(double lon, double lat) {MapPoint gps = transform(lon, lat);double lontitude = lon * 2 - gps.getLng();double latitude = lat * 2 - gps.getLat();return new MapPoint(lontitude, latitude);}/*** 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 将 GCJ-02 坐标转换成 BD-09 坐标** @param gg_lat* @param gg_lon*/public static MapPoint gcj02_To_Bd09(double gg_lon, double gg_lat) {double x = gg_lon, y = gg_lat;double z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * pi);double theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * pi);double bd_lon = z * Math.cos(theta) + 0.0065;double bd_lat = z * Math.sin(theta) + 0.006;return new MapPoint(bd_lon, bd_lat);}/*** * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 * * 将 BD-09 坐标转换成GCJ-02 坐标 * * @param* bd_lat * @param bd_lon * @return*/public static MapPoint bd09_To_Gcj02(double bd_lon, double bd_lat) {double x = bd_lon - 0.0065, y = bd_lat - 0.006;double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * pi);double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * pi);double gg_lon = z * Math.cos(theta);double gg_lat = z * Math.sin(theta);return new MapPoint(gg_lon, gg_lat);}/*** (BD-09)-->84* @param bd_lat* @param bd_lon* @return*/public static MapPoint bd09_To_Gps84(double bd_lon, double bd_lat) {MapPoint gcj02 = PositionUtil.bd09_To_Gcj02(bd_lat, bd_lon);MapPoint map84 = PositionUtil.gcj_To_Gps84(gcj02.getLat(),gcj02.getLng());return map84;}public static void main(String[] args) {MapPoint mapPoint = gps84_To_Gcj02(100.280295411,25.6047747599);System.out.println(JSONUtil.toJsonStr(mapPoint));}public static boolean outOfChina(double lon, double lat) {if (lon < 72.004 || lon > 137.8347) {return true;}if (lat < 0.8293 || lat > 55.8271) {return true;}return false;}public static MapPoint transform(double lon, double lat) {if (outOfChina(lon, lat)) {return new MapPoint(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 = 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);double mgLat = lat + dLat;double mgLon = lon + dLon;return new MapPoint(mgLon, mgLat);}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) + 320 * 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;}}

运行效果（转化有精度丢失，不影响大局）：

这篇关于WGS ddmm.mmmm 与度 互相转化的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/1057157。

相关文章

js对象和字符串互相转化

js对象和字符串互相转化

1、JSON.stringify() 例子： <!DOCTYPE html><html><head><meta charset="utf-8" /><script src='jquery-1.6.1.js'></script><title></title></head><body> </body><script>var clock={hour:12,minute:10,secon
阅读更多...
WGS_1984_UTM_Zone_50N全面深入介绍

WGS_1984_UTM_Zone_50N全面深入介绍

在地理信息系统（GIS）和遥感领域，WGS_1984_UTM_Zone_50N是一种广泛使用的投影坐标系，它结合了WGS 84地理坐标系和通用横轴墨卡托投影（UTM）的优点，为地图制作、空间分析、遥感数据处理等领域提供了强大的技术支持。本文将全面深入介绍WGS_1984_UTM_Zone_50N的定义、特点、应用及其相关概念。 一、WGS 84地理坐标系基础 WGS 84（World Geod
阅读更多...
Cell-在十字花科植物中年生和多次开花多年生开花行为的互相转化-文献精读21

Cell-在十字花科植物中年生和多次开花多年生开花行为的互相转化-文献精读21

Reciprocal conversion between annual and polycarpic perennial flowering behavior in the Brassicaceae 在十字花科植物中年生和多次开花多年生开花行为的互相转化 亮点 喜马拉雅须弥芥 和 内华达糖芥 是两个多年生植物模型 MADS-box 基因的剂量效应决定了一年生、二年生和多年生的生长习
阅读更多...
【java】Date String 互相转化

【java】Date String 互相转化

/*** Date型の日付をyyyy/MM/dd HH:mm:ss形式の文字列に変換する。** @param　value Date型の日付* @return 変換後の値*/public static String dateTimeToString(Date value, String format) {if (null == value) {return "";}DateFormat df
阅读更多...
char 、 unsigned char 互相转化

char 、 unsigned char 互相转化

1. 利用unsigned char （即uchar） 保存char 数据  ，直接赋值即可 unsigned char uc; char c=-33; uc= c; cout<<(int)((char)uc)<<endl; 2. 当移位运算和减法运算混合起来的时候记得加括号 int ss =  256 - (200*200>>8) ;
阅读更多...
java,int与integer互相转化

java,int与integer互相转化

int转Integer Java代码 int i = 0; Integer a = new Integer(i); Integer转int Java代码 Integer b = new Integer(0); int i = b.intValue();
阅读更多...
WES/WGS检测基因突变与对应表型汇总

WES/WGS检测基因突变与对应表型汇总

更新日期：2024年4月17日 1. CYP21A2 - 先天性肾上腺皮质增生 Cytochrome P450 Family 21 Subfamily A Member 2 细胞色素P450家族21亚家族A成员2 该基因编码酶的细胞色素P450超家族的成员。细胞色素P450蛋白是单加氧酶，其催化参与药物代谢和胆固醇、类固醇和其它脂质的合成的许多反应。该蛋白定位于内质网并在21位羟基化类固醇。
阅读更多...
【概率基础】生成式模型与判别式模型最大的区别是什么？两者可以互相转化吗?

【概率基础】生成式模型与判别式模型最大的区别是什么？两者可以互相转化吗?

1. 生成式模型与判别式模型最大的区别 生成式模型与判别式模型最大的区别在于它们对数据的建模方式和目标。 建模方式： 生成式模型（Generative models）尝试学习输入数据和输出标签的联合概率分布 P ( X , Y ) P(X, Y) P(X,Y)。它关注于如何根据这个分布生成数据，即模型试图理解数据是如何产生的，包括数据的特征以及特征与标签之间的关联。判别式模型（Discrim
阅读更多...
重测序（WGS）-群体遗传

重测序（WGS）-群体遗传

进化论 - 达尔文进化论 生物物种存在广泛变异自然选择，适者生存和繁殖，不适者淘汰生命是单起源的 - 中性学说 大多数突变是“中性”的进化主要是中性突变进行随机的遗传漂变 什么是群体？ 概念 遗传学上的群体，与生态学上的种群，是一个概念。 种群 population：指同一时间生活在一定自然区域内，同种生物的所有个体。种群中的个体并不是机械地集合在一起，而是彼此可以交配，并通过繁殖将
阅读更多...
百度坐标、火星坐标与WGS-84坐标转换（纠偏加偏）工具及其开发思路

百度坐标、火星坐标与WGS-84坐标转换（纠偏加偏）工具及其开发思路

一、杂言 百度坐标（BD-09）、火星坐标（GCJ-02）坐标与WGS-84坐标相关说明这里不再赘述，其他博客文章对于三类坐标系的说明已经比较详细，可以自行搜索了解，本篇文章主要对本人开发的各类转换工具进行说明（注：本人主要使用Python脚本进行开发）。文章末会将工具下载链接附上，工具为我在读研期间学习与实践的成果，可以提供数据服务，请尊重知识成果，欢迎提供意见或建议。转换源码为GitHub开
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2