initial sd-vbs

initial sd-vbs

add sd-vbs

sd-vbs

1 files changed, 41 insertions, 0 deletions

diff --git a/SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/imStats.m b/SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/imStats.m
new file mode 100755
index 0000000..3d79b4c
--- /dev/null
+++ b/SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/imStats.m
@@ -0,0 +1,41 @@
+% imStats(IM1,IM2)
+%
+% Report image (matrix) statistics.
+% When called on a single image IM1, report min, max, mean, stdev, 
+% and kurtosis.
+% When called on two images (IM1 and IM2), report min, max, mean, 
+% stdev of the difference, and also SNR (relative to IM1).
+
+% Eero Simoncelli, 6/96.
+
+function [] = imStats(im1,im2)
+
+if (~isreal(im1))
+  error('Args must be real-valued matrices');  
+end
+
+if (exist('im2') == 1)
+  difference = im1 - im2;
+  [mn,mx] = range2(difference);
+  mean = mean2(difference);
+  v = var2(difference,mean);
+  if (v < realmin) 
+    snr = Inf;
+  else
+    snr = 10 * log10(var2(im1)/v);
+  end
+  fprintf(1, 'Difference statistics:\n');
+  fprintf(1, '  Range: [%c, %c]\n',mn,mx);
+  fprintf(1, '  Mean: %f,  Stdev (rmse): %f,  SNR (dB): %f\n',...
+      mean,sqrt(v),snr);
+else
+  [mn,mx] = range2(im1);
+  mean = mean2(im1);
+  var = var2(im1);  
+  stdev = sqrt(real(var))+sqrt(imag(var));
+  kurt = kurt2(im1, mean, stdev^2);
+  fprintf(1, 'Image statistics:\n');
+  fprintf(1, '  Range: [%f, %f]\n',mn,mx);
+  fprintf(1, '  Mean: %f,  Stdev: %f,  Kurtosis: %f\n',mean,stdev,kurt);
+end
+