WGS ddmm.mmmm 与度 互相转化

已于 2024-06-13 10:09:15 修改
阅读量235 收藏
点赞数 3
文章标签： java 算法 开发语言
于 2024-06-13 09:21:10 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/jiejiegua/article/details/139642011
版权

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.mmmm
		System.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;
    }

}

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

carry杰
关注
解析经纬度，将度分转成度
清宵尚温的博客
11-10 1761
1. 问题描述 解析地图的经纬度，将GNSS协议的ddmm.mmmm转成dd.ddddd。 2. 工具类 目前将GNSS协议中的ddmm.mmmmm去点然后会回传到函数的的参数中。 /** * 解析转出经纬度，转成度的形式 * GNSS协议 dd.mm.mmmmm 转成dd.mm.sssss * * @param str GNSS的String表示 * @return dd.ddddd度,小数点后保留5位 */ public sta
GpsToUtil.java
01-21
GPS转高德 转百度 ddmm.mmmm转dd,dddd
GNSS NMEA ddmm.mmmmm格式转换
qlexcel的专栏
04-30 7157
ddmm.mmmm（度分）转换成常用坐标dd.ddddd GPS NMEA采集到的经度和纬度格式分别是 纬度ddmm.mmmmm（度分） 经度dddmm.mmmmm（度分） 其中dd是度数，mm.mmmmm是分数。 但实际如果想要在百度地图或Google地图中显示,则需要dd.mmmm和ddd.mmmm 其中dd.mmmm和ddd.mmmm全是度数。 因此，把ddmm.mmmmm中的分数转换为度数即可（1度=60分）。下面是转换算法: 转换后的纬度 = dd + (mm.mmmm/60) 转换后的经度 =
DDMM
图灵与对话
08-29 769
000300.SH    IF   沪深300 000905.SH    IC  中证 500 000016.SH    IH  上证50
GPS 原始坐标 (WGS-84) ddmm.mmmm 格式转高德（GCJ）百度（BD09）坐标
仙_
10-21 6411
本博客适用于ddmm.mmmm或者dd.dddd格式转高德（GCJ）百度（BD09）坐标 最近的项目中有个定位功能，定位是GPS原始坐标 (WGS-84)，数据格式是ddmm.mmmm类型的，前端大佬用的地图是百度的，了解过经纬度这一块的应该都知道经纬度偏移吧，这里就不多说明了，所以现在需要把WGS-84转换成BD09坐标。 先列举一下数据(ddmm.mmmm) Lat:223...
java ddmm.mmmm转ddd,GPS坐标转换ddmm.mmmm（度分）-> dd.dddd(度)
weixin_42177768的博客
03-12 1221
#include #include #include int ddmm2dd(const char *ddmm, char *dd){if (NULL == ddmm || NULL == dd){return -1;}int lenSrc = strlen(ddmm)+1;int lenMm = 0;int flag = 1;memcpy(dd,ddmm,lenSrc);char *pcMm;d...
GPS坐标转换ddmm.mmmm（度分）-> dd.dddd(度)
热门推荐
allon19的专栏
12-01 1万+
#include #include #include int ddmm2dd(const char *ddmm, char *dd) { if (NULL == ddmm || NULL == dd) { return -1; } int lenSrc = strlen(ddmm)+1; int lenMm = 0; int flag = 1; memcpy(dd,dd
GPS坐标转换ddmm.mmmm（度分）-＞ dd.dddd(度)
爱她就要努力
04-30 1207
#include <stdio.h>#include <string.h>#include <stdlib.h> int ddmm2dd(const char *ddmm, char *dd){ if (NULL == ddmm || NULL == dd) { return -1; } int lenSrc = strlen(ddmm)+1; int len...
GPS模块输出的NMEA数据ddmm.mmmm转换成dd.ddddd,在google Earth Pro中描点
qingfengxd1的博客
05-06 4752
GPS模块输出的数据是NMEA格式，其中GPGGA字段包含我们需要的经纬度信息。 例：$GPGGA,092204.999,4250.5589,S,14718.5084,E,1,04,24.4,12.2,M,19.7,M,,0000*1F 其中 4250.5589,S,14718.5084,E 就是经度和纬度数据，其格式是 纬度：ddmm.mmmmm（度分） 经度：dddmm.mmmmm（度...
将度分（ddmm.mmmm）的坐标转换为常用坐标（dd.ddddd）
10-18
GNSS NMEA 格式的坐标（ddmm.mmmmm）转换为原始 GPS 格式
PositionUtil.java
10-21
GPS 原始坐标 (WGS-84)转百度BD09坐标系，适用于ddmm.mmmm或者dd.dddd格式转高德（GCJ）百度（BD09）坐标
C++ GPS导航定位坐标转换程序
03-15
GPS卫星向地面发射包含时间信息的信号，接收机通过测量信号到达的时间来计算与各个卫星的距离，进而通过三角定位方法确定地理位置。然而，GPS系统中使用的坐标系是WGS84（World Geodetic System 1984），这并不是...
w042基于web的IT技术交流和分享平台的设计与实现
卓怡工作室的博客
11-18 710
我国科学技术的不断发展，计算机的应用日渐成熟，其强大的功能给人们留下深刻的印象，它已经应用到了人类社会的各个层次的领域，发挥着重要的不可替换的作用。信息管理作为计算机应用的一部分，使用计算机进行管理，具有非常明显的优点，利用网络的优势特开发了本基于SpringBoot的IT技术交流和分享平台。 本IT技术交流和分享平台是基于SpringBoot框架，采用Java技术，MYSQL数据库进行开发的。系统具有灵活的一体化设计方式，圆满完成了整个系统的界面设计。本系统实现了用户功能模块和管理员功能模块两大部分，
Mybatis
庸不易的博客
11-15 827
1、定义不同：#{} 预处理：而 ${} 是直接替换2、使用不同：#{} 适用于所有类型的参数匹配；但 ${} 只适用于数值类型。3、安全性不同：#{} 性能高，并且没有安全问题；但 $ {} 存在 SQL 注入的安全问题。4、使用场景不同：当传递的是一个 SQL 关键字 的时候，只能 使用 ${}。当传递的是一个字段，总之，就是需要获取到参数类型 与 内容，只能使用 #{}。（PS：数字类型，#{} 和 ${} 都是可以使用的！5、 ${} 不能用于 模糊匹配查询；
java加锁问题详解
qq_74056922的博客
11-14 799
java(JVAM)对加锁的API进行了封装,通过使用synchronized关键字来完成加锁操作。
什么是面向对象编程？什么是面向过程编程？
m0_70208894的博客
11-15 657
不同的编程范式适用于不同的场景，选择合适的范式可以提高代码的可维护性、可读性和效率。Java 和 Go 的选择反映了不同的设计哲学，Java 更加注重封装和对象的使用，而 Go 更加简洁，强调过程和函数的使用。
HashMap源码分析下
L25000的专栏
11-16 701
如果这个cap刚好是2的幂次方呢，即4,8,16,32,64等这样的，这样的数已经满足条件了，再减去1，把数变成3,7,31,63。，即hash值的第5为是1，这样在扩容到32后，取索引得出来的结果就和16的不一样了，索引位置比oldCap大了一个oldCap值；,求余数，n是一个2的幂次方数，减去1，即后面全是1，前面全是0，与上这个hash，只能保证经过运算，结果都会在。，即hash值的第5为是0，这样在扩容到32后，取索引得出来的结果还是和16的一样；cap=1，2的0次方就等于1，返回的是就是4。
阶乘递归javascript
最新发布
2401_83905000的博客
11-18 83
【代码】阶乘递归javascript。
BP新增页签增强
BinGeneral的博客
11-15 996
以上就是BP新增页签的介绍，希望对您有所帮助。在总结过程中，该篇文档给予了极大的帮助和借鉴。
Cesium.Ellipsoid.WGS84.cartAttitudeToFixedFrame 报不是一个方法
10-24
抱歉，我的描述有误。在Cesium.js中，`Cesium.Ellipsoid.WGS84.cartAttitudeToFixedFrame` 是一个静态方法，而不是实例方法。它用于从地球坐标系（地心坐标系）转换到固定的框架（如当地水平面）。正确的用法应该是直接通过 `Cesium.Ellipsoid.WGS84` 对象调用，不需要创建其实例。这是一个例子： ```javascript const cartographic = new Cesium.Cartographic(longitude, latitude); const attitude = { heading: userHeading * (Math.PI / 180), pitch: userPitch * (Math.PI / 180), roll: 0 // 简单示例，通常可以包含roll }; // 调用 static 方法计算旋转矩阵 const rotationMatrix = Cesium.Ellipsoid.WGS84.cartAttitudeToFixedFrame(cartographic, attitude); // 使用旋转矩阵来更新 Billboard 的方向 const direction = new Cesium.Cartesian3.fromRotationMatrix(rotationMatrix); billboard.orientation.direction = direction; ``` 这里的 `rotationMatrix` 可以用来处理 Billboard 的方向变换。请注意，实际应用中可能需要处理更复杂的场景，比如考虑地球曲率等因素。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

carry杰

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值