python射线法检测点是否位于区域外接矩形内检测方法

时间：2020-04-15 08:26:38
3932人关注

这篇文章主要为大家详细介绍了python射线法判断检测点是否位于区域外接矩形内，具有一定的参考价值，感兴趣的小伙伴们可以参考一下，另外这篇文章主要知识点是关于python、射线法、python的内容，如果大家想对相关知识点有系统深入的学习，可以参阅以下电子资料：

教程详情	电子资料
教程类别：python 编辑整理：菱飞星教程字数：1800字节阅读时间：大概11分钟下载本教程（DOC版）	《编程的乐趣：用Python解算法谜题》解题代码《给Python点颜色：青少年学编程》配套资源《中学生Python创意编程》源码课件《Python快乐编程:中学学科创意编程实例》课件,源码,视频《用图形学Python 3》源码资源

本文实例为大家分享了python射线法判断点是否位于区域内的具体代码，供大家参考，具体内容如下

#!/usr/bin/env python
# -*- coding: utf-8 -*-
# @Date : 2018-10-07 15:49:37
# @Author : Sheldon (thisisscret@qq.com)
# @Blog : 谢耳朵的派森笔记
# @Link : https://www.cnblogs.com/shld/
# @Version : 0.0.1

def isinpolygon(point,vertex_lst:list, contain_boundary=True):
 #检测点是否位于区域外接矩形内
 lngaxis, lataxis = zip(*vertex_lst)
 minlng, maxlng = min(lngaxis),max(lngaxis)
 minlat, maxlat = min(lataxis),max(lataxis)
 lng, lat = point
 if contain_boundary:  
  isin = (minlng<=lng<=maxlng) & (minlat<=lat<=maxlat)
 else:
  isin = (minlng<lng<maxlng) & (minlat<lat<maxlat)
 return isin

def isintersect(poi,spoi,epoi):
 #输入：判断点，边起点，边终点，都是[lng,lat]格式数组
 #射线为向东的纬线
 #可能存在的bug，当区域横跨本初子午线或180度经线的时候可能有问题
 lng, lat = poi
 slng, slat = spoi
 elng, elat = epoi
 if poi == spoi:
  #print("在顶点上")
  return None
 if slat==elat: #排除与射线平行、重合，线段首尾端点重合的情况
  return False
 if slat>lat and elat>lat: #线段在射线上边
  return False
 if slat<lat and elat<lat: #线段在射线下边
  return False
 if slat==lat and elat>lat: #交点为下端点，对应spoint
  return False
 if elat==lat and slat>lat: #交点为下端点，对应epoint
  return False
 if slng<lng and elat<lat: #线段在射线左边
  return False
 #求交点
 xseg=elng-(elng-slng)*(elat-lat)/(elat-slat)
 if xseg == lng:
  #print("点在多边形的边上")
  return None
 if xseg<lng: #交点在射线起点的左侧
  return False
 return True #排除上述情况之后

def isin_multipolygon(poi,vertex_lst, contain_boundary=True): 
 # 判断是否在外包矩形内，如果不在，直接返回false 
 if not isinpolygon(poi, vertex_lst, contain_boundary):
  return False
 sinsc = 0  
 for spoi, epoi in zip(vertex_lst[:-1],vertex_lst[1::]):
  intersect = isintersect(poi, spoi, epoi)
  if intersect is None:
   return (False, True)[contain_boundary]
  elif intersect:
   sinsc+=1   
 return sinsc%2==1


if __name__ == '__main__':
 vertex_lst = [[0,0],[1,1],[1,2],[0,2],[0,0]]
 poi = [0.82,0.75]
 print(isin_multipolygon(poi,vertex_lst, contain_boundary=True))

以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持码农之家。

上一篇：如何用python实现按比例随机切分数据

下一篇：Python操作配置文件方法分享

学习笔记

23小时27分钟前回答

如何使用python 判断矩阵中每行非零个数

如下所示： # -*- coding: utf-8 -*-# @Time : 2018/5/17 15:05# @Author : Sizer# @Site : # @File : test.py# @Software: PyCharmimport timeimport numpy as np# data = np.array([# [5.0, 3.0, 4.0, 4.0, 0.0],# [3.0, 1.0, 2.0, 3.0, 3.0],# [4.0, 3.0, 4.0, 3.0, 5.0],# [3.0, 3.0, 1.0, 5.0, 4.0],# [1.0, 5.0, 5.0, 2.0, 1.0]# ])data = np.random.random((1000, 1000))print(data.shape)start_time = time.time()# avg = [float(np.mean(data[i, :])) for i in range(data.shape[0])]# print(avg)start_time = time.time()avg = []for i in range(data.shape[0]): sum = 0 cnt = 0 for rx in data[i, :]: if rx 0: sum += rx cnt += 1 if cnt 0: avg.append(sum/cnt) else: avg.append(0)end_time = time.time()print(op 1:, end_time - start_time)start_time = time.time()avg = []isexist = (data 0) * 1for i in range(data.shape[0]): sum =……

30小时13分钟前回答

Python判断字符串是否为字母或者数字方法总结

str为字符串s为字符串 str.isalnum() 所有字符都是数字或者字母 str.isalpha() 所有字符都是字母 str.isdigit() 所有字符都是数字 str.isspace() 所有字符都是空白字符、t、n、r 检查字符串是数字/浮点数方法 float部分 float(Nan)nan float(Nan)nan float(nan)nan float(INF)inf float(inf)inf float(-INF)inf float(-inf)inf 第一种：最简单 def is_number(str): try: # 因为使用float有一个例外是NaN if str==NaN: return False float(str) return True except ValueError: return Falsefloat例外示例 float(NaN) nan 使用complex() def is_number(s): try: complex(s) # for int, long, float and complex except ValueError: return False return True 综合1 def is_number(s): try: float(s) # for int, long and float except ValueError: try: complex(s) #……

43小时33分钟前回答

python判断设备联网的代码

本文实例为大家分享了python判断设备是否联网的具体代码，供大家参考，具体内容如下直接上代码，就是用判断socket能不连上的方法来判断。 #!/usr/bin/env python# -*- coding: utf-8 -*-import socketdef isNetOK(testserver): s=socket.socket() s.settimeout(3) try: status = s.connect_ex(testserver) if status == 0: s.close() return True else: return False except Exception as e: return Falsedef isNetChainOK(testserver=(www.baidu.com,443)): isOK = isNetOK(testserver) return isOKdef isNetUSAOK(testserver=(www.google.com,443)): isOK = isNetOK(testserver) return isOKdef isNetYouTubeOK(testserver=(www.youtube.com,443)): isOK = isNetOK(testserver) return isOKdef main(): chinanet = isNetChainOK() print chinanet usanet = isNetUSAOK() print usanet youtubenet……

32小时14分钟前回答

python判断文件是否存在的方法

通常在读写文件之前，需要判断文件或目录是否存在，不然某些处理方法可能会使程序出错。所以最好在做任何操作之前，先判断文件是否存在。 python中如何判断文件是否存在？判断文件是否存在可使用os模块，使用Try语句，使用pathlib模块。这里将介绍三种判断文件或文件夹是否存在的方法，分别使用os模块、Try语句、pathlib模块。 1.使用os模块 os模块中的os.path.exists()方法用于检验文件是否存在。判断文件是否存在 import osos.path.exists(test_file.txt)#Trueos.path.exists(no_exist_file.txt)#False 或者 import osos.path.exists(test_dir)#Trueos.path.exists(no_exist_dir)#False 可以看出用os.path.exists()方法，判断文件和文件夹是一样。其实这……