Java计算两个坐标之间的距离

背景：

在做项目的时候需要做手机端页面，需要计算两个坐标的距离，由于之前没有接触过这样的计算，于是就找上度娘帮忙，顺便记录下来

1、添加JTS依赖

JTS是加拿大的 Vivid Solutions公司做的一套开放源码的 Java API。它提供了一套空间数据操作的核心算法。为在兼容OGC标准的空间对象模型中进行基础的几何操作提供2D空间谓词API。

<!--空间地理坐标算法-->
<dependency>
	<groupId>com.vividsolutions</groupId>
	<artifactId>jts-core</artifactId>
	<version>1.14.0</version>
</dependency>

2、坐标系转换

WGS-84：是国际标准，GPS坐标（Google Earth使用、或者GPS模块）
GCJ-02：中国坐标偏移标准，Google Map、高德、腾讯使用
BD-09：百度坐标偏移标准，Baidu Map使用

/**
 * @Description 坐标标转换工具类
 */
public class CoordinateUtil {
   
 
	/**
     * pi
     */
	public final static double pi = 3.1415926535897932384626;
	/**
     * a
     */
	public final static double a = 6378245.0;
	/**
     * ee
     */
	public final static double ee = 0.00669342162296594323;
 
	/**
	 * @Description WGS84 to 火星坐标系 (GCJ-02)
	 * @param lon
	 * @param lat
	 * @return
	 */
	public static double[] wgs84_To_Gcj02(double lon, double lat) {
   
		if (outOfChina(lat, lon)) {
   
			return null;
		}
		double dLat = transformLat(lon - 105.0, lat - 35.0);
		double dLon = transformLon(lon - 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);
		dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
		double mgLat = lat + dLat;
		double mgLon = lon + dLon;
		return new double[] {
    mgLon, mgLat };
	}
 
	/**
	 * @Description 火星坐标系 (GCJ-02) to WGS84
	 * @param lon
	 * @param lat
	 * @return
	 */
	public static double[] gcj02_To_Wgs84(double lon, double lat) {
   
		double[] gps = transform(lat, lon);
		double lontitude = lon * 2 - gps[1];
		double latitude = lat * 2 - gps[0];
		return new double[] {
    lontitude, latitude };
	}
 
	/**
	 * @Description 火星坐标系 (GCJ-02) to 百度坐标系 (BD-09)
	 * @param gg_lon
	 * @param gg_lat
	 * @return
	 */
	public static double[] gcj02_To_Bd09(double gg_lon, double gg_lat) {
   
		double x = gg_lon, y = gg_lat;
		double z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * pi);
		double theta = Math