经纬度相关公式及实现

最新推荐文章于 2024-03-30 08:40:46 发布
最新推荐文章于 2024-03-30 08:40:46 发布
阅读量795 收藏 3
点赞数
文章标签： java
原文链接：http://www.cnblogs.com/shenjie2017/articles/8467274.html
版权

需求

  1. 已知两点经纬度，求两点间的距离
  2. 已知两点经纬度，求两点间的方位角
  3. 已知一点经纬度为中心、r为半径，求出该经纬度区间
  4. 已知一点的经纬度、距离，方位角，求另一个点的经纬度

公式

  1. 已知两点经纬度，求两点间的距离:

    AB=R*arccos[sin(wA)sin(wB)+cos(wA)cos(wB)*cos(jA-jB)]

  2. 已知两点经纬度，求两点间的方位角

    r*arccos[cos(y1)cos(y2)cos(x1-x2)+sin(y1)*sin(y2)]
    这里得到的是弧度，我们需要将弧度转换成我们需要的方位角，这里需要考虑象限问题
    纬度是横轴，经度时纵轴
    B点在第一象限，Azimuth=A；
    B点在第四象限，Azimuth=360+A;
    B点在第二三象限，Azimuth=180-A

  3. 已知一点经纬度为中心、r为半径，求出该经纬度区间

  4. 已知一点的经纬度、距离，方位角，求另一个点的经纬度

实现代码

import java.util.HashMap;
import java.util.Map;

public class DistanceUtil {     

    /*
    * 大地坐标系资料WGS-84 长半径a=6378137 短半径b=6356752.3142 扁率f=1/298.2572236
    */
    /** 长半径a=6378137 */
    private static double EARTH_RADIUS = 6378137;
    /** 短半径b=6356752.3142 */
    private static double b = 6356752.3142;
    /** 扁率f=1/298.2572236 */
    private static double f = 1 / 298.2572236;
    
    private static double rad(double d) {
           return d * Math.PI / 180.0;
    }
    
    /**
     * 弧度换成度
     *
     * @param x
     *弧度
     * @return 度
     */
    private static double deg(double x) {
         return x * 180 / Math.PI;
    }
    
