fm_guifun.m

介绍了傅立叶变换的基本原理

data.input2 = rotsclIm; data.input2_windowed = window2d(size(rotsclIm,1),size(rotsclIm,2),data.windowType).*rotsclIm; 
elseif (size_rs(1)==size_in1(1))&amt;(size_rs(2)==size_in1(2)) > then no need to crop anything 
data.input2 = rotsclIm; data.input2_windowed = window2d(size(rotsclIm,1),size(rotsclIm,2),data.windowType).*rotsclIm; 
data.input1 = in1ref; data.input1_windowed = window2d(size(in1ref,1),size(in1ref,2),data.windowType).*in1ref; 
else 
disp('fm_guifun.m (290ish): AUTOCROP does not currently work on these types of images (ie, certain rectangle shapes)........') 
end 
end 

end 
end 

set(gcbf,'CurrentAxes',findobj(gcbf,'Tag','Axes2')); 
data.h2 = findobj(gcbf,'Tag','Axes2'); 

cla; 
set(gcbf,'Userdata',data); 

updateImage(1,0); 
updateImage(2,0); 
data = get(gcbf,'Userdata'); 

vel1 = size(data.input1_windowed); 
vel2 = size(data.input2_windowed); 

set(findobj(gcbf,'Tag','Ref_im_c'),'String',[data.filename1 ', ' int2str(vel1(1)) ' x ' int2str(vel1(2))],'ForegroundColor','k'); 
set(findobj(gcbf,'Tag','Sens_im_c'),'String',[data.filename2 ', ' int2str(vel2(1)) ' x ' int2str(vel2(2))],'ForegroundColor','k'); 
set(gcbf,'Userdata',data); 

dispTag('Sens_im','k'); 
dispText('','b'); 

else 
disp('Cannot rotate/scale a non-existent image') 
end 

if data.performanceLevel==1 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peak)']); 
else 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peaks)']); 
end 
set(gcbf,'Userdata',data); 


case('dispText') 
data = get(gcbf,'Userdata'); 
val = get(findobj(gcbf,'Tag','cb_dispText'),'Value'); 
if val==1 
data.dispText = 1; 
else 
data.dispText = 0; 
end 
set(gcbf,'Userdata',data); 

case('register') 
data = get(gcbf,'Userdata'); 
set(findobj(gcbf,'Tag','Pushbutton1'),'String','please wait...'); 

if ((size(data.input2reference,1)>0)&amt;(size(data.input2reference,1)>0)&amt;(size(data.input1,1)>0)&amt;(size(data.input1,2)>0)) 
[data.registered,data.input1registered,data.input2registered,reg_output,data.pc_rs,data.pc_trans] = fourier_mellin(data); 

regImDims = size(data.registered); 
data.regInfo = reg_output; 
set(findobj(gcbf,'Tag','TransOut'),'String',[num2str(data.regInfo.translation(1)) 'x' num2str(data.regInfo.translation(2))],'ForegroundColor','k'); 
set(findobj(gcbf,'Tag','RotOut'),'String',num2str(data.regInfo.rotation),'ForegroundColor','k'); 
set(findobj(gcbf,'Tag','ScaleOut'),'String',num2str(data.regInfo.scale),'ForegroundColor','k'); 
set(findobj(gcbf,'Tag','TransPeakOut'),'String',num2str(data.regInfo.trans_peak),'ForegroundColor','k'); 
set(findobj(gcbf,'Tag','RSPeakOut'),'String',num2str(data.regInfo.rs_peak),'ForegroundColor','k'); 

set(gcbf,'CurrentAxes',findobj(gcbf,'Tag','Axes3')); 
data.h2 = findobj(gcbf,'Tag','Axes3'); 
end 
if data.performanceLevel==1 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peak)']); 
else 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peaks)']); 
end 

set(gcbf,'Userdata',data); 
dispText('','b'); 

case('setPerformanceLevel') 
data = get(gcbf,'Userdata'); 
data.performanceLevel = floor(get(findobj(gcbf,'Tag','performLevel'),'Value')); 
if data.performanceLevel==1 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peak)']); 
else 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peaks)']); 
end 
set(gcbf,'Userdata',data); 


> invoke help window 
case('help') 
data = get(gcbf,'Userdata'); 

