Java坐标标转换工具类

最新推荐文章于 2024-06-25 08:07:00 发布
最新推荐文章于 2024-06-25 08:07:00 发布
阅读量1.4k 收藏 5
点赞数 2
分类专栏： Java 文章标签： java
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_37024898/article/details/126399017
版权
本文介绍了如何在Java中使用JTS库进行坐标系统的转换，包括WGS84到GCJ-02（火星坐标）、GCJ-02到BD-09（百度坐标）以及BD-09回GCJ-02的转换方法，同时提供了经纬度范围判断和距离计算功能。
摘要由CSDN通过智能技术生成

1、pom.xml添加依赖

<dependency>
    <groupId>com.vividsolutions</groupId>
    <artifactId>jts-core</artifactId>
    <version>1.14.0</version>
</dependency>

2、创建坐标标转换工具类CoordinateUtil

import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.GeometryFactory;
import com.vividsolutions.jts.geom.Point;

import java.math.BigDecimal;
import java.math.RoundingMode;
import java.text.DecimalFormat;

/**
 * 坐标标转换工具类
 */
public class CoordinateUtil {

    /**
     * 元周率
     */
    public final static double pi = Math.PI;
    /**
     * 卫星椭球坐标投影到平面地图坐标系的投影因子
     */
    public final static double a = 6378245.0;

    /**
     * 椭球的偏心率 (a^2 - b^2) / a^2
     */
    public final static double ee = 0.00669342162296594323;

    /**
     * @param lon
     * @param lat
     * @return
     * @Description WGS84 to 火星坐标系 (GCJ-02)
     */
    public static double[] wgs84_To_Gcj02(double lon, double lat) {
        if (outOfChina(lat, lon)) {
            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 double[]{mgLon, mgLat};
    }

    /**
     * @param lon
     * @param lat
     * @return
     * @Description 火星坐标系 (GCJ-02) to WGS84
     */
    public static double[] gcj02_To_Wgs84(double lon, double lat) {
        double[] gps = transform(lat, lon);
        double lontitude = lon * 2 - gps[1];
        double latitude = lat * 2 - gps[0];
        return new double[]{lontitude, latitude};
    }

    /**
     * @param gg_lon
     * @param gg_lat
     * @return
     * @Description 火星坐标系 (GCJ-02) to 百度坐标系 (BD-09)
     */
    public static double[] 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 double[]{bd_lon, bd_lat};
    }

