1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
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
|
