【文章內(nèi)容簡介】 if q0 p=p。 end q=abs(q)。 etemp=e。 e=d。 if not (abs(p)=abs(*q*etemp)||p=q*(ax)||p=q*(bx)) d=p/q。 u=x+d。 if uaEpsilon*2||buEpsilon*2 d=MySign(Epsilon,xmx)。 end else if x=xm e=ax。 else e=bx。 end d=cgold*e。 end else if x=xm。 e=ax。 else e=bx。 end d=cgold*e。 end if abs(d)=Epsilon u=x+d。 else u=x+MySign(Epsilon,d)。 end fu=Fx(u,X,direction)。 if fu=fx if u=x a=x。 else b=x。 end v=w。 fv=fw。 w=x。 fw=fx。 x=u。 fx=fu。 else if ux a=u。 else b=u。 end if fu=fw||w==x v=w。 fv=fw。 w=u。 fw=fu。 else if fu=fv||v==x||v==w v=u。 fv=fu。 end end end end landa=x。 Y=X+x*direction。 fY=Fx(x,X,direction)。 function[mySign]=MySign(a,b) if b0 Result=abs(a)。 else Resulr=abs(a)。 end mySign=Result。 function [fx]=Fx(x,X,direction,handles) fx=PV(X(1)+direction(1)*x,X(2)+direction(2)*x,(X(3)+direction(3)*x),handles)。 function [mi]=PV(x,y,ang,handles) a=。 a=double(a)。 b=。 b=double(b)。 [M,N]=size(a)。 hab=zeros(256,256)。 ha=zeros(1,256)。 hb=zeros(1,256)。 if max(max(a))~=min(min(a)) a=(amin(min(a)))/(max(max(a))min(min(a)))。 else a=zeros(M,N)。 end if max(max(b))~=min(min(b)) b=(bmin(min(b)))/(max(max(b))min(min(b)))。 else b=zeros(M,N)。 end a=double(int16(a*255))+1。 b=double(int16(b*255))+1。 [width,height]=size(b)。 u=(width1)/2。 v=(height1)/2。 rad=pi/180*ang。 t1=[1 0 0。0 1 0。x y 1]。 t2=[1 0 0。0 1 0。u v 1]。 t3=[cos(rad) sin(rad) 0。sin(rad) cos(rad) 0。0 0 1]。 t4=[1 0 0。0 1 0。u v 1]。 T=t2*t3*t4*t1。 tform=maketform(39。affine39。,T)。 coordinate_x=zeros(width,height)。 coordinate_y=zeros(width,height)。 for i=1:width for j=1:height coordinate_x(i,j)=i。 end end for i=1:width for j=1:height coordinate_y(i,j)=j。 end end [w z]=tforminv(tform,coordinate_x,coordinate_y)。 for i=1:width for j=1:height source_x=w(i,j)。 source_y=z(i,j)。 if(source_xwidth1||source_yheight1||double(uint16(source_x))=1||double(uint16(source_y))=1) hab(a(1,1),a(1,1))=hab(a(1,1),a(1,1))+1。 else m=fix(source_x)。 n=fix(source_y)。 index_b=b(i,j)。 index_a0=a(m,n)。 index_a1=a(m+1,n)。 index_a2=a(m,n+1)。 index_a3=a(m+1,n+1)。 dx=source_xm。 dy=source_yn。 hab(index_a0,index_b)=hab(index_a0,index_b)+(1dx)*(1dy)。 hab(index_a1,index_b)=hab(index_a1,index_b)+dx*(1dy)。 hab(index_a2,index_b)=hab(index_a2,index_b)+(1dx)*dy。 hab(index_a3,index_b)=hab(index_a3,index_b)+dx*dy。 end end end habsum=sum(sum(hab))。 index=find(hab~=0)。 pab=hab/habsum。 Hab=sum(sum(pab(index).*log2(pab(index))))。 pa=sum(pab,2)。 index=find(pa~=0)。 Ha=sum(sum(pa(index).*log2(pa(index))))。 pb=sum(pab,1)。 index=find(pb~=0)。 Hb=sum(sum(pb(index).*log2(pb(index))))。 mi=Ha+HbHab。 function[RegistrationImage]=Register(handles) I=。 J=。 y=(1)。 x=(2)。 ang=(3)。 [nrows,ncols]=size(J)。 width=nrows。 height=ncols。 new_J=uint8(zeros(width,height))。 a=(width1)/2。 c=a。 b=(height1)/2。 d=b。 rad=pi/180*ang。 t1=[1 0 0。0 1 0。x y 1]。 t2=[1 0 0。0 1 0。a b 1]。 t3=[cos(rad) sin(rad) 0。sin(rad) cos(rad) 0。0 0 1]。 t4=[1 0 0。0 1 0。c d 1]。 T=t2*t3*t4*t1。 tform=maketform(39。affine39。,T)。 tx=zeros(width,height)。 ty=zeros(width,height)。 for i=1:width for j=1:height tx(i,j)=i。 end end for i=1:width for j=1:height ty(i,j)=j。 end end [w z]=tforminv(tform,tx,ty)。 for i=1:width for j=1:height source_x=w(i,j)。 source_y=z(i,j)。 if(source_xwidth1||source_yheight1||double(uint16(source_x))=0||double(uint16(source_y))=0) new_J(i,j)=J(1,1)。 else if(source_x/double(uint16(source_x))==amp。amp。(source_y/double(uint16(source_y)==))) ne