【正文】 算法依據(jù)最大互信息理論判斷所得參數(shù)是否最優(yōu)，若不是，則繼續(xù)搜索較優(yōu)參數(shù)，在搜索時會不斷重復“空間幾何變換（ affine）統(tǒng)計聯(lián)合直方圖（ PV插值法） 計算互信息值 最優(yōu)化判斷”的過程，直至搜索到滿足精度要求的參數(shù)；最后輸出配準參數(shù)。 M文件及其源代碼 ． function varargout = ImageRegistration(varargin) gui_Singleton = 1。 gui_State = struct(39。gui_Name39。, mfilename, ... 39。gui_Singleton39。, gui_Singleton, ... 39。gui_OpeningF39。, @ImageRegistration_OpeningF, ... 39。gui_OutputF39。, @ImageRegistration_OutputF, ... 39。gui_LayoutF39。, [], ... 39。gui_Callback39。, [])。 if nargin amp。amp。 ischar(varargin{1}) 輸入?yún)⒖紙D像 輸入浮動圖像 設置初始點和初始搜素方向 空間幾何變換 計算互信息值 最優(yōu)化 否 是 輸出配準參數(shù) = str2func(varargin{1})。 end if nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。 else gui_mainf(gui_State, varargin{:})。 end addpath(pwd)。 function ImageRegistration_OpeningF(hObject, eventdata, handles, varargin) = hObject。 guidata(hObject, handles)。 function varargout = ImageRegistration_OutputF(hObject, eventdata, handles) varargout{1} = 。 function pushbutton1_Callback(hObject, eventdata, handles) global I。 %%%調用 、圖像大小 %%% [filename ,pathname]=uigetfile({39。*.jpg39。39。*.bmp39。39。*.bmp39。},39。209。161。212。241。205。188。198。172。39。)。 str=[pathname filename]。 I=imread(str)。 axes()。 imshow(I)。 =I。 guidata(hObject,handles)。 figure(1)。 imshow()。 function pushbutton3_Callback(hObject, eventdata, handles) =。 =。 [I,J]=GLPF(handles)。 =I。 =J。 guidata(hObject,handles)。 tic RegistrationParameters=Powell(handles)。 toc ElapsedTime=toc。 =RegistrationParameters。 y=RegistrationParameters(1)。 x=RegistrationParameters(2)。 ang=RegistrationParameters(3)。 MI_Value=RegistrationParameters(4)。 RegistrationResult=sprintf(39。X,Y,Angle=[%.5f][%.5f][%.5f]39。,x,y,ang)。 MI_Value=sprintf(39。MI_Value=[%.4f]39。,MI_Value)。 ElapsedTime=sprintf(39。Elapsed Time=[%.3f]39。,ElapsedTime)。 axes() [RegistrationImage]=Register(handles)。 imshow(RegistrationImage) set(,39。string39。,RegistrationResult)。 set(,39。string39。,MI_Value)。 set(,39。string39。,ElapsedTime)。 function pushbutton2_Callback(hObject, eventdata, handles) global J。 [filename ,pathname]=uigetfile({39。*.jpg39。39。*.bmp39。39。*.bmp39。},39。209。161。212。241。205。188。198。172。39。)。 str=[pathname filename]。 J=imread(str)。 axes()。 imshow(J)。 =J。 guidata(hObject,handles)。 figure(2)。 imshow()。 function [RegistrationParameters]=Powell(handles) e=。 X0=[0 0 0]。 D=[1 0 0。0 1 0。0 0 1]。 num=0。 while(num200) num=num+1。 d1=D(1,:)。%d1206。170。190。216。213。243。D181。196。181。218。210。187。208。208。163。172。179。245。202。188。203。209。203。247。189。207。242。 [X1,fX1]=OneDimSearch(X0,d1,handles)。 d2=D(2,:)。%d2206。170。190。216。213。243。D181。196。181。218。182。254。208。208。163。172。179。245。202。188。203。209。203。247。189。207。242。 [X2,fX2]=OneDimSearch(X1,d2,handles)。 d3=D(3,:)。%d3206。170。190。216。213。243。D181。196。181。218。200。253。208。208。163。172。179。245。202。188。203。209。203。247。189。207。242。163。172。200。253。206。172。203。209。203。247。185。202。211。208。200。253。184。246。189。207。242。 [X3,fX3]=OneDimSearch(X2,d3,handles)。 fX0=PV(X0(1),X0(2),X0(3),handles)。 Diff=[fX1fX0 fX2fX1 fX3fX2]。 [maxDiff,m]=max(Diff)。%max186。175。202。253。181。196。211。195。168。163。172。181。187。216。maxdiff206。170。207。242。193。191。Diff181。196。238。180。243。212。170。203。216。163。172。m206。170。238。180。243。212。170。203。216。181。196。208。242。186。197。 d4=X3X0。 temp1=X3X0。 Conditon1=sum(temp1.*temp1)。 if Conditon1=e break end [X4,fX4,landa]=OneDimSearch(X0,d4,handles)。 X0=X4。 temp2=X4X3。 Conditon2=sum(temp2.*temp2)。 if Conditon2=e X3=X4。 break end temp3=sqrt((fX4fX0)/(maxDiff+eps))。 if(abs(landa)temp3) D(4,:)=d4。 for i=m:3 D(i,:)=D(i+1,:)。 end end end RegistrationParameters(1)=X3(1)。 RegistrationParameters(2)=X3(2)。 RegistrationParameters(3)=X3(3)。 RegistrationParameters(4)=fX3。 function [Y,fY,landa]=OneDimSearch(X,direction,handles) %210。187。206。172。203。209。203。247。178。201。211。195。brent189。168。 a=5。b=5。 Epsilon=。 cgold=。 IterTimes=200。 if ab temp=a。 a=b。 b=temp。 end v=a+cgold*(ba)。 w=v。 x=v。 e=。 fx=Fx(x,X,direction)。 fv=fx。 fw=fx。 for iter=1:IterTimes xm=*(a+b)。 if abs(xxm)=Epsilon**(ba) break end if abs(e)Epsilon r=(xw)*(fxfv)。 q=(xv)*(fxfw)。 p=(xv)*(q(xw)*r。 q=2*(qr)。