Python計算任意多邊形間的重疊面積的示例代碼
跟某人討論一個排樣問題。他說,算法搜索速度很慢,每兩個物體間的重疊面積計算時間若按1s來算,300個物體需要計算將近9萬次。我說,這用計算機視覺難道不是幾句話解決的嘛!(小小的嘚瑟一把,雖然做了這么久的CV,一直覺得自己一無所成,但是沒想到默默的就能解決別人的問題了哈哈哈~~)
本文檔目的為:給定的數(shù)據(jù)為多邊形的各個頂點,為N*2的矩陣,N 為多邊形的頂點個數(shù),計算任意兩個多邊形重疊面積計算的工具介紹及程序。注意,并不涉及IOU的計算(雖然只是一句話的事哈哈哈)等,只是要重疊面積。
1. shapely工具箱判斷任意兩個多邊形的面積是否有交集有函數(shù)A.intersects(B)實現(xiàn),若A和B有交集,返回為TRUE,若沒有交集,返回為FALSE;計算任意兩個多邊形的面積有函數(shù)A.intersection(B).area,輸出直接為A和B的交集的面積。
那么問題在于,如何將多邊形頂點的坐標(biāo),換為組成多邊形的所有內(nèi)部點的坐標(biāo)。python有個工具箱shapely。用于解決多邊形有關(guān)問題。其中有個多邊形填充函數(shù)Polygon.
shapely的安裝方法為:
pip install Shapely
但是不知道為什么,用該方法安裝時,一直報錯,后來直接從https://www.lfd.uci.edu/~gohlke/pythonlibs/#shapely (直接ctrl+F搜索shapely快速轉(zhuǎn)到下載位置)下載了該文件,用地址名安裝。
安裝成功后,即可用下面示意的程序計算面積。
2. 程序import numpy as npimport timefrom shapely.geometry import Polygon # 多邊形import scipy.io as iodef Cal_area_2poly(data1,data2): ''' 任意兩個圖形的相交面積的計算 :param data1: 當(dāng)前物體 :param data2: 待比較的物體 :return: 當(dāng)前物體與待比較的物體的面積交集 ''' poly1 = Polygon(data1).convex_hull # Polygon:多邊形對象 poly2 = Polygon(data2).convex_hull if not poly1.intersects(poly2):inter_area = 0 # 如果兩四邊形不相交 else:inter_area = poly1.intersection(poly2).area # 相交面積 return inter_area data1 = [] # 帶比較的第一個物體的頂點坐標(biāo)data2 = [] #待比較的第二個物體的頂點坐標(biāo)area = Cal_area_2poly(data1,data2)
到此這篇關(guān)于Python計算任意多邊形間的重疊面積的示例代碼的文章就介紹到這了,更多相關(guān)Python 計算重疊面積內(nèi)容請搜索好吧啦網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持好吧啦網(wǎng)!
相關(guān)文章:
1. asp讀取xml文件和記數(shù)2. asp批量添加修改刪除操作示例代碼3. 使用Spry輕松將XML數(shù)據(jù)顯示到HTML頁的方法4. ASP中解決“對象關(guān)閉時,不允許操作?!钡脑幃悊栴}……5. CSS可以做的幾個令你嘆為觀止的實例分享6. CSS Hack大全-教你如何區(qū)分出IE6-IE10、FireFox、Chrome、Opera7. 低版本IE正常運行HTML5+CSS3網(wǎng)站的3種解決方案8. CSS3中Transition屬性詳解以及示例分享9. html小技巧之td,div標(biāo)簽里內(nèi)容不換行10. 告別AJAX實現(xiàn)無刷新提交表單
