中国航天空间站最新消息,灰色行业推广平台网站,wordpress 做的网站,WordPress去掉文章摘要一张图片中两个任意多边形相交的面积计算方法。本文参考https://blog.csdn.net/PanYHHH/article/details/110940428#xff1b;加了一个简单的示例#xff0c;也对代码做了一点清淅化。原博客中还有其他链接#xff0c;是C代码#xff0c;没有看原理#xff0c;但以下代码…一张图片中两个任意多边形相交的面积计算方法。本文参考https://blog.csdn.net/PanYHHH/article/details/110940428加了一个简单的示例也对代码做了一点清淅化。原博客中还有其他链接是C代码没有看原理但以下代码也可以很容易转成C代码用opencv 来实现。
import cv2
import numpy as np
ImH200
ImW200
Polygon1np.array([[0,0],[100,0],[100,100],[0,100]],dtypenp.int32)
Polygon2np.array([[50,50],[150,50],[150,150],[50,150]],dtypenp.int32)
ImShape(ImH,ImW,3)
def DrawPolygon(ImShape,Polygon,Color):Im np.zeros(ImShape, np.uint8) cv2.fillPoly(Im, [Polygon], Color) # 这个函数可以画凹凸多边形所以这个更稳妥#cv2.fillConvexPoly(Im,Polygon,Color) #这行代码也是可以的因为是凸多边形return Im
Im1 DrawPolygon(ImShape, Polygon1, (255, 0, 0))
Im2 DrawPolygon(ImShape, Polygon2, (0, 255, 0))
plt.subplot(311)
plt.imshow(Im1)
plt.subplot(312)
plt.imshow(Im2)Im1 DrawPolygon(ImShape[:-1],Polygon1,122)#多边形1区域填充为122
Im2 DrawPolygon(ImShape[:-1], Polygon2, 133)#多边形2区域填充为133
Im Im1 Im2
ret, OverlapIm cv2.threshold(Im, 200, 255, cv2.THRESH_BINARY)#根据上面的填充值因此新图像中的像素值为255就为重叠地方
plt.subplot(313)
plt.imshow(OverlapIm,cmapgray)
IntersectAreanp.sum(np.greater(OverlapIm, 0))#求取两个多边形交叠区域面积
print(cumstom calcuate area:{}\n.format(IntersectArea))#下面使用opencv自带的函数求取一下最为对比
contours, hierarchy cv2.findContours(OverlapIm,cv2.RETR_TREE,cv2.CHAIN_APPROX_SIMPLE)
contourAreacv2.contourArea(contours[0])
print(contourArea{}\n.format(contourArea))
perimeter cv2.arcLength(contours[0], True)
print(contourPerimeter{}\n.format(perimeter))
RealContourAreacontourAreaperimeter
print(RealContourArea{}\n.format(RealContourArea))结果如下图
