1 I = vl_impattern('roofs1') ;
2 image(I) ;
3 %vl_sift函数的输入是一个单精度的灰度图像，灰度值区间归一化到[0, 255]。
4 %因此图像 I 需要通过下面的函数转成相应的格式：
5 I = single(rgb2gray(I)) ;
6 %vl_sift( )的输出矩阵 这里f是一个[x,y,s,th]的向量矩阵，
7 %其中一列的前两个x,y表示兴趣点的中心位置，s表示兴趣点的尺度大小，th是其梯度方向。
8 %d表示128维的特征向量
9 [f,d] = vl_sift(I) ;
10 perm = randperm(size(f,2)) ;
11 sel = perm(1:50) ;
12 h1 = vl_plotframe(f(:,sel)) ;
13 h2 = vl_plotframe(f(:,sel)) ;
14 set(h1,'color','k','linewidth',3) ;
15 set(h2,'color','y','linewidth',2) ;