经纬度坐标工具

最新推荐文章于 2024-09-12 05:48:32 发布
最新推荐文章于 2024-09-12 05:48:32 发布
阅读量456 收藏
点赞数
分类专栏： Utils 文章标签： 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yiqiu1959/article/details/132139883
版权

LngLatUtil :用于计算里程数


import cn.hutool.core.util.ArrayUtil;
import com.alibaba.fastjson.JSON;
import com.alibaba.fastjson.JSONObject;
import lombok.Getter;
import lombok.Setter;

import java.io.FileInputStream;
import java.io.Serializable;
import java.text.DecimalFormat;
import java.util.ArrayList;
import java.util.List;
import java.util.Objects;

public class LngLatUtil {
    /**
     * 地球半径，m
     */
    public static final double EARTH_RADIUS = 6371393d;
    /**
     * 偏离参考弧度大小
     */
    public static final double OFFSET_RADIANS = Math.PI * 0.25d;

    public static double getRand(double min, double max) {
        double v = Math.random() * (max - min) + min;
        DecimalFormat df = new DecimalFormat("#0.00");
        return Double.parseDouble(df.format(v));
    }

    public static double getOffset(Point p1, Point p2) {
        // 纬度之差
        double a = Math.toRadians(p1.lat) - Math.toRadians(p2.lat);
        // 经度之差
        double b = Math.toRadians(p1.lng) - Math.toRadians(p2.lng);
        // 计算两点距离的公式
        double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) +
                Math.cos(p1.lat) * Math.cos(p2.lat) * Math.pow(Math.sin(b / 2), 2)));
        // 弧长乘地球半径, 返回单位: 米
        return s * EARTH_RADIUS;
    }

    public static Point offsetPoint(double lng, double lat, double offset, double radians) {
        Point point = new Point();
        if (offset > 0) {
            point.setLng(lng + offset * Math.sin(radians) * 180d / (Math.PI * EARTH_RADIUS * Math.cos(lat * Math.PI / 180d)));
            point.setLat(lat + offset * Math.cos(radians) / (Math.PI * EARTH_RADIUS / 180d));
        } else {
            point.setLng(lng);
            point.setLat(lng);
        }
        return point;
    }

    public static Point getNewPoint(Point samplePrev, Point sampleCur,int index) {
        double radians = 0d;
        if (Objects.nonNull(samplePrev)) {
            radians = getRadians(sampleCur, samplePrev);
        }
        radians = offsetRadians(radians);
        double offset;
        if(index > 0){
            offset = getRand(10, 50);
        }else{
            offset = getRand(150, 200);
        }
        return offsetPoint(sampleCur.lng, sampleCur.lat, offset, radians);
    }

    public static Point getNewPoint(Point prev, Point samplePrev, Point sampleCur) {
        double radians = 0d;
        if (Objects.nonNull(samplePrev)) {
            radians = getRadians(samplePrev, sampleCur);
        }
        radians = offsetRadians(radians);
        double offset = getRand(50, 100);
        if (Objects.isNull(prev)) {
            prev = sampleCur;
        }
        return offsetPoint(prev.lng, prev.lat, offset, radians);
    }

    public static double getRadians(Point p1, Point p2) {
        double y = Math.sin(p2.getLng() - p1.getLng()) * Math.cos(p2.getLat());
        double x = Math.cos(p1.getLat()) * Math.sin(p2.getLat()) - Math.sin(p1.getLat()) * Math.cos(p2.getLat()) * Math.cos(p2.getLng() - p1.getLng());
        return Math.atan2(y, x);
    }


    public static double offsetRadians(double radians) {
        double min = radians - OFFSET_RADIANS;
        double max = radians + OFFSET_RADIANS;
        return getRand(min, max);
    }

    @Getter
    @Setter
    public static class Point implements Serializable {
        private double lng;
        private double lat;
    }


    public static void main(String[] args) throws Exception {
        String path = "D:\\Desktop\\imei-data\\";
        FileInputStream fileInputStream = new FileInputStream(path + "4.json");
        byte[] bytes = new byte[1024];//每一次读取四个字节
        int readCount = 0;
        byte[] data = new byte[0];
        while ((readCount = fileInputStream.read(bytes)) != -1) {
            byte[] read = new byte[readCount];
            System.arraycopy(bytes, 0, read, 0, readCount);
            data = ArrayUtil.addAll(data, read);
        }
        String s = new String(data);
        List<Point> testData = JSONObject.parseArray(s, Point.class);
        Point samplePrev = null;
        List<Double[]> lines = new ArrayList<>();
        List<Double[]> origin = new ArrayList<>();
        for (int i = 0; i < testData.size(); i++) {
            Point sampleCur = testData.get(i);
            origin.add(new Double[]{sampleCur.lng, sampleCur.lat});
            Point p = getNewPoint(samplePrev, sampleCur, i);
            lines.add(new Double[]{p.lng, p.lat});
            samplePrev = sampleCur;
        }
//        System.out.println(JSON.toJSONString(origin));
        System.out.println();
        System.out.println(JSON.toJSONString(lines));
    }
}

