i=imread('picture.tif');%將要濾波的影象讀到i中,當然要換成你的名字了
imshow(i);
f=double(i);
g=fft2(f);
g=fftshift(g);
[n1,n2]=size(g);
n=2;
d0=50; %此處d0為截止頻率,你可以依次換成60、70、80進行濾波比較
n1=fix(n1/2);
n2=fix(n2/2);
for i=1:n1
for j=1:n2
d=sqrt((i-n1)^2+(j-n2)^2);
h=1/(1+0.414*(d/d0)^(2*n));
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
x2=ifft2(result);
j=uint8(real(x2));
figure,imshow(j);
高通 低通濾波
一 關於影象高頻和低頻 1.1 對影象高頻訊號和低頻訊號的理解 影象中的低頻訊號和高頻訊號也叫做低頻分量和高頻分量。簡單一點說,影象中的高頻分量,指的是影象強度 亮度 灰度 變化劇烈的地方,也就是 我們常說的邊緣 輪廓 影象中的低頻分量,指的是影象強度 亮度 灰度 變換平緩的地方,也就是大片色塊的地...
低通濾波opencv
理論基礎知識,高通部分已經講過,這裡不重複了。本部分做低通濾波,保留細節。即頻譜中間部分保留,與高通相反。o cv2.imread image lena256.bmp 0 讀入 dft cv2.dft np.float32 o flags cv2.dft complex output dshift ...
高通濾波與低通濾波公式
from 設 x n 為取樣資料 y n 為濾波結果 初始y 0 x 0 高通濾波 i從1到n 1迴圈操作 y i y i 1 x i x i 1 a 1,迴圈乘就不斷弱化上次濾波的值,使取樣間的差值相對強化,a 1時沒有濾波效果。在不斷變弱上次濾波值的基礎上,增加本次取樣值與上次取樣值之間的變化。...