Matlab实现 理想低通、巴特沃斯低通、高斯低通、理想高通、巴特沃斯高通、高斯高通（d=10,50,150）

最新推荐文章于 2024-08-12 21:32:34 发布
最新推荐文章于 2024-08-12 21:32:34 发布
阅读量1.3w 收藏 104
点赞数 18
分类专栏： matlab算法 文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/dooonald/article/details/78525631
版权

数字图像处理第二次编程课后作业

理想低通（d=10,50,150）:

close all;
clear all;

%% ---------Ideal Lowpass Filters (Fre. Domain)------------
f = imread('D:/testData/cameraman.tif');
f = mat2gray(f,[0 255]);

[M,N] = size(f);
P = 2*M;
Q = 2*N;
fc = zeros(M,N);

for x = 1:1:M
    for y = 1:1:N
        fc(x,y) = f(x,y) * (-1)^(x+y);
    end
end

F = fft2(fc,P,Q);

H_0 = zeros(P,Q);
H_1 = zeros(P,Q);
H_2 = zeros(P,Q);

for x = (-P/2):1:(P/2)-1
     for y = (-Q/2):1:(Q/2)-1
        D = (x^2 + y^2)^(0.5);
        if(D <= 10)  H_0(x+(P/2)+1,y+(Q/2)+1) = 1; end    
        if(D <= 50)  H_1(x+(P/2)+1,y+(Q/2)+1) = 1; end    
        if(D <= 150)  H_2(x+(P/2)+1,y+(Q/2)+1) = 1; end
     end
end

G_0 = H_0 .* F;
G_1 = H_1 .* F;
G_2 = H_2 .* F;

g_0 = real(ifft2(G_0));
g_0 = g_0(1:1:M,1:1:N);

g_1 = real(ifft2(G_1));
g_1 = g_1(1:1:M,1:1:N);

g_2 = real(ifft2(G_2));
g_2 = g_2(1:1:M,1:1:N);         

for x = 1:1:M
    for y = 1:1:N
        g_0(x,y) = g_0(x,y) * (-1)^(x+y);
        g_1(x,y) = g_1(x,y) * (-1)^(x+y);
        g_2(x,y) = g_2(x,y) * (-1)^(x+y);
    end
end

%% -----show-------
figure();
subplot(1,2,1);
imshow(f,[0 1]);
xlabel('1).原图');

subplot(1,2,2);
imshow(log(1 + abs(F)),[ ]);
xlabel('2).1图的傅里叶光谱');



figure();
subplot(1,2,1);
imshow(H_0,[0 1]);
xlabel('3).理想低通滤波(D=10)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_0(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_0)),[ ]);
xlabel('4).3得到的滤波结果');

subplot(1,2,2);
imshow(g_0,[0 1]);
xlabel('5).效果图(D=10)');



figure();
subplot(1,2,1);
imshow(H_1,[0 1]);
xlabel('6).理想低通滤波(D=50)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_1(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_1)),[ ]);
xlabel('7).6得到的滤波结果');

subplot(1,2,2);
imshow(g_1,[0 1]);
xlabel('8).效果图(D=50)');



figure();
subplot(1,2,1);
imshow(H_2,[0 1]);
xlabel('9).理想低通滤波(D=150)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_2(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_2)),[ ]);
xlabel('10).9得到的滤波结果');

subplot(1,2,2);
imshow(g_2,[0 1]);
xlabel('11).效果图(D=150)');

巴特沃斯低通（d=10,50,150）:

close all;
clear all;

%% ---------Butterworth Lowpass Filters (Fre. Domain)------------
f = imread('D:/testData/cameraman.tif');
f = mat2gray(f,[0 255]);

[M,N] = size(f);
P = 2*M;
Q = 2*N;
fc = zeros(M,N);

for x = 1:1:M
    for y = 1:1:N
        fc(x,y) = f(x,y) * (-1)^(x+y);
    end
end

F = fft2(fc,P,Q);

H_0 = zeros(P,Q);
H_1 = zeros(P,Q);
H_2 = zeros(P,Q);

for x = (-P/2):1:(P/2)-1
     for y = (-Q/2):1:(Q/2)-1
        D = (x^2 + y^2)^(0.5);
        %n=1
        D_0 = 10;
        H_0(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D/D_0)^2);
        D_0 = 50;
        H_1(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D/D_0)^2);   
        D_0 = 150;
        H_2(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D/D_0)^2);
        %H_2(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D/D_0)^6);
     end
