• 根据经纬度以及距离求周边正方形区间


    1、概念

    纬线、经线、维度、经度

    2、根据两点的经纬度,计算距离

    /**
         * <p>Description: 计算两点之间距离</p>
         * @param slng 起始经度
         * @param slat 起始纬度
         * @param elng 结束经度
         * @param elat 结束纬度
         * @return*/
        public static double getDistanceByGaoDe(double slng, double slat, double elng, double elat) {
            double lng1 = (Math.PI / 180) * slng;
            double lng2 = (Math.PI / 180) * elng;
            double lat1 = (Math.PI / 180) * slat;
            double lat2 = (Math.PI / 180) * elat;
    
            // 两点间距离 km,如果想要米的话,结果*1000就可以了
            double d = Math.acos(
                    Math.sin(lat1) * Math.sin(lat2) 
                    + Math.cos(lat1) * Math.cos(lat2) * Math.cos(lng2 - lng1)
                ) * EARTH_RADIUS;
            // 精度1位小数
            d = Math.round(d * 10000d) / 10000d;
            
            return d * 1000;
        }

    3、根据当前经纬度以及距离,求周边正方形区间(通过求出的区间范围作为SQL入参过滤数据)

    /**
         * <p>Description: 计算两点之间距离</p>
         * @param dis 距离
         * @param lng 起始经度
         * @param lat 起始纬度
         * @return 正方形范围
         */
        public static double[] getPoiRange(double dis, double lng, double lat) {
            double[] poi = new double[4];
            
            double lng1 = (Math.PI / 180) * lng;
            double lng2 = (Math.PI / 180) * lng;
            double lat1 = (Math.PI / 180) * lat;
            double lat2 = (Math.PI / 180) * lat;
            
            // 角度
            double theta = dis / EARTH_RADIUS;
            
            // 维度相同,反算经度
            double lngVal = Math.cos(theta) - Math.sin(lat1) * Math.sin(lat2);
            double lngThetaVal = lngVal / (Math.cos(lat1) * Math.cos(lat2));
            double elng1 = (-Math.acos(lngThetaVal) + lng1) / (Math.PI / 180);
            double elng2 = (Math.acos(lngThetaVal) + lng1) / (Math.PI / 180);
            
            LogLIB.LogRun.info("elng1=" + UtilTools.format(elng1) + " elng2=" + UtilTools.format(elng2));
            
            // 经度相同,反算维度(根据辅助角公式计算 asinx + bcosx)
            // theta = a * Math.sin(lat2) + b * Math.cos(lat2) = Math.sqrt(a平方 + b平方)sin(x + arctan(b/a))
            double a = Math.sin(lat1);
            double b = Math.cos(lat1) * Math.cos(lng2 - lng1);
            double sqrt = Math.sqrt(Math.pow(a, 2) + Math.pow(b, 2));
            double latTheteVal = Math.cos(theta) / sqrt;
            double elat1 = (Math.asin(latTheteVal) - Math.atan(b/a)) / (Math.PI / 180);
            double elat2 = (Math.asin(-latTheteVal) - Math.atan(b/a)) / (Math.PI / 180) + 180;
            
            LogLIB.LogRun.info("elat1=" + UtilTools.format(elat1) + " elat2=" + UtilTools.format(elat2));
            
            poi[0] = elng1;
            poi[1] = elng2;
            poi[2] = elat1;
            poi[3] = elat2;
            
            return poi;
        }
  • 相关阅读:
    安装VC6.0遇到的问题
    开发、测试环境
    OPENGL绘制文字
    C++实现文件关联
    MFC多国语言——配置文件
    MFC 资源记录
    如何解决——汉化英文界面出现乱码
    项目配置——添加第三方资源
    队列&生产者消费者模型
    抢票小程序
  • 原文地址:https://www.cnblogs.com/feong/p/7251403.html
Copyright © 2020-2023  润新知