dispText('','b'); 
global FM_PATH; 
web(['file:///' FM_PATH 'help/help_fm.html']); 

> BELOW IS THE SELECTIONS FOR THE ANALYSIS WINDOW 
case('input1Select') 
data = get(gcbf,'Userdata'); 
set(findobj(gcbf,'Tag','Pushbutton1'),'String','please wait...'); 
sel = get(findobj(gcbf,'Tag','input1analysis'),'Value'); 

if (sel==1) > input image 'Input|Magnitude Spectrum|Phase Spectrum|Log-Polar' 
mx = data.input1; 
elseif (sel==2) 
mx = data.input1_windowed; 
elseif(sel==3) > magnitude spectrum 
mx = data.input1_magSpec; 
elseif(sel==4) > phase spectrum 
mx = angle(data.input1_freq); 
elseif(sel==5) > log-polar 
mx = data.input1_freq_lp; 
elseif(sel==6) > windowed log-polar 
mx = data.windowed_input1_freq_lp; 
elseif(sel==7) > invariant 
mx = data.logmagspec_lp_ms1; 
else 
disp('input1Select: Should never get here') 
end 

set(gcbf,'CurrentAxes',findobj(gcbf,'Tag','Axes1')); > move to the correct axes 

cla; > clear the image corresponding to these axes 
imagesc(mx); > display image 
axis image; > fit axis box tightly around image 
axis ij; > puts MATLAB into its "matrix" axes mode. The coordinate system origin is at the upper left corner. The i axis is vertical and is numbered from top to bottom. The j axis is horizontal and is numbered from left to right. 
axis off; > turns off labelling 
colormap('gray'); > not sure where 'm' is coming from 
vel = size(data.input1); 
set(findobj(gcbf,'Tag','Ref_im_c'),'Sting',[data.filename1 ', ' int2str(vel(1)) ' x ' int2str(vel(2))],'ForegroundColor','k'); 

if data.performanceLevel==1 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peak)']); 
else 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peaks)']); 
end 

set(gcbf,'Userdata',data); 

case('input2Select') 
data = get(gcbf,'Userdata'); 
set(findobj(gcbf,'Tag','Pushbutton1'),'String','please wait...'); 
sel = get(findobj(gcbf,'Tag','input2analysis'),'Value'); 

if (sel==1) > input image 'Input|Magnitude Spectrum|Phase Spectrum|Log-Polar|Invariant' 
mx = data.input2; 
elseif(sel==2) 
mx = data.input2_windowed; 
elseif(sel==3) > magnitude spectrum 
mx = data.input2_magSpec; 
elseif(sel==4) > phase spectrum 
mx = angle(data.input2_freq); 
elseif(sel==5) > log-polar 
mx = data.input2_freq_lp; 
elseif(sel==6) > windowed log-polar 
mx = data.windowed_input2_freq_lp; 
elseif(sel==7) > invariant 
mx = data.logmagspec_lp_ms2; 
else 
disp('input2Select: Should never get here') 
end 

set(gcbf,'CurrentAxes',findobj(gcbf,'Tag','Axes2')); > move to the correct axes 

cla; > clear the image corresponding to these axes 
imagesc(mx); > display image 
axis image; > fit axis box tightly around image 
axis ij; > puts MATLAB into its "matrix" axes mode. The coordinate system origin is at the upper left corner. The i axis is vertical and is numbered from top to bottom. The j axis is horizontal and is numbered from left to right. 
axis off; > turns off labelling 
colormap('gray'); > not sure where 'm' is coming from 
if data.performanceLevel==1 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peak)']); 
else 
set(findobj(gcbf,'Tag','Pushbutton1'),'String',['REGISTER (' num2str(data.performanceLevel) 'peaks)']); 
end 

set(gcbf,'Userdata',data); 

case('regSelect') 
data = get(gcbf,'Userdata'); 
sel = get(findobj(gcbf,'Tag','reganalysis'),'Value'); 

if (sel==1) > regiistered image 'Registered|Registered Image 2|Log-Polar PC|Spatial PC' 
mx = data.registered; 
elseif(sel==2) 
mx = data.input1registered; 
elseif(sel==3) > magnitude spectrum 
mx = data.input2registered; 
elseif(sel==4) > phase spectrum 
mx = data.pc_rs; 
elseif(sel==5) > log-polar 
mx = data.pc_trans; 
else 
disp('input2Select: Should never get here') 
end 