    /**
     * 根据两个位置的经纬度，来计算两地的距离（单位为KM）
     * 参数为String类型
     * @param lat1 用户经度
     * @param lon1 用户纬度
     * @param lat2 商家经度
     * @param lon2 商家纬度
     * @return
     */
    public static double getDistance(int lat1, int lon1, int lat2, int lon2) {
        double tmplat1 = lat1/1000000.0;
        double tmplon1 = lon1/1000000.0;
        double tmplat2 = lat2/1000000.0;
        double tmplon2 = lon2/1000000.0;;
        
        double radLat1 = rad(tmplat1);
        double radLat2 = rad(tmplat2);
        double difference = radLat1 - radLat2;
        double mdifference = rad(tmplon1) - rad(tmplon2);
        double distance = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(difference / 2), 2)
        + Math.cos(radLat1) * Math.cos(radLat2)
        * Math.pow(Math.sin(mdifference / 2), 2)));
        distance = distance * EARTH_RADIUS;
        distance = Math.round(distance * 10000) / 10000;
        
        return distance;
    }
    
    /**
     * 获取当前用户一定距离以内的经纬度值,距离单位是米
     * 单位米 return minLat
     * 最小经度 minLon
     * 最小纬度 maxLat
     * 最大经度 maxLon
     * 最大纬度 minLat
     */
    public static Map<String,Integer> getAround(int lat, int lon, double raidus) {
        Map<String,Integer> map = new HashMap();
        
        double latitude = lat/1000000.0;// 传值给经度
        double longitude = lon/1000000.0;// 传值给纬度
        
        Double degree = (24901 * 1609) / 360.0; // 获取每度
        double raidusMile = raidus;
        
        Double mpdLon = Double.parseDouble((degree * Math.cos(latitude * (Math.PI / 180))+"").replace("-", ""));
        Double dpmLon = 1 / mpdLon;
        Double radiusLon = dpmLon * raidusMile;
        //获取最小经度
        double minLat = latitude - radiusLon;
        // 获取最大经度
        double maxLat = latitude + radiusLon;
        
        Double dpmLat = 1 / degree;
        Double radiusLat = dpmLat * raidusMile;
        // 获取最小纬度
        double minLon = longitude - radiusLat;
        // 获取最大纬度
        double maxLon = longitude + radiusLat;
        
        map.put("minLat", new Double(minLat*1000000).intValue());
        map.put("maxLat", new Double(maxLat*1000000).intValue());
        map.put("minLon", new Double(minLon*1000000).intValue());
        map.put("maxLon", new Double(maxLon*1000000).intValue());
        
        return map;
    }
    
    /**
     * 已知两个点的坐标，求两个的方位角
     * @param lat1 起始点
     * @param lon1
     * @param lat2 结束点
     * @param lon2
     * @return
     */
    public static double getAngle(int lat1, int lon1, int lat2, int lon2) {
        double tmplat1 = lat1/1000000.0;
        double tmplon1 = lon1/1000000.0;
        double tmplat2 = lat2/1000000.0;
        double tmplon2 = lon2/1000000.0;
        double dRotateAngle = Math.atan2(Math.abs(tmplat1 - tmplat2), Math.abs(tmplon1 - tmplon2));
        
        
        /*
        纬度是横轴，经度时纵轴
         */
        if (tmplat2 >= tmplat1){
            if (tmplon2 >=tmplon1){
                  //第一象限
            } else {
                  // 第四象限
                  dRotateAngle =2 * Math.PI - dRotateAngle;
            }
        } else{
            if (tmplon2 >= tmplon1) {
                 //第二象限
                  dRotateAngle = Math.PI - dRotateAngle;
            } else {
                 //第三象限
                 dRotateAngle = Math.PI + dRotateAngle;
            }
        }
        
        dRotateAngle = dRotateAngle * 180 / Math.PI;
        return dRotateAngle;
    }
    
    /**
     * 求两经纬度方向角
     *
     *@param lat1
     *第一点的纬度
     * @param lon1
     *第一点的经度
     * @param lat2
     *第二点的纬度
      @param lon2
     *第二点的经度
     * @return 方位角，角度（单位：°）
     * */
    public static double getAzimuth(int lat1,int lon1,int lat2,int lon2) {
        double tmplat1 = lat1/1000000.0;
        double tmplon1 = lon1/1000000.0;
        double tmplat2 = lat2/1000000.0;
        double tmplon2 = lon2/1000000.0;
        tmplat1 = rad(tmplat1);
        tmplat2 = rad(tmplat2);
        tmplon1 = rad(tmplon1);
        tmplon2 = rad(tmplon2);
        double azimuth = Math.sin(tmplat1) * Math.sin(tmplat2) + Math.cos(tmplat1)
        * Math.cos(tmplat2) * Math.cos(tmplon2 - tmplon1);
        azimuth = Math.sqrt(1 - azimuth * azimuth);
        azimuth = Math.cos(tmplat2) * Math.sin(tmplon2 - tmplon1) / azimuth;
        azimuth = Math.asin(azimuth) * 180 / Math.PI;
        
        if (Double.isNaN(azimuth)) {
            if (lon1 < lon2) {
                azimuth = 90.0;
            } else {
                azimuth = 270.0;
            }
        }
        
        /*
        纬度是横轴，经度时纵轴
        B点在第一象限，Azimuth=A；
        B在第四象限，Azimuth=360+A;
        B在第二三象限，Azimuth=180-A
         */
        
        if (tmplat2 >= tmplat1){
            if (tmplon2 >=tmplon1){
                //第一象限
            } else {
                // 第四象限
                azimuth =360 + azimuth;
            }
        } else{
            if (tmplon2 >= tmplon1) {
                 //第二象限
                azimuth = 180 - azimuth;
            } else {
                //第三象限
                 azimuth = 180 - azimuth;
            }
        }
        
        return azimuth;
    }
    
    /**
     * 计算另一点经纬度
     *
     * @param lon
     *经度
     * @param lat
     *维度
     * @param brng
     *方位角
     * @param dist
     *距离（米）
     */
    public static Map<String,Integer> computerThatLonLat(int lat,int lon, double brng, double dist) {
    
        double tmplon = lon / 1000000.0;
        double tmplat = lat / 1000000.0;
        
        double alpha1 = rad(brng);
        double sinAlpha1 = Math.sin(alpha1);
        double cosAlpha1 = Math.cos(alpha1);
        
        double tanU1 = (1 - f) * Math.tan(rad(tmplat));
        double cosU1 = 1 / Math.sqrt((1 + tanU1 * tanU1));
        double sinU1 = tanU1 * cosU1;
        double sigma1 = Math.atan2(tanU1, cosAlpha1);
        double sinAlpha = cosU1 * sinAlpha1;
        double cosSqAlpha = 1 - sinAlpha * sinAlpha;
        double uSq = cosSqAlpha * (EARTH_RADIUS * EARTH_RADIUS - b * b) / (b * b);
        double A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq)));
        double B = uSq / 1024 * (256 + uSq * (-128 + uSq * (74 - 47 * uSq)));
        
        double cos2SigmaM=0;
        double sinSigma=0;
        double cosSigma=0;
        double sigma = dist / (b * A), sigmaP = 2 * Math.PI;
        while (Math.abs(sigma - sigmaP) > 1e-12) {
        cos2SigmaM = Math.cos(2 * sigma1 + sigma);
        sinSigma = Math.sin(sigma);
        cosSigma = Math.cos(sigma);
        double deltaSigma = B * sinSigma * (cos2SigmaM + B / 4 * (cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM)
        - B / 6 * cos2SigmaM * (-3 + 4 * sinSigma * sinSigma) * (-3 + 4 * cos2SigmaM * cos2SigmaM)));
        sigmaP = sigma;
        sigma = dist / (b * A) + deltaSigma;
        }
        
        double tmp = sinU1 * sinSigma - cosU1 * cosSigma * cosAlpha1;
        double lat2 = Math.atan2(sinU1 * cosSigma + cosU1 * sinSigma * cosAlpha1,
        (1 - f) * Math.sqrt(sinAlpha * sinAlpha + tmp * tmp));
        double lambda = Math.atan2(sinSigma * sinAlpha1, cosU1 * cosSigma - sinU1 * sinSigma * cosAlpha1);
        double C = f / 16 * cosSqAlpha * (4 + f * (4 - 3 * cosSqAlpha));
        double L = lambda - (1 - C) * f * sinAlpha
        * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM)));
        
        double revAz = Math.atan2(sinAlpha, -tmp); // final bearing
        
        Map<String,Integer> map = new HashMap<String, Integer>();
        map.put("lon",+new Double((tmplon+deg(L))*1000000).intValue());
        map.put("lat",+new Double(deg(lat2)*1000000).intValue());
        return map;
    }
    
    /**
     * 判断是否漂移
     * @return
     */
    /**
     * 判断是否漂移，如果漂移则返回true,否则返回false
     * @param lat1 单位是10^-6
     * @param lon1 单位是10^-6
     * @param raidus 半径单位是米
     * @param lat2
     * @param lon2
     * @return
     */
    public static boolean isDrift(int lat1,int lon1,int raidus,int lat2,int lon2){
    
        //中国纬度范围:4°~53°
        if(lat1<4000000 || lat1 >53000000 ){
        return true;
        }
        //中国经度范围:73°~135°
        if(lon1<73000000 || lon1 >135000000){
        return true;
        }
        
        Map<String,Integer> map = getAround(lat1,lon1,raidus);
        
        if(lat2>map.get("maxLat") || lat2<map.get("minLat")){
            return true;
        }
        if(lon2>map.get("maxLon") || lon2<map.get("minLon")){
            return true;
        }

        return false;
    }

}

