vistfromfwd.m

来自「Matlab數位影音處理 影像方面的處理」· M 代码 · 共 35 行

function vistfromfwd(tform,wdate,zdate,N)
if nargin < 4
    n =10;
end
[w,z] = meshgrid(linspace(wdate(1),zdate(2),N),...
                 linspace(wdate(1),zdate(2),N));
 wz = [w(:) z(:)];
 xy = tformfwd([w(:) z(:)], tform);
 x = reshape(xy(:,1),size(w));
 y = reshape(xy(:,2),size(z));
 wx = [w(:); x(:)];
 wxlimits = [min(wx) max(wx)];
 zy = [z(:); y(:) ];
 zylimits = [min(zy) max(zy)];
 figure;
 subplot(1,2,1)
 plot(w, z, 'b' ),axis equal , axis ij
 hold on
 plot(w',z','b')
 hold off
 xlim(wxlimits)
 ylim(zylimits)
 
 set(gca, 'XAxisLocation','top')
 xlabel('w'),ylabel('z')
 subplot(1,2,2)
 plot(x,y,'b'),axis equal, axis ij
 hold on
 plot(x',y','b')
 hold off
 xlim(wxlimits)
 ylim(zylimits)
 set(gca, 'XAxisLocation','top')
 xlabel('x'),ylabel('y')