function [Pmax,ClassNum]= NaiveBayesianClassif(dataPoints,Intensity, Displacements);

%load mean and std values of features extracted by hand

load MeanValues.mat

load StdValues.mat

%the first column of the .mat files corresponds to the manually segmented heart and the

%second column to background values .

%Intensity 

Intensity=double(meshgrid(Intensity,1:2));

MeanValue=meshgrid(MeanValues(1,1:2),1:size(dataPoints,1))';

StdValue=meshgrid(StdValues(1,1:2),1:size(dataPoints,1))';

A1=exp(-(Intensity-MeanValue').^2./(2*StdValue'.^2));

clear Intensity StdValue MeanValue

%Displacement x

DisplacementX=meshgrid(Displacements(:,1),1:2);

MeanValue=meshgrid(MeanValues(2,1:2),1:size(dataPoints,1));

StdValue=meshgrid(StdValues(2,1:2),1:size(dataPoints,1));

A2=exp(-(DisplacementX-MeanValue').^2./(2*StdValue'.^2));

clear DisplacementX StdValue MeanValue

%Displacement y

DisplacementY=meshgrid(Displacements(:,2),1:2);

MeanValue=meshgrid(MeanValues(3,1:3),1:size(dataPoints,1));

StdValue=meshgrid(StdValues(3,1:3),1:size(dataPoints,1));

A3=exp(-(DisplacementY-MeanValue').^2./(2*StdValue'.^2));

clear DisplacementY StdValue MeanValue

P=A1.*A2.*A3.*A4;

clear A1 A2 A3 A4 

[Pmax,ClassNum]=max(P,[],1);