%butterworth高通滤波
p=ones(100);
q=ones(size(p));
d=0.5; %决定漩涡的大小
n=2; %决定漩涡的尖点
[f1,f2]=freqspace(size(p),'meshgrid');
r=sqrt(f1.^2+f2.^2);
for i=1:size(p,1)
for j=1:size(p,2)
t=r(i,j)/(d*d);
q(i,j)=t^n/(t^n+1);
end
end
f=p.*q;
p=ones(100);
q=ones(size(p));
d=0.5; %决定漩涡的大小
n=2; %决定漩涡的尖点
[f1,f2]=freqspace(size(p),'meshgrid');
r=sqrt(f1.^2+f2.^2);
for i=1:size(p,1)
for j=1:size(p,2)
t=r(i,j)/(d*d);
q(i,j)=t^n/(t^n+1);
end
end
f=p.*q;
mesh(f1,f2,f);
模拟图像代码:
%butterworth高通滤波
p=ones(100);
q=ones(size(p));
d=0.5; %决定漩涡的大小
n=2; %决定漩涡的尖点
[f1,f2]=freqspace(size(p),'meshgrid');
r=sqrt(f1.^2+f2.^2);
for i=1:size(p,1)
for j=1:size(p,2)
t=r(i,j)/(d*d);
q(i,j)=t^n/(t^n+1);
end
end
f=p.*q;
mesh(f1,f2,f);