CoordinateTransformUtils ： 坐标转换工具 BD09 、GCJ02、WGS84


/**
 * 百度坐标（BD09）、国测局坐标（火星坐标，GCJ02）、和WGS84坐标系之间的转换的工具
 *
 * @see 参考https://github.com/wandergis/coordtransform实现的Java版本
 * @author geosmart
 */
public class CoordinateTransformUtils {
    static double x_pi = 3.14159265358979324 * 3000.0 / 180.0;
    // π
    static double pi = 3.1415926535897932384626;
    // 长半轴
    static double a = 6378245.0;
    // 扁率
    static double ee = 0.00669342162296594323;

    /**
     * 百度坐标系(BD-09)转WGS坐标
     *
     * @param lng 百度坐标纬度
     * @param lat 百度坐标经度
     * @return WGS84坐标数组
     */
    public static double[] bd09towgs84(double lng, double lat) {
        double[] gcj = bd09togcj02(lng, lat);
        double[] wgs84 = gcj02towgs84(gcj[0], gcj[1]);
        return wgs84;
    }

    /**
     * WGS坐标转百度坐标系(BD-09)
     *
     * @param lng WGS84坐标系的经度
     * @param lat WGS84坐标系的纬度
     * @return 百度坐标数组
     */
    public static double[] wgs84tobd09(double lng, double lat) {
        double[] gcj = wgs84togcj02(lng, lat);
        double[] bd09 = gcj02tobd09(gcj[0], gcj[1]);
        return bd09;
    }

