Files
sf-op/script/test/CamtoolCheckProfile.py
2017-07-04 10:19:29 +02:00

32 lines
1.2 KiB
Python

import ch.psi.pshell.imaging.Data as Data
path = "2017/06/16/20170616_122305_snapshot.h5|/data"
path = "2017/06/19/20170619_113633_snapshot.h5|/data"
path = "2017/06/19/20170619_114839_snapshot.h5|/data"
path = "2017/06/19/20170619_115340_snapshot.h5|/data"
path = "2017/06/19/20170619_120203_snapshot.h5|/data"
path = "2017/06/19/20170619_141032_snapshot.h5|/data"
path = "2017/06/19/20170619_141203_snapshot.h5|/data"
path = "2017/06/19/20170619_142222_snapshot.h5|/data"
path = "2017/06/19/20170619_142310_snapshot.h5|/data"
path = "2017/06/19/20170619_142826_snapshot.h5|/data"
path = "2017/06/19/20170619_142836_snapshot.h5|/data"
path = "2017/06/19/20170619_144625_snapshot.h5|/data"
d = load_data(path)
atts = get_attributes(path)
data = Data(d)
iv = data.integrateVertically(False)
ih = data.integrateHorizontally(False)
xp = atts["x_profile"]
yp = atts["y_profile"]
p1, p2 = plot([xp, yp], ["X profile", "Y profile"])
p1.addSeries(LinePlotSeries("data X"))
p2.addSeries(LinePlotSeries("data Y"))
p1.getSeries(1).setData(iv)
p2.getSeries(1).setData(ih)
p1.addSeries(LinePlotSeries("error X", None, 2))
p2.addSeries(LinePlotSeries("error Y", None, 2))
p1.getSeries(2).setData(arrsub(xp,iv))
p2.getSeries(2).setData(arrsub(yp,ih))