转载于:https://www.cnblogs.com/shenjie2017/articles/8467274.html

Niki20034655
关注
米转经纬度_经纬度换算米(经纬度精度换算米数)
weixin_39777464的博客
02-12 1万+
经度113度55分24秒。纬度22度29分3秒。经纬度的距离都是一度=111公里。一公里=1000米，经纬度分度、分、秒都是60进制。记忆中好象是每一度是1852米，也就是一海里.那么既然知道了一度有多长，你也知道度分秒的换算关系，分秒的长度就很好算了可以自己算的嘛，用地球周长除以360就是.jing'wei'du经纬度怎么转换成米呢这个需要具体的经纬度来计算的，运算方法很复杂。惯例上，经度1°等...
根据经纬度算出两点距离以及方位角
03-21
在网上找个好几种方法来根据经纬度算距离和方位角，但误差都比较多，这个是我个人优化的，通过测试误差比较像小，希望能帮到大家
高德转百度经纬度excel公式.xlsx
11-19
网上相关代码已经汗毛充栋了，但对于临时使用一次性使用来讲，写代码的效率显然是不高的。 自己根据代码整理了一个excel文档，见附件，可批量实现高德坐标点向百度坐标点的转换。 您值得下载！
经纬度相关计算公式
li10081325的博客
02-23 7473
一、 1、计算两点之间球面距离 matlab版 function [dist] = computeDistanceLatLon(lon1,lat1,lon2,lat2) pi = 3.14159265358979324; Rc = 6378137.0; Rj = 6356725.0; DEG2RAD = pi / 180.0; r_lat_1 = lat1 * DEG2RAD; ...
根据两点的经纬度求方位角和距离
热门推荐
WJsuperrunner
04-17 5万+
转自：http://blog.sina.com.cn/s/blog_5e7960620101vi0d.html 拿过来学习下；这里主要解决四个问题：1、已知两点经纬度，求两点间距离；2、已知两点经纬度，求一点相对于另一点航向；3、已知一点经纬度及与另一点距离和航向，求另一点经纬度；总声明：因为地球偏心率极低，所以此处将地球看做球体，以下所有公式在适用范围内使用可以保证计算结果与WGS84坐标系统...
计算两个经纬度点之间的方位角（初始方位角）
weixin_43578304的博客
03-30 2221
最近做地图上轨迹计算方位角的时候，需要计算两个经纬度点之间的方位角，这里把计算过程记录分享一下。
通过经纬度求解方位角
ControlLearner的博客
06-08 2万+
拟解决四个问题 1、已知两点经纬度，求一点相对于另一点方位角； 2、已知两点经纬度，求两点间距离； 3、已知一点经纬度及与另一点距离和方位角，求另一点经纬度； 4、问题1与问题2的简化算法。 注：简化算法的运算量和对系统的运算精度要求都大大降低，但只在短距离内（高纬地区建议10km以下）可以保证精度，除简化算法之外的算法可适用于地球上任意两点。这里只是出于便于理解的目的来解释“原理”，具体到不同的编程环境还要自己做化简和注意单位。 符号和单位约定 此处设定求B相对于A的方位角，即A为当前位置，B为目标位置
经纬度计算距离公式经纬度计算距离公式
最新发布
05-17
了解了经纬度计算距离的基本原理和公式，我们可以更有效地处理与地理定位相关的任务，如地图应用、导航系统、物流配送路线规划等。在实际使用时，还需要考虑其他因素，比如地球的非球形性质、地球自转、地形障碍等，...
java实现经纬度计算距离
03-23
这里我们将深入探讨如何利用Java实现这一功能，以及涉及到的相关数学原理。 首先，地球并非完美的球体，而是略微扁平的椭球体，但在小范围内，我们可以使用球面模型简化计算，即假设地球是一个完美的球体。在这种...
经纬度计算距离,经纬度计算距离公式,matlab源码.zip
10-15
本资源“经纬度计算距离,经纬度计算距离公式,matlab源码.zip”提供了一个基于MATLAB的实现，用于解决这个问题。下面将详细介绍经纬度计算距离的基本原理、常用的计算公式以及MATLAB源码的相关知识。 首先，我们了解...
根据经纬度距离和方位角计算新的经纬度
04-03
根据经纬度距离和方位角计算新的经纬度，计算经纬度两点间距离
地理经纬度弧度角度转换
04-29
支持弧度和角度的互相转化，可以弧度转角度，也可以角度转弧度
计算经纬度之间的距离（百度经纬度）
10-30
计算百度经纬度的公式，已经写在Excel里面了，直接输入两个点的经纬度就可以自动计算出公司
经纬度，弧度角度转换软件
02-26
度分秒转换器，可以实现经纬度不同单位进制的相互转换
tems经纬度转换公式
12-28
文件"tems 经纬度转换公式.xls"可能包含了具体的转换值或步骤，这些数据可以输入到GIS软件如QGIS、ArcGIS或MapInfo中，通过坐标变换工具实现批量转换。在实际应用中，用户需要确保输入数据的精度和单位与软件要求...
根据经纬度计算地球两个坐标之间的方位角
qiannianlaoyao2010的博客
10-29 1万+
近年来随着人工智能机器人、无人机、自动驾驶等领域技术的不断发展，卫星精准定位，自动导航成为了这些领域中的核心关键。随着人造卫星和互联网的发展，我们很容易就能通过卫星或基站信号获取个人设备在地球上所在的坐标，在利用姿态传感器，获取自身的航向角便可实时矫正自身航向往目的方向前进。而自动导航在人工智能领域的主要应用就是主动的寻找目的地坐标，并通过反馈执行相应运动抵达目的坐标。 要想实现自动导航，机器人首...
【GIS - 地理信息系统】经纬度计算 ( 经度、纬度概念 | 地球周长计算 | 地球经线周长计算 | 经纬度相关计算 | 经纬度坐标距离计算公式 | 经纬度与实际距离换算 | 相关代码实现 )
让 学习 成为一种 习惯 ( 韩曙亮 の 技术博客 )
08-04 1万+
一、经度、纬度概念 二、地球周长计算 1、地球半径、周长计算 2、地球经线周长计算 3、地球纬线周长计算 三、经纬度相关计算 1、经纬度坐标距离计算公式 2、经纬度与实际距离换算 1 米对应经度 1 米对应纬度 3、实际距离与经纬度换算 1 度经度对应东西距离 1 度纬度对应南北距离 四、相关代码实现 1、计算两个经纬度之间的距离 2、距离与经纬度之间换算 3、弧度与角度计算
通过地址获取经纬度的三种方法; 通过经纬度获取省市的方法
weixin_30412013的博客
04-22 690
# coding: utf-8 #通过地址获取经纬度的三种方法 #方法一 http://www.cnblogs.com/liangto/p/6287957.html中的方法二 import urllib2 import json address ='西安市高新区锦业路69号创业研发园' url = 'http://api.map.baidu.com/geocoder/v2/?...
已知两点经纬度坐标（大地坐标）计算方位角
gushangu的专栏
04-22 6107
使用经纬度坐标计算方位角
echarts经纬度计算距离公式
06-11
Echarts经纬度计算距离公式可以使用Haversine公式来实现。Haversine公式根据两个经纬度点的经度、纬度坐标，计算出这两个点之间的距离。下面是Haversine公式的实现： ``` function getDistance(lat1, lng1, lat2, lng2) { var radLat1 = lat1 * Math.PI / 180.0; var radLat2 = lat2 * Math.PI / 180.0; var a = radLat1 - radLat2; var b = lng1 * Math.PI / 180.0 - lng2 * Math.PI / 180.0; var s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a/2),2) + Math.cos(radLat1)*Math.cos(radLat2)*Math.pow(Math.sin(b/2),2))); s = s * 6378.137; // 地球半径，单位为公里 s = Math.round(s * 10000) / 10000; return s; } ``` 其中，lat1、lng1和lat2、lng2分别表示两个经纬度点的经度和纬度坐标，计算出的距离单位为公里。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值