光流计算

matlab实现的光流计算实例，实现图像中光流的计算，可用于计算人群量，物体跟踪等

function psidashFS = computePsidashFS_brox（ u v ）
% This function computes the value psi（i~j）（kl） as mentioned in the PhD
% thesis of Dr. Brox.

% Author: Visesh Chari 
% Centre for Visual Information Technology
% International Institute of Information Technology
% http://cvit.iiit.ac.in/
% http://research.iiit.ac.in/~visesh
%
% The Software is provided “as is“ without warranty of any kind.

[h w] = size（ u ） ;
psidashFS = zeros（ 2 * h + 1 2 * w + 1 ） ;

ux = convn（ u [1 -1] ） ; % x and y derivatives of u.
uy = convn（ u [1 -1]‘ ） ;

vx = convn（ v [1 -1] ） ; % x and y derivatives of v.
vy = convn（ v [1 -1]‘ ） ;

uxd = convn（ ux [1 1] / 2 ‘valid‘ ） ; % averaging as per the numerics section
vxd = convn（ vx [1 1] / 2 ‘valid‘ ） ; % of the second chapter.

uyd = convn（ uy [1 1]‘ / 2 ‘valid‘ ） ;
vyd = convn（ vy [1 1]‘ / 2 ‘valid‘ ） ;


t = convn（ uyd [1 1] / 2 ） ; % Computes the delta u（i+1/2 j） and delta u（i-1/2 j）.
uxpd = ux .^ 2 + t .^ 2 ;
t = convn（ uxd [1 1]‘ / 2 ） ; % Computes the delta u（i j+1/2） and delta u（i j-1/2）.
uypd = uy .^ 2 + t .^ 2 ;

t = convn（ vyd [1 1] / 2 ） ; % Computes the delta v（i+1/2 j） and delta v（i-1/2 j）.
vxpd = vx .^ 2 + t .^ 2 ;
t = convn（ vxd [1 1]‘ / 2 ） ; % Computes the delta v（i+1/2 j） and delta v（i-1/2 j）.
vypd = vy .^ 2 + t .^ 2 ;

% Finally g（ delta u ） （i+1/2 j） and （i-1/2 j） and （i j+1/2） and （i j-1/2）.
psidashFS（ 1:2:end 2:2:end ） = psiDerivative（ uypd + vypd ） ;
psidashFS（ 2:2:end 1:2:end ） = psiDerivative（ uxpd + vxpd ） ;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2009-06-23 16:27  Flow\
     目录           0  2009-06-23 16:23  Flow\brox_zip\
     文件        1937  2007-11-06 14:44  Flow\brox_zip\sor.m
     文件         835  2007-11-05 21:32  Flow\brox_zip\image_warp.m
     文件        2417  2009-06-23 13:50  Flow\brox_zip\resolutionProcess_brox.m
     文件        2415  2009-06-23 13:49  Flow\brox_zip\optic_flow_brox.m
     文件        1485  2009-06-23 13:48  Flow\brox_zip\gaussianSmooth.m
     文件         601  2009-06-23 13:49  Flow\brox_zip\psiDerivative.m
     文件         683  2009-06-23 13:48  Flow\brox_zip\gaussianRescaling.m
     文件        3788  2009-06-23 13:47  Flow\brox_zip\constructMatrix_brox.m
     文件         691  2009-06-23 13:47  Flow\brox_zip\gaussDeriv.m
     文件        1027  2007-11-14 21:28  Flow\brox_zip\split.m
     文件         804  2009-06-23 13:48  Flow\brox_zip\mywarp_rgb.m
     文件        1547  2009-06-23 13:47  Flow\brox_zip\computePsidashFS_brox.m
     文件         799  2009-06-23 13:49  Flow\brox_zip\imgGrad.m
     目录           0  2009-06-23 16:24  Flow\sand_zip\
     文件         600  2009-06-23 13:54  Flow\sand_zip\psiDerivative.m
     文件        2085  2009-06-23 13:54  Flow\sand_zip\resolutionProcess_sand.m
     文件        2279  2009-06-23 13:54  Flow\sand_zip\optic_flow_sand.m
     文件         804  2009-06-23 13:53  Flow\sand_zip\mywarp_rgb.m
     文件         799  2009-06-23 13:53  Flow\sand_zip\imgGrad.m
     文件         571  2009-06-23 13:53  Flow\sand_zip\getalphaImg.m
     文件        1486  2009-06-23 13:53  Flow\sand_zip\gaussianSmooth.m
     文件         684  2009-06-23 13:52  Flow\sand_zip\gaussianRescaling.m
     文件         692  2009-06-23 13:52  Flow\sand_zip\gaussDeriv.m
     文件        2981  2009-06-23 13:52  Flow\sand_zip\constructMatrix_sand.m
     文件        1074  2009-06-23 13:51  Flow\sand_zip\channels_new.m
     文件        1027  2007-11-14 21:28  Flow\sand_zip\split.m
     文件        1132  2009-06-23 13:52  Flow\sand_zip\computePsidashFS_sand.m
     文件        1965  2007-11-05 21:36  Flow\sand_zip\sor.m
     文件         835  2007-11-05 21:36  Flow\sand_zip\image_warp.m
............此处省略9个文件信息