    /**
     * 火星坐标系(GCJ-02)转百度坐标系(BD-09)
     * 谷歌、高德——>百度
     *
     * @param lng 火星坐标经度
     * @param lat 火星坐标纬度
     * @return 百度坐标数组
     */
    public static double[] gcj02tobd09(double lng, double lat) {
        double z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * x_pi);
        double theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * x_pi);
        double bd_lng = z * Math.cos(theta) + 0.0065;
        double bd_lat = z * Math.sin(theta) + 0.006;
        return new double[]{bd_lng, bd_lat};
    }

    /**
     * 百度坐标系(BD-09)转火星坐标系(GCJ-02)
     * 百度——>谷歌、高德
     *
     * @param bd_lon 百度坐标纬度
     * @param bd_lat 百度坐标经度
     * @return 火星坐标数组
     */
    public static double[] bd09togcj02(double bd_lon, double bd_lat) {
        double x = bd_lon - 0.0065;
        double y = bd_lat - 0.006;
        double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
        double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
        double gg_lng = z * Math.cos(theta);
        double gg_lat = z * Math.sin(theta);
        return new double[]{gg_lng, gg_lat};
    }

    /**
     * WGS84转GCJ02(火星坐标系)
     *
     * @param lng WGS84坐标系的经度
     * @param lat WGS84坐标系的纬度
     * @return 火星坐标数组
     */
    public static double[] wgs84togcj02(double lng, double lat) {
        if (out_of_china(lng, lat)) {
            return new double[]{lng, lat};
        }
        double dlat = transformlat(lng - 105.0, lat - 35.0);
        double dlng = transformlng(lng - 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);
        dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * pi);
        double mglat = lat + dlat;
        double mglng = lng + dlng;
        return new double[]{mglng, mglat};
    }

    /**
     * GCJ02(火星坐标系)转GPS84
     *
     * @param lng 火星坐标系的经度
     * @param lat 火星坐标系纬度
     * @return WGS84坐标数组
     */
    public static double[] gcj02towgs84(double lng, double lat) {
        if (out_of_china(lng, lat)) {
            return new double[]{lng, lat};
        }
        double dlat = transformlat(lng - 105.0, lat - 35.0);
        double dlng = transformlng(lng - 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);
        dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * pi);
        double mglat = lat + dlat;
        double mglng = lng + dlng;
        return new double[]{lng * 2 - mglng, lat * 2 - mglat};
    }

    /**
     * 纬度转换
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double transformlat(double lng, double lat) {
        double ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng));
        ret += (20.0 * Math.sin(6.0 * lng * pi) + 20.0 * Math.sin(2.0 * lng * 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 += (160.0 * Math.sin(lat / 12.0 * pi) + 320 * Math.sin(lat * pi / 30.0)) * 2.0 / 3.0;
        return ret;
    }

    /**
     * 经度转换
     *
     * @param lng
     * @param lat
     * @return
     */
    public static double transformlng(double lng, double lat) {
        double ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));
        ret += (20.0 * Math.sin(6.0 * lng * pi) + 20.0 * Math.sin(2.0 * lng * pi)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(lng * pi) + 40.0 * Math.sin(lng / 3.0 * pi)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(lng / 12.0 * pi) + 300.0 * Math.sin(lng / 30.0 * pi)) * 2.0 / 3.0;
        return ret;
    }

    /**
     * 判断是否在国内，不在国内不做偏移
     *
     * @param lng
     * @param lat
     * @return
     */
    public static boolean out_of_china(double lng, double lat) {
        if (lng < 72.004 || lng > 137.8347) {
            return true;
        } else if (lat < 0.8293 || lat > 55.8271) {
            return true;
        }
        return false;
    }
}
一湫1959
关注
专栏目录
计算经纬度之间的距离、通过指定经纬度查询一定范围内的数据
故事码的博客
07-17 2316
计算经纬度之间的距离、通过指定经纬度查询一定范围内的数据 1.1 计算两个经纬度之间的距离 此方法需要导入一个依赖包，专门用于计算经纬度的数据的 <dependency> <groupId>org.gavaghan</groupId> <artifactId>geodesy</artifactId> <version>1.1.3</vers
经纬度计算距离经纬度计算距离工具
06-15
经纬度计算距离经纬度计算距离工具, 挺不错的。经纬度计算距离经纬度计算距离工具, 挺不错的。
计算经纬度、距离、方位角
zybiao2的专栏
07-26 2851
经纬度、距离、方位角计算
高德GIS地图动态显示打点数据
最新发布
2401_86358825的博客
09-12 651
12: 1,13: 1zooms: [5, 16], // 设置当前打点在地图的5-16级别区间内执行styles: [{ //第一组样式（本例中，4-9级用 第一个样式（第0项）10-13级用第二个样式（第1项））icon: { // 打点的iconsize: [300, 200], //可见区域的大小ancher: [40, 36], //锚点，图标原始大小下锚点所处的位置，相对左上角//可缺省，当使用精灵图时可用，标示图标在整个图片中左上角的位置。
计算距离的工具类
04-09
计算距离的工具类计算距离的工具类计算距离的工具类计算距离的工具类计算距离的工具类
java工具类(六)根据经纬度计算距离
热门推荐
IT全栈 华强工作室
08-14 5万+
Java实现根据经纬度计算距离 在项目开发过程中，需要根据两地经纬度坐标计算两地间距离，所用的工具类如下: Demo1: public static double getDistatce(double lat1, double lat2, double lon1, double lon2) { double R = 6371; double distanc
JAVA计算两个经纬度之间的距离
Feel Word
09-08 2万+
1、只不过是封装好的一个计算方法，计算的方式也就是上学时就学过的定理，老规矩直接贴代码，主要是方便实用，计算结果单位：米 public class MapUtils { //private static double EARTH_RADIUS = 6378.137; private static double EARTH_RADIUS = 6371.393; private static
批量经纬度坐标转换工具
08-11
标题中的“批量经纬度坐标转换工具”是一种软件或服务，其主要功能是处理大量地理坐标数据，将经纬度坐标转换为对应的省份或城市信息。在地理信息系统（GIS）和地图应用中，这种工具非常实用，因为它们能帮助用户...
经纬度坐标和xy坐标转换程序.rar
05-09
"经纬度坐标和xy坐标转换程序.rar" 提供了一个实用工具，帮助我们处理这两种坐标系统之间的转换。接下来，我们将深入探讨经纬度坐标系统和XY坐标系统的概念，以及它们之间的转换方法。 经纬度坐标系统是地球表面...
中国省份地图边界经纬度坐标点集合.zip
08-18
中国省份地图边界经纬度坐标点集合是一个非常实用的数据资源，主要包含了中国所有省份的地理边界信息，以经纬度坐标的形式呈现。这些数据在各种地图应用、地理信息系统（GIS）、数据分析和可视化项目中都有着广泛的...
经纬度坐标工具类
李卫康的博客
07-30 1256
在开发地图的时候难免会遇到经纬度的一些算法，这里我分享一些经纬度的算法，主要包括， 判断两个经纬度坐标之间的距离 判断一个经纬度坐标是否在一个圆形面中（如共享单车的电子围栏为圆形） 判断一个经纬度坐标是否在一个多边形的中（电子围栏类型为多边形） 具体算法: import java.awt.geom.Point2D; import java.util.ArrayList; import ...
经纬度坐标转换成空间直角坐标
10-15
基于matlab软件，可将经纬度坐标转换成空间直角坐标。内含WGS84参考椭球参数，可通过输入自己的经纬度坐标，实现三维空间直角坐标的转换。
经纬度距离计算小工具
04-02
经纬度距离计算小工具，经纬度距离计算小工具
根据两个经纬度坐标计算两个坐标间的距离
weixin_39297798的博客
10-15 3264
/**      * 根据两点间的经纬度计算距离      * @param $lat1 //纬度值      * @param $lng1 //经度值      * @param $lat2      * @param $lng2      * @return float 单位为(km)      */     function getDistance($lat1, $lng1,
经纬度和距离之间的计算
weixin_34129696的博客
07-22 292
不是我写的, 也看不太懂, 备忘 package com.iwhere.easy.travel.tool; import org.junit.Test; import org.springframework.util.Assert; /** */ public class GeoUtils { public final static float EARTH_RADIUS = 63...
根据经纬度，计算两点间的距离
MR.Zhang的博客
12-18 653
根据经纬度，计算两点间的距离第一步：声明一个全局变量（地球的半径）第二步：创建一个封装经纬度的VO类第三步：计算两个经纬度之间的距离 第一步：声明一个全局变量（地球的半径） /** * @Author: Mr.Zhang * 地球半径,单位 km(如果距离单位为M，即将数值*1000) */ private static final double EARTH_RADIUS = 6378.137; 第二步：创建一个封装经纬度的VO类 /** * @Author:
【百度地图API】自行获取区域经纬度的工具
weixin_30301183的博客
06-20 80
行政区域工具 转载于:https://www.cnblogs.com/bober/p/3146101.html
坐标转换工具(地理坐标经纬度)
09-19
北京54，西安80，WGS等，3度，6度带转换
天地图连续拾取经纬度坐标
05-19
天地图连续拾取经纬度坐标
根据经纬度计算两点间距离
midpzc的博客
10-24 445
SELECT *, ROUND(6378.138*2*ASIN(SQRT(POW(SIN((36.072785*PI()/180-lat*PI()/180)/2),2)+COS(36.072785*PI()/180)*COS(lat*PI()/180)*POW(SIN((103.79748*PI()/180-lon*P
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值