对一幅灰度图像:
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 梁某对强某的死亡是...