图形图像通信与处理技术在林火监测中的应用 第18页
．7．
e(b+l，C，：)=a(b+l，c，：)；
end
end
figure；imshow(q)；
figure；imshow(w)；
figure；imshow(e)；
3、db4小波分解后高频系数采用文纳(自适应)滤波去噪程序：
Y--imread(’d：、3．bmp‘)；
[X,mapl=rgb2ind(Y,256)；
SrIVC。It'Xmap；
figure；image(X)；eolormap(map)；
axis square；
％loadIt
％因为图像不光滑，因此必须索引图像进一步处理
Y--double(X)；
R--m；~p(Y+l，1)；R=rcshape(R,size(Y))；
G=map(Y+I，2)；‰h印e(G’size(Y))；
B=map(Y+l，3)；B--reshape(B，size(Y))；
Xgray=O．2990‘R-卜o．5870‘G加．1 140‘B；
n=256；
,~-,ound(XFay‘(n·1))+l；
map2=gray(n)；
figure；image(z)；colormap(map2)；
axis square；
tptxl=z；
map---map2；
SaVCtptxl tptxl map；
【c,s]---wavedec2(tptx l，2,'db4’)；
siz=s(size(s，1)，：)；
ca l=appcoef2(c，S，’db4',l，siz)；
ehd 1=deteoef2(W,C，s，1)；
cvdl---deteoef2(V，c，S，1)；
edd l=deteoef2(’d’，c，S，1)；
figure；
subplot(22 1)；image(ca 1)；colormap(map)；
subplot(222)；image(chd 1)；colormap(map)；
subplot(223)；image(cvd 1)；colormap(map)；
subplot(224)；image(edd 1)；colormap(map)；
％h=【l，9 1／9 1／9；1／9 1／9 1／9；1／9 1／9 1／91；
％cal-~lter2(h,cal)；
chdl---wiener2(chdl，【3 4】)；
cvdl=wiener2(cvdl，【3 4】)；
cddl---wiener2(eddl，【3 4】)；
figure；
subplot(221)；image(cal)；colormap(map)；
subplot(222)；image(chd 1)；colormap(map)；
subplot(223)；image(cvdl)；colormap(map)；
subplot(224)；image(cdd 1)；colormap(map)；
figure；
dd=idwt2(cal，chdl，cvdl，cddl,'db4’，siz)；
image(dd)；coiormap(m印)；
4、最大值灰度法后采用阈值化提取火焰区域程序：
i--imread(’d：＼4．bmp’)；imshow(i)；
info--imfmfo(’d：＼4．bmp’)； ’
width=info．Width；height=info．Height；
form=1：height ，
forn=l：width
j(m,n)--max(i(m,n,1)，i(m'Il'2))；
j(m,n)--max(j(m,n)，i(m'II，3))；
end
end
figure；imshow(j)；
level=graythrosh(j)；％[~动获取阈值
bw=im2bw(j，level)；％24位真彩色转换为2值图像(黑白)
figure；imshow(bw)；
％对图像取反
form=1：height
forn=l：width
if bw(m，n)==O
bw(m，n)=l；
else
bw(m'11)=O；
end
end
end
figure；imshow(bw)； 。
％提取火焰区域代码
for m=l：height
forn=l：width
if bw(m，n)一l
i(m，^：)=255；
end
end
end
figure；imshow(i)；
<< 上一页 [11] [12] [13] [14] [15] [16] [17] [18]
图形图像通信与处理技术在林火监测中的应用 第18页下载如图片无法显示或论文不完整，请联系qq752018766
