根据经纬度,判断点位是否在范围内(多边形内)

根据经纬度,判断点位是否在范围内(多边形内)

问题描述:

  • 有个需求,需要判断一个点位,在多边形内,并展示出来,研究了许多方法,还是这个射线法比较简单好用

  • 数据库的 coordinate多边形字段存储如下:

    129.26203608256768,26.085271312376 129.26801230123188,26.085807124332913 129.26858068565312,26.0231714249 129.26524138157777,26.0780112349326908 129.2601231480449,26.08158773429283 129.26203608256768,26.081231310453276
    

    由多个经纬度,组成一个多边形的形状,逗号区分经纬度,空格隔开多个坐标的存储方式(不同地方不一样)

最终效果如下:

在这里插入图片描述

代码部分:

  • 代码思路:1、写好判断点位的方法 2、切割多边形,并存入List集合

1、判断是否在点位的方法( 用于判断点 (x0, y0) 是否在给定多边形 ps 内部的方法 )

    /**
     * 判断点是否在范围内 , (point类的代码在下面)
     * @param x0
     * @param y0
     * @param ps - > ps集合 存放coordiante切割的多个经纬度点
     * @return
     */
    public  boolean contations(double x0, double y0, List<Point> ps) {
        int crossings = 0;
        for (int i = 0; i < ps.size() - 1; i++) {
            double slope = (ps.get(i + 1).getY() - ps.get(i).getY()) / (ps.get(i + 1).getX() - ps.get(i).getX());
            boolean cond1 = (ps.get(i).getX() <= x0) && (x0 < ps.get(i + 1).getX());
            boolean cond2 = (ps.get(i + 1).getX() <= x0) && (x0 < ps.get(i).getX());
            boolean above = (y0 < slope * (x0 - ps.get(i).getX()) + ps.get(i).getY());
            if ((cond1 || cond2) && above) {
                crossings++;
            }
        }
        return (crossings % 2 != 0);
    }
    

2、Point类(简易版)

@data
public class Point {

    public Point(double x0, double y0){
        this.x = x0;
        this.y = y0;
    }
    /**
     * 
     */
    private double x;
    
    /**
     * 
     */
    private double y;
}

3、切割coordinate 并存入List集合去判断

 List<Point> points = new ArrayList<Point>();

 if (StringUtils.isNotBlank(coor)) {
            String[] cs = coor.split(" ");
            for (int i = 0; i < cs.length; i++) {
                String s = cs[i];
                if(StringUtils.isBlank(s)){
                    continue;
                }
                String[] xy = s.split(",");
                Double x = new Double(xy[0]);
                Double y = null;
                try {
                    y = new Double(xy[1]);
                } catch (Exception e) {
                    //System.out.println(xy[1]);
                    //System.out.println(l.getStr("areacode"));
                    throw new RuntimeException(e.getMessage(), e);
                }
                //所有的点 都存入这里
                points.add(new Point(x, y));
            }
        }

4、最后的判断 是否在点位内


if( contations( lon,lat,points ) ){
	sout("存在多变形内");
}else{
    sout("不存在多变形内");
}
可以使用 Python 的 Shapely 库来实现判断某个经纬度是否在某个多边形范围内的功能。具体步骤如下: 1. 安装 Shapely 库:使用 pip 命令安装 Shapely 库,命令如下: ``` pip install shapely ``` 2. 导入 Shapely 库:在 Python 代码中导入 Shapely 库,命令如下: ``` from shapely.geometry import Point, Polygon ``` 3. 创建多边形对象:使用 Shapely 库的 Polygon 类创建多边形对象,多边形对象需要传入多个点的坐标,格式为 (x, y)。 ``` polygon = Polygon([(x1, y1), (x2, y2), (x3, y3), ...]) ``` 4. 创建点对象:使用 Shapely 库的 Point 类创建点对象,点对象需要传入点的坐标,格式为 (x, y)。 ``` point = Point(x, y) ``` 5. 判断是否多边形内:使用 Shapely 库的 contains 方法判断是否多边形内,在多边形内返回 True,不在多边形内返回 False。 ``` if polygon.contains(point): print("点在多边形内") else: print("点不在多边形内") ``` 完整的示例代码如下: ``` from shapely.geometry import Point, Polygon # 创建多边形对象 polygon = Polygon([(116.397, 39.908), (116.410, 39.908), (116.410, 39.920), (116.397, 39.920)]) # 创建点对象 point = Point(116.403, 39.914) # 判断是否多边形内 if polygon.contains(point): print("点在多边形内") else: print("点不在多边形内") ``` 注意:在使用 Shapely 库时,经度和纬度的顺序需要注意,通常情况下经度在前,纬度在后,例如 (116.397, 39.908)
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值