matlab运动估计代码.里面有源图像和运行结果.对写论文应该有帮助.rar

matlab运动估计代码.里面有源图像和运行结果.对写论文应该有帮助

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%  Name: Musawir Ali Shah
%%%  Assignment #3
%%%  Project title: Lucas Kanade Motion Estimation Using Pyramids （Level 4）
%%%  Note: The project specification says use density of 10 in plotting
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%% Usage: Lucas_Kanade（‘1.bmp‘‘2.bmp‘10）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function Lucas_Kanade（file1file2density）;
%% Read Images %%
img1 = im2double （imread （file1））;
%% Take alternating rows and columns %%
[odd1 even1] = split （img1）;
img2 = im2double （imread （file2））;
[odd2 even2] = split （img2）;
%% Run Lucas Kanade %%
[Dx Dy] = Estimate （odd1 odd2）;
%% Plot %%
figure; 
[maxImaxJ]=size（Dx）;
Dx=Dx（1:density:maxI1:density:maxJ）;
Dy=Dy（1:density:maxI1:density:maxJ）;
quiver（1:density:maxJ（maxI）:（-density）:1Dx-Dy1）;
axis square;


%% Run Lucas Kanade on all levels and interpolate %%
function [Dx Dy] = Estimate （img1 img2）
level = 4;
half_window_size=2;
[m n] = size （img1）;
G00 = img1; G10 = img2;
if （level>0）
    G01 = reduce （G00）; G11 = reduce （G10）;
end
if （level>1）
    G02 = reduce （G01）; G12 = reduce （G11）;
end
if （level>2）
    G03 = reduce （G02）; G13 = reduce （G12）;
end
if （level>3）
    G04 = reduce （G03）; G14 = reduce （G13）;
end
l = level;

for i=level:-1:0
    if （l == level）
        switch （l）
        case 4 Dx = zeros （size （G04））; Dy = zeros （size （G04））;
        case 3 Dx = zeros （size （G03））; Dy = zeros （size （G03））;
        case 2 Dx = zeros （size （G02））; Dy = zeros （size （G02））;
        case 1 Dx = zeros （size （G01））; Dy = zeros （size （G01））;
        case 0 Dx = zeros （size （G00））; Dy = zeros （size （G00））;
        end
    else
        Dx = expand （Dx）; Dy = expand （Dy）;
        Dx = Dx .* 2; Dy = Dy .* 2;
    end
    switch （l）
    case 4 
        W = warp （G04 Dx Dy）; 
        [Vx Vy] = EstimateMotion （W G14 half_window_size）;
    case 3 
        W = warp （G03 Dx Dy）; 
        [Vx Vy] = EstimateMotion （W G13 half_window_size）;
    case 2 
        W = warp （G02 Dx Dy）; 
        [Vx Vy] = EstimateMotion （W G12 half_window_size）;
    case 1 
        W = warp （G01 Dx Dy）; 
        [Vx Vy] = EstimateMotion （W G11 half_window_size）;
    case 0 
        W = warp （G00 Dx Dy）; 
        [Vx Vy] = EstimateMotion （W G10 half_window_size）;
    end
    [m n] = size （W）;
    Dx（1:m 1:n） = Dx（1:m1:n） + Vx; Dy（1:m 1:n） = Dy（1:m 1:n） + Vy;
    smooth （Dx）; smooth （Dy）;
    l = l - 1;
end


%% Lucas Kanade on the image sequence at pyramid step %%
function [Vx Vy] = EstimateMotion （W G1 half_window_size）
[m n] = size （W）;
Vx = zeros （size （W））; Vy = zeros （size （W））;
N = zeros （2*half_window_size+1 5）;
for i = 1:m
    l = 0;
    for j = 1-half_window_size:1+half_window_size
        l = l + 1;
        N （l:） = getSlice （W G1 i j half_window_size）;
    end
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     目录          0  2006-10-05 17:57  lk

     文件      64078  2003-04-08 18:03  lk\1.bmp

     文件      64078  2003-04-08 18:03  lk\2.bmp

     文件      21958  2003-04-07 14:03  lk\3.bmp

     文件      21958  2003-04-07 14:03  lk\4.bmp

     文件     169984  2004-04-18 03:57  lk\Lucas.doc

     文件       6932  2004-04-18 03:07  lk\Lucas_Kanade.m

     文件        298  2004-04-18 04:01  lk\tests.m

    ..A.SH.     20480  2006-10-05 18:07  lk\Thumbs.db

----------- ---------  ---------- -----  ----

               369984                    10