四、地图、java方式实现-根据距离范围查找该坐标内的其他坐标

实现思路：先定义PI和R地球半径

准备阶段：

public class CoordinateConversionUtils {

    private static final double PI = 3.14159265358979323;

    private static final double R = 6371229;

    /**
     * 计算两点距离
     * @param latLonReq1
     * @param latLonReq2
     * @return
     */
    public static double getDistance(LatLonReq latLonReq1, LatLonReq latLonReq2) {
        GlobalCoordinates source = new GlobalCoordinates(latLonReq1.getLatitude(), latLonReq1.getLongitude());
        GlobalCoordinates target = new GlobalCoordinates(latLonReq2.getLatitude(), latLonReq2.getLongitude());
        return new GeodeticCalculator().calculateGeodeticCurve(Ellipsoid.Sphere, source, target).getEllipsoidalDistance();
    }
}

计算两点的距离，这里使用了geodesy-1.1.3.jar这个jar库，需要的话，可以直接前往maven官网搜索添加对应的pom依赖。此插件的误差在0.1米内

我使用的百度坐标，百度坐标系 (BD-09) ，根据官网javascript的api进行测试，误差在0.1米内

1.根据纬度，先算出1纬度的距离，不同经度下的1纬度的距离都是相等的

double oneLatitudeDistance = PI * R / 180;

 

2.根据经度，先算出1经度的距离，这里先把纬度固定，然后再加上1经度，算出不同纬度下的1经度距离是多少米

　　　　　double latitude = latLon.getLatitude();//纬度
        double longitude = latLon.getLongitude();//经度
        double longitudeEnd = longitude + 1.0;//经度+1
        double oneLongitudeDistance = getDistance(latLon, new LatLonReq(latitude, longitudeEnd));//经度加1后，计算当前纬度下的1经度的距离

 

3.根据1纬度和1经度的距离，反推距离内，经度和纬度的最大值和最小值

    /**
     * latitude范围
     * @param latLon
     * @param range
     * @return
     */
    public static double[] getLatitudeRange(LatLonReq latLon, double range){
        double[] latitudeRange = {0.0,0.0};
        double oneLatitudeDistance = PI * R / 180;
        double latitude = latLon.getLatitude();
        double calculateLatitude = range / oneLatitudeDistance;
        latitudeRange[0] = latitude - calculateLatitude;
        latitudeRange[1] = latitude + calculateLatitude;
        return latitudeRange;
    }

    /**
     * longitude范围
     * @param latLon
     * @param range
     * @return
     */
    public static double[] getLongitudeRange(LatLonReq latLon, double range){
        double[] longitudeRange = {0.0,0.0};
        double latitude = latLon.getLatitude();
        double longitude = latLon.getLongitude();
        double longitudeEnd = longitude + 1.0;
        double oneLongitudeDistance = getDistance(latLon, new LatLonReq(latitude, longitudeEnd));
        double calculateLongitude = range / oneLongitudeDistance;
        longitudeRange[0] = longitude - calculateLongitude;
        longitudeRange[1] = longitude + calculateLongitude;
        return longitudeRange;
    }

 

4.根据上面的三个步骤，可推算出经纬度内的距离内的最大经度和纬度的范围，然后根据经纬度的范围可以形成一个类似正方形的经纬度范围

　数据库存储了多个点的坐标，我们需要查找某个点，范围内的其他点，此处可以直接通过sql进行查询

 SELECT qlp.cgi,qlp.lonc,qlp.latc FROM qj_lte_param qlp
        WHERE cast(qlp.lonc as decimal(9,6) ) BETWEEN #{minLon} AND #{maxLon} AND cast(qlp.latc as decimal(9,6)) BETWEEN #{minLat} AND #{maxLat}
          AND qlp.project_id = #{projectId}

具体代码步骤：先算出对应的经纬度范围，然后使用sql查询范围内的其他点

double[] latitudeRange = CoordinateConversionUtils.getLatitudeRangeByLatLon(latc, lonc, lteRange);
double[] longitudeRange = CoordinateConversionUtils.getLongitudeRangeByLatLon(latc, lonc, lteRange);
List<QjLteParam> plots = qjLteParamMapper.searchPlot(longitudeRange[0], longitudeRange[1],
                    latitudeRange[0], latitudeRange[1], projectId);

如果还需要计算距离的话，再使用getDistance方法，计算即可算出两点的距离

double distance = CoordinateConversionUtils.getDistance(new LatLonReq(Double.parseDouble(plot.getLatc()),
                        Double.parseDouble(plot.getLonc())), new LatLonReq(Double.parseDouble(latc), Double.parseDouble(lonc)));