WIP: compute centroid of CT and Expanded brain mask

expandBrainMask.m

@@ -9,4 +9,14 @@ brainMaskImg = Brain_mask.img;
 structEl = strel('sphere', 3);
 expndBrainMask = imdilate(brainMaskImg, structEl);
 
-isosurface(expndBrainMask);
\ No newline at end of file
+isosurface(expndBrainMask);
+
+% compute the centroid of the CT to later map the electrode with the
+% expanded brain Mask and get rid of the artefact far away from brain
+stats_CT_post = regionprops(maskedImg);
+centroid_CT_post = stats_CT_post.Centroid;
+
+% compute the centroid of the expanded brain mask  
+stats_expndBrainMask = regionprops(expndBrainMask);
+centroid_expndBrainMask = stats_expndBrainMask.Centroid; 
+