clc
clear
%圖像拼接
f=imread（‘1.jpg‘）;
g=imread（‘21.jpg‘）;%加載匹配圖像
f=double（f）;
g=double（g）;


????B=f（::1）;

????G=f（::2）;

????R=f（::3）;

????B1=g（::1）;

????G1=g（::2）;

????R1=g（::3）;

????
Y=R*0.299+G*0.587+B*0.114;?%YUV法??彩圖轉化為灰度圖
Y1=R1*0.299+G1*0.587+B1*0.114;??%YUV法??彩圖轉化為灰度圖

[mn]=size（Y）;

%~~~~~~~~~~~~~小波分析取概貌~~~~~~~~~~~~~~~~~~~
ori_im=xbbh（Ymn）;
[mn]=size（Y1）;
ori_im1=xbbh（Y1mn）;
ori_im=uint8（ori_im）;
ori_im1=uint8（ori_im1）;

fx?=?[-2?-1?0?1?2];?????????????????%?x方向梯度算子（用于Harris角點提取算法）?
Ix?=?filter2（fxori_im）;??%?x方向濾波?
Ix1?=?filter2（fxori_im1）;
fy?=?[-2;-1;0;1;2];?????????????????%?y方向梯度算子（用于Harris角點提取算法）?
Iy?=?filter2（fyori_im）;?%?y方向濾波?
Iy1?=?filter2（fyori_im1）;?
Ix2?=?Ix.^2;?
Iy2?=?Iy.^2;?
Ixy?=?Ix.*Iy;?
Ix21?=?Ix1

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

?????文件??????33280??2010-05-25?21:30??我的綜述.doc

?????文件?????109620??2009-06-01?14:46??2.jpg

?????文件??????83751??2009-06-01?14:46??1.jpg

?????文件???????4714??2010-05-31?17:39??pinjie.m

?????文件????????687??2010-05-13?10:09??xbbh.m

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

???????????????232052????????????????????5