%%?Copyright?（c）?2007?Christopher?Wellons?
%%
%%?Permission?to?use?copy?modify?and?distribute?this?software?for?any
%%?purpose?with?or?without?fee?is?hereby?granted?provided?that?the
%%?above?copyright?notice?and?this?permission?notice?appear?in?all
%%?copies.
%%
%%?THE?SOFTWARE?IS?PROVIDED?“AS?IS“?AND?THE?AUTHOR?DISCLAIMS?ALL
%%?WARRANTIES?WITH?REGARD?TO?THIS?SOFTWARE?INCLUDING?ALL?IMPLIED
%%?WARRANTIES?OF?MERCHANTABILITY?AND?FITNESS.?IN?NO?EVENT?SHALL?THE
%%?AUTHOR?BE?LIABLE?FOR?ANY?SPECIAL?DIRECT?INDIRECT?OR?CONSEQUENTIAL
%%?DAMAGES?OR?ANY?DAMAGES?WHATSOEVER?RESULTING?FROM?LOSS?OF?USE?DATA?OR
%%?PROFITS?WHETHER?IN?AN?ACTION?OF?CONTRACT?NEGLIGENCE?OR?OTHER
%%?TORTIOUS?ACTION?ARISING?OUT?OF?OR?IN?CONNECTION?WITH?THE?USE?OR
%%?PERFORMANCE?OF?THIS?SOFTWARE.

function?project3（folder?alpha）
%?Chris?Wellons
%?Ron?Ledwich
%?Ben?Hvostal
%?CSE?486?Project?3

%?Set?number?of?frames?（images）?based?on?folder?entered
frames?=?length（dir（fullfile（folder?‘f*jpg‘）））;

%?Set?alpha
if?~exist（‘alpha‘?‘var‘）
	alpha?=?.25;
end

%?object?tracking?colors
color_list?=?[];
color_list?=?[color_list?{[1?0?0]}];?%?Red
color_list?=?[color_list?{[0?1?0]}];?%?Green
color_list?=?[color_list?{[0?0?1]}];?%?Blue
color_list?=?[color_list?{[0?1?1]}];?%?Cyan
color_list?=?[color_list?{[1?0?1]}];?%?Fuchsia?（sp?）
color_list?=?[color_list?{[1?1?1]}];?%?White
color_list?=?[color_list?{[1?.5?0]}];?%?Orange
color_list?=?[color_list?{[.914?.5882?.4784]}];?%?Dark?Salmon
color_list?=?[color_list?{[0?1?.5]}];?%?Spring?green
color_list?=?[color_list?{[.42?.557?.137]}];?%?Olive?drab
image_list?=?[{}?{}?{}?{}?{}];

%?Initialize?some?variables?and?stuff
tic;?figure;
boxes?=?[];
SE?=?strel（‘disk‘?3）;
outm?=?[];
impad?=?5;?%?Padding?around?ncc?box
ncc_thresh?=?0.8;?%?Threshold?for?good?ncc?match
timeave?=?0;

%?Adaptive?Background?Subtraction
for?i=1:frames
????
????%?Figure?out?time?left
????timeave?=?（toc+6*timeave）/7;
????timeleft?=?timeave?*?（frames?-?i）;
????minutes?=?floor（timeleft/60）;
????seconds?=?floor（mod（timeleft60））;
????tic;
	disp（[‘Working?on?frame?‘?num2str（i）?‘?（‘?num2str（minutes）?‘:‘?num2str（seconds）?‘?left）‘]）;
	[I?cI]?=?getimage（folder?i）;

	%?Initialize?B
	if?~exist（‘B‘?‘var‘）
		B?=?I;
	end

	%?Get?difference?data
	diff?=?abs（B?-?I）?*?255;
	M?=?threshold（diff?15）;
	
	%?Prepare?next?B
	B?=?alpha*I?+?（1-alpha）*B;
	
	%?Initalize?current?frame?box?data
	curbox?=?[];
	ncc_pass?=?0;
	ncc_fail?=?0;
	
	%?Dilate?the?difference?image
	dim?=?imdilate（M?SE）;
	
	%?Find?each?object
	for?j?=?1:size（dim?1）
		for?k?=?1:size（dim?2）
			if?dim（jk）?==?1
				doit?=?0;
				
				%?Trace?out?the?object
				T?=?bwtraceboundary（dim?[jk]?‘E‘）;
				
				%?Filter?small?objects
				if?length（T）?>?100
					doit?=?1;
				end
				
				%?Convert?trace?into?a?blob
				T?=?full（spconvert（[T?ones（length（T）?1）]））;
				[x1?y1?x2?y2]?=?find_bounds（T）;
				T（y1:y2x1:x2）?=?ones（y2-y1+1x2-x1+1）;
				
				%?Find?original?parts?and?draw?a?box
				if?doit
					C?=?and（T