vis3d/tign2ros.m

   1 function r=tign2ros(p,i)
   2 dxlon=(p.fxlong(3:end,2:end-1,i)-p.fxlong(1:end-2,2:end-1,i));
   3 dxlat=(p.fxlat(3:end,2:end-1,i)-p.fxlat(1:end-2,2:end-1,i));
   4 dx=sqrt(dxlon.^2+dxlat.^2);
   5 gx=(p.tign_g(3:end,2:end-1,i)-p.tign_g(1:end-2,2:end-1,i))./dx;
   6 dylon=(p.fxlong(2:end-1,3:end,i)-p.fxlong(2:end-1,1:end-2,i));
   7 dylat=(p.fxlat(2:end-1,3:end,i)-p.fxlat(2:end-1,1:end-2,i));
   8 dy=sqrt(dylon.^2+dylat.^2);
   9 gy=(p.tign_g(2:end-1,3:end,i)-p.tign_g(2:end-1,1:end-2,i))./dy;
  10 gn=sqrt(gx.^2+gy.^2);
  11 r = 1./gn;
  12 end