initial sd-vbs

initial sd-vbs

add sd-vbs

sd-vbs

14 files changed, 790 insertions, 0 deletions

diff --git a/SD-VBS/common/toolbox/toolbox_basic/stella/afromncut.m b/SD-VBS/common/toolbox/toolbox_basic/stella/afromncut.m
new file mode 100755
index 0000000..ec014d0
--- /dev/null
+++ b/SD-VBS/common/toolbox/toolbox_basic/stella/afromncut.m
@@ -0,0 +1,73 @@
+% function a = afromncut(v,s,d,visimg,no_rep,pixel_loc)
+% Input:
+%    v = eigenvectors of d*a*d, starting from the second.
+%      (the first is all one over some constant determined by d)
+%    s = eigenvalues
+%    d = normalization matrix 1/sqrt(rowsum(abs(a)))
+%    visimg = 1/0 if each eigenvector is/not 2D (so v is 3D)
+%    no_rep = 1 (default), affinity has attraction only
+%        if 1, the first column of v is the second eigenvector
+%        if 0, the first column of v is the first eigenvector.
+%    pixel_loc = nx1 matrix, each is a pixel location
+% Output:
+%    a = diag(1/d) * na * diag(1/d);
+%    If pixel_loc = []; a is returned, if not out of memory
+%    otherwise, only rows of a at pixel_loc are returned.
+%
+% This routine is used to estimate the original affinity matrix
+% through the first few eigenvectors and its normalization matrix.
+
+% A test sequence includes:
+% a = randsym(5);
+% [na,d] = normalize(a);
+% [v,s] = ncut(a,5);
+% v = v(:,2:end); s = s(2:end);
+% aa = afromncut(v,s,d);
+% max(abs(aa(:) - a(:)))
+
+% Stella X. Yu, 2000.
+
+function a = afromncut(v,s,d,visimg,no_rep,pixel_loc)
+
+[nr,nc,nv] = size(v);
+if nargin<4 | isempty(visimg),
+   visimg = (nv>1);
+end
+
+if nargin<5 | isempty(no_rep),
+   no_rep = 1;
+end
+
+if visimg,
+   nr = nr * nc;
+else
+   nv = nc;
+end
+
+if nargin<6 | isempty(pixel_loc),
+   pixel_loc = 1:nr;
+end
+
+% D^(1/2)
+d = 1./(d(:)+eps);
+
+% first recover the first eigenvector
+if no_rep,
+    u = (1/norm(d)) + zeros(nr,1);
+    s = [1;s(:)];
+    nv = nv + 1;
+else
+    u = [];
+end
+
+% the full set of generalized eigenvectors
+v = [u, reshape(v,[nr,nv-no_rep])];
+
+% This is the real D, row sum
+d = d.^2;
+
+% an equivalent way to compute v = diag(d) * v;
+v = v .* d(:,ones(nv,1)); % to avoid using a big matrix diag(d)
+
+% synthesis
+a = v(pixel_loc,:)*diag(s)*v';
diff --git a/SD-VBS/common/toolbox/toolbox_basic/stella/dispimg.m b/SD-VBS/common/toolbox/toolbox_basic/stella/dispimg.m
new file mode 100755
index 0000000..4e419a0
--- /dev/null
+++ b/SD-VBS/common/toolbox/toolbox_basic/stella/dispimg.m
@@ -0,0 +1,65 @@
+% function dispimg(g,fmt,lgd,cmap) display multiple images in one figure. 
+% Input:
+%    g = a cell and fmt is a 1x2 vector specifying the layout.
+%    lgd = a string cell for the title of each image. 
+%    cmap = the colormap (default is the gray, -1 for the inverted gray).
+%    ishori = a vector of 1/0 to display real and imag parts horizontally / vertically
+
+% Stella X. Yu, 2000.
+
+function dispimg(g,fmt,lgd,cmap,ishori);
+
+cellg = iscell(g);
+if cellg,
+    num_fig = length(g);
+else
+    num_fig = size(g,3);
+end;
+
+if nargin<2 | isempty(fmt),
+    m = ceil(sqrt(num_fig));
+    n = ceil(num_fig / m);
+else
+    m = fmt(1);
+    n = fmt(2);
+end
+
+if nargin<3 | isempty(lgd),
+    lgd = 1:num_fig;
+end 
+if isnumeric(lgd),
+   lgd = cellstr(num2str(lgd(:),3)); 
+end
+i = size(lgd);
+if i(1)==1,
+   lgd = [lgd, cell(1,num_fig-i(2))];
+else
+   lgd = [lgd; cell(num_fig-i(1),1)];
+end
+
+if nargin<5 | isempty(ishori),
+   ishori = ones(num_fig,1);
+end
+ishori(end+1:num_fig) = ishori(end);
+
+for k=1:num_fig,
+    subplot(m,n,k);
+    if cellg,
+        showim(g{k},[],ishori(k));
+    else
+        showim(g(:,:,k),[],ishori(k));
+    end
+    title(lgd{k});
+end
+
+if nargin<4 | isempty(cmap),
+   cmap = gray;
+end
+if length(cmap)==1,
+   if cmap==1,
+      cmap = gray;
+   else
+      cmap = flipud(gray);
+   end
+end
+colormap(cmap);
diff --git a/SD-VBS/common/toolbox/toolbox_basic/stella/firstncut.m b/SD-VBS/common/toolbox/toolbox_basic/stella/firstncut.m
new file mode 100755
index 0000000..a22077d
--- /dev/null
+++ b/SD-VBS/common/toolbox/toolbox_basic/stella/firstncut.m
@@ -0,0 +1,67 @@
+% function [v,s,d] = firstncut(base_name,rec_num)
+% Input:
+%    base_name = image name
+%    rec_num = parameter record number 
+% Output:
+%    v = eigenvectors
+%    s = eigenvalues
+%    d = normalization matrix d = 1/sqrt(rowsum(abs(a)))
+% Convert Jianbo Shi's Ncut Ccode results from images to matlab matrices.
+
+% Stella X. Yu, 2000.
+
+function [v,s,d] = firstncut(base_name,rec_num);
+
+if nargin<2 | isempty(rec_num),
+   rec_num = 1;
+end
+
+cur_dir = pwd;
+globalenvar;
+cd(IMAGE_DIR);
+cd(base_name);
+   feval([base_name,'_par']);
+   j = length(p);
+   if rec_num>j,