java判断是否是经纬度_java实现判断一个经纬度坐标是否在一个多边形内(经自己亲测)...

public classTest01 {public static voidmain(String[] args) {//114.331951,30.64091#114.341049,30.610185#114.331436,30.588058#114.312038,30.56393#114.293498,30.558609#114.267922,30.563784#114.231185,30.57945#114.212303,30.601616#114.235649,30.626878#114.280624,30.646818#

Map [] map=newMap[]{};

Point[] ps= new Point[] { new Point(114.309914,30.599556),//114.309914,30.599556

new Point(114.295688,30.592879),//114.295688,30.592879

new Point(114.292812,30.587726), //114.292812,30.587726

new Point(114.292812,30.587726), //114.292812,30.587726

new Point(114.30058,30.580318),//114.30058,30.580318

new Point(114.303606,30.586959),//114.303606,30.586959

new Point(114.304534,30.594751),//114.304534,30.594751

new Point(114.30838,30.590131),//114.30838,30.590131

new Point(114.308651,30.584182),//114.308651,30.584182

new Point(114.304495,30.584015),//114.304495,30.584015

new Point(114.301301,30.578759),//114.301301,30.578759

new Point(114.309437,30.578528),//114.309437,30.578528

new Point(114.323282,30.592786)};//114.323282,30.592786

Point n1 = new Point(114.303217,30.583553);

Point n2= new Point(114.307336,30.597592);

Point n3= new Point(114.286565,30.590056);

Point y1= new Point(114.227342,30.587987);

Point y2= new Point(120.1866 , 30.2672);

Point y4= new Point(120.1869 , 30.2718);

System.out.println("n1:" +isPtInPoly(n1.getX() , n1.getY() , ps));

System.out.println("n2:" +isPtInPoly(n2.getX() , n2.getY() , ps));

System.out.println("n3:" +isPtInPoly(n3.getX() , n3.getY() , ps));

System.out.println("y1:" +isPtInPoly(y1.getX() , y1.getY() , ps));

System.out.println("y2:" +isPtInPoly(y2.getX() , y2.getY() , ps));

System.out.println("y4:" +isPtInPoly(y4.getX() , y4.getY() , ps));

}public static boolean isPtInPoly (double ALon , doubleALat , Point[] ps) {intiSum, iCount, iIndex;double dLon1 = 0, dLon2 = 0, dLat1 = 0, dLat2 = 0, dLon;if (ps.length < 3) {return false;

}

iSum= 0;

iCount=ps.length;for (iIndex = 0; iIndex

dLon1=ps[iIndex].getX();

dLat1=ps[iIndex].getY();

dLon2= ps[0].getX();

dLat2= ps[0].getY();

}else{

dLon1=ps[iIndex].getX();

dLat1=ps[iIndex].getY();

dLon2= ps[iIndex + 1].getX();

dLat2= ps[iIndex + 1].getY();

}//以下语句判断A点是否在边的两端点的水平平行线之间,在则可能有交点,开始判断交点是否在左射线上

if (((ALat >= dLat1) && (ALat < dLat2)) || ((ALat >= dLat2) && (ALat 0) {//得到 A点向左射线与边的交点的x坐标:

dLon = dLon1 - ((dLon1 - dLon2) * (dLat1 - ALat) ) / (dLat1 -dLat2);//如果交点在A点左侧(说明是做射线与 边的交点),则射线与边的全部交点数加一:

if (dLon

iSum++;

}

}

}

}if ((iSum % 2) != 0) {return true;

}return false;

}

}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
判断一个经纬度坐标是否在一个多边形内,可以使用Java的GIS相关的库,如GeoTools等。下面是一个基于GeoTools的实现示例: ```java import java.io.IOException; import java.util.ArrayList; import java.util.List; import org.geotools.geometry.jts.JTSFactoryFinder; import org.geotools.referencing.CRS; import org.locationtech.jts.geom.Coordinate; import org.locationtech.jts.geom.Geometry; import org.locationtech.jts.geom.GeometryFactory; import org.locationtech.jts.geom.Point; import org.locationtech.jts.geom.Polygon; import org.opengis.referencing.FactoryException; import org.opengis.referencing.crs.CoordinateReferenceSystem; import org.opengis.referencing.operation.MathTransform; import org.opengis.referencing.operation.TransformException; public class PolygonUtil { /** * 判断是否多边形内 * * @param lon 经度 * @param lat 纬度 * @param polygonCoords 多边形顶点坐标列表,顺序为顺时针或逆时针 * @param polygonSRID 多边形坐标系的SRID * @return true表示点在多边形内,false表示不在 */ public static boolean isPointInPolygon(double lon, double lat, List<double[]> polygonCoords, int polygonSRID) { try { // 创建多边形的JTS几何对象 GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory(null); Coordinate[] coords = new Coordinate[polygonCoords.size()]; for (int i = 0; i < polygonCoords.size(); i++) { double[] coord = polygonCoords.get(i); coords[i] = new Coordinate(coord[0], coord[1]); } Polygon polygon = geometryFactory.createPolygon(coords); // 创建点的JTS几何对象 Point point = geometryFactory.createPoint(new Coordinate(lon, lat)); // 获取多边形和点的坐标系 CoordinateReferenceSystem polygonCRS = CRS.decode("EPSG:" + polygonSRID); CoordinateReferenceSystem pointCRS = CRS.decode("EPSG:4326"); // WGS84经纬度坐标系 MathTransform transform = CRS.findMathTransform(pointCRS, polygonCRS); // 将点的坐标系转换为多边形坐标系 Point transformedPoint = (Point) org.geotools.geometry.jts.JTS.transform(point, transform); // 判断是否多边形内 return transformedPoint.within(polygon); } catch (FactoryException | TransformException e) { e.printStackTrace(); return false; } } public static void main(String[] args) throws IOException { List<double[]> polygonCoords = new ArrayList<>(); polygonCoords.add(new double[] { 116.35886, 39.89931 }); polygonCoords.add(new double[] { 116.36315, 39.89354 }); polygonCoords.add(new double[] { 116.37324, 39.89805 }); polygonCoords.add(new double[] { 116.36795, 39.90383 }); polygonCoords.add(new double[] { 116.35886, 39.89931 }); boolean inPolygon = PolygonUtil.isPointInPolygon(116.364, 39.896, polygonCoords, 4326); System.out.println(inPolygon); // true inPolygon = PolygonUtil.isPointInPolygon(116.373, 39.896, polygonCoords, 4326); System.out.println(inPolygon); // false } } ``` 注意事项: - 多边形顶点坐标列表的顺序要按照顺时针或逆时针方向排列,否则会得到错误的结果。 - 多边形坐标系的SRID要与经纬度坐标系的SRID不同,否则可能会得到错误的结果。 - 由于GeoTools的性能较低,处理大量数据时可能会较慢。如果需要高性能的GIS库,可以考虑使用JTS或Spatial4j等。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值