%融合圖像
function?[FusionImageRegistrationImage]=Fusion（handles）
%%%%%%%%%%%%%%%%%%%獲取原圖像及配準參數%%%%%%%%%%%%%%%%%%%
I=handles.Old_I;
J=handles.Old_J;
x=handles.RegistrationParameters（1）;
y=handles.RegistrationParameters（2）;
ang=-handles.RegistrationParameters（3）;
%%%%%%%%%%%%%%%%對浮動圖像空間變換及插值，以便得到配準后的輸出圖像%%%%%%%%%
[nrowsncols]=size（J）;
width=nrows;
height=ncols;
new_J=uint8（zeros（widthheight））;

a=（width-1）/2;
c=a;
b=（height-1）/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=makeform（‘affine‘T）;

tx=zeros（widthheight）;
ty=zeros（widthheight）;
for?i=1:width
????for?j=1:height
????????tx（ij）=i;
????end
end
for?i=1:width
????for?j=1:height
????????ty（ij）=j;
????end
end
[w?z]=tforminv（tformtxty）;

for?i=1:width
????for?j=1:height
????????source_x=w（ij）;
????????source_y=z（ij）;
????????if?（source_x>=width-1||source_y>=height-1||...
????????????????double（uint16（source_x））<=0||...
????????????????double（uint16（source_y））<=0）
????????????new_J（ij）=J（11）;
????????else
????????????if?（source_x/double（uint16（source_x））==1.0）&...
????????????????????（source_y/double（uint16（source_y））==1.0）
????????????????new_J（ij）=J（int16（source_x）int16（source_y））;
????????????else
????????????????%a=double（round（source_x））;
????????????????%b=double（round（source_y））;
????????????????%new_J（des_xdes_y）=J（ab）;
????????????????a=double（uint16（source_x））;
????????????????b=double（uint16（source_y））;
????????????????x11=double（J（ab））;
????????????????x12=double（J（ab+1））;
????????????????x21=double（J（a+1b））;
????????????????x22=double（J（a+1b+1））;
????????????????new_J（ij）=uint8（（b+1-source_y）*（（sourec_x-a）*x21+...
????????????????????（a+1-sourec_x）*x11）+（source_y-b）*（（source_x-a）*x22...
????????????????????+（a+1-source_x）*x12））;
????????????end
????????end
????end
end
J=new_J;

I=uint8（I）;
J=uint8（J）;
RegistrationImage=uint8（J）;

%%%%%%%%%%%%%%%%%%%%產生“融合”后的效果圖%%%%%%%%%%%%%%%%%%%%%%
I=double（I）/255;
J=double（J）/255;
IJ=double（zeros（widthheight））;
for?m=1:width
????for?n=1:height
????????if?I（mn）>0.999||J（mn）>0.999
????????????IJ（mn）=0.8;
????????elseif?I（mn）==0||J（mn）==0
????????????IJ（mn）=0.01;
????????else
????????????IJ（mn）=（I（mn）*0.3+J（mn）*0.7）;
????????end
????end
end
IJ=IJ*255;
IJ=uint8（IJ）;
FusionImage=IJ;

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

?????文件???????2521??2011-08-15?00:55??8?圖像配準實例\Fusion.m

?????文件????????516??2011-08-15?00:55??8?圖像配準實例\GLPF.m

?????文件???????4235??2011-08-15?00:55??8?圖像配準實例\ImageRegistration.m

?????文件????????707??2011-08-15?00:55??8?圖像配準實例\OneDimensionSearch.m

?????文件???????4002??2011-08-15?00:55??8?圖像配準實例\OpenImage.m

?????文件???????1008??2011-08-15?00:55??8?圖像配準實例\Powell.m

?????文件???????2159??2011-08-15?00:55??8?圖像配準實例\PV.m

?????文件?????????87??2011-08-15?00:55??8?圖像配準實例\Test.m

?????目錄??????????0??2016-06-02?23:23??8?圖像配準實例

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

????????????????15235????????????????????9