initial sd-vbs

initial sd-vbs

add sd-vbs

sd-vbs

1 files changed, 38 insertions, 0 deletions

diff --git a/SD-VBS/common/toolbox/toolbox_basic/stella/test_ncutm.m b/SD-VBS/common/toolbox/toolbox_basic/stella/test_ncutm.m
new file mode 100755
index 0000000..c9b46ab
--- /dev/null
+++ b/SD-VBS/common/toolbox/toolbox_basic/stella/test_ncutm.m
@@ -0,0 +1,38 @@
+fn = 'walk1';
+
+repulsion_test = 1;
+
+if 1,
+   f = getimage2(fn);
+   par = jshincutdefpar;
+   par.fname_base = fn;
+   par.spatial_neighborhood_x = 10;
+   par.sigma_x = 3 * par.spatial_neighborhood_x;
+   par.sig_IC = 0.03;
+   par.num_eigvecs = 10;
+   par.offset = 0.00;
+   par.sig_filter = 1.0;
+   par.elong_filter = 3.0;
+   [par,rec_num] = jshincut(par);
+   [g,lgd,v,s,dd] = showncut(fn,rec_num);
+
+if repulsion_test,
+   par.offset = 0.1;
+   [par,rec_num] = jshincut(par);
+   figure;
+   [g,lgd,v,s,dd] = showncut(fn,rec_num);
+end
+end
+
+if 0,
+   x = v(:,:,1);
+   y = v(:,:,2);
+   figure;
+   subplot(2,1,1); plot(x(:),y(:),'ro');
+   r = sqrt(x.^2+y.^2);
+   x = x./r;
+   y = y./r;
+   subplot(2,1,2); im([x,y]*[x,y]');
+ %  mask = (x>0) & y>0;
+ %  showmask(f,mask);
+end