Python DynamicWarpingRの例

コード例 #1

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)

コード例 #2

ファイル: sdwfigures.py プロジェクト: cgraziano/idh

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)

コード例 #3

ファイル: sino2r.py プロジェクト: BKJackson/idh

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()

コード例 #4

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")

コード例 #5

ファイル: sino2s.py プロジェクト: xuewengeophysics/idh

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)

コード例 #6

ファイル: sdwfigures.py プロジェクト: cgraziano/idh

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")

コード例 #7

ファイル: sino2s.py プロジェクト: xuewengeophysics/idh

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")
    """

コード例 #8

ファイル: sino2s.py プロジェクト: BKJackson/idh

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)

コード例 #9

ファイル: sino2s.py プロジェクト: xuewengeophysics/idh

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")

コード例 #10

ファイル: sino2s.py プロジェクト: scompton/idh

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")
    """

コード例 #11

ファイル: sdwfigures.py プロジェクト: cgraziano/idh

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")

コード例 #12

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()