def goSino():
na, ka = 81, -20 # sampling for inverse wavelet A
nh, kh = 181, -90 # sampling for wavelet H
nt, dt, ft = 501, 0.004, 0.000 # used for plotting only
nx, dx, fx = 721, 0.015, 0.000
f, g, u = getSinoImages()
st = Sampling(nt, dt, ft)
sx = Sampling(nx, dx, fx)
tmin, tmax = 88, 363
fmin, fmax = 0.0, 0.5
sfac = 1.000
ww = WaveletWarping()
ww.setFrequencyRange(fmin, fmax)
ww.setTimeRange(tmin, tmax)
ww.setStabilityFactor(sfac)
aw = ww.getInverseA(na, ka, u, f, g) # estimated inverse wavelet
hw = ww.getWaveletH(na, ka, aw, nh, kh) # estimated wavelet
nhw = normalize(hw)
af = ww.applyA(na, ka, aw, f)
ag = ww.applyA(na, ka, aw, g)
slag = ww.applyS(u, ww.applyL(u, ag))
slg = ww.applyS(u, ww.applyL(u, g))
cmin, cmax = -0.5, 0.5
#plotImage(st,sx,haslg,fmin=cmin,fmax=cmax,zoom=True)
#plotImage(st,sx,hslag,fmin=cmin,fmax=cmax,zoom=True)
plotImage(st, sx, f, fmin=cmin, fmax=cmax, zoom=True)
plotImage(st, sx, slg, fmin=cmin, fmax=cmax, zoom=True)
plotImage(st, sx, af, fmin=cmin, fmax=cmax, zoom=True)
plotImage(st, sx, slag, fmin=cmin, fmax=cmax, zoom=True)
f0 = ww.rms(f)
f1 = ww.rms(af)
g0 = ww.rms(slg)
g1 = ww.rms(slag)
e0 = ww.rms(sub(f, slg))
e1 = ww.rms(sub(af, slag))
print "f0 =", f0, " f1 =", f1
print "g0 =", g0, " g1 =", g1
print "e0 =", e0, " e1 =", e1
def goSino(): na,ka = 81,-20 # sampling for inverse wavelet A nh,kh = 181,-90 # sampling for wavelet H nt,dt,ft = 501,0.004,0.000 # used for plotting only nx,dx,fx = 721,0.015,0.000 f,g,u = getSinoImages() st = Sampling(nt,dt,ft) sx = Sampling(nx,dx,fx) tmin,tmax = 88,363 fmin,fmax = 0.0,0.5 sfac = 1.000 ww = WaveletWarping() ww.setFrequencyRange(fmin,fmax) ww.setTimeRange(tmin,tmax) ww.setStabilityFactor(sfac) aw = ww.getInverseA(na,ka,u,f,g) # estimated inverse wavelet hw = ww.getWaveletH(na,ka,aw,nh,kh) # estimated wavelet nhw = normalize(hw) af = ww.applyA(na,ka,aw,f) ag = ww.applyA(na,ka,aw,g) slag = ww.applyS(u,ww.applyL(u,ag)) slg = ww.applyS(u,ww.applyL(u,g)) cmin,cmax = -0.5,0.5 #plotImage(st,sx,haslg,fmin=cmin,fmax=cmax,zoom=True) #plotImage(st,sx,hslag,fmin=cmin,fmax=cmax,zoom=True) plotImage(st,sx,f,fmin=cmin,fmax=cmax,zoom=True) plotImage(st,sx,slg,fmin=cmin,fmax=cmax,zoom=True) plotImage(st,sx,af,fmin=cmin,fmax=cmax,zoom=True) plotImage(st,sx,slag,fmin=cmin,fmax=cmax,zoom=True) f0 = ww.rms(f) f1 = ww.rms(af) g0 = ww.rms(slg) g1 = ww.rms(slag) e0 = ww.rms(sub(f,slg)) e1 = ww.rms(sub(af,slag)) print "f0 =",f0," f1 =",f1 print "g0 =",g0," g1 =",g1 print "e0 =",e0," e1 =",e1
