图像处理第二次作业 2

发布 2022-07-12 13:09:28 阅读 1463

对一幅灰度图像:

1)计算并画出此图像的中心化频率谱。

2)分别用高斯低通和高斯高通滤波器对图像进行频域处理,并对结果进行分析。

3)用频域拉普拉斯算子对此图像进行锐化处理,并对结果进行分析。

1)程序。clear;

clc; a=imread(''jpg');

b=rgb2gray(a);

figure;

subplot(2,2,1);imshow(b);title('原图像');

b_d=double(b);

m,n] =size(b);

fft2函数实现。

if ismatrix(b_d)%ismatrix:确定输入的是否为矩阵。

if nargin==1%nargin用来判断输入变量的个数。

f = fftn(b_d);%fftn:n维离散傅里叶变换。

elsef = fftn(b_d,[m n]);

endelse

ifnargin==1

f = fft(fft(b_d,2),[1);

elsef = fft(fft(b_d,n,2),m,1);

endend

c3=log(1+abs(f));

subplot(2,2,2);imshow(c3,title('原图像频谱');

计算并画出此图像的中心化频率谱。

m,n] =size(b);

fori=1:m

for j=1:n

c1(i,j)=(1)^(i+j))*b_d(i,j));

endend

c2=log(1+abs(fft2(c1)))

c4=fftshift(log(1+abs(fft2(b_d)))

subplot(2,2,3);imshow(c2,

title('编写函数实现频谱中心化')

subplot(2,2,4);imshow(c4,

title('matlab函数实现频谱中心化')

分别用低通滤波和高通滤波对此图像进行频域处理。

c5=fftshift(fft2(b_d));

p1=round(m/2);

q1=round(n/2);

d0=30;%半径取30,随着d0的增加,低通滤波图像越来越模糊。

b_lf=c5;

b_hf=c5;

低通滤波器。

fori=1:m

for j=1:n

d=sqrt((i-p1)^2+(j-q1)^2);

if d<=d0

flag=1;

elseflag=0;

endb_lf(i,j)=flag*c5(i,j);

endend

b_lf=uint8(real(ifft2(ifftshift(b_lf)))

figure;

高通滤波器。

fori=1:m

for j=1:n

d=sqrt((i-p1)^2+(j-q1)^2);

if d>d0

flag=1;

elseflag=0;

endb_hf(i,j)=flag*c5(i,j);

endend

b_hf=uint8(real(ifft2(ifftshift(b_hf)))

subplot(1,3,1);imshow(b);title('原图像');

subplot(1,3,2);imshow(b_lf);

title('高斯低通滤波');

subplot(1,3,3);imshow(b_hf);

title('高斯高通滤波');

用频域拉普拉斯算子对此图像进行锐化处理。

fori=1:m

for j=1:n

d1=(i-p1)^2+(j-q1)^2;

h(i,j)=(4)*(pi)^2)*d1;

h1(i,j)=h(i,j)*c5(i,j);

endend

h2=real(ifft2(ifftshift(h1)))

h2max=max(h2(:)

h3=h2(1:m,1:n);

k=1;h4=b_d-k*h3/h2max;

figure;

subplot(1,2,1);imshow(b);title('原图像');

subplot(1,2,2); imshow(h4,

title('拉普拉斯图像');

使用moon图像。

task3=imread(''tif');

task3_1=double(task3);

task3_2=fftshift(fft2(task3_1));

m1,n1]=size(task3);

p2=round(m1/2);

q2=round(n1/2);

fori=1:m1

for j=1:n1

d1=(i-p2)^2+(j-q2)^2;

h(i,j)=(4)*(pi)^2)*d1;

h1(i,j)=h(i,j)*task3_2(i,j);

endend

h2=real(ifft2(ifftshift(h1)))

h2max=max(h2(:)

h3=h2(1:m1,1:n1);

k=1;h4=task3_1-k*h3/h2max;

figure;

subplot(1,2,1);imshow(task3);title('原图像');

subplot(1,2,2); imshow(h4,

title('拉普拉斯图像');

2)实验结果。

2 第二次作业 答案

2011 2012 2 经济法第二次作业。姓名班级学号得分 一 单选题。1 下列合同中不属于要式合同的有 d a 中外合资经营企业合同 b 中外合作经营企业合同。c 技术转让合同。d 国际 合同。2 下列合同中,应当采取书面形式订立的有 a a 技术转让合同b 租赁合同。c 房屋租赁合同d 保管合同...

第二次作业

第 次作业日期姓名 得分 一 给下列词语注音。废墟调遣履行凄凉。荒无人烟潸然泪下井然有序。生死攸关绘声绘色薄雾笼罩。二 查字典填空。负 字用音序查字法,应查音序 音节 负 在字典中的解释有 背 担任 遭受 享有 背弃 违背 败 输。请为下面词语中的 负 选择合适的解释。序号填空 负荆请罪 久负盛名 ...

第二次作业

三 案例题。梁某与好友强某深夜在酒吧喝酒。强某醉酒后,钱包从裤袋里掉到地上,梁某拾后见钱包里有5000元现金就将其隐匿。强某要梁某送其回家,梁某怕钱包之事被发现,托辞拒绝。强某在回家途中醉倒在地,被人发现时已经冻死。关于本案,请问 1 梁某占有财物的行为构成何种犯罪并说明理由。2 梁某对强某的死亡是...