[問題] 請問有關分析圖片的方式已回收
我想要用這個分析出來的傅利葉的頻譜
來找出我其他做contrast的處理
看看有哪些cpd範圍有變化
哪些不變
對不起
我現在遇到了一個要分析圖片中的空間頻率上的問題
我不知道如何定義多少cpd
但我知道人類的閾值是從0.1~70cpd左右
但我現在想要把每個範圍算出來 到底有多少
我想分成
0~0.1、
0.1~0.5、
0.5~2、
2~10、
10~30、
30~50、
50~70、
70~100,
這幾種cpd範圍
找出其中的關連性
請問一下 要用什麼function 才可以寫出這樣的需求
還有要怎麼定義分析出來的圖片的fourier中的t
可以請大大幫忙嗎?
對不起
以下是我現在用的分析出來的程式
img = imread('D:\Documents and Settings\soiipopoalan.NONE-6EC144771A\桌面\新資
料夾 (3)\5\5.tif');
subplot(1,2,1), imshow(img)
f1 = fft2(img);
f2 = fftshift(f1); %平移的2D傅利葉頻譜
P = log(1+abs(f2));
maxP = max(P(:));
subplot(1,2,2), imshow(P/maxP)
[x,y]=meshgrid(-128:127,-128:127);
z=sqrt(x.^2 + y.^2);
c=(z<15)
cf=fft2shift(fft2(z));
fftshow(cf,'log')
cm=imread('D:\Documents and Settings\soiipopoalan.NONE-6EC144771A\桌面\新資料
夾 (3)\5\5.tif');
cf=fftshift(fft2(cm));
figure,fftshow(cf,'log')
cf1=cf.*c;
figure,fftshow(cf1,'log')
cfli=ifft2(cf1);
figure,fftshow(cfli,'abs')
[x,y]=meshgrid(-128:127,-128:127);
z=sqrt(x.^2 + y.^2);
c=(z>15)
cfh=cf.*c;
figure,fftshow(cfh,'log')
cfhi=ifft2(cfh);
figure,fftshow(cfhi,'abs')
謝謝你們
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.137.51.66
※ 編輯: soiialan 來自: 140.137.51.66 (03/09 13:49)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章