1 2 2 (POC)

Size: px

Start display at page:

Download "1 2 2 (POC)"

えりかかくはり
7 years ago
Views:

1 2007 Taichi Kawasaki Kinki University School of Engineerings

2 1 2 2 (POC)

3 1 [1]. 1 [2] P.Chazal [3] 2 3 2

4 1: 3

5 2: 4

6 3: 5

7 2 (POC) f(n 1, n 2 ) g(n 1, n 2 ) 2 F (k 1, k 2 ) G(k 1, k 2 ) F (k 1, k 2 ) G(k 1, k 2 ) R F G (k 1, k 2 ) R F G (k 1, k 2 ) = F (k 1, k 2 )G(k 1, k 2 ) F (k 1, k 2 )G(k 1, k 2 ) (1) G(k 1, k 2 ) G(k 1, k 2 ) r fg (n 1, n 2 ) R F G (k 1, k 2 ) 2 r fg (n 1, n 2 ) = IDF T [R F G (k 1, k 2 )] (2) IDFT 6

8 3 4 Lenna 5 MANDRILL Lenna 6 2 POC r fg (n 1, n 2 ) Lenna MANDRILL 7 2 POC r fg (n 1, n 2 ) 2 7

9 4: Lenna 8

10 5: MANDRILL 9

11 6: 2 POC 10

12 7: POC 11

13 4 MAT- LAB MATLAB (GUI) [4, 5] 12

14 5 8 9 POC ( ) 13

15 8: 14

16 9: 15

17 10: 16

18 11: ( ) 17

19 6 12: 18

20 13: 19

21 14: 20

22 15: 21

23 16: 22

24 17: 23

25 7 MATLAB 24

26 [1] [2] [3] P.Chazal, J.Flynn, and B.Rcilly Automatcd Processing of Images Based on the Fourier Transform for Use in Forensic Science IEEE Trans., Vol. PAMI-27, 3, pp (Mar. 2005). [4] MATLAB [5] MATLAB

27 8 26

28 9 function varargout = lab2(varargin) gui Singleton = 1; gui State = struct( gui Name, mfilename,... gui Singleton, gui Singleton,... gui OpeningFcn,... gui OutputFcn,... gui LayoutFcn, [],... gui Callback, []); if nargin && ischar(varargin1) gui State.gui Callback = str2func(varargin1); end if nargout [varargout1:nargout] = gui mainfcn(gui State, varargin:); else gui mainfcn(gui State, varargin:); end function lab2 OpeningFcn(hObject, eventdata, handles, varargin) handles.output = hobject; guidata(hobject, handles); function varargout = lab2 OutputFcn(hObject, eventdata, handles) varargout1 = handles.output; // function Untitled 1 Callback(hObject, eventdata, handles) // function Untitled 2 Callback(hObject, eventdata, handles) axes(handles.axes1); [fname, dpath] = uigetfile( *.bmp,, C:\Documents and Settings\ \My Documents\MATLAB\ ); IMG=imread(fname); handles.img=img; imagesc(img); 27

29 colormap(gray); guidata(hobject,handles); // function Untitled 3 Callback(hObject, eventdata, handles) axes(handles.axes2); [fname, dpath] = uigetfile( *.bmp,, C:\Documents and Settings\ \My Documents\MATLAB\ ); IMG2=imread(fname); handles.img2=img2; imagesc(img2); colormap(gray); guidata(hobject,handles); // function Untitled 5 Callback(hObject, eventdata, handles) axes(handles.axes3); x=fft2(double(handles.img)); y=fft2(double(handles.img2)); p=(y.*conj(x))./abs(y.*conj(x)); IMG=fftshift(ifft2(p)); imagesc([-256,256],[-256,256],img); // function Untitled 7 Callback(hObject, eventdata, handles) axes(handles.axes4); x=fft2(double(handles.img)); y=fft2(double(handles.img2)); p=(y.*conj(x))./abs(y.*conj(x)); IMG=fftshift(ifft2(p)); xx = 1; yy = 1; max=img(1,1); [y,x]=size(img); for y1=1:y for x1=1:x if (IMG(y1,x1) max) 28

30 max = IMG(y1,x1); // yy = y1; xx = x1; end end end x2=xx-257; y2=yy-257; mesh(img); // xlabel( x ); ylabel( y ); // HT=uicontrol(lab2, style, text ); set(ht, position,[ ]); set(ht, string,max); set(ht, BackgroundColor,[1,1,1]); set(ht, fontsize,15); set(ht, ForegroundColor,[0,0,0]); // (x ) HT=uicontrol(lab2, style, text ); set(ht, position,[ ]); set(ht, string,xx); set(ht, BackgroundColor,[1,1,1]); set(ht, fontsize,15); set(ht, ForegroundColor,[0,0,0]); // (y ) HT=uicontrol(lab2, style, text ); set(ht, position,[ ]); set(ht, string,yy); set(ht, BackgroundColor,[1,1,1]); set(ht, fontsize,15); set(ht, ForegroundColor,[0,0,0]); 29

31 //x HT=uicontrol(lab2, style, text ); set(ht, position,[ ]); set(ht, string,x2); set(ht, BackgroundColor,[1,1,1]); set(ht, fontsize,15); set(ht, ForegroundColor,[0,0,0]); //y HT=uicontrol(lab2, style, text ); set(ht, position,[ ]); set(ht, string,y2); set(ht, BackgroundColor,[1,1,1]); set(ht, fontsize,15); set(ht, ForegroundColor,[0,0,0]); // function Untitled 4 Callback(hObject, eventdata, handles) // function Untitled 9 Callback(hObject, eventdata, handles) exit 30

32 10 31

Ver.1 1/17/2003 2

Ver.1 1/17/2003 1 Ver.1 1/17/2003 2 Ver.1 1/17/2003 3 Ver.1 1/17/2003 4 Ver.1 1/17/2003 5 Ver.1 1/17/2003 6 Ver.1 1/17/2003 MALTAB M GUI figure >> guide GUI GUI OK 7 Ver.1 1/17/2003 8 Ver.1 1/17/2003 Callback