logo资料库

Python计算不规则图形面积算法实现解析.pdf

第1页 / 共6页
第2页 / 共6页
第3页 / 共6页
第4页 / 共6页
第5页 / 共6页
第6页 / 共6页
资料共6页,全文预览结束
Python计算不规则图形面积算法实现解析 计算不规则图形面积算法实现解析 主要介绍了Python计算不规则图形面积算法实现解析,文中通过示例代码介绍的非常详细,对大家的学习或者工 作具有一定的参考学习价值,需要的朋友可以参考下 这篇文章主要介绍了Python计算不规则图形面积算法实现解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具 有一定的参考学习价值,需要的朋友可以参考下 介绍:大三上做一个医学影像识别的项目,医生在原图上用红笔标记病灶点,通过记录红色的坐标位置可以得到病灶点的外接 矩形,但是后续会涉及到红圈内的面积在外接矩形下的占比问题,有些外接矩形内有多个红色标记,在使用网上的opencv的 fillPoly填充效果非常不理想,还有类似python计算任意多边形方法也不理想的情况下,自己探索出的一种效果还不错的计算多 圈及不规则图形的面积的算法。 能较为准确的计算出不规则图形的面积 能较为准确的计算出不规则图形的面积 正文:算法的思想很简单,遍历图片每一列,通过色差判断是否遇到标记圈,将坐标全部记录,对每一列的坐标都进行最小行 和最大行记录,确定每一列的最小和最大的坐标,然后上色(类似opencv的fillPoly的实现,但是细节有些区别),只是这样 效果并不好,将图片旋转90度,再做一边,将两个图片的结果放在一起做与操作,得到结果就能很好的处理多圈的标记问题 和多算面积的问题(比如上面的08-LM), 算法实现 算法实现 全程只用pillow库 首先先用屏幕拾色器获取目标颜色的rgb值,我这种情况下就是(237,28,36),前期截取外接矩形也是要这一步的,颜色 也一致
def pixel_wanted(pix): return pix==(237,28, 36) 每一列都设定翻转位初始为False,如果上一个像素点不是目标色,当前是目标色则开始记录,一旦不是目标色,停止检测 top_Pixel都设定为黑色(0,0,0)因为有图片最上方就是目标色,导致判定出问题,直接让最上面的像素初始化是黑色 coordinate_List记录了所有符合的点坐标 coordinate_List = [] top_Pixel = (0,0,0) for x in range(im.size[0]): flag = False #初始化每一列翻转位为False for y in range(im.size[1]): current_pixel = im.getpixel((x,y)) last_pixel = im.getpixel((x,y-1)) if y>0 else top_Pixel #翻转判定 if pixel_wanted(current_pixel) and \ not pixel_wanted(last_pixel): flag = True if flag and not pixel_wanted(current_pixel): flag = False if(flag): coordinate_List.append((x,y)) coordinate_List中的点如下图 然后就是将上面获得coordinate列表进行处理 将coordinate列表中每一列的最小坐标和最大坐标进行记录 因为每一列记录的数量并不确定(应该可以在上一步改进一下),所以需要遍历多次 首先找到第一个列出现的坐标,将它的行信息记录(行信息最小确定), 然后遍历出全部的同列的坐标,比较行坐标,如果大的就将最大的代替(行信息最大确定),用一个新的列表记录数据 coordinate_Min_Max_List = [] #找最小最大 for i in range(im.size[0]): min=-1 max=-1 for coordinate in coordinate_List: if coordinate[0] == i: min = coordinate[1] max = coordinate[1] break for coordinate in coordinate_List: if coordinate[0] == i: if coordinate[1]>max: max = coordinate[1] coordinate_Min_Max_List.append(min) coordinate_Min_Max_List.append(max) 其中要将min和max都初始化为一个坐标不存在的值比如-1,为了在下一步多圈且有空隙情况下,不会出现残影现象,如下图
上一步的最后得到一个列表,第n列的最小行和最大行分别是第2n和2n+1元素,结果中的-1,为了让下一步不会画进去 然后就是绘制图片了,每一列将列表中对应的最小行到最大行涂满 #上色 for x in range(im.size[0]): for y in range(im.size[1]): min = coordinate_Min_Max_List[x*2] max = coordinate_Min_Max_List[x*2+1] if min
可以看到效果非常不错,但是依旧有个别图像有问题,比如十字分布的, 但现在的话误差已经降低非常多了,这些极其个别的十字现象可以手动把原图切割一下,或者干脆不处理了
所有代码,画出绿图片为了方便直观的查看,函数中可以把图片顺便保存一下,总体看一下效果 from PIL import Image def pixel_wanted(pix): return pix==(237,28, 36) def fillPoly(im): coordinate_List = [] top_Pixel = (0,0,0) for x in range(im.size[0]): flag = False #初始化每一列翻转位为False for y in range(im.size[1]): current_pixel = im.getpixel((x,y)) last_pixel = im.getpixel((x,y-1)) if y>0 else top_Pixel #翻转判定 if pixel_wanted(current_pixel) and \ not pixel_wanted(last_pixel): flag = True if flag and not pixel_wanted(current_pixel): flag = False if(flag): coordinate_List.append((x,y)) coordinate_Min_Max_List = [] #找最小最大 for i in range(im.size[0]): min=-1 max=-1 for coordinate in coordinate_List: if coordinate[0] == i: min = coordinate[1] max = coordinate[1] break for coordinate in coordinate_List: if coordinate[0] == i: if coordinate[1]>max: max = coordinate[1] coordinate_Min_Max_List.append(min) coordinate_Min_Max_List.append(max)
#上色 for x in range(im.size[0]): for y in range(im.size[1]): min = coordinate_Min_Max_List[x*2] max = coordinate_Min_Max_List[x*2+1] if min
分享到:
收藏