// Sorts input date/time wave by the time of day. NaNs (blanks) go to the end.
// NOTE: This overwrites the input wave.
Function SortByTimeOfDay(w)
	Wave w					// Date/time wave to be sorted
 
	Duplicate/FREE w, timeOfDay
	timeOfDay = mod(w, 24*60*60)		// timeOfDay now contains just time of day
	Sort {timeOfDay,w}, w			// Sort w by timeOfDay, then by date
End

Function test()
	make/n=(100,100, 100)/free w
	w = gauss(x, 10, 3, y, 20, 3, z, 30, 3)
	wave temp = com3D(w)
	print temp
End
 
Function/wave com3D(w)
	Wave w
	variable ii, t1, t2, t3
	make/n=(wavedims(w))/free/d COM
	switch(wavedims(w))
		case 1:
			COM[0] = Pla_peakcentroid(w)
			break
		case 2:
			imagetransform sumallrows w
			COM[0] = Pla_peakcentroid(W_sumrows)
			imagetransform sumallcols w
			COM[1] = Pla_peakcentroid(W_sumcols)
			break
		case 3:
			imagetransform sumplanes w
			Wave temp = com3D(M_SumPlanes)
			COM[0] = temp[0]
			COM[1] = temp[1]

// Interpolate to create a path with equal spacing between successive points in the XY plane.
// the path described by xwave and ywave as inputs needs to be quite smooth for this to work.
// play with step size (step), number of nodes for interpolation (numNodes) and linear vs. cubic spline 
// interpolation (flagT = 1, 2) to optimise the path output as xout and yout.
// step size must be rather small with respect to spacing of path points in xwave, ywave
 
function StepInterpolate(xwave, ywave, step, flagT, numNodes)
	wave xwave, ywave
	variable step, flagT, numNodes