java高德点到ian距离,高德地图（点到线段的最短距离算法）不调用高德API

package com.sage.dss.util;

import java.math.BigDecimal;

/**

* 提取自 AMapUtils(版本 com.amap.api:map2d:2.9.1)。

* 修改了参数，解除对高德 API 的依赖。

*/

public class MapUtil {

/**

* 此方法与高德API  distance方法算法一致

* 计算两点间的距离，单位：米

*/

public static BigDecimal calculateLineDistance(double x1, double y1, double x2, double y2) {

double lon1 = (Math.PI / 180) * x1;

double lon2 = (Math.PI / 180) * x2;

double lat1 = (Math.PI / 180) * y1;

double lat2 = (Math.PI / 180) * y2;

if (lon2 == lon1) {

return new BigDecimal(0);

}

// 地球半径

double R = 6378;

// 两点间距离 km，如果想要米的话，结果*1000就可以了

double d = Math.acos(Math.sin(lat1) * Math.sin(lat2) + Math.cos(lat1) * Math.cos(lat2) * Math.cos(lon2 - lon1)) * R;

return new BigDecimal(d);

}

double PointToSegDist(double x, double y, double x1, double y1, double x2, double y2){

double cross = (x2 - x1) * (x - x1) + (y2 - y1) * (y - y1);

if (cross <= 0){

return calculateLineDistance(x, y, x1, y1).doubleValue();

}

double d2 = (x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1);

if (cross >= d2){

return calculateLineDistance(x, y, x2, y2).doubleValue();

}

BigDecimal l1 = calculateLineDistance(x, y, x1, y1);

BigDecimal l2 = calculateLineDistance(x, y, x2, y2);

BigDecimal l12 = calculateLineDistance(x1, y1, x2, y2);

BigDecimal l = l1.add(l2).add(l12).divide(new BigDecimal(2), 6); // 半周长

double s = Math.sqrt(l.multiply(l.subtract(l12).multiply(l.subtract(l1).multiply(l.subtract(l2)))).doubleValue());

return new BigDecimal(s).multiply(new BigDecimal(2)).divide(l12, 6).doubleValue();

}

}