数字图像处理-第三周 实践课

创建一个10x10的矩阵，$A(i,j)=i+j$
创建一个10x10的矩阵，$A(i,j)=\sqrt{(i-1{)}^2+(j-1{)}^2}$

clear;

A=1:10;
B=ones(10,1);
J=B*A;
I=J';
AAA=I+J;
BBB=sqrt((I-1).^2+(J-1).^2);

创建一个100x100的图像。以A(20,50)为圆心，由内向外颜色由黑渐变为白色。

clear;
A=1:100;
B=ones(100,1);
J=B*A;
I=J';
BBB=sqrt((I-20).^2+(J-50).^2);
g=BBB./max(BBB);
imshow(g);

图像相加降噪 1、给定一个原始图像 2、通过添加随机噪声产生50张图像 3、比较 10，30，50 张图像叠加(求平均)后的效果 4、并用直方图显示直方图随叠加张数的改变

clear; 
img=imread('d:\dsp\flower.jpg'); 
img=rgb2gray(img); 
for i=1:50 
    B(:,:, i) = uint8(double(img)+normrnd(0,50,size(img))); 
end
C3_10 = uint8(mean( B(:,:,1:10), 3)); 
C3_30 = uint8(mean( B(:,:,1:30), 3)); 
C3_50 = uint8(mean( B(:,:,1:50), 3)); 
subplot(2,5, 1 ); 
imshow(img); 
subplot(2,5, 5+1); 
histogram(img); 
subplot(2,5, 2 ); 
imshow(B(:,:,5)); 
subplot(2,5, 5+2); 
histogram(B(:,:,5)); 
subplot(2,5, 3 ); 
imshow( C3_10 ); 
subplot(2,5, 5+3); 
histogram(C3_10); 
subplot(2,5, 4 ); 
imshow( C3_30 ); 
subplot(2,5, 5+4); 
histogram(C3_30); 
subplot(2,5, 5 ); 
imshow( C3_50 ); 
subplot(2,5, 5+5); 
histogram(C3_50);

幂次变换1、给定一张图像 S0 2、对图像进行大于1的幂次变换 S1 3、对图像进行小于1的幂次变换 S2 4、分别对S0,S1,S2进行直方图均衡化，S0E,S1E,S2E 5、显示三张图 1） S0,S0直方图，S0E, S0E直方图 2） S1,S1直方图，S1E, S1E直方图 3） S2,S2直方图，S2E, S2E直方图 比较三张图的直方图均衡化后的效果差异。

clear;
S0=rgb2gray(imread('d:\dsp\flower.jpg'));
S1=im2uint8(power((double(S0)./255),0.3));
S2=im2uint8(power((double(S0)./255),3));
S0E=histeq(S0);
S1E=histeq(S1);
S2E=histeq(S2);
figure(1);
subplot(221);imshow(S0);subplot(223);histogram(S0); 
subplot(222);imshow(S0E);subplot(224);histogram(S0E); 
figure(2);
subplot(221);imshow(S1);subplot(223);histogram(S1); 
subplot(222);imshow(S1E);subplot(224);histogram(S1E); 
figure(3);
subplot(221);imshow(S2);subplot(223);histogram(S2); 
subplot(222);imshow(S2E);subplot(224);histogram(S2E); 

生成浮雕效果的图像,模拟运动模糊图像

clear
img = imread('d:\dsp\flower.jpg'); 
img = rgb2gray(img); 
H = [1 0 0; 0 0 0; 0 0 -1]*5; 
G = uint8(filter2(H, img)+128); %最后统一加上128用于控制灰度,显示出类似浮雕的灰色
figure;
subplot(1,2,1); imshow(img); 
subplot(1,2,2); imshow(G); 

clear 
img = imread('d:\dsp\flower.jpg'); 
img = rgb2gray(img); 
H = eye(30,30)./30.0; 
G = uint8(filter2(H, img)); 
figure; 
subplot(1,2,1); imshow(img); 
subplot(1,2,2); imshow(G);