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) ;