import ch.psi.pshell.imaging.Data as Data

path = "2017_06/20170619/20170619_153743_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))