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;