From f618466c25d43f3bae9e40920273bf77de1e1149 Mon Sep 17 00:00:00 2001
From: leochanj105 <leochanj@live.unc.edu>
Date: Mon, 19 Oct 2020 23:09:30 -0400
Subject: initial sd-vbs

initial sd-vbs

add sd-vbs

sd-vbs
---
 .../texture_synthesis/src/matlab/mkAngularSine.m   | 42 ++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100755 SD-VBS/benchmarks/texture_synthesis/src/matlab/mkAngularSine.m

(limited to 'SD-VBS/benchmarks/texture_synthesis/src/matlab/mkAngularSine.m')

diff --git a/SD-VBS/benchmarks/texture_synthesis/src/matlab/mkAngularSine.m b/SD-VBS/benchmarks/texture_synthesis/src/matlab/mkAngularSine.m
new file mode 100755
index 0000000..f5238cc
--- /dev/null
+++ b/SD-VBS/benchmarks/texture_synthesis/src/matlab/mkAngularSine.m
@@ -0,0 +1,42 @@
+% IM = mkAngularSine(SIZE, HARMONIC, AMPL, PHASE, ORIGIN)
+%
+% Make an angular sinusoidal image:
+%     AMPL * sin( HARMONIC*theta + PHASE),
+% where theta is the angle about the origin.
+% SIZE specifies the matrix size, as for zeros().  
+% AMPL (default = 1) and PHASE (default = 0) are optional.
+
+% Eero Simoncelli, 2/97.
+
+function [res] = mkAngularSine(sz, harmonic, ampl, ph, origin)
+
+sz = sz(:);
+if (size(sz,1) == 1)
+  sz = [sz,sz];
+end
+
+mxsz = max(sz(1),sz(2));
+
+%------------------------------------------------------------
+%% OPTIONAL ARGS:
+
+if (exist('harmonic') ~= 1)
+  harmonic = 1;
+end
+
+if (exist('ampl') ~= 1)
+  ampl = 1;
+end
+
+if (exist('ph') ~= 1)
+  ph = 0;
+end
+
+if (exist('origin') ~= 1)
+  origin = (sz+1)/2;
+end
+
+%------------------------------------------------------------
+
+res = ampl * sin(harmonic*mkAngle(sz,ph,origin) + ph);
+
-- 
cgit v1.2.2