--- a
+++ b/Thoracic Organs Segmentation code/MAP/NaiveBayesianClassifier.asv
@@ -0,0 +1,34 @@
+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);