bzr branch
http://darksoft.org/webbzr/normxcorr/trunk
41
by Suren A. Chilingaryan
Synchronize with Chris version, set precision to 1, initial implementation of labview wrapper |
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
|