python基于opencv实现数籽粒

千粒重是一个重要的农艺性状,通过对其的测量和研究,我们可以更好地理解作物的生长状况,优化农业生产,提高作物产量和品质。但数籽粒数目是一个很繁琐和痛苦的过程,我们现在用一个简单的python程序来数水稻籽粒。代码的大致流程如下:

  1. 对照片进行预处理,包括灰度化、二值化等操作,以便更好地识别籽粒。
  2. 使用轮廓检测算法来找到照片中的籽粒,并计算轮廓面积;
  3. 对轮廓面积进行判断,以便区分重叠的籽粒。
  4. 计算籽粒数目,将结果输出到屏幕上或保存到文件中。

具体代码如下:

1. 加载需要的包

python 复制代码
import cv2
from matplotlib import pyplot as plt
import numpy as np
from PIL import Image, ImageDraw, ImageFont

2. 导入图片并展示

python 复制代码
img = cv2.imread('pic2.jpg',0)
img_rgb = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)
# 显示图像
plt.imshow(img_rgb)
plt.show()

3. 轮廓检测算法来找到照片中的籽粒

python 复制代码
ret, thresh = cv2.threshold(img, 180, 255, cv2.THRESH_BINARY)
contours, hierarchy = cv2.findContours(thresh, cv2.RETR_TREE, cv2.CHAIN_APPROX_SIMPLE)
area = np.zeros(len(contours))
for i in np.arange(0, len(contours)):
    area[i] = cv2.contourArea(contours[i])

4. 对轮廓面积进行判断,以便区分重叠的籽粒

python 复制代码
med = np.median(area)
j = 0
i = 0
a = area/med
draw = Image.fromarray(img)
for c in contours:
    if round(a[j]) == 1:   
        i = i+1
        draw1 = ImageDraw.Draw(draw)
        font = ImageFont.truetype('arial.ttf', 30)
        draw1.text(tuple(c[0][0]), str(i), fill = 'black', font = font)
    if round(a[j]) > 1:   
        i = i+round(a[j])
        draw1 = ImageDraw.Draw(draw)
        font = ImageFont.truetype('arial.ttf', 30)
        draw1.text(tuple(c[0][0]), str(i+1)+"~"+str(i+round(a[j])), fill = 'black', font = font)
    j = j + 1
print("籽粒数目:", i)

5. 展示标记情况

python 复制代码
img_rgb = cv2.cvtColor(np.array(draw), cv2.COLOR_BGR2RGB)
# 显示图像
plt.imshow(img_rgb)
plt.savefig("rice.pdf")
plt.show()
相关推荐
qq_12498707533 分钟前
基于Flask的穷游网酒店数据分析系统(源码+论文+部署+安装)
后端·python·flask·毕业设计
夜晚中的人海7 分钟前
【C++】位运算算法习题
开发语言·c++·算法
Brianna Home11 分钟前
PyTorch实战:CV模型搭建全指南
人工智能·pytorch·经验分享·python·神经网络·结对编程
裸奔在上海16 分钟前
使用Java做URL短连接还原长链接获取参数
java·开发语言·程序人生·spring
爱吃烤鸡翅的酸菜鱼16 分钟前
如何掌握【Java】 IO/NIO设计模式?工厂/适配器/装饰器/观察者模式全解析
java·开发语言·后端·设计模式·nio
喜欢吃豆25 分钟前
spec-kit深度解析:AI驱动的规范驱动开发(SDD)的架构、哲学与实践启示
人工智能·驱动开发·python·架构·大模型
rengang6627 分钟前
28-Keras:探讨高层神经网络API及其快速原型设计能力
人工智能·python·神经网络·机器学习·keras
椰羊sqrt27 分钟前
MetaTwo靶机实战:SQL注入到权限提升全解析
python·学习·网络安全
程序员大雄学编程34 分钟前
「用Python来学微积分」17. 导数与导函数
开发语言·python·数学·微积分
扶尔魔ocy36 分钟前
【QT常用技术讲解】可拖拽文件的Widget--QListWidget
开发语言·qt