// usage : ImageJ -macro svc.txt TA.tif Dz

A=split(getArgument," ,");
if (A.length==2) {
  open(A[0]); Nv=getWidth; Ny=getHeight; oID=getImageID;
  Dz=parseFloat(A[1]); Nz=Ny+abs(Dz);
  newImage("Dz = "+A[1],"32-bit",Nv,Nz,1); nID=getImageID;
  setBatchMode(true);
  Dz_Nv1=Dz/(Nv-1); if (Dz<0) y0=Dz+0.5; else y0=0.5;
  for (z=0; z<Nz; z++) {
    showProgress(z/Nz);
    for (v=0; v<Nv; v++) {
      y=y0-Dz_Nv1*v;
      if (y<0 || y>=Ny)
        TA=0;
      else {
        selectImage(oID); TA=getPixel(v,floor(y));
      }
      selectImage(nID); setPixel(v,z,TA);
    }
    y0++;
  }
  setBatchMode(false);
  resetMinAndMax();
}