Python XyDiscreteFun 예제들

예제 #1

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)

예제 #2

파일: imgshift_test.py 프로젝트: XuGroup/Pydao

#####################################################
# 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);

예제 #3

파일: fin_spect_shift.py 프로젝트: caoshi/Pydao

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

예제 #4

파일: fin_spect_shift.py 프로젝트: caoshi/Pydao

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;