判斷坐標是否在某一個區域內（圓形，多邊形）

import java.awt.geom.Point2D;
import java.util.ArrayList;
import java.util.List;

public class regionUtil {
     private static double EARTH_RADIUS = 6378137;
   //private static double EARTH_SEA = 1.852;海里
private static double rad(double d) {
        return d * Math.PI / 180.0;
   }

    /**
     * 判斷是否在多邊形區域內
     *
     * @param pointLon 要判斷的點的橫坐標 經度
     * @param pointLat 要判斷的點的縱坐標 維度
     * @param lon      區域各頂點的橫坐標數組
     * @param lat      區域各頂點的縱坐標數組
     * @return
*/
public static boolean isInPolygon(double pointLon, double pointLat, Double[] lon,
Double[] lat) {
        // 將要判斷的橫縱坐標組成一個點
Point2D.Double point = new Point2D.Double(pointLon, pointLat);
// 將區域各頂點的橫縱坐標放到一個點集合里面
List<Point2D.Double> pointList = new ArrayList<Point2D.Double>();
        double polygonPoint_x = 0.0, polygonPoint_y = 0.0;
        for (int i = 0; i < lon.length; i++) {
            polygonPoint_x = lon[i];
polygonPoint_y = lat[i];
Point2D.Double polygonPoint = new Point2D.Double(polygonPoint_x, polygonPoint_y);
pointList.add(polygonPoint);
}
        return check(point, pointList);
}

    /**
     * 多邊形判斷
     *
     * @param point   要判斷的點的橫縱坐標
     * @param polygon 組成的頂點坐標集合
     * @return
*/
private static boolean check(Point2D.Double point, List<Point2D.Double> polygon) {
        java.awt.geom.GeneralPath peneralPath = new java.awt.geom.GeneralPath();

Point2D.Double first = polygon.get(0);
// 通過移動到指定坐標（以雙精度指定），將一個點添加到路徑中
peneralPath.moveTo(first.x, first.y);
polygon.remove(0);
        for (Point2D.Double d : polygon) {
            // 通過繪制一條從當前坐標到新指定坐標（以雙精度指定）的直線，將一個點添加到路徑中。
peneralPath.lineTo(d.x, d.y);
}
        // 將幾何多邊形封閉
peneralPath.lineTo(first.x, first.y);
peneralPath.closePath();
// 測試指定的 Point2D 是否在 Shape 的邊界內。
return peneralPath.contains(point);
}

    /**
     * 通過經緯度獲取距離(單位：米)
     *
     * @param lat1 緯度1
     * @param lng1 經度1
     * @param lat2 緯度2
     * @param lng2 經度2
     * @return 距離
     */
public static double getDistance(double lat1, double lng1, double lat2,
                                     double lng2) {
        double radLat1 = rad(lat1);
        double radLat2 = rad(lat2);
        double a = radLat1 - radLat2;
        double b = rad(lng1) - rad(lng2);
        double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) +
                Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(b / 2), 2)));
s = s * EARTH_RADIUS;
s = Math.round(s * 10000d) / 10000d;
//double len=s/EARTH_SEA;
        //s = s / 1000;
return s;
}

    /**
     *  圓形區域判斷
     * @param lat1   緯度
     * @param lat2   緯度
     * @param lng1   經度
     * @param lng2   經度
     * @param radius 判斷一個點是否在圓形區域內,比較坐標點與圓心的距離是否小於半徑
     */
public static boolean isInCircle(double lng1, double lat1, double lng2, double lat2, double radius) {
        double distance = getDistance(lat1, lng1, lat2, lng2);
        if (distance > radius) {
            return false;
} else {
return true;
}
    }

}