xlabel('row')
#####################################################
# linear shift fit (consider factor and background)
#####################################################
from pydao.math import XyDiscreteFun
#print img0.nthmoment(1);
#print img1.nthmoment(1);
#print new_mimg.nthmoment(1);
N, M = data0.shape
x = arange(M)
y0 = data0.sum(0)
y1 = data1.sum(0)
sp0 = XyDiscreteFun()
sp0.set('x', x)
sp0.set('y', y0)
sp1 = XyDiscreteFun(x)
sp1.set('x', x)
sp1.set('y', y1)
paras = sp1.find_xshift_renormbg(sp0)
delta = paras[0]
factor = paras[1]
bg = paras[2]
print "shift in columns", paras
x1fit = x + delta
y1fit = y0 * factor + bg
subplot(2, 3, 4)
plot(x1fit, y1fit, linewidth=4)
from pylab import *;
from pydao.math import XyDiscreteFun;
data = loadtxt(r'C:\Home\Dropbox\Programming\data_analysis\2014_11_10_find_spect_shift\S003_004_1.4mj_t2t+2.txt');
x = data[:,0];
y1 = data[:,1];
y2 = data[:,2];
sp1 = XyDiscreteFun();
sp1.set('x',x);
sp1.set('y',y1);
sp1.pick(xmin=50,xmax=130);
sp2 = XyDiscreteFun();
sp2.set('x',x);
sp2.set('y',y2);
sp2.pick(xmin=50,xmax=130);
subplot(1,2,1);
plot(x,y1)
plot(x,y2);
paras = sp1.find_xshift_renorm(sp2);
delta = paras[0];
factor = paras[1];
print paras;
x2fit = x+delta;
y2fit = y2*factor;
