本文实例为大家分享了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))
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
华山资源网 Design By www.eoogi.com
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
华山资源网 Design By www.eoogi.com
暂无评论...
更新日志
2024年09月25日
2024年09月25日
- 纯音入心系列纯音乐《天籁古筝》1CD[MP3][331MB]
- 男女对唱典藏天碟《发烧对唱·那个季节里的歌DSD》10CD[WAV]
- 群星2010-歌林精选辑[歌林][WAV+CUE]
- TrioZimbalist-PianoTriosofWeinberg,Auerbach,Dvorak(2024)[24-44,1]wav
- 群星.1992-滚石第一流台湾歌·12王牌大车拼【滚石】【WAV+CUE】
- 林子祥.1984-创作歌集【华纳】【WAV+CUE】
- 周汤豪.2010-周汤豪【华纳】【FLAC分轨】
- Mozart-TheStringQuintets-AmadeusQuartet,CecilAronowitz(2017)[24-44,1][WAV+CUE]
- JamesWilliamsDennisIrwin-Focus(2024,Red)[24-48]FLAC
- 藤泽麻衣《空みあげて》[WAV]
- 纯音入心系列纯音乐《中国古筝经典名曲》1CD[MP3][1.3GB]
- 纯音入心系列纯音乐《古筝新奏:岁月静好与筝语,细水流年与筝同》1CD[MP3][846.9MB]
- 670《脆弱敏感小女生》[320K/MP3][27.53MB]
- 曾琳.2011-好好爱我【南方】【WAV+CUE】
- 范晓萱.1995-RAIN【福茂】【WAV+CUE】