数字图像相关算法

matlab数字图像相关算法，已打包好，希望能帮到大家！！

function [validx1validy1]=automate_image_session1（grid_x1grid_y1filenamelist）;

% Code to start actual image correlation
% Programmed by Chris and Rob
% Last revision: 10/1/06

% The automation function is the central function and processes all markers and 
% images by the use of the matlab function cpcorr.m. 
% 
% automate_image_session1 and automate_image_session2 are called by the
% master function automate_image_twothreads and open grid_x1.dat or
% grid_x2.dat as well as filenamelist.m and processes the job. Since both
% threads are running independent from each other computers with dual core
% processors or two processor computers can run jobs twice as fast as
% single core/singleprocessor computers.

clear 

if exist（‘workdirectory.dat‘‘file‘）==2
load（‘workdirectory.dat‘）;
cd（char（workdirectory））;
end

% Load necessary files
if exist（‘grid_x1‘）==0
load（‘grid_x1.dat‘）              % file with x position created by grid_generator.m
end
if exist（‘grid_y1‘）==0
load（‘grid_y1.dat‘）              % file with y position created by grid_generator.m
end
if exist（‘filenamelist‘）==0
load（‘filenamelist‘）            % file with the list of filenames to be processed
end

grid_x=grid_x1;
grid_y=grid_y1;
% Initialize variables

input_points_x=grid_x;
base_points_x=grid_x;

input_points_y=grid_y;
base_points_y=grid_y;

[rowcol]=size（base_points_x）;      % this will determine the number of rasterpoints we have to run through
[rc]=size（filenamelist）;                   % this will determine the number of images we have to loop through

% Open new figure so previous ones （if open） are not overwritten

h=figure;
imshow（filenamelist（1:））           % show the first image
title（‘Initial Grid For Image Correlation （Note green crosses）‘）        % put a title
hold on
plot（grid_xgrid_y‘g+‘）            % plot the grid onto the image
hold off

% Start image correlation using cpcorr.m
g = waitbar（0sprintf（‘Processing images‘））;        % initialize the waitbar
set（g‘Position‘[2755027550]）                               % set the position of teh waitbar [left bottom width height]

for i=1:（r-1）               % run through all images
    
  
    tic             % start the timer
    base = uint8（mean（double（imread（filenamelist（1:）））3））;            % read in the base image （ which is always  image number one. You might want to change that to improve correlation results in case the light conditions are changing during the experiment
    input = uint8（mean（double（imread（filenamelist（（i+1）:）））3））;       % read in the image which has to be correlated
    
    input_points_for（:1）=reshape（input_points_x[]1）;         % we reshape the input points to one row of values since this is the shape cpcorr will accept
    input_points_for（:2）=reshape（input_points_y[]1）;
    base_points_for（:1）=reshape（base_points_x[]1）;
    base_points_for（:2）=reshape（base_points_y[]1）;
    input_correl（::

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        5343  2006-12-28 14:58  automate_image_session1.m
     文件        5628  2006-12-28 14:58  automate_image_session2.m
     文件        5296  2006-12-28 14:38  automate_image_twothreads.m
     文件        1267  2006-12-30 12:57  jobskript_twoinstances.m
     文件         170  2006-10-01 22:50  matlabsessions.au3
     文件      382850  2006-10-01 22:51  matlabsessions.exe
     文件        1308  2014-02-12 12:26  license.txt