Synopsis
		sfresamplextnd < in.rsf > out.rsf o1=o[0] d1=d[0] o2=o[1] d2=d[1] n1=(d[0]*(n[0]-1)+o[0]-o1)/d1 n2=(d[1]*(n[1]-1)+o[1]-o2)/d2

Parameters
		float d1=d[0] sample interval on 1st axis   float d2=d[1] sample interval on 2nd axis   int n1=(d[0]*(n[0]-1)+o[0]-o1)/d1 number of samples on first axis   int n2=(d[1]*(n[1]-1)+o[1]-o2)/d2 number of samples on second axis   float o1=o[0] first sample sample on 1st axis   float o2=o[1] first sample on 2nd axis