From 770e6624b950d949baf3e28ad5c743f30fcede29 Mon Sep 17 00:00:00 2001
From: Florian Vollrath <florian.vollrath@brain.mpg.de>
Date: Wed, 21 Sep 2016 15:43:50 +0200
Subject: [PATCH] added the centerofgravity function

---
 centerofgravity.m | 21 +++++++++++++++++++++
 1 file changed, 21 insertions(+)
 create mode 100644 centerofgravity.m

diff --git a/centerofgravity.m b/centerofgravity.m
new file mode 100644
index 0000000..036c90b
--- /dev/null
+++ b/centerofgravity.m
@@ -0,0 +1,21 @@
+function cog = centerofgravity(d,thr)
+%Determine center of mass of a matrix
+% cog = centerofgravity(d,thr)
+
+if nargin == 2
+    d(d<thr) = 0;
+end
+d=double(d);
+
+
+
+for dim = 1:numel(size(d))
+    dd = d;
+    ind = 1:numel(size(d));
+    ind(dim) = [];
+    dd = permute(dd,[dim ind]);
+    dd = reshape(dd,size(d,dim),[]);
+    dsum=nansum(dd,2);
+    distance = 1:size(d,dim);
+    cog(dim) = nansum(dsum(:).*distance(:)) / nansum(dsum);
+end
\ No newline at end of file