clear all
rmin =0; %min row value for search window
rmax = 0; %max row value for search window
cmin =0; %min col value for search window
cmax = 0; %max col value for search window
numofframes = 0; %number of frames in the avi
threshold = 1; %threshold for convergence
centerold = [0 0]; %for convergence... previous center of window
centernew = [0 0]; %for convergence... new center of window
% ******************************************************************
% Pre code... load movie and select initial frame
% ******************************************************************
% prompt user for avi file name
% user_entry = input('Please enter an avi filename: ','s');
% load the avi file... handle is M
% M =mmreader(user_entry);
M=VideoReader('E:\tp\car1.avi','tag','myreader1');
% get number of frames
% [dontneed numberofframes] = size(M);
vidFrames=read(M);
numberofframes=get(M,'numberofFrames');
% initialize to hold center coordinates
imagecenters = zeros(numberofframes, 2);
% extract the first frame from the avi
% Frame1 = getframe(1,1);
% Frame1 =M(1,1);
% Image1= frame2im(Frame1);
Image1=vidFrames(:,:,:,45);
% Image1=read(M,20);
%%% ********** images(:, :, numberofframes) = G(:,:);
% get search window for first frame
[ cmin, cmax, rmin, rmax ] = select( Image1 );
cmin = round(cmin);
cmax = round(cmax);
rmin = round(rmin);
rmax = round(rmax);
% 计算目标尺寸wsize:长与高
wsize(1) = abs(rmax - rmin);
wsize(2) = abs(cmax - cmin);
% create histogram
% translate to hsv
hsvimage = rgb2hsv(Image1);
% pull out the h
huenorm = hsvimage(:,:,1);
% scale to 0 to 255
hue = huenorm*255;
% set unit type
hue=uint8(hue);
% Getting Histogram of Image:
histogram = zeros(256);
% 计算目标矩形内的色调直方图histogram