function?[fx?fy?ft]?=?computeDerivatives（im1?im2）

if?size（im21）==0
????im2=zeros（size（im1））;
end

%?Horn-Schunck?original?method
fx?=?conv2（im10.25*?[-1?1;?-1?1]‘same‘）?+?conv2（im2?0.25*[-1?1;?-1?1]‘same‘）;
fy?=?conv2（im1?0.25*[-1?-1;?1?1]?‘same‘）?+?conv2（im2?0.25*[-1?-1;?1?1]?‘same‘）;
ft?=?conv2（im1?0.25*ones（2）‘same‘）?+?conv2（im2?-0.25*ones（2）‘same‘）;

%?derivatives?as?in?Barron
%?fx=?conv2（im1（1/12）*[-1?8?0?-8?1]‘same‘）;
%?fy=?conv2（im1（1/12）*[-1?8?0?-8?1]‘‘same‘）;
%?ft?=?conv2（im1?0.25*ones（2）‘same‘）?+?conv2（im2?-0.25*ones（2）‘same‘）;
%?fx=-fx;fy=-fy;

%?An?alternative?way?to?compute?the?spatiotemporal?derivatives?is?to?use?simple?finite?difference?masks.
%?fx?=?conv2（im1[1?-1]）;
%?fy?=?conv2（im1[1;?-1]）;
%?ft=?im2-im1;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-11-11?14:50??光流法matlab程序\
?????文件?????????759??2011-11-11?09:15??光流法matlab程序\computeDerivatives.m
?????文件??????526106??2011-10-12?09:48??光流法matlab程序\CT0.dcm
?????文件??????526110??2011-10-12?09:48??光流法matlab程序\CT5.dcm
?????文件?????????921??2011-11-11?09:15??光流法matlab程序\ffd.m
?????文件????????1080??2011-11-11?09:15??光流法matlab程序\gaussFilter.m
?????文件????????2818??2011-11-11?09:15??光流法matlab程序\HS.asv
?????文件????????2788??2011-11-11?14:40??光流法matlab程序\HS.m
?????文件??????236079??2011-11-11?09:15??光流法matlab程序\HSsnapshot.JPG
?????文件???????34678??2011-11-11?09:15??光流法matlab程序\lenag1.png
?????文件???????36184??2011-11-11?09:15??光流法matlab程序\lenag2.png
?????文件????????1382??2011-11-11?09:15??光流法matlab程序\plotFlow.asv
?????文件????????1382??2011-11-11?09:15??光流法matlab程序\plotFlow.m
?????文件?????????944??2011-11-11?09:15??光流法matlab程序\smoothImg.m
?????文件???????78886??2011-11-11?09:15??光流法matlab程序\yos10.tif
?????文件???????78908??2011-11-11?09:15??光流法matlab程序\yos9.tif