创建一个10x10的矩阵,
A
(
i
,
j
)
=
i
+
j
A(i,j)=i+j
A(i,j)=i+j
创建一个10x10的矩阵,
A
(
i
,
j
)
=
(
i
−
1
)
2
+
(
j
−
1
)
2
A(i,j)=\sqrt{(i-1)^{2}+(j-1)^{2}}
A(i,j)=(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);