在Shape Context中画出直方图和极坐标图

Serge Belongie一干人等放出的Shape Context代码[1 ]没有显示极坐标图，和统计的二维直方图，我加了一些代码，调算法的时候看着舒服多了，也方便分析了~

 

代码如下（不太规范）：

% % sc_compute.m最后更改如下 by visionfans @ 2011.05.03 % for n=1:nsamp fzn=fz(n,:)&in_vec; Sn=sparse(theta_array_q(n,fzn),r_array_q(n,fzn),1,nbins_theta,nbins_r); BH(n,:)=Sn(:)'; scTest(); end

里面调用了一个scTest.m的脚本文件，内容如下：

% % 画极坐标图和二维直方图 by visionfans @ 2011.05.03 % %% draw polar diagram h = figure(); axs = subplot(2,1,1); axis equal; point = [Bsamp(1,n) Bsamp(2,n)]; scDrawPolar(Bsamp,point,r_inner*mean_dist,r_outer*mean_dist,nbins_theta,nbins_r); %% draw 2d histogram subplot(2,1,2); axis equal; hold on; %% temp = flipud(full(Sn)'); imagesc(temp);colormap(gray); axis image; xlabel('{/theta}'); ylabel('log(r)'); for i=1:size(temp,2) for j=1:size(temp,1) if temp(j,i) text(i,j,sprintf('%d',temp(j,i)),'Color','r'); end end end hold off; %% pause; close(h);

 

我在这里面调用了一个函数scDrawPolar来画极坐标图：

function scDrawPolar(samp,point,r_min,r_max,nbins_theta,nbins_r) %SCDRAWPOLAR draw a polar on the center point % point - the center point % r_min - min radius % r_max - max radius % nbins_theta - theta divide % nbins_r - r divide % fig_handle - draw the diagram on which figure % by visionfans @ 2011.05.03 gca; hold on; plot(samp(1,:)',samp(2,:)','r.'); plot(point(1),point(2),'ko'); r_bin_edges=logspace(log10(r_min),log10(r_max),nbins_r); % draw circles th = 0 : pi / 50 : 2 * pi; xunit = cos(th); yunit = sin(th); for i=1:length(r_bin_edges) line(xunit * r_bin_edges(i) + point(1), yunit * r_bin_edges(i) + point(2),... 'LineStyle', ':', 'Color', 'k', 'LineWidth', 1); end % draw spokes th = (1:nbins_theta) * 2*pi / nbins_theta; cs = [cos(th);zeros(1,size(th,2))]; sn = [sin(th);zeros(1,size(th,2))]; line(r_max*cs + point(1), r_max*sn + point(2),'LineStyle', ':', ... 'Color', 'k', 'LineWidth', 1); axis equal; axis off; hold off; end

 

Just thest, make it!

References:

[1]Matching with Shape Contexts，Serge Belongie, Jitendra Malik and Jan Puzicha，http://www.eecs.berkeley.edu/Research/Projects/CS/vision/shape/sc_digits.html