matlab基于SAD的区域立体匹配方法

%function SADMatch


 %im1=imread（‘Test_left.jpg‘）;
 %im2=imread（‘Test_right.jpg‘）;
 % im1=imread（‘im0.jpg‘）;
 % im2=imread（‘im1.jpg‘）;
 im1=imread（‘left.BMP‘）;
 im2=imread（‘right.BMP‘）;


%if isrgb（im1）
    %im1=rgb2gray（im1）;
%end
%imshow（im1）;
im1=double（im1）;

%if isrgb（im2）
    %im2=rgb2gray（im2）;
%end
%figure
%imshow（im2）;
im2=double（im2）;

D=20; %最大视差
N=5; %窗口大小的一半
[HW]=size（im1）;

%计算右图减去左图，相减产生D个矩阵放到imgDiff中 
imgDiff=zeros（HWD）;
e=zeros（HW）;
for i=1:D
    fprintf（‘%g\n‘i）
    e（:1:（W-i））=abs（im2（:1:（W-i））- im1（: