function [feature,feature_value1_eq]=Algor1_new()
%UNTITLED4 此处显示有关此函数的摘要
% 此处显示详细说明
%视频为1秒25帧
video= VideoReader('small2.avi');
tic
vFrames = video.NumberOfFrames; %视频总帧数
vHeight = video.Height;
vWidth = video.Width;
vDuration=video.Duration; %视频时长
frame_per_sec=video.FrameRate; %每秒帧数
%1秒间隔边缘重叠超过多少有预警
% frame_per_sec=fix(vFrames/vDuration); %每秒帧数
% feature=zeros(vHeight,vWidth,ceil(vFrames/frame_per_sec));
thresD=1000; %需先知道没车的时候边缘值为1的点是多少
time_interval=1; %用于分析的采集图像时间间隔
frame_interval=max(fix(time_interval*vFrames/vDuration),1); %向下取整，避免为0，硬性设最小为1
time_begin=0; %分析
time_end=vDuration;
stat_frames=fix(((time_end-time_begin)*frame_per_sec)/frame_interval); %统计的视频帧数
basic_edge=read(video,10);
basic_edge=rgb2gray(basic_edge);
basic_edge=edge(basic_edge,'sobel');
% for i=1:25:vFrames,
% for i=1:1:vFrames/frame_per_sec,
% for i=1:1:24*25/15,
for i=1:1:stat_frames,
fprintf('processing Frame = %d,Video time=%f, processing time=%f\n',(i-1)*frame_interval+1,(i-1)*time_interval,toc);
im=read(video,(i-1)*frame_interval+1);
im=rgb2gray(im);
% feature(:,:,i)=edge(im,'sobel');
temp_feature=edge(im,'sobel');
temp_feature=temp_feature-basic_edge; %帧边缘减底图边缘
% figure()
% imshow(temp_feature)
% feature(:,:,i)=temp_feature(170:end-70,:); %small视频去掉上下有水印的重复部分（边缘特征重复多）
feature(:,:,i)=temp_feature(60:end-70,:); %big视频去掉上下有水印的重复部分（边缘特征重复多）
if i>=2,
feature_value1_1=find(feature(:,:,i-1)==1);
feature_value1_2=find(feature(:,:,i)==1);
eq_num=intersect(feature_value1_1,feature_value1_2);
eq_num=length(eq_num);
% if isempty(eq_num),
% eq_num=0;
% end
feature_value1_eq(i)=eq_num;
end
end
end