    /**
     * @param bd_lon
     * @param bd_lat
     * @return
     * @Description 百度坐标系 (BD-09) to 火星坐标系 (GCJ-02)
     */
    public static double[] 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 double[]{gg_lon, gg_lat};
    }

    /**
     * @param bd_lat
     * @param bd_lon
     * @return
     * @Description 百度坐标系 (BD-09) to WGS84
     */
    public static double[] bd09_To_Wgs84(double bd_lon, double bd_lat) {
        double[] gcj02 = CoordinateUtil.bd09_To_Gcj02(bd_lon, bd_lat);
        double[] map84 = CoordinateUtil.gcj02_To_Wgs84(gcj02[0], gcj02[1]);
        return map84;

    }

    /**
     * @param lat
     * @param lon
     * @return
     * @Description 判断是否在中国范围内
     */
    public static boolean outOfChina(double lat, double lon) {
        if (lon < 72.004 || lon > 137.8347)
            return true;
        if (lat < 0.8293 || lat > 55.8271)
            return true;
        return false;
    }

    /**
     * @param lat
     * @param lon
     * @return
     * @Description transform
     */
    private static double[] transform(double lat, double lon) {
        if (outOfChina(lat, lon)) {
            return new double[]{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 double[]{mgLat, mgLon};
    }

    /**
     * @param x
     * @param y
     * @return
     * @Description transformLat
     */
    private 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;
    }

    /**
     * @param x
     * @param y
     * @return
     * @Description transformLon
     */
    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;
    }


    /**
     * @param longitude 经度
     * @param latitude  纬度
     * @return double[] x y
     * @Description WGS84 to 高斯投影(6度分带)
     */
    public static double[] wgs84_To_Gauss6(double longitude, double latitude) {
        int ProjNo = 0;
        int ZoneWide; // //带宽
        double[] output = new double[2];
        double longitude1, latitude1, longitude0, X0, Y0, xval, yval;
        double a, f, e2, ee, NN, T, C, A, M, iPI;
        iPI = 0.0174532925199433; // //3.1415926535898/180.0;
        ZoneWide = 6; //6度带宽
        a = 6378137.0;
        f = 1.0 / 298.257223563; //WGS84坐标系参数
        //a = 6378245.0;f = 1.0 / 298.3; // 54年北京坐标系参数
        // //a=6378140.0; f=1/298.257; //80年西安坐标系参数
        ProjNo = (int) (longitude / ZoneWide);
        longitude0 = (double) (ProjNo * ZoneWide + ZoneWide / 2);
        longitude0 = longitude0 * iPI;
        longitude1 = longitude * iPI; // 经度转换为弧度
        latitude1 = latitude * iPI; // 纬度转换为弧度
        e2 = 2 * f - f * f;
        ee = e2 / (1.0 - e2);
        NN = a
                / Math.sqrt(1.0 - e2 * Math.sin(latitude1)
                * Math.sin(latitude1));
        T = Math.tan(latitude1) * Math.tan(latitude1);
        C = ee * Math.cos(latitude1) * Math.cos(latitude1);
        A = (longitude1 - longitude0) * Math.cos(latitude1);
        M = a
                * ((1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256)
                * latitude1
                - (3 * e2 / 8 + 3 * e2 * e2 / 32 + 45 * e2 * e2 * e2
                / 1024) * Math.sin(2 * latitude1)
                + (15 * e2 * e2 / 256 + 45 * e2 * e2 * e2 / 1024)
                * Math.sin(4 * latitude1) - (35 * e2 * e2 * e2 / 3072)
                * Math.sin(6 * latitude1));
        // 因为是以赤道为Y轴的，与我们南北为Y轴是相反的，所以xy与高斯投影的标准xy正好相反;
        xval = NN
                * (A + (1 - T + C) * A * A * A / 6 + (5 - 18 * T + T * T + 14
                * C - 58 * ee)
                * A * A * A * A * A / 120);
        yval = M
                + NN
                * Math.tan(latitude1)
                * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * A * A * A * A / 24 + (61
                - 58 * T + T * T + 270 * C - 330 * ee)
                * A * A * A * A * A * A / 720);
        X0 = 1000000L * (ProjNo + 1) + 500000L;
        Y0 = 0;
        xval = xval + X0;
        yval = yval + Y0;
        output[0] = xval;
        output[1] = yval;
        return output;
    }


    /**
     * 两经纬度之间的距离
     *
     * @param lat_a
     * @param lng_a
     * @param lat_b
     * @param lng_b
     * @return
     */
    public static double distanceToPoint(double lat_a, double lng_a, double lat_b, double lng_b) {
        //jts提供的几何要素工厂类
        GeometryFactory geometryFactory = new GeometryFactory();

        //火星坐标(gcj02)转GPS坐标(WGS84)
        double[] wgsPntA = CoordinateUtil.gcj02_To_Wgs84(lng_a, lat_a);
        double[] wgsPntB = CoordinateUtil.gcj02_To_Wgs84(lng_b, lat_b);

        //WGS84->高斯6度分带投影
        double[] gaussPntA = wgs84_To_Gauss6(wgsPntA[0], wgsPntA[1]);
        double[] gaussPntB = wgs84_To_Gauss6(wgsPntB[0], wgsPntB[1]);

        //通过几何要素工厂得到point实体
        Point pntA = geometryFactory.createPoint(new Coordinate(gaussPntA[0], gaussPntA[1]));
        Point pntB = geometryFactory.createPoint(new Coordinate(gaussPntB[0], gaussPntB[1]));
        // 两点距离
        return pntA.distance(pntB);
    }

    /**
     * 保留6位精度
     *
     * @return
     */
    public static String format(String resource) {
        //保留两位小数
        DecimalFormat decimalFormat = new DecimalFormat("###.000000");
        decimalFormat.setRoundingMode(RoundingMode.DOWN);
        String target = decimalFormat.format(new BigDecimal(resource));
        return target;
    }
}
我有一台FC
关注
  • 2
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
java实现坐标系转换
09-25
java实现2000坐标转wgs84，WGS84和cgcs2000互转，使用String字符串解析和三参数算法实现，只需坐标就能实现
Gps坐标系转换Java工具类WGS坐标与Google和百度坐标互转
12-28
Gps坐标系转换Java工具类WGS坐标与Google和百度坐标互转，偏差很小，与百度Api调用转换几乎相差无几，程序为Java程序，工具类直接传经纬度调用相应转换方法即可获取转换后的返回值
Java精准地图坐标转换 高德 百度 谷歌 腾讯 高德 批量转换工具类达到项目生产精度
最新发布
2401_85699160的博客
06-25 1062
大型分布式系统犹如一个生命，系统中各个服务犹如骨骼，其中的数据犹如血液，而Kafka犹如经络，串联整个系统。这份Kafka源码笔记通过大量的设计图展示、代码分析、示例分享，把Kafka的实现脉络展示在读者面前，帮助读者更好地研读Kafka代码。麻烦帮忙转发一下这篇文章+关注我分 如：112°10.25’’//计算并转换为stringreturn jwd;/**经纬度转换，度分转度度.如：112°30.4128 = 112.50688。
坐标转高德工具
Fog_cloud的博客
11-18 293
坐标转高德工具 package org.tongzhou.unit; /** * Created by Wonshine on 2018-11-14. */ public class GPSUtil { public static double pi = 3.1415926535897932384626; public static double x_pi = 3.14159...
java坐标系转换
06-09
ARCGIS，坐标系转换以及纠偏（GS84以及GCJ02B），欢迎大家下载
GIS java坐标转换工具类
成长的笨熊
04-10 397
【代码】GIS java坐标转换工具类。
java# 坐标系转换工具类
myshare2022的专栏
11-20 1007
或者简单地来说，WGS84就是一堆参数，用于建立球体或者类球体，来近似地球。这里讲的坐标系统，相对于笛卡尔坐标系统，要复杂许多，但作用却都是一样，主要用于定位，也就是精确地定位地表上的一点。：地理坐标系统，Google Earth和中国外的Google Map使用，另外，目前基本上所有定位空间位置的设备都使用这种坐标系统，如手机的GPS系统。：投影坐标系统，也就是我们平常所说的火星坐标系，Google Map中国、高德和腾讯，这个是中国自己在WGS84基础上加密而成，目的显而易见。
java坐标转换工具类.docx
10-24
根据提供的文档信息，本文将详细解析“java坐标转换工具类”中的关键知识点，特别是与WGS84坐标系、火星坐标（GCJ-02）以及百度坐标（BD-09）之间的转换方法。 ### 一、WGS84坐标系简介 WGS84坐标系（World ...
geojson文件坐标转换工具
05-24
本项目提供了一个用C#编写的GeoJSON文件坐标转换工具，支持WGS84、GCJ02（国测局加密坐标系，也称火星坐标）以及百度坐标之间的相互转换。 1. **WGS84坐标系**：这是一个全球通用的大地坐标系统，广泛应用于GPS设备...
java_adnroid_GPS坐标转换工具类的标准代码_百度坐标系 GCJ-02 WGS84 经纬度转换
05-15
一些工具类代码块的标准代码,包括但不限于： 百度坐标系(BD-09)转WGS坐标(百度坐标纬度,百度坐标经度),WGS84坐标数组 WGS坐标转百度坐标系(BD-09)(WGS84坐标系的经度,WGS84坐标系的纬度),百度坐标数组 火星坐标系...
坐标转换工具
11-29
集合了各个坐标系之间的转换，以及计算坐标直接的距离等多个方法
GIS坐标转换 JAVA版
12-12
用于GIS各坐标的转换！ 我是为了将EPSG:4326坐标转换成EPSG:900913坐标，找到的工具，然后添加了两者之间的转换方式！ 使用方式： private static Projection proj = ProjectionFactory .getNamedPROJ4CoordinateSystem("epsg:900913"); public static double[] cast(double lon, double lat) { Point2D.Double src = new Point2D.Double(lon, lat); Point2D.Double dst = new Point2D.Double(0, 0); proj.transform(src, dst); return new double[] { dst.x, dst.y }; }
收集的java坐标转换WGS84、GCJ02B、BD09及纠偏
12-13
收集的java坐标转换工具，包括WGS84、GCJ02B、BD09，及国内纠偏。 收集的java坐标转换工具，包括WGS84、GCJ02B、BD09，及国内纠偏。
【Java】 # java实现坐标转换工具类
LRcoding
05-24 3288
1. 主流坐标系简介 WGS84坐标系 地球坐标系，国际上通用的坐标系 使用GPS芯片或者北斗芯片的设备，获取的经纬度就是WGS84地理坐标系。 地图API：谷歌地图使用的是WGS84坐标系，但是中国范围使用的是GCJ02坐标系 GCJ02坐标系 火星坐标系，WGS84坐标系经过加密后的坐标系 出于国家安全考虑，国内的所有导航电子地图必须使用国家测绘局制定的加密坐标系，即将一个真实的经纬度通过加密变成一个对应位置不正确的经纬度 地图API：高德MapABC地图、腾讯搜搜地图、阿里云地图 BD09
java常用地图坐标系转换工具类,支持谷歌，腾讯，百度等主流的地理坐标转换
qq_45926602的博客
02-16 2552
java常用地图坐标系转换工具类,支持谷歌，腾讯，百度等主流的地理坐标转换
Java - 坐标转换工具类 WGS84 ⇋ GCJ02
ego_rainki
12-29 477
原JS代码 double计算小数点后7位后有误差 public class CoordinateUtils { private static final int a = 6378245; private static final double ee = 0.006693421622965823; private static final double PI = Math.PI; private static double sin(double d) { r
Java实现WGS84、百度、腾讯、高德等主流的地理坐标转换
分享DotNet技术和日常开发笔记，DotNet，Python为主。
09-22 7341
主要是WGS84(大地坐标系)、GCJ02(国测局坐标系，如高德地图、腾讯地图、谷歌中国范围地图)、BD09(百度坐标系)三种主流坐标系之间的转换。Java实现WGS84、百度坐标系、腾讯左边西、高德坐标系转化工具类。
坐标系转换工具类Java版&Js版
CSDN已弃用
07-29 1014
Java版本： public class GPSUtils { private static double x_pi = 3.14159265358979324 * 3000.0 / 180.0; // π private static double pi = 3.1415926535897932384626; // 长半轴 private static double a = 6378245.0; // 扁率 private static double
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值