1、掌握圖像加減運(yùn)算的運(yùn)用。

2、掌握分段線性點(diǎn)運(yùn)算的運(yùn)用。

實(shí)驗(yàn)內(nèi)容及要求：

1、把一幅圖像加上高斯噪聲，再通過100次相加求平均的方法去除噪聲。

參考代碼如下：（找到實(shí)驗(yàn)圖像并完成實(shí)驗(yàn)）

2、已知一幅受“椒鹽”噪聲干擾的圖像，通過減法運(yùn)算提取出噪聲。

3、將胸腔的斷層掃描圖像灰度級(jí)開窗，分別顯示骨骼、肌肉、肺部。

參考代碼如下：

（1）請補(bǔ)充代碼，讀入實(shí)驗(yàn)圖片并轉(zhuǎn)換成灰度圖像。

（2）給出實(shí)驗(yàn)結(jié)果。

（3）列出分段函數(shù)式，并結(jié)合實(shí)驗(yàn)結(jié)果畫圖分析。

實(shí)驗(yàn)過程（記錄）：

I=imread('yb.tif');

J=imnoise(I,'gaussian',0,0.02);

subplot(1,2,1),imshow(I);title('原始圖像');

subplot(1,2,2),imshow(J);title('加噪聲的圖像');

K=zeros(135,152);

for?i=1:100

????J=imnoise(I,'gaussian',0,0.02);

????J1=im2double(J);

????K=K+J1;

end

K=K/100;

figure;imshow(K);title('求平均后的圖像');

I=imread('lenna1.jpg');

J=imread('lenna.jpg');

K=imsubtract(I,J);

K1=255-K;

figure;imshow(I);title('有噪聲的圖像');

figure;imshow(J);title('原始圖像');

figure;imshow(K1);title('提取的圖像');

f1= imread('xq.jpg');

f = rgb2gray(f1);

subplot(2,2,1),imshow(f),title('(a)Original Image');

f=double(f);

[m,n]=size(f);

a=min(f(:));

b=max(f(:));

for?i=1:m

?????for?j=1:n

?????????if(f(i,j)<210)

?????????????g(i,j)=0;

?????????else

?????????????g(i,j)=255/(255-210)*(f(i,j)-210)+0;

?????????end

?????end

end

?????g=uint8(g);

?????subplot(2,2,2),imshow(g),title('(b)Bone');

?????clear g;

??????for?i=1:m

?????????for?j=1:n

?????????????if(f(i,j)<160)

?????????????????g(i,j)=0;

?????????????else?if(f(i,j)>210)

?????????????????????g(i,j)=0;

?????????????????else

?????????????????????g(i,j)=255/(210-160)*(f(i,j)-160)+0;

?????????????????end

?????????????end

?????????end

??????end

??????g=uint8(g);

?????subplot(2,2,3),imshow(g),title('(c)Muscle');

?????clear g;

?????for?i=1:m

?????????for?j=1:n

?????????????if(f(i,j)<45)

?????????????????g(i,j)=0;

?????????????else?if(f(i,j)>160)

?????????????????????g(i,j)=0;

?????????????????else

?????????????????????g(i,j)=255/(160-45)*(f(i,j)-45)+0;

?????????????????end

?????????????end

?????????end

?????end

?????g=uint8(g);

?????subplot(2,2,4),imshow(g),title('(d)Lung');

imnoise 函數(shù)說明

　　在MATLAB中，可以使用imnoise函數(shù)給一幅圖像添加噪聲。

　　例子：g = imnoise(f, type, parameters)

　　參數(shù)說明：

　　輸出：g是添加噪聲之后的圖像。

　　輸入：f是原圖像，type是加入的噪聲類型，parameters是噪聲的一些參數(shù)。