end

G_0 = H_0 .* F;
G_1 = H_1 .* F;
G_2 = H_2 .* F;

g_0 = real(ifft2(G_0));
g_0 = g_0(1:1:M,1:1:N);

g_1 = real(ifft2(G_1));
g_1 = g_1(1:1:M,1:1:N);

g_2 = real(ifft2(G_2));
g_2 = g_2(1:1:M,1:1:N);         

for x = 1:1:M
    for y = 1:1:N
        g_0(x,y) = g_0(x,y) * (-1)^(x+y);
        g_1(x,y) = g_1(x,y) * (-1)^(x+y);
        g_2(x,y) = g_2(x,y) * (-1)^(x+y);
    end
end

%% -----show-------
figure();
subplot(1,2,1);
imshow(f,[0 1]);
xlabel('1).原图');

subplot(1,2,2);
imshow(log(1 + abs(F)),[ ]);
xlabel('2).1图的傅里叶光谱');



figure();
subplot(1,2,1);
imshow(H_0,[0 1]);
xlabel('3).巴特沃斯低通滤波(D=10)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_0(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_0)),[ ]);
xlabel('4).3得到的滤波结果');

subplot(1,2,2);
imshow(g_0,[0 1]);
xlabel('5).效果图(D=10)');



figure();
subplot(1,2,1);
imshow(H_1,[0 1]);
xlabel('6).巴特沃斯低通滤波(D=50)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_1(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_1)),[ ]);
xlabel('7).6得到的滤波结果');

subplot(1,2,2);
imshow(g_1,[0 1]);
xlabel('8).效果图(D=50)');



figure();
subplot(1,2,1);
imshow(H_2,[0 1]);
xlabel('9).巴特沃斯低通滤波(D=150)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_2(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_2)),[ ]);
xlabel('10).9得到的滤波结果');

subplot(1,2,2);
imshow(g_2,[0 1]);
xlabel('11).效果图(D=150)');

高斯低通（d=10,50,150）:

close all;
clear all;
clc;
%% ---------Gaussian Lowpass Filters (Fre. Domain)------------
f = imread('D:/testData/cameraman.tif');
f = mat2gray(f,[0 255]);

[M,N] = size(f);
P = 2*M;
Q = 2*N;
fc = zeros(M,N);

for x = 1:1:M
    for y = 1:1:N
        fc(x,y) = f(x,y) * (-1)^(x+y);
    end
end

F = fft2(fc,P,Q);

H_0 = zeros(P,Q);
H_1 = zeros(P,Q);
H_2 = zeros(P,Q);

for x = (-P/2):1:(P/2)-1
     for y = (-Q/2):1:(Q/2)-1
        D = (x^2 + y^2)^(0.5);
        D_0 = 10;
        H_0(x+(P/2)+1,y+(Q/2)+1) = exp(-(D*D)/(2*D_0*D_0)); 
        D_0 = 50;
        H_1(x+(P/2)+1,y+(Q/2)+1) = exp(-(D*D)/(2*D_0*D_0));   
        D_0 = 150;
        H_2(x+(P/2)+1,y+(Q/2)+1) = exp(-(D*D)/(2*D_0*D_0));
     end
end

G_0 = H_0 .* F;
G_1 = H_1 .* F;
G_2 = H_2 .* F;

g_0 = real(ifft2(G_0));
g_0 = g_0(1:1:M,1:1:N);

g_1 = real(ifft2(G_1));
g_1 = g_1(1:1:M,1:1:N);

g_2 = real(ifft2(G_2));
g_2 = g_2(1:1:M,1:1:N);         

for x = 1:1:M
    for y = 1:1:N
        g_0(x,y) = g_0(x,y) * (-1)^(x+y);
        g_1(x,y) = g_1(x,y) * (-1)^(x+y);
        g_2(x,y) = g_2(x,y) * (-1)^(x+y);
    end
end


%% -----show-------
close all;

figure();
subplot(1,2,1);
imshow(f,[0 1]);
xlabel('1).原图');

subplot(1,2,2);
imshow(log(1 + abs(F)),[ ]);
xlabel('2).1图的傅里叶光谱');



figure();
subplot(1,2,1);
imshow(H_0,[0 1]);
xlabel('3).高斯低通滤波(D=10)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_0(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_0)),[ ]);
xlabel('4).3得到的滤波结果');

subplot(1,2,2);
imshow(g_0,[0 1]);
xlabel('5).效果图(D=10)');



figure();
subplot(1,2,1);
imshow(H_1,[0 1]);
xlabel('6).高斯低通滤波(D=50)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_1(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_1)),[ ]);
xlabel('7).6得到的滤波结果');

subplot(1,2,2);
imshow(g_1,[0 1]);
xlabel('8).效果图(D=50)');



figure();
subplot(1,2,1);
imshow(H_2,[0 1]);
xlabel('9).高斯低通滤波(D=150)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_2(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_2)),[ ]);
xlabel('10).9得到的滤波结果');

subplot(1,2,2);
imshow(g_2,[0 1]);
xlabel('11).效果图(D=150)');

理想高通（d=10,50,150）:

close all;
clear all;

%% ---------Ideal Lowpass Filters (Fre. Domain)------------
f = imread('D:/testData/cameraman.tif');
f = mat2gray(f,[0 255]);

[M,N] = size(f);
P = 2*M;
Q = 2*N;
fc = zeros(M,N);

for x = 1:1:M
    for y = 1:1:N
        fc(x,y) = f(x,y) * (-1)^(x+y);
    end
end

F = fft2(fc,P,Q);

H_0 = zeros(P,Q);
H_1 = zeros(P,Q);
H_2 = zeros(P,Q);

for x = (-P/2):1:(P/2)-1
     for y = (-Q/2):1:(Q/2)-1
        D = (x^2 + y^2)^(0.5);
        if(D > 10)  H_0(x+(P/2)+1,y+(Q/2)+1) = 1; end    
        if(D > 50)  H_1(x+(P/2)+1,y+(Q/2)+1) = 1; end    
        if(D > 150)  H_2(x+(P/2)+1,y+(Q/2)+1) = 1; end
     end
end

G_0 = H_0 .* F;
G_1 = H_1 .* F;
G_2 = H_2 .* F;

g_0 = real(ifft2(G_0));
g_0 = g_0(1:1:M,1:1:N);

g_1 = real(ifft2(G_1));
g_1 = g_1(1:1:M,1:1:N);

g_2 = real(ifft2(G_2));
g_2 = g_2(1:1:M,1:1:N);         

for x = 1:1:M
    for y = 1:1:N
        g_0(x,y) = g_0(x,y) * (-1)^(x+y);
        g_1(x,y) = g_1(x,y) * (-1)^(x+y);
        g_2(x,y) = g_2(x,y) * (-1)^(x+y);
    end
end

%% -----show-------
figure();
subplot(1,2,1);
imshow(f,[0 1]);
xlabel('1).原图');

subplot(1,2,2);
imshow(log(1 + abs(F)),[ ]);
xlabel('2).1图的傅里叶光谱');



figure();
subplot(1,2,1);
imshow(H_0,[0 1]);
xlabel('3).理想高通滤波(D=10)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_0(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_0)),[ ]);
xlabel('4).3得到的滤波结果');

subplot(1,2,2);
imshow(g_0,[0 1]);
xlabel('5).效果图(D=10)');



figure();
subplot(1,2,1);
imshow(H_1,[0 1]);
xlabel('6).理想高通滤波(D=50)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_1(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_1)),[ ]);
xlabel('7).6得到的滤波结果');

subplot(1,2,2);
imshow(g_1,[0 1]);
xlabel('8).效果图(D=50)');



figure();
subplot(1,2,1);
imshow(H_2,[0 1]);
xlabel('9).理想高通滤波(D=150)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_2(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_2)),[ ]);
xlabel('10).9得到的滤波结果');

subplot(1,2,2);
imshow(g_2,[0 1]);
xlabel('11).效果图(D=150)');

巴特沃斯高通（d=10,50,150）:

close all;
clear all;

%% ---------Butterworth Lowpass Filters (Fre. Domain)------------
f = imread('D:/testData/cameraman.tif');
f = mat2gray(f,[0 255]);

[M,N] = size(f);
P = 2*M;
Q = 2*N;
fc = zeros(M,N);

for x = 1:1:M
    for y = 1:1:N
        fc(x,y) = f(x,y) * (-1)^(x+y);
    end
end

F = fft2(fc,P,Q);

H_0 = zeros(P,Q);
H_1 = zeros(P,Q);
H_2 = zeros(P,Q);

for x = (-P/2):1:(P/2)-1
     for y = (-Q/2):1:(Q/2)-1
        D = (x^2 + y^2)^(0.5);
        %n=1
        D_0 = 10;
        H_0(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D_0/D)^2);
        D_0 = 50;
        H_1(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D_0/D)^2);   
        D_0 = 150;
        H_2(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D_0/D)^2);
        %H_2(x+(P/2)+1,y+(Q/2)+1) = 1/(1+(D/D_0)^6);
     end
end

G_0 = H_0 .* F;
G_1 = H_1 .* F;
G_2 = H_2 .* F;

g_0 = real(ifft2(G_0));
g_0 = g_0(1:1:M,1:1:N);

g_1 = real(ifft2(G_1));
g_1 = g_1(1:1:M,1:1:N);

g_2 = real(ifft2(G_2));
g_2 = g_2(1:1:M,1:1:N);         

for x = 1:1:M
    for y = 1:1:N
        g_0(x,y) = g_0(x,y) * (-1)^(x+y);
        g_1(x,y) = g_1(x,y) * (-1)^(x+y);
        g_2(x,y) = g_2(x,y) * (-1)^(x+y);
    end
end

%% -----show-------
figure();
subplot(1,2,1);
imshow(f,[0 1]);
xlabel('1).原图');

subplot(1,2,2);
imshow(log(1 + abs(F)),[ ]);
xlabel('2).1图的傅里叶光谱');



figure();
subplot(1,2,1);
imshow(H_0,[0 1]);
xlabel('3).巴特沃斯高通滤波(D=10)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_0(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_0)),[ ]);
xlabel('4).3得到的滤波结果');

subplot(1,2,2);
imshow(g_0,[0 1]);
xlabel('5).效果图(D=10)');



figure();
subplot(1,2,1);
imshow(H_1,[0 1]);
xlabel('6).巴特沃斯高通滤波(D=50)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_1(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_1)),[ ]);
xlabel('7).6得到的滤波结果');

subplot(1,2,2);
imshow(g_1,[0 1]);
xlabel('8).效果图(D=50)');



figure();
subplot(1,2,1);
imshow(H_2,[0 1]);
xlabel('9).巴特沃斯高通滤波(D=150)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_2(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_2)),[ ]);
xlabel('10).9得到的滤波结果');

subplot(1,2,2);
imshow(g_2,[0 1]);
xlabel('11).效果图(D=150)');

高斯高通（d=10,50,150）:

close all;
clear all;
clc;
%% ---------Gaussian Lowpass Filters (Fre. Domain)------------
f = imread('D:/testData/cameraman.tif');
f = mat2gray(f,[0 255]);

[M,N] = size(f);
P = 2*M;
Q = 2*N;
fc = zeros(M,N);

for x = 1:1:M
    for y = 1:1:N
        fc(x,y) = f(x,y) * (-1)^(x+y);
    end
end

F = fft2(fc,P,Q);

H_0 = zeros(P,Q);
H_1 = zeros(P,Q);
H_2 = zeros(P,Q);

for x = (-P/2):1:(P/2)-1
     for y = (-Q/2):1:(Q/2)-1
        D = (x^2 + y^2)^(0.5);
        D_0 = 10;
        H_0(x+(P/2)+1,y+(Q/2)+1) = 1 - exp(-(D*D)/(2*D_0*D_0)); 
        D_0 = 50;
        H_1(x+(P/2)+1,y+(Q/2)+1) = 1 - exp(-(D*D)/(2*D_0*D_0));   
        D_0 = 150;
        H_2(x+(P/2)+1,y+(Q/2)+1) = 1 - exp(-(D*D)/(2*D_0*D_0));
     end
end

G_0 = H_0 .* F;
G_1 = H_1 .* F;
G_2 = H_2 .* F;

g_0 = real(ifft2(G_0));
g_0 = g_0(1:1:M,1:1:N);

g_1 = real(ifft2(G_1));
g_1 = g_1(1:1:M,1:1:N);

g_2 = real(ifft2(G_2));
g_2 = g_2(1:1:M,1:1:N);         

for x = 1:1:M
    for y = 1:1:N
        g_0(x,y) = g_0(x,y) * (-1)^(x+y);
        g_1(x,y) = g_1(x,y) * (-1)^(x+y);
        g_2(x,y) = g_2(x,y) * (-1)^(x+y);
    end
end


%% -----show-------
close all;

figure();
subplot(1,2,1);
imshow(f,[0 1]);
xlabel('1).原图');

subplot(1,2,2);
imshow(log(1 + abs(F)),[ ]);
xlabel('2).1图的傅里叶光谱');



figure();
subplot(1,2,1);
imshow(H_0,[0 1]);
xlabel('3).高斯高通滤波(D=10)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_0(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_0)),[ ]);
xlabel('4).3得到的滤波结果');

subplot(1,2,2);
imshow(g_0,[0 1]);
xlabel('5).效果图(D=10)');



figure();
subplot(1,2,1);
imshow(H_1,[0 1]);
xlabel('6).高斯高通滤波(D=50)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_1(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_1)),[ ]);
xlabel('7).6得到的滤波结果');

subplot(1,2,2);
imshow(g_1,[0 1]);
xlabel('8).效果图(D=50)');



figure();
subplot(1,2,1);
imshow(H_2,[0 1]);
xlabel('9).高斯高通滤波(D=150)');

subplot(1,2,2);
h = mesh(1:20:P,1:20:Q,H_2(1:20:P,1:20:Q));
set(h,'EdgeColor','k');
axis([0 P 0 Q 0 1]);
xlabel('u');ylabel('v');
zlabel('|H(u,v)|');

figure();
subplot(1,2,1);
imshow(log(1 + abs(G_2)),[ ]);
xlabel('10).9得到的滤波结果');

subplot(1,2,2);
imshow(g_2,[0 1]);
xlabel('11).效果图(D=150)');
bllddee
关注
专栏目录
matlab设计理想高斯巴特沃斯低通滤波器_完整二阶有源带通滤波器设计！（下载：教程+原理图+视频+代码）...
weixin_39929465的博客
11-21 1982
1、背景对于微弱的信号的处理方式一般是：放大和滤波，这个过程中就涉及到放大电路的选取、滤波器的选择以及偏置电路的设计。本例以实例的方式讲解并附带参数计算、仿真、实物测试三个环节。假设需要处理一个20mV的正弦信号，该信号的频率范围是15~35Hz，经过处理后幅值不超过3.3V，且需要经过带通滤波器滤除杂波。2、滤波器定义滤波电路又称为滤波器，是一种选频电路，能够使特定频率范围的信号通过，...
matalb 图像处理 低通滤波和高通滤波 （理想，巴特沃斯，高斯 含代码）
joreng的博客
11-24 1万+
主要分为理想低通滤波，巴特沃斯低通滤波，高斯低通滤波。
matlab实现理想滤波、Butterworth、Gaussian滤波器
05-12
用MATLAB在频域实现理想滤波、Butterworth、Gaussian滤波器，
matlab设计理想高斯巴特沃斯低通滤波器_巴特沃斯切比雪夫贝塞尔
weixin_39831170的博客
11-24 1152
看到这三个人的名字，学过模电的同学应该有点面熟，他们出现在有源滤波那一部分。要是学过《信号与系统》和《数字信号处理》课程的，对这三个名字就更熟了。这三个人是谁？是干啥的呢？斯蒂芬·巴特沃斯(Stephen Butterworth)，是个英国工程师，他1930年在英国《无线电工程》期刊发表了一篇有关滤波器的论文，后来以他的名字命名这种滤波器。巴特沃斯滤波器也被称作最大平坦滤波器, 其特点是...
c++ opencv数字图像处理：频率域滤波--同态滤波
最新发布
生活需要深度
08-12 126
文章目录 前言 一、同态滤波原理 1.处理原理 二、同态滤波器模板及MATLAB代码 1.同态滤波器 2.代码 前言 数字图像处理c++ opencv（VS2019 opencv4.53）持续更新一、同态滤波原理 1.处理原理 （1）认为图像f(x,y)由两部分组成：照射分量i(x,y)，反射分量r(x,y): f ( x , y ) = i ( x , y ) ∗ r ( x , y ) f(x,y)=i(x,y)*r(x,y) f(x,y)=i(x,y)∗r(x,y)（2）但上式不能直接用于对两个分量在
一维高斯低通滤波器：返回一维高斯低通滤波器的系数-matlab开发
06-01
此函数返回高斯低通滤波器的系数。 高斯滤波器的优点：时域无振铃或过冲。 缺点：频域滚降慢。 将 SR=sampling rate，fco=cutoff freq（均以 Hz 为单位）传递给函数。 计算长度为 L（L 总是奇数）的 FIR 滤波器的系数。 这个长度为 L=2N+1 的对称 FIR 滤波器具有延迟 N/SR 秒。 使用示例： 计算 SR=1000,fco=50 Hz 时的高斯滤波器频率响应： freqz(gaussfiltcoef(1000,50),1,256,1000); 使用 fco=500 的高斯滤波器以 5kHz 采样的滤波信号 X： y=filter(gaussfiltcoef(5000,500),1,X); SR,fco 没有经过健全性检查。 WCR 2006-10-11。
ideal low pass filter
04-21
matlab ideal low pass filter
计算机视觉（二）-matlab之理想低通滤波器，布特沃斯低通、高斯低通，理想高通、布特沃斯。。。.pdf
04-21
理想低通滤波器是一种最简单的滤波器，它可以将频谱中的高频分量完全去除，仅保留低频分量。其频率响应函数为： H(u,v) = {1, D(u,v) ≤ D0 {0, otherwise 其中，D(u,v)是频率坐标，D0是截止频率。 2. 布特沃斯...
巴特沃斯低通滤波matlab实现代码-DIP-Filter:DIP过滤器
05-20
巴特沃斯低通滤波matlab实现代码 DIP-Filter 1.概况 项目：实现一个通用的高通、低通、带通和带阻滤波器函数。其中又分别实现理想、巴特沃思和指数等滤波形式。用实现的函数对图1（lena_noise.bmp）进行低通处理，图...
matlab设计理想高斯巴特沃斯低通滤波器_方程推导：二阶有源带通滤波器设计！（下载：教程+原理图+视频+代码）...
weixin_39739170的博客
11-21 994
1、背景对于微弱的信号的处理方式一般是：放大和滤波，这个过程中就涉及到放大电路的选取、滤波器的选择以及偏置电路的设计。本例以实例的方式讲解并附带参数计算、仿真、实物测试三个环节。假设需要处理一个20mV的正弦信号，该信号的频率范围是15~35Hz，经过处理后幅值不超过3.3V，且需要经过带通滤波器滤除杂波。2、滤波器定义滤波电路又称为滤波器，是一种选频电路，能够使特定频率范围的信号通过，...
OpenCV C++实现频域滤波器汇总 理想低通、理想高通、巴特沃思低通、巴特沃思高通、高斯低通、高斯高通
cyf15238622067的博客
02-26 6401
根据各个滤波公式及前几章内容，可实现理想低通、理想高通、巴特沃思低通、巴特沃思高通、高斯低通、高斯高通滤波器； 效果图： 高斯低通 D0=30 巴特沃思低通 D0=30,n=2 理想低通D0=30： 高斯高通D0=80： 巴特沃思高通 D0=80,n=2 理想高通D0=80： 代码实现： #include&lt;opencv2/opencv.hpp&gt;...
verilog 浮点转定点_在FPGA上优化实现复数浮点计算
weixin_39608301的博客
12-08 571
电子万花筒中国最活跃的射频微波天线雷达技术群欢迎您加入高性能浮点处理一直与高性能CPU相关联。在过去几年中，GPU也成为功能强大的浮点处理平台，超越了图形，称为GP-GPU(通用图形处理单 元)。新创新是在苛刻的应用中实现基于FPGA的浮点处理。本文的重点是FPGA及其浮点性能和设计流程，以及OpenCL的使用，这是高性能浮点计算前沿的编程语言。各种处理平台的GFLOP指标在不断提高，现...
计算机视觉（二）-matlab之理想低通滤波器，布特沃斯低通、高斯低通，理想高通、布特沃斯高通、高斯高通滤波器
博观而约取，厚积而薄发
04-08 2万+
在滤波器之前，先讲解傅里叶变换 理想低通滤波器 f = imread('Fig0441.tif'); f = im2double(f); % 计算填充图像大小 [M,N] = size(f); M2 = 2*M; N2 = 2*N; % 傅里叶变换 F = fftshift(fft2(f,M2,N2)); figure;imshow(mat2gray(log(1+abs(F))));ti...
MATLAB图像处理--高斯低通滤波、高斯高通滤波（代码及示例）
热门推荐
留缘的博客
04-28 7万+
常用的滤波器有以下几种： 频域滤波器 频域平滑滤波器 理想低通滤波器 巴特沃斯低通滤波器 高斯低通滤波器 频域锐化滤波器 理想高通滤波器 巴特沃斯高通滤波器 高斯高通滤波器 频域拉普拉斯算子 高频提升滤波 同态滤波器 低通滤波器通过过滤掉频域的高频成分达到平滑目的。 滤波的基本模型G(u,v) = H(u,v)F(u,v)，其中 F(u,v) 是图像的傅立叶变换结果， H(u,v) 称为滤波器传输......
IIR数字滤波器设计（数字信号处理）
Five—菜鸟级的博客
12-18 7万+
~~ 如果有什么问题可以在我的个人博客留言 ，我会及时回复。欢迎来访交流 ~~ IIR数字滤波器设计（数字信号处理） 一、实验目的 1.熟悉双线性变换法设计IIR数字滤波器的原理与方法。 2.掌握IIR数字滤波器的MATLAB实现方法，会调用ellipord()和ellip() 函数设计各种滤波器。 3.观察分析滤波器输入输出数据波形，理解数字滤波的概念...
巴特沃斯数字低通滤波器
诚朴勇毅
10-08 1万+
matlab实现信号滤波有很多种方法，本文采用filter函数，先由butter函数生成数字滤波器分子、分母多项式，然后将这些参数传给filter函数，便可实现滤波的功能。 知识点如下： 1、[b,a] = butter(N,wc); 其中，b和a分别为滤波器分子、分母多项式，N为滤波器阶数，wc为3dB截止频率，巴特沃斯低通滤波器的特点为：整个频带内单调下降，且非常光滑，并且阶数越高，阻带内下降的越陡峭、越干脆。无论阶数多高，都在经过（wc，3dB)一点。 2、y = filter(b,a,x);
图像处理之理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器的matlab简单实现
qq_44111805的博客
08-13 3万+
经理想高频滤波后的图像把信息丰富的低频去掉了，丢失了许多必要的信息**。一般情况下，高通滤波对噪声没有任何抑制作用**，若简单的使用高通滤波，图像质量可能由于噪声严重而难以达到满意的改善效果。这种滤波器实际上是由一个高通滤波器和一个全通滤波器构成的，这样便能在高通滤波的基础上保留低频信息。高通滤波可以分为理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器。其中D0为巴特沃斯高通滤波器的截止频率，其中D0为理想高通滤波器的截止频率，其中D0为高斯高通滤波器的截止频率，三、巴特沃斯高通滤波器(BHPF)....
数字图像处理实验十理想、高斯、巴特沃斯低通和高通滤波器
氮气蜗牛的博客
10-10 3055
利用理想低通、高斯低通、巴特沃斯低通滤波器进行滤波 利用理想高通、高斯高通、巴特沃斯高通滤波器进行滤波
数字图像处理---低高通滤波实验（MATLAB实现）
OpenCV学习笔记
04-29 1万+
数字图像处理MATLAB实现高低通滤波器代码以及实验原理和实验效果图
matlab均值滤波中值滤波理想低通滤波巴特沃斯低通滤波高斯高通滤波代码
07-14
以上是MATLAB中实现均值滤波、中值滤波、理想低通滤波、巴特沃斯低通滤波和高斯高通滤波的代码。参数说明：image为输入图像，windowSize为滤波窗口的大小，D0为截止频率，n为巴特沃斯滤波器的阶数。输出结果为滤波后...
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值