1、电子围栏

男娘i 2021-11-04 10:46 306阅读 0赞

电子围栏是用来判定点是否在一定的范围之内的程序  
1、其中电子围栏分为不规则电子围栏和规则的电子围栏，规则的电子围栏是指，圆形的电子围栏，只要设置了中心点和半径即可划分一个电子围栏，其余的还有不规则的电子围栏，比如矩形、长方形、不规则图形的电子围栏，例如这种电子围栏主要采用描边的形式进行电子围栏的划分。

1、不规则电子围栏代码

def ispointinpolygon(point, rangelist):
        """
        电子围栏程序
        point:目标点坐标：如point=[1.5,3]
        rangelist:电子围栏的数组：如rangelist=[[118.804556,32.257551],[118.811950,32.251226],[118.814639,32.2490114]
        [118.793418,32.247962],[118.804556,32.257551]].注意：电子围栏的数组第一个坐标必须和最后一个坐标相似     
        return：点是否在围栏内,是的话返回True,否则返回False
        """
        ##经度列表 维度列表
        lnglist=[]
        latlist=[]
        for i in range(len(rangelist)-1):
            lnglist.append(rangelist[i][0])
            latlist.append(rangelist[i][1])
        #print(lnglist, latlist)
        count = 0 
        point1 = rangelist[0]
        for i in range(1, len(rangelist)):
            point2 = rangelist[i]
           # 判断线段两端点是否在射线两侧 不在肯定不相交 射线（-∞，lat）（lng,lat）
            if(point1[0]<=point[0]<=point2[0]) or (point2[0]<=point[0]<=point1[0]):
                target_y=(point2[1]-point1[1])/(point2[0]-point1[0])*(point[0]-point2[0])+point2[1]
                if target_y<=point[1]:
                    count=count+1
            point1=point2
        #print(count)
        if count % 2 == 1:
           #print("在圈内")
            return True
        else:
            if ((point[0] == point1[0] and point[1] == point1[1]) or (point[0] == point2[0] and point[1] == point2[1])):
                return True
            else:
                return False
                ##所有点
                
    ##测试数据
    #1、广东省公安厅交通管理局
    rangelist=[[113.296193,23.229597],[113.296193,23.229597],[113.29582,23.230884],
        [113.296553,23.230722],[113.296193,23.229597]]
    point1=[113.29622,23.230394]
    point2=[113.296831,23.230581]   
    ispointinpolygon(point1, rangelist)
    ispointinpolygon(point2, rangelist)

\#其中的rangelist为电子围栏的外围经纬度，point为标记点。

2、规则电子围栏代码  
主要是圆形的电子围栏

3、坐标系的转换  
坐标主要分为gps坐标，高德坐标、百度坐标

import numpy as np
    x_pi = 3.14159265358979324 * 3000.0 / 180.0
    pi = 3.1415926535897932384626 # π
    a = 6378245.0 # 长半轴
    ee = 0.00669342162296594323 # 偏心率平方
    #转换维度
    def _transformlat(lng, lat):
      ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + \
         0.1 * lng * lat + 0.2 *np.sqrt(np.fabs(lng))
      ret += (20.0 *np.sin(6.0 * lng * pi) + 20.0 *
         np.sin(2.0 * lng * pi)) * 2.0 / 3.0
      ret += (20.0 *np.sin(lat * pi) + 40.0 *
         np.sin(lat / 3.0 * pi)) * 2.0 / 3.0
      ret += (160.0 *np.sin(lat / 12.0 * pi) + 320 *
         np.sin(lat * pi / 30.0)) * 2.0 / 3.0
      return ret
    # 转换经度
    def _transformlng(lng, lat):
      ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \
         0.1 * lng * lat + 0.1 *np.sqrt(np.fabs(lng))
      ret += (20.0 *np.sin(6.0 * lng * pi) + 20.0 *
         np.sin(2.0 * lng * pi)) * 2.0 / 3.0
      ret += (20.0 *np.sin(lng * pi) + 40.0 *
         np.sin(lng / 3.0 * pi)) * 2.0 / 3.0
      ret += (150.0 *np.sin(lng / 12.0 * pi) + 300.0 *
         np.sin(lng / 30.0 * pi)) * 2.0 / 3.0
      return ret
    
    ###wgs84（gps坐标）转高德坐标
    def wgs84_to_gcj02(lng, lat):
      """
      WGS84转GCJ02(火星坐标系)
      :param lng:WGS84坐标系的经度
      :param lat:WGS84坐标系的纬度
      :return:高德地图的坐标
      经检验有效
      """
      dlng=np.array(lng)
      dlat=np.array(lat)
      dlat = _transformlat(lng - 105.0, lat - 35.0)
      dlng = _transformlng(lng - 105.0, lat - 35.0)
      radlat = lat / 180.0 * pi
      magic =np.sin(radlat)
      magic = 1 - ee * magic * magic
      sqrtmagic =np.sqrt(magic)
      dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
      dlng = (dlng * 180.0) / (a / sqrtmagic *np.cos(radlat) * pi)
      mglat = lat + dlat
      mglng = lng + dlng
      return [mglng, mglat]
    
    ###百度坐标转换为高德坐标
    def bd09_to_gcj02(bd_lon, bd_lat):
      """
      百度坐标系(BD-09)转火星坐标系(GCJ-02)
      百度——>谷歌、高德
      :param bd_lat:百度坐标纬度
      :param bd_lon:百度坐标经度
      :return:转换后的坐标列表形式
      """
      bd_lon=np.array(bd_lon)
      bd_lat=np.array(bd_lat)
      x = bd_lon - 0.0065
      y = bd_lat - 0.006
      z =np.sqrt(x * x + y * y) - 0.00002 *np.sin(y * x_pi)
      theta =np.arctan2(y, x) - 0.000003 *np.cos(x * x_pi)
      gg_lng = z *np.cos(theta)
      gg_lat = z *np.sin(theta)
      return [gg_lng, gg_lat]
    ##高德坐标转换为gps坐标
    def gcj02_to_wgs84(lng, lat):
      """
      GCJ02(火星坐标系)转GPS84
      :param lng:火星坐标系的经度
      :param lat:火星坐标系纬度
      :return:
      """
      dlat = _transformlat(lng - 105.0, lat - 35.0)
      dlng = _transformlng(lng - 105.0, lat - 35.0)
      radlat = lat / 180.0 * pi
      magic = np.sin(radlat)
      magic = 1 - ee * magic * magic
      sqrtmagic = np.sqrt(magic)
      dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * pi)
      dlng = (dlng * 180.0) / (a / sqrtmagic * np.cos(radlat) * pi)
      mglat = lat + dlat
      mglng = lng + dlng
      return [lng * 2 - mglng, lat * 2 - mglat]
    
    ##百度坐标转化为gps坐标
    def bd09_to_wgs84(bd_lon, bd_lat):
      lon, lat = bd09_to_gcj02(bd_lon, bd