%?“高速公路汽車陰影”視頻的陰影去除的第一種方法
clc;
clear?all;
close?all;
%?打開視頻文件對話框
h1=figure（‘toolbar‘‘none‘‘name‘‘Shadow?Removal‘）;??%?讀取受損圖像
[FileNamePathName]=uigetfile（??{?‘*.avi‘‘a(chǎn)vi?（*.avi）‘;‘*.*‘‘All?Files?（*.*）‘}‘Open?video?sequence‘??）;
if?isequal（[FileNamePathName][00]）
???return;
else
???pic_path=fullfile（PathNameFileName）;
???mov=aviread（‘highwayII_raw.avi‘）;
end

n=10;?????????????%?讀取視頻幀數(shù)
s=1;???????????????%?起始為1
e=n;???????????????%?結束為n

for?i=s:e??
????mov_gray（i）.image=mov（i）.cdata;???????%mov_gray（i）.image為第i幀圖像
end?????????????????????????????????????

[ihiwid]=size（mov_gray（1）.image）;???????%?取得視頻幀的大小

%?計算1到第n幀的各個像素的均值
for?i=s:e
????u=mov_gray（i）.image;
????h=u（:）;
????t（i:）=h（:）;
end
tmp=median（t）;????????????????????????????%?對t數(shù)組各列取均值?????????????????????????????????????
mean=reshape（tmpihiwid）;???????????????%?將n幀均值圖像恢復成h*w*d的圖像矩陣得到背景
figureimshow（uint8（mean））;title（‘背景圖像‘）;

fn=10;

ok=mean;
tmp（1:3）=0;
di=20;
th=ones（113）*di;??????????????????????????%?預設閥值
diff=double（ok）-double（mov_gray（fn）.image）;?%?計算第fn幀圖像與背景的差值
diff1=abs（diff）;
for?ai=1:ih
????for?bi=1:iw
????????if?abs（diff（aibi1:3））>=th（11:）
???????????ok（aibi1:3）=mov_gray（fn）.image（aibi1:3）;
????????else
???????????ok（aibi1:3）=0;
???????????diff1（aibi1:3）=0;
????????end
????end
end
figureimshow（uint8（ok））;title（‘前景圖像‘）;

impixelinfo;

g=rgb2hsv（mean）;?????????????????????????????%?將背景圖像由RGB格式轉(zhuǎn)為HSV格式
figureimshow（g）;title（‘背景圖像的hsv格式‘）;
impixelinfo;
figureimtool（g）;

g1=rgb2hsv（ok）;??????????????????????????????%?將前景圖像由RGB格式轉(zhuǎn)為HSV格式
figureimshow（g1）;title（‘前景圖像的hsv格式‘）;
impixelinfo;
figureimtool（g1）;

%?閥值分割
for?ci=1:ih
????for?di=1:iw
????????if??（g1（cidi3）./g（cidi3）>=0.4）...?????????????????????????????????????%這個地方的閾值選取是手動的，根據(jù)經(jīng)驗得到！
????????????&&（g1（cidi3）./g（cidi3）<=0.5）...
????????????&&abs（g1（cidi2）-g（cidi2））<=0.2?&&?abs（g1（cidi1）-g（cidi1））<=0.2;???
????????????h1（cidi1:3）=1;
????????else
????????????h1（cidi1:3）=0;
????????end
????end
end
h=hsv2rgb（h1）;
figureimshow（h）;title（‘陰影去除結果‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2328??2013-05-24?23:33??hsv_file.m