python - 如何实时平滑阈值处理？

Question

首先我想用相机找到水滴入织物的时间变化。算法检测到运动后，用户将水滴到织物上，直到水完全吸收并绘制图表变化时间并显示吸收时间，面积等。

为了检测运动，我使用了具有恒定变化率的 absdif 函数。我将检测帧的开始时间计算到结束 like this image .这里没有问题。但是为了计算水的吸收率，我对帧进行了阈值处理，并使用 countNonZero 函数来计算黑色像素的数量。但是这里有一个问题，黑色像素显示红线thresholded images不断变化(如晃动、振动等)。因此绘图过程失败。

尝试

1. 我尝试更改网络摄像头设备(通过 İpcam 使用手机摄像头)
2. 我尝试使用自适应阈值方法(otsu 等)来找到最佳阈值
3. 平滑闪电条件并在没有背景的情况下捕获

成功

1. 当我使用手机摄像头拍摄的视频作为输入时，抖动和振动效果会降低，我可以取得成功this graph正如预期的那样

问题

• 如何实时平滑阈值图像
• 另一种方法

代码

import cv2
from datetime import datetime
import numpy as np
import matplotlib.pyplot as plt
import operator

def pixelHesaplayici(x):
    siyaholmayanpixel=cv2.countNonZero(x)
    height,width=x.shape
    toplampixel=height*width
    siyahpixelsayisi=toplampixel-siyaholmayanpixel
    return siyahpixelsayisi

def grafikciz(sure,newblackpixlist,maxValue,index,totaltime,cm):
    plt.figure(figsize=(15,15))
    plt.plot(sure,newblackpixlist)
    line,=plt.plot(sure,newblackpixlist)
    plt.setp(line,color='r')
    plt.text(totaltime/2,maxValue/2, r'$Max- 
    Pixel=%d$'%maxValue,fontsize=14,color='r')
    plt.text(totaltime/2,maxValue/2.5, r'$Max-emilim- 
    zamanı=%f$'%sure[index],fontsize=14,color='b')
    plt.text(totaltime/2,maxValue/3, r'$Max- 
    Alan=%fcm^2$'%cm,fontsize=14,color='g')
    plt.ylabel('Black Pixels')
    plt.xlabel('Time(s)')
    plt.grid(True)
    plt.show()


static_back=None
i=0
blackpixlist=[]
newblackpixlist=[]

t=[]
video=cv2.VideoCapture("kumas1.mp4")

while(True):
    ret,frame=video.read()

    if ret==True:
        gray=cv2.cvtColor(frame,cv2.COLOR_BGR2GRAY)
        gray=cv2.GaussianBlur(gray,(5,5),0)
        _,threshforgraph=cv2.threshold(gray,0,255,
       cv2.THRESH_BINARY+cv2.THRESH_OTSU)
        if static_back is None:
            static_back=gray
            continue
        diff_frame=cv2.absdiff(static_back,gray)

        threshfortime=cv2.threshold(diff_frame,127,255,cv2.THRESH_BINARY)[1]
        #threshfortime=cv2.dilate(threshfortime,None,iterations=2)
        (_,cnts,_)=cv2.findContours(threshfortime.copy(),
                               cv2.RETR_EXTERNAL,cv2.CHAIN_APPROX_SIMPLE)

        for contour in cnts:
            if cv2.contourArea(contour)<450:
                continue
            an=datetime.now()
            t.append(an.minute*60+an.second+(an.microsecond/1000000))
            cv2.fillPoly(frame,contour, (255,255,255), 8,0)
            cv2.imwrite("samples/frame%d.jpg"%i,threshforgraph)


            i+=1


        cv2.imshow("org2",frame)
        #cv2.imshow("Difference Frame",diff_frame)
        #cv2.imshow("Threshold Frame",threshfortime)
        #cv2.imshow("Threshforgraph",threshforgraph)

        if cv2.waitKey(1) & 0xFF == ord('q'):
            break
    else:
        break
ti=t[1::3]

lasttime=ti[-1]
firsttime=ti[-len(ti)]
totaltime=lasttime-firsttime

for i in range(0,i):
        img=cv2.imread('samples/frame%d.jpg'%i,0)
        blackpixlist.append(pixelHesaplayici(img))
ilkpix=blackpixlist[0]

for a in blackpixlist:
    newblackpixlist.append(a-ilkpix)
newblackpixlisti=newblackpixlist[1::3]  
index , maxValue=max(enumerate(newblackpixlisti),
key=operator.itemgetter(1))
sure=np.linspace(0,totaltime,len(newblackpixlisti))
cm=0.0007*maxValue # For 96 dpi

grafikciz(sure,newblackpixlisti,maxValue,index,totaltime,cm)

Answer 1

从下一帧中减去第一帧怎么样？如果你可以知道或可以检测到什么时候没有掉落并减去它，那么差值只会给你掉落的结果。

如果您在不同的地方有多个掉落并且想要丢弃前一个掉落，这种方法也可能很有趣。请注意，您可以在阈值化之前和之后进行减法。我会在设置阈值之前推荐。

如果你知道你在你的过程中有很多晃动，你可能需要应用数字稳定，在这种情况下我建议看看这个教程: https://www.learnopencv.com/video-stabilization-using-point-feature-matching-in-opencv/

当然，稳定应该在减法之前完成。

一般来说，对于您的问题，我不会使用自适应方法。所有帧的阈值都应该相同，如果它根据图像进行调整，您可能会得到无效的结果。

希望我正确理解了您的问题!