Python一行代码处理地理围栏

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了Python一行代码处理地理围栏，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含2115字，纯文字阅读大概需要4分钟。

内容图文

最近在工作中遇到了这个一个需求，用户设定地理围栏，后台获取到实时位置信息后通过与围栏比较，判断是否越界等。

这个过程需要用到数据协议为GEOjson，通过查阅资料后，发现python的shapely库可以非常简单的解决这个问题,接下来演示一下我处理这个问题的过程。

测试数据：

通过http://geojson.io/来获得测试数据，如下图，在地图上绘制了一个多边形设为地理围栏，分别取了围栏内外两个点来进行测试。

得到GEOjson数据如下：

{
  "type": "FeatureCollection",
  "features": [
    {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "type": "Polygon",
        "coordinates": [
          [
            [
              114.3458104133606,
              30.476167529462785
            ],
            [
              114.34512376785278,
              30.475575748963195
            ],
            [
              114.34576749801636,
              30.474540124433936
            ],
            [
              114.3467652797699,
              30.475363076967565
            ],
            [
              114.34693694114685,
              30.476102803645833
            ],
            [
              114.3458104133606,
              30.476167529462785
            ]
          ]
        ]
      }
    },
    {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "type": "Point",
        "coordinates": [
          114.34605717658997,
          30.475584995561178
        ]
      }
    },
    {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "type": "Point",
        "coordinates": [
          114.346604347229,
          30.476518897432545
        ]
      }
    }
  ]
}

安装shapely

本测试基于python——python3.6

$ pip install shapely

windows安装shapely会报错

shapely解析地理围栏

话不多说直接上代码

from shapely.geometry import Point
from shapely.geometry.polygon import Polygon

point = Point(0.5, 0.5)
polygon = Polygon([(0, 0), (0, 1), (1, 1), (1, 0)])
print(polygon.contains(point))

下面是实际的实例：

from shapely.geometry import Point
from shapely.geometry.polygon import Polygon

polygon_data= [
            [
              114.3458104133606,
              30.476167529462785
            ],
            [
              114.34512376785278,
              30.475575748963195
            ],
            [
              114.34576749801636,
              30.474540124433936
            ],
            [
              114.3467652797699,
              30.475363076967565
            ],
            [
              114.34693694114685,
              30.476102803645833
            ],
            [
              114.3458104133606,
              30.476167529462785
            ]
          ]
          
point1 = Point([114.34605717658997,30.475584995561178])
point2 = Point([114.346604347229,30.476518897432545])
polygon = Polygon(polygon_data)
print(polygon.contains(point1))
print(polygon.contains(point2))

输出结果：

True
False

这样一来我们就快速的实现了，目标点是否在地理围栏内的判断。

总结

Python还是挺好用的：）

参考资料：

https://stackoverflow.com/questions/36399381/whats-the-fastest-way-of-checking-if-a-point-is-inside-a-polygon-in-python

内容总结

以上是互联网集市为您收集整理的Python一行代码处理地理围栏全部内容，希望文章能够帮你解决Python一行代码处理地理围栏所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端