matlab处理视频代码,Matlab视频处理心脏跳动。代码补充

我正在努力编写一段代码，帮助我进行生物学工作。

代码的概念是分析在组织中收缩细胞的视频文件

Example 1

示例2：youtube.com/watch?v=uG_WOdGw6Rk

并绘出以下内容：

每分钟跳动次数

击败力量

跳动的规律

所以我写了一个Matlab代码，它可以循环播放一段视频，并比较每一帧和随后的帧，然后查看帧中是否有任何变化，并将这些变化绘制在曲线上。

我的代码结果示例

当前代码的核心我写道：

for i=2:totalframes

compared=read(vidObj,i);

ref=rgb2gray(compared);%% convert to gray

level=graythresh(ref);%% calculate threshold

compared=im2bw(compared,level);%% convert to binary

differ=sum(sum(imabsdiff(vid,compared))); %% get sum of difference between 2 frames

if (differ ~=0) && (any(amp==differ)==0) %%0 is = no change happened so i dont wana record that !

amp(end+1)=differ; % save difference to array amp wi

time(end+1)=i/framerate; %save to time array with sec's, used another array so i can filter both later.

vid=compared; %% save current frame as refrence to compare the next frame against.

end

end

figure,plot(amp,time);=====================

这就是我的代码，但有没有办法改进它，这样我可以得到更好的结果？

因为我担心imabsdiff并不完全是我应该使用的，因为我的视频包含很多噪音，并且会影响我的结果，而且我认为我所有的放大器数据实际上都是假的！

另外，我实际上只能通过计算峰值来提取跳动率，但是如何提高我的代码以便能够获取所有必需的数据？

也非常感谢你的帮助，这是代码的一小部分，如果你需要更多的信息，请让我知道。

谢谢