我正在尝试编写代码,这有助于我的生物学工作。代码的概念是分析组织中收缩细胞的视频文件
示例 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 不完全是我应该使用的,因为我的视频包含很多噪音,这会影响我的结果,而且我认为我所有的 amp 数据实际上都是伪造的!
此外,我实际上只能通过计算峰值来从中提取跳动率,但是我如何改进我的代码以便能够从中获取所有必需的数据?
也非常感谢您的帮助,这是一小部分代码,如果您需要更多信息,请告诉我。谢谢