计算已知经纬度的两点距离（两种方法GeoTools和Haversine公式）

计算已知经纬度的两点距离（两种方法）

法一：GeoTools

要使用GeoTools，你需要在Maven项目中添加以下依赖坐标：

<dependency>
    <groupId>org.geotoolsgroupId>
    <artifactId>gt-mainartifactId>
    <version>24.0version>
dependency>

这个依赖坐标指定了GeoTools的主要模块gt-main以及对应的版本号24.0。请确保在你的Maven项目的pom.xml文件中添加了上述坐标，并根据需要修改版本号。

在添加了依赖坐标后，Maven会自动下载相应的GeoTools库，并将其包含在你的项目中。这样，你就可以使用GeoTools提供的类和方法来进行经纬度计算等地理信息处理任务。

使用GeoTools来计算经纬度之间的距离非常方便。下面是一个使用GeoTools的示例代码来计算经纬度之间的距离：

import org.geotools.referencing.GeodeticCalculator;
import org.opengis.referencing.crs.CoordinateReferenceSystem;

public class CoordinateDistanceCalculator {

    public static void main(String[] args) {
        double lat1 = 40.7128; // 起始点纬度
        double lon1 = -74.0060; // 起始点经度
        double lat2 = 34.0522; // 目标点纬度
        double lon2 = -118.2437; // 目标点经度

        // 创建CoordinateReferenceSystem对象，表示经纬度坐标系（通常是WGS84坐标系）
        CoordinateReferenceSystem crs = org.geotools.referencing.crs.DefaultGeographicCRS.WGS84;

        // 创建GeodeticCalculator对象
        GeodeticCalculator calculator = new GeodeticCalculator(crs);

        // 设置起始点
        calculator.setStartingGeographicPoint(lon1, lat1);

        // 设置目标点
        calculator.setDestinationGeographicPoint(lon2, lat2);

        // 计算两点之间的距离（返回结果的单位为米）
        double distance = calculator.getOrthodromicDistance();

        System.out.println("距离： " + distance + " 米");
    }
}

上述示例代码演示了如何使用GeoTools的GeodeticCalculator来计算两个经纬度坐标之间的距离。首先，我们创建了一个CoordinateReferenceSystem（坐标参考系统）对象，表示使用的经纬度坐标系（通常是WGS84坐标系）。然后，我们创建了一个GeodeticCalculator对象，并设置起始点和目标点的经纬度坐标。最后，调用getOrthodromicDistance方法计算两点之间的距离，结果以米为单位。

请确保在项目中包含GeoTools库的依赖，以便成功编译和运行上述示例代码。

法二：Haversine公式

使用Haversine公式：Haversine公式是一种常用的方法，可以在不考虑地球的不规则性的情况下，通过经纬度来计算两点之间的直线距离。公式如下：

a = sin²(Δlat/2) + cos(lat1) * cos(lat2) * sin²(Δlon/2)
c = 2 * atan2(√a, √(1-a))
d = R * c

其中，lat1和lon1表示起始点的纬度和经度，lat2和lon2表示目标点的纬度和经度，Δlat = lat2 - lat1，Δlon = lon2 - lon1，R为地球半径（一般取6371公里）。

public class HaversineFormula {
    private static final double EARTH_RADIUS = 6371; // 地球半径，单位：千米

    public static double calculateDistance(double lat1, double lon1, double lat2, double lon2) {
        double dLat = Math.toRadians(lat2 - lat1);
        double dLon = Math.toRadians(lon2 - lon1);

        double a = Math.sin(dLat / 2) * Math.sin(dLat / 2)
                + Math.cos(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2))
                * Math.sin(dLon / 2) * Math.sin(dLon / 2);

        double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));

        double distance = EARTH_RADIUS * c;
        return distance;
    }

    public static void main(String[] args) {
        double lat1 = 40.7128; // 纬度1
        double lon1 = -74.0060; // 经度1
        double lat2 = 34.0522; // 纬度2
        double lon2 = -118.2437; // 经度2

        double distance = calculateDistance(lat1, lon1, lat2, lon2);
        System.out.println("Distance between the two locations: " + distance + " km");
    }
}