test_features.m

   1 tic
   2 [Visuals_train Boundary1_train Boundary2_train Filenames_train]=read_data('/home/tsix/tsix/high_jump/09sem/gmm/OCT/Lernstichprobe/');
   3
   4 [Visuals_test Boundary1_test Boundary2_test Filenames_test]=read_data('/home/tsix/tsix/high_jump/09sem/gmm/OCT/ttt1/');
   5 [s1, S1, s2, S2] = get_probabilities(Visuals_train, Boundary1_train,Boundary2_train);
   6 sigma1=10;
   7 sigma2=40;
   8
   9 %for i=1:length(Visuals_test)
  10 %for i=1:9
  11  for i=[34 24 9 8]
  12     tic
  13   f1=extract_features1(Visuals_test{i}, 3);
  14   %toc
  15   %tic
  16
  17   X=get_emission_probabilities(f1, s1,sigma1);
  18   %toc
  19   %tic
  20   fb_mean1=hmm_forward_backward(s1, S1, X,'mean');
  21   %toc
  22   %tic
  23
  24   f2=extract_features2(Visuals_test{i}, fb_mean1,s2,sigma2);
  25   %toc
  26   %tic
  27   s=compose_priors(s2,fb_mean1);
  28   %toc
  29   X=get_emission_probabilities(f2, s,sigma2);
  30   %toc
  31   %tic
  32   fb_mean2=hmm_forward_backward(s(:,1), S2(1:size(s,1),1:size(s,1)), X,'mean');
  33   toc
  34
  35   figure(i);
  36   show_result(Visuals_test{i}, fb_mean1,Boundary1_test{i},fb_mean2,Boundary2_test{i});
  37
  38   %figure(i+3);
  39  %show_features(Visuals_test{i},f1,f2);
  40   [mse1 bias1 var1 count1 mse2 bias2 var2 count2]=validate(fb_mean1,Boundary1_test{i},fb_mean2,Boundary2_test{i});
  41   fprintf('Results for %s:\n',Filenames_test{i});
  42   fprintf('  Boundary1:\n');
  43   fprintf('    Bias=%f, MSE=%f, var=%f, size of ground truth: %d\n',bias1, mse1, var1, count1);
  44   fprintf('  Boundary2:\n');
  45   fprintf('    Bias=%f, MSE=%f, var=%f, size of ground truth: %d\n',bias2, mse2, var2, count2);
  46   fprintf('\n');
  47
  48 end
  49
  50
  51