code/get_emission_probabilities.m

   1 % get_emission_probabilities - computes p(x,s) for a single image
   2 %
   3 % usage
   4 %        X = get_emission_probabilities(features, s,sigma)
   5 % input
   6 %       features - a matrix of size(width(image),2*?) containing feature
   7 %       vectors for each column of a given image
   8 %       s - prior probabilities p(s)
   9 %       sigma - variance of a single part of the underlying distribution
  10 %
  11 % output
  12 %       X - emission probabilities (rows correspond to
  13 %       observations at given times, columns to hidden states)
  14 function [X] = get_emission_probabilities(features, s,sigma)
  15 if(size(s,2)==1); s=repmat(s,1,size(features,1)); end
  16 X = zeros(size(features,1),size(s,1));
  17
  18 for i=1:size(X,1)
  19
  20   for j=1:size(X,2)
  21     if(s(j,i)==0)
  22       X(i,j)=0;
  23     else
  24       X(i,j)=get_prob(j,features(i,:),sigma);
  25     end
  26
  27   end
  28 end
  29
  30
  31