set(gcbf,'CurrentAxes',findobj(gcbf,'Tag','Axes3')); > move to the correct axes 

cla; > clear the image corresponding to these axes 
imagesc(mx); > display image 
axis image; > fit axis box tightly around image 
axis ij; > puts MATLAB into its "matrix" axes mode. The coordinate system origin is at the upper left corner. The i axis is vertical and is numbered from top to bottom. The j axis is horizontal and is numbered from left to right. 
axis off; > turns off labelling 
> colormap(m); > not sure where 'm' is coming from 


set(gcbf,'Userdata',data); 
> dispTag('Ref_im','k'); 
> dispText('','b'); 

otherwise 
disp(['fm_guifun.m: trying to call non-existent switch...',action]) 
end 

> ------------------------------------------------------------------------------------ 
> *************************** auxilliary functions 

function dispText(txt,colr); 

set(findobj(gcbf,'Tag','MessText'),'String',txt,'ForegroundColor',colr); 

> ------------------------------------------------------------------------------------ 
function vr=vrat(ktery); 

idx = get(findobj(gcbf,'Tag',ktery),'value'); 
val = get(findobj(gcbf,'Tag',ktery),'String'); 
vr=str2num(val(idx)); 


> ------------------------------------------------------------------------------------ 
function ramek(kde,barva); 

set(findobj(gcbf,'Tag',kde),'XColor',barva); 
set(findobj(gcbf,'Tag',kde),'YColor',barva); 
set(findobj(gcbf,'Tag',kde),'ZColor',barva); 

> ------------------------------------------------------------------------------------ 
function dispTag(txt,colr); 

set(findobj(gcbf,'Tag',txt),'ForegroundColor',colr); 

> ------------------------------------------------------------------------------------ 
function updateImage(im,LP_ONLY) 
> USAGE: updateImage(im,LP_ONLY) A.I.Wilmer, 2002 
> 
> function to update magnitude spectra, log-polar plots etc of image 'im' 
> LP_ONLY : only update the log-polar plot and stuff dependent on it 

data = get(gcbf,'Userdata'); 

if (im==1) > then update image 1 information 
if (~LP_ONLY) > if log-polar setting is changed then don't need to do the next couple of lines 
data.input1_freq = fftshift(fft2(data.input1_windowed)); 
data.input1_magSpec = hipass_filter(size(data.input1_freq,1),size(data.input1_freq,2)).*abs(data.input1_freq); 
> data.input1_magSpec = log10(abs(data.input1_freq)); 
end 
data.input1_freq_lp = imlogpolar(data.input1_magSpec,data.logpolarScaleRes,data.logpolarAngleRes,data.LogInterp); 
data.windowed_input1_freq_lp = repmat(window1d(size(data.input1_freq_lp,1),data.windowType),1,size(data.input1_freq_lp,2)).*data.input1_freq_lp; 
data.logmagspec_lp_ms1 = hipass_filter(size(data.input1_freq_lp,1),size(data.input1_freq_lp,2)).*abs(fftshift(fft2(data.input1_freq_lp))); 
elseif (im==2) > update image 2 plots 
if (~LP_ONLY) > if log-polar setting is changed then don't need to do the next couple of lines 
data.input2_freq = fftshift(fft2(data.input2_windowed)); 
data.input2_magSpec = hipass_filter(size(data.input2_freq,1),size(data.input2_freq,2)).*abs(data.input2_freq); 
> data.input2_magSpec = log10(abs(data.input2_freq)); 
end 
data.input2_freq_lp = imlogpolar(data.input2_magSpec,data.logpolarScaleRes,data.logpolarAngleRes,data.LogInterp); 
data.windowed_input2_freq_lp = repmat(window1d(size(data.input2_freq_lp,1),data.windowType),1,size(data.input2_freq_lp,2)).*data.input2_freq_lp; 
data.logmagspec_lp_ms2 = hipass_filter(size(data.input2_freq_lp,1),size(data.input2_freq_lp,2)).*abs(fftshift(fft2(data.input2_freq_lp))); 
else 
disp('updateImage(): incorrect image number used.') 
end 
set(gcbf,'Userdata',data);