def warp1(f, g):
usmooth = 8.0
rsmooth = 101
strainMax1 = 0.125
#strainMax1 = 1.0
dstrainMax1 = 0.1
shiftMin = 0
shiftMax = 160
#shiftMax = 250
nr = rsmooth
dr = 2.0 * strainMax1 / (nr - 1)
fr = -strainMax1
sr = Sampling(nr, dr, fr)
dw = DynamicWarping(shiftMin, shiftMax)
dw.setErrorExtrapolation(DynamicWarping.ErrorExtrapolation.REFLECT)
dw.setStrainMax(strainMax1)
dw.setShiftSmoothing(usmooth, usmooth)
e1 = dw.computeErrors1(f, g)
d1 = dw.accumulateForward(e1)
u1 = dw.backtrackReverse(d1, e1)
nl = len(e1[0])
def plotShifts():
nl = len(e1[0])
sp = SimplePlot()
sp.setSize(1400, 500)
sl = Sampling(nl, s1f.delta, shiftMin * s1f.delta)
pv = sp.addPixels(s1f, sl, pow(transpose(e1), 1.0))
pv.setInterpolation(PixelsView.Interpolation.NEAREST)
pv.setColorModel(ColorMap.JET)
pv.setPercentiles(2, 98)
pv = sp.addPoints(s1f, mul(s1f.delta, u1))
plotShifts()
#u1 = dw.smoothShifts(sr,u1)
#u1 = dw.smoothShifts(u1)
n1, n2 = len(f[0]), len(f)
e1 = dw.computeErrorsRaw(f[n2 / 2], g[n2 / 2])
dw.findShiftsSmooth(dstrainMax1, u1[0], e1, u1)
plotShifts()
#u1 = dw.findShifts1(f,g)
n1, n2 = len(f[0]), len(f)
h = zerofloat(n1, n2)
u = zerofloat(n1, n2)
for i2 in range(n2):
copy(u1, u[i2])
dw.applyShifts(u1, g[i2], h[i2])
print "warp1: u min =", min(u), " max =", max(u)
return u, h
def warp1(f,g):
usmooth = 8.0
rsmooth = 101
strainMax1 = 0.125
#strainMax1 = 1.0
dstrainMax1 = 0.1
shiftMin = 0
shiftMax = 160
#shiftMax = 250
nr = rsmooth; dr = 2.0*strainMax1/(nr-1); fr = -strainMax1
sr = Sampling(nr,dr,fr)
dw = DynamicWarping(shiftMin,shiftMax)
dw.setErrorExtrapolation(DynamicWarping.ErrorExtrapolation.REFLECT)
dw.setStrainMax(strainMax1)
dw.setShiftSmoothing(usmooth,usmooth)
e1 = dw.computeErrors1(f,g)
d1 = dw.accumulateForward(e1)
u1 = dw.backtrackReverse(d1,e1)
nl = len(e1[0])
def plotShifts():
nl = len(e1[0])
sp = SimplePlot()
sp.setSize(1400,500)
sl = Sampling(nl,s1f.delta,shiftMin*s1f.delta)
pv = sp.addPixels(s1f,sl,pow(transpose(e1),1.0))
pv.setInterpolation(PixelsView.Interpolation.NEAREST)
pv.setColorModel(ColorMap.JET)
pv.setPercentiles(2,98)
pv = sp.addPoints(s1f,mul(s1f.delta,u1))
plotShifts()
#u1 = dw.smoothShifts(sr,u1)
#u1 = dw.smoothShifts(u1)
n1,n2 = len(f[0]),len(f)
e1 = dw.computeErrorsRaw(f[n2/2],g[n2/2]);
dw.findShiftsSmooth(dstrainMax1,u1[0],e1,u1)
plotShifts()
#u1 = dw.findShifts1(f,g)
n1,n2 = len(f[0]),len(f)
h = zerofloat(n1,n2)
u = zerofloat(n1,n2)
for i2 in range(n2):
copy(u1,u[i2])
dw.applyShifts(u1,g[i2],h[i2])
print "warp1: u min =",min(u)," max =",max(u)
return u,h
def goFakeErrors():
dw = DynamicWarping(shiftMin, shiftMax)
u = zerofloat(n1)
for cs in [1.0, 0.5]:
e = dw.fakeErrors(cs, nl, n1)
plotShifts(e)
for ds in [1.0, 0.1, 0.01]:
u[0] = 0.0
esum = dw.findShiftsSmooth(ds, e, u)
print "cs =", cs, " ds=", ds, " esum =", esum
plotShifts(e, u)
def goFakeErrors():
dw = DynamicWarping(shiftMin,shiftMax)
u = zerofloat(n1)
for cs in [1.0,0.5]:
e = dw.fakeErrors(cs,nl,n1)
plotShifts(e)
for ds in [1.0,0.1,0.01]:
u[0] = 0.0
esum = dw.findShiftsSmooth(ds,e,u)
print "cs =",cs," ds=",ds," esum =",esum
plotShifts(e,u)
