java坐标代码_Java实现地理坐标判断

Java实现地理坐标判断

在做一个海洋相关的项目时候，遇到这个需求，判断船只是否在雷达监控的区域内，或者是否在某个保护区内。抽象为在地理位置中，一个坐标点是否在由一个圆心和半径的圆区域内，一个点是否在由一组点构成的区域内。

分析

多边形和圆分开写，首先简单的就是判断是否在圆里面，如何判断一个坐标是否在圆形区域内，相信不用我说都知道，计算这个坐标点和圆心之间的距离，然后跟圆的半径进行比较，如果比半径大，就不在圆形区域内，如果小于等于圆的半径，则该坐标点在圆形区域内。

点在圆区域内

数学上的计算公式是这样的：

代码采用谷歌地图计算距离的方式，应该算是比较精确。

private static double EARTH_RADIUS = 6378.137;

private static double rad(double d) {

return d * Math.PI / 180.0;

}

/**

* 通过经纬度获取距离(单位：米)

*

* @param lat1

* @param lng1

* @param lat2

* @param lng2

* @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;

return s;

}

/**

* 判断一个点是否在圆形区域内

*/

public static boolean isInCircle(double lng1, double lat1, double lng2, double lat2, String radius) {

return getDistance(lat1, lng1, lat2, lng2) > Double.parseDouble(radius);

}

如果使用Math.hypot()方法，计算(经纬度距离时)结果会有偏差)：

double x = (lon1 - lon2) * PI * R * Math.cos(((lat1 + lat2) / 2) * PI / 180) / 180;

double y = (lat1 - lat2) * PI * R / 180;

double distance = Math.hypot(x, y);

点在多边形区域内

这里用到JAVA的关于坐标系和几何图形的一个类GeneralPath，使用这个类，结合传入的各顶点参数，画一个几何图形，并通过它自身的contains方法，判断一点是否在这个几何图形内。

也就是，通过JAVA已经封装好的方法，画一个几何多边形，判断一点是否在这个几何多边形里面。

/**

* 判断是否在多边形区域内

*

* @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 pointList = new ArrayList();

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 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 pointLon

* 要判断的点的纵坐标

* @param pointLat

* 要判断的点的横坐标

* @param lon

* 指定区域的纵坐标组成的数组

* @param lat

* 指定区域的横坐标组成的数组

* @return

*/

private static boolean isInMaxArea(double pointLon, double pointLat, double[] lon,

double[] lat) {

// 获取区域横纵坐标最大值和最小值

double temp = 0.0;

for (int i = 0; i < lon.length; i++) {

for (int j = 0; j < lon.length - i - 1; j++) {

if (lon[j] > lon[j + 1]) {

temp = lon[j];

lon[j] = lon[j + 1];

lon[j + 1] = temp;

}

}

}

for (int i = 0; i < lat.length; i++) {

for (int j = 0; j < lat.length - i - 1; j++) {

if (lat[j] > lat[j + 1]) {

temp = lat[j];

lat[j] = lat[j + 1];

lat[j + 1] = temp;

}

}

}

// 如果在最值组成的区域外，那肯定不在重点区域内

return (pointLon < lon[0] || pointLon > lon[lon.length - 1] || pointLat < lat[0]

|| pointLat > lat[lat.length - 1]);

}

如果通过了上面的判断，可以进行接下来的算法判断了

用到了两点间的斜率公式

这个方法就是，通过一点，画一条线，这条线与多边形相交，如果相交点数位奇数，就在区域内，如果为偶数，就不在区域内

代码：

/**

* 判断坐标是否在重点区域内

*

* @param pointLon

* 要判断的点的纵坐标

* @param pointLat

* 要判断的点的横坐标

* @param lon

* 指定区域的纵坐标组成的数组

* @param lat

* 指定区域的横坐标组成的数组

* @return

*/

private static boolean isInAccurateArea(double pointLon, double pointLat, double[] lon,

double[] lat) {

// 代表有几个点

int vertexNum = lon.length;

boolean result = false;

for (int i = 0, j = vertexNum - 1; i < vertexNum; j = i++) {

// 满足条件，与多边形相交一次，result布尔值取反一次，奇数个则在区域内

if ((lon[i] > pointLon) != (lon[j] > pointLon)

&& (pointLat < (lat[j] - lat[i]) * (pointLon - lon[i]) / (lon[j] - lon[i])

+ lat[i])) {

result = !result;

}

}

return result;

}