def goShifts2():
f, g = getSinoImages()
f = copy(ni, nx, f)
g = copy(ni + nl, nx, g)
h = copy(f)
sf = Sampling(ni)
sg = Sampling(ni + nl)
dw = DynamicWarpingR(sl.first, sl.last, si, sx)
dw.setStrainLimits(0.0, 5.0, -0.2, 0.2)
dw.setSmoothness(50, 50)
u = dw.findShifts(sf, f, sg, g)
v = vpvs(u)
h = dw.applyShifts(sg, g, u)
d = sub(h, f)
print "nrms(h,f) =", nrms2(h, f)
zoom = True
plotImage(f, fmax=5, pp=True, zoom=zoom)
plotImage(h, fmax=5, pp=True, zoom=zoom, png="pw2z")
plotImage(d, fmax=5, pp=True, zoom=zoom, png="pd2z")
if zoom:
plotImage(v, fmin=1.9, fmax=2.3, pp=True, zoom=zoom)
plotImage(u, fmin=150, fmax=250, pp=True, zoom=zoom, cv=True)
else:
plotImage(v, fmin=1.5, fmax=2.5, pp=True, zoom=zoom)
plotImage(u, fmin=50, fmax=350, pp=True, zoom=zoom, cv=True)
def goShifts2():
f,g = getSinoImages()
f = copy(ni,nx,f)
g = copy(ni+nl,nx,g)
h = copy(f)
sf = Sampling(ni)
sg = Sampling(ni+nl)
dw = DynamicWarpingR(sl.first,sl.last,si,sx)
dw.setStrainLimits(0.0,5.0,-0.2,0.2)
dw.setSmoothness(50,50)
u = dw.findShifts(sf,f,sg,g)
v = vpvs(u)
h = dw.applyShifts(sg,g,u)
d = sub(h,f)
print "nrms(h,f) =",nrms2(h,f)
zoom = True
plotImage(f,fmax=5,pp=True,zoom=zoom)
plotImage(h,fmax=5,pp=True,zoom=zoom,png="pw2z")
plotImage(d,fmax=5,pp=True,zoom=zoom,png="pd2z")
if zoom:
plotImage(v,fmin=1.9,fmax=2.3,pp=True,zoom=zoom)
plotImage(u,fmin=150,fmax=250,pp=True,zoom=zoom,cv=True)
else:
plotImage(v,fmin=1.5,fmax=2.5,pp=True,zoom=zoom)
plotImage(u,fmin=50,fmax=350,pp=True,zoom=zoom,cv=True)
def warp1(f,g):
smin,smax = 0.0,3.0
dw = DynamicWarpingR(smin,smax,s1,s2)
ss = dw.getSamplingS()
ns,n1,n2 = ss.count,s1.count,s2.count
e = zerofloat(ns,n1)
for i2 in range(n2/2-11,n2/2+12):
ei = dw.computeErrors(sf,f[i2],sg,g[i2])
add(ei,e,e)
dw.normalizeErrors(e)
def plotShifts():
sp = SimplePlot()
sp.setSize(1400,500)
pv = sp.addPixels(s1,ss,pow(transpose(e),1.0))
pv.setInterpolation(PixelsView.Interpolation.NEAREST)
pv.setColorModel(ColorMap.GRAY)
pv.setPercentiles(2,98)
#pv = sp.addPoints(sf,mul(sf.delta,u1))
plotShifts()
def goTeaser():
f, g = getSinoImages()
f = copy(ni, nx, f)
g = copy(ni + nl, nx, g)
h = copy(f)
sf = Sampling(ni)
sg = Sampling(ni + nl)
dw = DynamicWarpingR(sl.first, sl.last, si, sx)
dw.setStrainLimits(0.0, 5.0, -0.2, 0.2)
dw.setSmoothness(50, 50)
u = dw.findShifts(sf, f, sg, g)
h = dw.applyShifts(sg, g, u)
plotTeaser(f, fmax=5, png="ppt")
plotTeaser(h, fmax=5, png="pst")
plotTeaser(u, fmin=150, fmax=250, cv=True, png="ut")
def goWrite():
f, g = getSinoImages()
f = copy(ni, nx, f)
g = copy(ni + nl, nx, g)
h = copy(f)
sf = Sampling(ni)
sg = Sampling(ni + nl)
dw = DynamicWarpingR(sl.first, sl.last, si, sx)
dw.setStrainLimits(0.0, 2.0, -0.1, 0.1)
dw.setSmoothness(50, 20)
u = dw.findShifts(sf, f, sg, g)
h = dw.applyShifts(sg, g, u)
print "max Tps =", max(add(rampfloat(0.0, 1.0, 0.0, ni, nx), u))
writeImage("pp.dat", f)
writeImage("ps.dat", g)
writeImage("pswarped.dat", h)
writeImage("shifts.dat", u)
def goTeaser(): f,g = getSinoImages() f = copy(ni,nx,f) g = copy(ni+nl,nx,g) h = copy(f) sf = Sampling(ni) sg = Sampling(ni+nl) dw = DynamicWarpingR(sl.first,sl.last,si,sx) dw.setStrainLimits(0.0,5.0,-0.2,0.2) dw.setSmoothness(50,50) u = dw.findShifts(sf,f,sg,g) h = dw.applyShifts(sg,g,u) plotTeaser(f,fmax=5,png="ppt") plotTeaser(h,fmax=5,png="pst") plotTeaser(u,fmin=150,fmax=250,cv=True,png="ut")
def goShifts2():
f, g = getSinoImages()
f = copy(ni, nx, f)
g = copy(ni + nl, nx, g)
h = copy(f)
sf = Sampling(ni)
sg = Sampling(ni + nl)
dw = DynamicWarpingR(sl.first, sl.last, si, sx)
dw.setStrainLimits(0.0, 2.0, -0.1, 0.1)
dw.setSmoothness(50, 20)
u = dw.findShifts(sf, f, sg, g)
v = vpvs(u)
h = dw.applyShifts(sg, g, u)
d = sub(h, f)
print "nrms(h,f) =", nrms2(h, f)
zoom = True
#plotImage(f,fmax=5,pp=True,zoom=zoom)
plotImage(h, fmax=5, pp=True, zoom=zoom, png="pw2z")
plotImage(d, fmax=5, pp=True, zoom=zoom, png="pd2z")
"""
def goWrite():
f,g = getSinoImages()
f = copy(ni,nx,f)
g = copy(ni+nl,nx,g)
h = copy(f)
sf = Sampling(ni)
sg = Sampling(ni+nl)
dw = DynamicWarpingR(sl.first,sl.last,si,sx)
dw.setStrainLimits(0.0,2.0,-0.1,0.1)
dw.setSmoothness(50,20)
u = dw.findShifts(sf,f,sg,g)
h = dw.applyShifts(sg,g,u)
print "max Tps =",max(add(rampfloat(0.0,1.0,0.0,ni,nx),u))
writeImage("pp.dat",f)
writeImage("ps.dat",g)
writeImage("pswarped.dat",h)
writeImage("shifts.dat",u)
def goShifts1():
f, g = getSinoImages()
f = copy(ni, nx, f)
g = copy(ni + nl, nx, g)
h = copy(f)
e = computeErrors(nx, f, g)
sf = Sampling(ni)
sg = Sampling(ni + nl)
ix = nx / 2
dw = DynamicWarpingR(sl.first, sl.last, si)
dw.setStrainLimits(0.0, 2.0)
plotImage(f, fmax=5, pp=True, zoom=True)
for hs in [50]:
dw.setSmoothness(hs)
u = dw.findShifts(e)
for ix in range(nx):
h[ix] = dw.applyShifts(sg, g[ix], u)
d = sub(h, f)
print "nrms(h,f) =", nrms2(h, f)
plotImage(h, fmax=5, pp=True, zoom=True, png="pw1z")
plotImage(d, fmax=5, pp=True, zoom=True, png="pd1z")
def goShifts2():
f, g = getSinoImages()
f = copy(ni, nx, f)
g = copy(ni + nl, nx, g)
h = copy(f)
sf = Sampling(ni)
sg = Sampling(ni + nl)
dw = DynamicWarpingR(sl.first, sl.last, si, sx)
dw.setStrainLimits(0.0, 2.0, -0.1, 0.1)
dw.setSmoothness(50, 20)
u = dw.findShifts(sf, f, sg, g)
v = vpvs(u)
h = dw.applyShifts(sg, g, u)
d = sub(h, f)
print "nrms(h,f) =", nrms2(h, f)
zoom = True
# plotImage(f,fmax=5,pp=True,zoom=zoom)
plotImage(h, fmax=5, pp=True, zoom=zoom, png="pw2z")
plotImage(d, fmax=5, pp=True, zoom=zoom, png="pd2z")
"""
def goShifts1():
f,g = getSinoImages()
f = copy(ni,nx,f)
g = copy(ni+nl,nx,g)
h = copy(f)
e = computeErrors(nx,f,g)
sf = Sampling(ni)
sg = Sampling(ni+nl)
ix = nx/2
dw = DynamicWarpingR(sl.first,sl.last,si)
dw.setStrainLimits(0.0,2.0)
plotImage(f,fmax=5,pp=True,zoom=True)
for hs in [50]:
dw.setSmoothness(hs)
u = dw.findShifts(e)
for ix in range(nx):
h[ix] = dw.applyShifts(sg,g[ix],u)
d = sub(h,f)
print "nrms(h,f) =",nrms2(h,f)
plotImage(h,fmax=5,pp=True,zoom=True,png="pw1z")
plotImage(d,fmax=5,pp=True,zoom=True,png="pd1z")
def warp1(f, g):
smin, smax = 0.0, 3.0
dw = DynamicWarpingR(smin, smax, s1, s2)
ss = dw.getSamplingS()
ns, n1, n2 = ss.count, s1.count, s2.count
e = zerofloat(ns, n1)
for i2 in range(n2 / 2 - 11, n2 / 2 + 12):
ei = dw.computeErrors(sf, f[i2], sg, g[i2])
add(ei, e, e)
dw.normalizeErrors(e)
def plotShifts():
sp = SimplePlot()
sp.setSize(1400, 500)
pv = sp.addPixels(s1, ss, pow(transpose(e), 1.0))
pv.setInterpolation(PixelsView.Interpolation.NEAREST)
pv.setColorModel(ColorMap.GRAY)
pv.setPercentiles(2, 98)
#pv = sp.addPoints(sf,mul(sf.delta,u1))
plotShifts()
