tic;
clear?all;
close?all;
clc;
A?=?aviinfo（‘smoke.avi‘）;
LEN?=?A.Numframes;%得到視頻的幀數
for?x?=?1?:?LEN
pic?=?imread（strcat（‘smoke‘?int2str（x）?‘.jpg‘））;%讀入一張圖片
pic1?=?pic;%pic1中存放的是處理過后的彩色圖像
pic_hsv?=?rgb2hsv（pic）;
s?=?size（pic1）;
bw_pic?=?zeros（s（1）?s（2））;
for?i?=?1?:?s（1）
????for?j?=?1?:?s（2）
????????if?pic_hsv（i?j?1）?>?0.35?&&?pic_hsv（i?j?1）????????????pic_hsv（i?j?2）?>?0.06?&&?pic_hsv（i?j?2）????????????pic_hsv（i?j?3）?>?0.5?&&?pic_hsv（i?j?3）?<=?1;
???????????pic1（i?j?1）?=?255;
???????????pic1（i?j?2?:?3）?=?0;
????????end
????????if?pic1（i?j?1）?==?255?&&?pic1（i?j?2）?==?0
????????????bw_pic（i?j）?=?1;
????????end
????end
end
cov_bw_pic?=?1?-?bw_pic;
fill_pic?=?imfill（cov_bw_pic?‘holes‘）;
fill_pic1?=?1?-?fill_pic;
B_pic?=?bwboundaries（fill_pic1）;%得到的是一個P×1的單元
ss?=?size（B_pic）;%ss的值就是連通分量的數量
ss1?=?size（B_pic{1?1}）;%ss1是一個大小為m×2的矩陣
max?=?ss1（1）;
for?i?=?2?:?ss（1）%找到最大的連通區域的邊界總個數
????ss2?=?size（B_pic{i?1}）;
????if?ss2（1）>?max
????????max?=?ss2（1）;
????end
end
%⑤…………………………將小于最大連通區域的連通分量邊界賦零
for?k?=?1?:?ss（1）%將小于最大連通區域的連通分量邊界賦零，
%其作用就是將其余小的連通分量與邊界隔開
????Mid?=?B_pic{k?1};
????ss3?=?size（Mid）;
????if?ss3（1）?????????for?i?=?1?:?ss3（1）
????????????fill_pic1（Mid（i?1）?Mid（i?2））?=?0;??
????????end
????end
end
%反轉填充
cov_fill_pic1?=?1?-?fill_pic1;
A_fill_pic?=?imfill（cov_fill_pic1?‘holes‘）;%進行內部填充，再反變換回來，即得到所要求的區域
A_fill_pic1?=?1?-?A_fill_pic;
B_pic1?=?bwboundaries（A_fill_pic1）;
sss?=?size（B_pic1）;%sss的值就是連通分量的數量
for?k?=?1?:?sss（1）%該循環作用是找到最大區域的邊界，然后在原圖像的對應位置上賦紅
????Mid?=?B_pic1{k};%Mid是一個m×2的矩陣，m就是邊界點的數量
????%ss4（1）表示邊界點的數量，Mid（i?1）和Mid（i?2）表示邊界點的坐標
????ss4?=?size（Mid）;
????if?ss4（1）?==?max
????????for?i?=?1?:?ss4（1）
????????????pic（Mid（i?1）?Mid（i?2）?1）?=?255;??
????????????pic（Mid（i?1）?Mid（i?2）?2?:?3）?=?0;??
????????end
????end
end
imwrite（pic?strcat（‘lou_smoking‘?int2str（x）?‘.jpg‘）?‘jpg‘）;
end
aviobj?=?avifile（‘lou_smoking.avi‘）;%將處理后的圖片連成視頻文件
aviobj.fps?=?25;%初始化視頻文件的屬性
aviobj.Quality?=?100;
aviobj.compression?=?‘None‘;
for?i?=?1?:?LEN
????frame?=?imread（strcat（‘lou_smoking‘?int2str（i）?‘.jpg‘））;?
????aviobj?=?addframe（aviobj?uint8（frame））;?
end
aviobj=close（aviobj）;%關閉aviobj的對象
toc;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2534??2013-05-21?14:49??圖像處理中的煙霧檢測\smoke_inspect_two.m

?????目錄??????????0??2013-05-21?14:50??圖像處理中的煙霧檢測

-----------?---------??----------?-----??----

?????????????????2534????????????????????2