计算圆和多边形
package com.htstar.ovms.device.util;
import java.awt.geom.Point2D;
import java.util.ArrayList;
import java.util.List;
public class GeographyUtil {
// private static double EARTH_RADIUS = 6378.137;
private static double EARTH_RADIUS = 6371.393; // 地图半径
private static double rad(double d) {
return d * Math.PI / 180.0;
}
/**
* 是否在园内
*
* @param lat1
* @param lng1
* @param lat2
* @param lng2
* @return
*/
public static boolean getOutsideStatus(double lat1, double lng1, double lat2,
double lng2, Integer radius) {
double s = distance(lat1, lng1, lat2, lng2);
// 在圆外
if (s > radius) {
return false;
}
return true;
}
public static double distance(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 * 1000);
return s;
}
/**
* 判断是否在多边形区域内(正方形 测试通过)
*
* @param pointLon
* 要判断的点的纵坐标 (经度)
* @param pointLat
* 要判断的点的横坐标 (纬度)
* @param lon
* 区域各顶点的纵坐标数组
* @param lat
* 区域各顶点的横坐标数组
* @return
*/
public static boolean isInPolygon(double pointLon, double pointLat, String[] lon,
String[] 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++) {
for (int i1 = 0; i1 < lat.length; i1++) {
polygonPoint_x = Double.parseDouble(lon[i]);
polygonPoint_y = Double.parseDouble(lat[i1]);
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);
}
public static String getDistanceTest(double lat1, double lng1, double lat2,
double lng2, Integer radius) {
double s = distance(lat1, lng1, lat2, lng2);
return "围栏半径:" + radius + "_距离:" + s;
}
}