/normxcorr/trunk

To get this branch, use:
bzr branch http://suren.me/webbzr/normxcorr/trunk

« back to all changes in this revision

Viewing changes to validxy_mean.m

  • Committer: Suren A. Chilingaryan
  • Date: 2010-08-17 01:41:15 UTC
  • Revision ID: csa@dside.dyndns.org-20100817014115-xbh0h6086nz5sj2o
Synchronize with Chris version, set precision to 1, initial implementation of labview wrapper

Show diffs side-by-side

added added

removed removed

Lines of Context:
validxy_mean.m
 
1
%% This function allows to take the mean of a couple of analyszed markers.
 
2
% Furthermore it allows you to take also the mean of the time. Every
 
3
% experimentalist knows the issue, vibrations. The only way of reducing the
 
4
% influence of vibrations which you have in your images is to take a mean
 
5
% value through time or images.
 
6
 
 
7
function [validx_mean,validy_mean,time_image_mean]=validxy_mean(validx,validy, time_image);
 
8
 
 
9
 
 
10
if exist('validx')==0
 
11
    [validxname,Pathvalidx] = uigetfile('*.dat','Open validx.dat');
 
12
    if validxname==0
 
13
        disp('You did not select a file!')
 
14
        return
 
15
    end
 
16
    cd(Pathvalidx);
 
17
    validx=importdata(validxname,'\t');
 
18
end
 
19
if exist('validy')==0
 
20
    [validyname,Pathvalidy] = uigetfile('*.dat','Open validy.dat');
 
21
    if validyname==0
 
22
        disp('You did not select a file!')
 
23
        return
 
24
    end
 
25
    cd(Pathvalidy);
 
26
    validy=importdata(validyname,'\t');
 
27
end
 
28
if exist('time_image')==0
 
29
    [timename,Pathtime] = uigetfile('*.txt','Open time_image.txt');
 
30
    if timename==0
 
31
        time=0;
 
32
    else
 
33
        cd(Pathtime);
 
34
        time_image=importdata(timename,'\t');
 
35
    end
 
36
end
 
37
 
 
38
 
 
39
 
 
40
[nummarker numimages]=size(validx);
 
41
numaverage=3;
 
42
for j=1:numaverage
 
43
    validxtemp(:,:,j)=validx(:,j:numaverage:(numimages-numaverage)+j);
 
44
    validytemp(:,:,j)=validy(:,j:numaverage:(numimages-numaverage)+j);
 
45
    if time~=0
 
46
        time_image_temp(:,:,j)=time_image(j:numaverage:(numimages-numaverage+1)+j,:);
 
47
    end
 
48
end
 
49
validx_mean=mean(validxtemp,3);
 
50
validy_mean=mean(validytemp,3);
 
51
save validx_mean.dat validx_mean
 
52
save validy_mean.dat validy_mean
 
53
if time~=0
 
54
    time_image_mean=mean(time_image_temp,3);
 
55
    save time_image.dat time_image
 
56
end
 
 
b'\\ No newline at end of file'