www.gusucode.com > 无线传感器网络的能量控制分簇发送数据包 > 图像配准源码程序/code/demons3d/showgrid.m
%% Display vector field % Changed: Dec 11th, 2011 % function showgrid(ux,uy,uz,downsample,lim) if nargin<3; downsample = 1; end; sizex = size(ux,1); sizey = size(ux,2); sizez = size(ux,3); % Use mid slice as image slice = ceil(sizez/2); ux = ux(:,:,slice); uy = uy(:,:,slice); ux = ux(1:downsample:end, 1:downsample:end); uy = uy(1:downsample:end, 1:downsample:end); if nargin<5; scale = 3; end; % Scale vector to show small ones if nargin<6; lim = [0 sizex-1 0 sizey-1 0 sizez-1]; end; % Display whole image [y,x] = ndgrid((1:downsample:sizex)+downsample/2, (1:downsample:sizey)+downsample/2); % coordinate image z = zeros(size(x)); mesh(x+ux,y+uy,z); view(2); daspect([1 1 1]); axis([lim(3) lim(4) lim(1) lim(2)] + .5 + [downsample 0 downsample 0]/2); % which vector to show set(gca,'YDir','reverse'); end