flow-error/compute_std.m at master · chunmingchen/flow-error · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
base_list = '/data/flow/isabel_all/all.list';
skipped_list = '/data/flow/isabel_all/all8.list';
fitted_list = '/data/flow/isabel_all/fitted_quad/all.list';
fitted_endpoint_list = '/data/flow/isabel_all/fitted_quad_endpoint/all.list';
global data_w data_d data_t data_h

base_files = load_list(base_list);
skipped_files = load_list(skipped_list);
fitted_list = load_list(fitted_list);
fitted_endpoint_list = load_list(fitted_endpoint_list);

data_t = 9;

% fitted
err = zeros(3, data_w, data_h, 1, data_t);
for t=1:data_t
   vec = load_vec(base_files{t});
   vec_fitted = load_vec(fitted_list{t});

   % trim
   vec = vec(:,:,:,50);
   vec_fitted = vec_fitted(:,:,:,50);

   err(:,:,:,:,t) = vec-vec_fitted;

end
err_std = std(err, 1, 5);
err_std_overtime = std(reshape(err(1, :,:,:,:), data_w*data_h, data_t))

img = squeeze(err_std(1,:,:));
imagesc(img');
axis image;
axis xy;
colormap(jet);


% fitted_endpoint
err = zeros(3, data_w, data_h, 1, data_t);
for t=1:data_t
   vec = load_vec(base_files{t});
   vec_fitted_endpoint = load_vec(fitted_endpoint_list{t});

   % trim
   vec = vec(:,:,:,50);
   vec_fitted_endpoint = vec_fitted_endpoint(:,:,:,50);

   err(:,:,:,:,t) = vec-vec_fitted_endpoint;

end
err_std = std(err, 1, 5);
err_std_overtime = std(reshape(err(1, :,:,:,:), data_w*data_h, data_t))

img = squeeze(err_std(1,:,:));
imagesc(img');
axis image;
axis xy;
colormap(jet);


% gen interp
vec1=load_vec(skipped_files{1});
vec2=load_vec(skipped_files{2});
for t=1:data_t
    vec = load_vec(base_files{t});
    vec_interp = vec1*(data_t-t)/(data_t-1) + vec2*(t-1)/(data_t-1);

    % trim
    vec = vec(:,:,:,50);
    vec_interp = vec_interp(:,:,:,50);

    err(:,:,:,:,t) = vec-vec_fitted;
end
err_std = std(err, 1, 5);

figure
img = squeeze(err_std(1,:,:));
imagesc(img');
axis image;
axis xy;
colormap(jet);

% gen polyfit
slice_polyfit =
for x=1:data_w
    for y=1:data_h

vec_polyfit =
for t=1:data_t
    vec = load_vec(base_files{t});