initial sd-vbs

initial sd-vbs

add sd-vbs

sd-vbs

1 files changed, 33 insertions, 0 deletions

diff --git a/SD-VBS/benchmarks/texture_synthesis/src/matlab/mkZonePlate.m b/SD-VBS/benchmarks/texture_synthesis/src/matlab/mkZonePlate.m
new file mode 100755
index 0000000..2b6a4a5
--- /dev/null
+++ b/SD-VBS/benchmarks/texture_synthesis/src/matlab/mkZonePlate.m
@@ -0,0 +1,33 @@
+% IM = mkZonePlate(SIZE, AMPL, PHASE)
+%
+% Make a "zone plate" image:
+%     AMPL * cos( r^2 + PHASE)
+% SIZE specifies the matrix size, as for zeros().  
+% AMPL (default = 1) and PHASE (default = 0) are optional.
+
+% Eero Simoncelli, 6/96.
+
+function [res] = mkZonePlate(sz, ampl, ph)
+
+sz = sz(:);
+if (size(sz,1) == 1)
+  sz = [sz,sz];
+end
+
+mxsz = max(sz(1),sz(2));
+
+%------------------------------------------------------------
+%% OPTIONAL ARGS:
+
+if (exist('ampl') ~= 1)
+  ampl = 1;
+end
+
+if (exist('ph') ~= 1)
+  ph = 0;
+end
+
+%------------------------------------------------------------
+
+res = ampl * cos( (pi/mxsz) * mkR(sz,2) + ph );
+