1 files changed, 24 insertions, 0 deletions

diff --git a/SD-VBS/benchmarks/texture_synthesis/src/matlab/kurt2.m b/SD-VBS/benchmarks/texture_synthesis/src/matlab/kurt2.m
new file mode 100755
index 0000000..8351eb4
--- /dev/null
+++ b/SD-VBS/benchmarks/texture_synthesis/src/matlab/kurt2.m
@@ -0,0 +1,24 @@
+% K = KURT2(MTX,MEAN,VAR)
+%
+% Sample kurtosis (fourth moment divided by squared variance) 
+% of a matrix.  Kurtosis of a Gaussian distribution is 3.
+%  MEAN (optional) and VAR (optional) make the computation faster.
+
+% Eero Simoncelli, 6/96.
+
+function res = kurt2(mtx, mn, v)
+
+if (exist('mn') ~= 1)
+        mn =  mean(mean(mtx));
+end
+
+if (exist('v') ~= 1)
+        v =  var2(mtx,mn);
+end
+
+if (isreal(mtx))
+  res = mean(mean(abs(mtx-mn).^4)) / (v^2);
+else
+  res = mean(mean(real(mtx-mn).^4)) / (real(v)^2) + ...
+      i*mean(mean(imag(mtx-mn).^4)) / (imag(v)^2);
+end