clc
clear
Rmuban=imread（‘Rmuban.jpg‘）;%%原圖像
Rmuban=rgb2gray（Rmuban）;
%?imshow（Rmuban）
J=Rmuban;
%?J=medfilt2（Rmuban）;
%?figure;imshow（J）???%%中值濾波獲取目標圖像


muban=uint8（zeros（256256））;
for?i=94:160
????for?j=58:128
????????muban（ij）=Rmuban（i-93j-57）;
????end
end
muban1=uint8（zeros（256256））;
for?m=94:160
????????for?n=129:199
????????????muban1（mn）=J（m-93n-128）;
????????end
end
zonghe=muban+muban1;
imshow（zonghe）
%%生成輸入圖像

?Forizonghe=fft2（double（zonghe）256256）;
?Foricenterzonghe=fftshift（Forizonghe）;
?Forizongheabs=abs（Foricenterzonghe）;
?Forizongheabs=double（Forizongheabs）;
?Forizonghe2=Forizongheabs.*Forizongheabs;
?figure;mesh（double（Forizonghe2））;axis?tight;title（‘綜合功率譜‘）;%%計算聯合變換功率譜
?
?Iforizero=ifft2（Forizongheabs）;
?Iforizerocenter=fftshift（Iforizero）;
?figure;mesh（Iforizerocenter）;?axis?tight;title（‘相關峰輸出的三維顯示‘）;?%%實現相關峰的三維顯示
?Iforizerocenterim=uint8（（abs（Iforizerocenter）/max（max（abs（Iforizerocenter））））*255）;
?figure;imshow（Iforizerocenterim）;title（‘去零級相關輸出‘）;%%去零級強干擾項

?figure;?mesh（abs（Iforizerocenter））;axis?tight;title（‘去零級相關輸出的三維顯示‘）;%%實現去零級聯合相關的三維圖像顯示

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????44046??2011-05-28?23:41??圖像聯合變換器\2.jpg

?????文件???????1132??2009-04-22?13:54??圖像聯合變換器\Rmuban.JPG

?????文件???????1219??2011-06-01?08:41??圖像聯合變換器\xgglv.m

?????文件?????352392??2011-05-30?13:38??圖像聯合變換器\光學聯合變換相關系統的設計和研究_（1）.kdh

?????文件?????933223??2011-05-30?00:22??圖像聯合變換器\基于Matlab仿真的聯合變換相關器實驗.pdf

?????文件?????320437??2011-05-28?21:08??圖像聯合變換器\基于Matlab的圖像聯合變換相關識別的實現.pdf

?????目錄??????????0??2011-06-24?09:55??圖像聯合變換器

-----------?---------??----------?-----??----

??????????????1652449????????????????????7