Exemplo n.º 1
0
def run_script(args):

    
    matplotlib.interactive(False)
        
    left = (args.left[0], args.left[1])
    right = (args.right[0], args.right[1])


    l1 = (150,110,110)
    l2 = (110,150,110)
    l3 = (110,110,150)
    layers= [l1,l2]

    # This is a hack to conserve colors
    colourmap = { rgb(l1[0],l1[1],l1[2]): args.Rock0,
                  rgb(l2[0],l2[1],l2[2]): args.Rock1 }
    
    if not isinstance(args.Rock2, str):
        colourmap[rgb( l3[0],l3[1],l3[2])] = args.Rock2
        layers.append( l3 )
    
    model = mb.body( traces = args.ntraces,
                     pad = args.pad,
                     margin=args.margin,
                     left = left,
                     right = right,
                     layers = layers
                   )

    return modelr_plot(model, colourmap, args)
Exemplo n.º 2
0
def run_script(args):

    matplotlib.interactive(False)

    left = (args.left[0], args.left[1])
    right = (args.right[0], args.right[1])

    l1 = (150, 110, 110)
    l2 = (110, 150, 110)
    l3 = (110, 110, 150)
    layers = [l1, l2]

    # This is a hack to conserve colors
    colourmap = {
        rgb(l1[0], l1[1], l1[2]): args.Rock0,
        rgb(l2[0], l2[1], l2[2]): args.Rock1
    }

    if not isinstance(args.Rock2, str):
        colourmap[rgb(l3[0], l3[1], l3[2])] = args.Rock2
        layers.append(l3)

    model = mb.body(traces=args.ntraces,
                    pad=args.pad,
                    margin=args.margin,
                    left=left,
                    right=right,
                    layers=layers)

    return modelr_plot(model, colourmap, args)
Exemplo n.º 3
0
def run_script(args):
    
    matplotlib.interactive(False)
    
    """if args.transparent == 'False' or args.transparent == 'No':
        transparent = False
    else:
        transparent = True"""

    args.ntraces = 300
    args.pad = 150
    args.reflectivity_method = zoeppritz
    args.title = 'Wedge - AVA gather'
    args.theta = (0.0,50,.5)
    args.wavelet = ricker
    args.wiggle_skips = 10
    args.aspect_ratio = 1
    args.left = (0,0)
    args.right = (0,50)
    args.margin=1 
    args.slice='angle'
   
    
    transparent = False
    # This is a hack to conserve colors
    l1 = (150,110,110)
    l2 = (110,150,110)
    l3 = (110,110,150)
    layers= [l1,l2]
    colourmap = { rgb(150,110,110): args.Rock0,
                  rgb(110,150,110): args.Rock1 }
    
    if not isinstance(args.Rock2, str):
        colourmap[rgb( 110,110,150)] = args.Rock2
        layers.append( l3 )
    # Get the physical model (an array of rocks)    
    model = mb.body( traces = args.ntraces,
                     pad = args.pad,
                     margin=args.margin,
                     left = args.left,
                     right = args.right,
                     layers = layers
                   )

    
    return modelr_plot( model, colourmap, args )
Exemplo n.º 4
0
def run_script(args):
    from modelr.constants import dt, duration
    
    matplotlib.interactive(False)

    args.ntraces = 300
    args.pad = 150
    args.reflectivity_method = zoeppritz
    args.title = 'Wedge Model - RGB Blended Cross Section'
    #args.colourmap = 'Greys' #might need to overwrite this (it may not be used)
    args.wavelet = ricker
    args.wiggle_skips = 10
    args.aspect_ratio = 1
    args.margin=1
    args.left = (0,0)
    args.right = (0,50)
    args.slice='spatial'
    args.trace = 0
    
    left = (args.left[0], args.left[1])
    right = (args.right[0], args.right[1])


    l1 = (150,110,110)
    l2 = (110,150,110)
    l3 = (110,110,150)
    layers= [l1,l2]

    # This is a hack to conserve colors
    colourmap = { rgb(l1[0],l1[1],l1[2]): args.Rock0,
                  rgb(l2[0],l2[1],l2[2]): args.Rock1 }
    
    if not isinstance(args.Rock2, str):
        colourmap[rgb( l3[0],l3[1],l3[2])] = args.Rock2
        layers.append( l3 )
    
    model = mb.body( traces = args.ntraces,
                     pad = args.pad,
                     margin=args.margin,
                     left = left,
                     right = right,
                     layers = layers
                   )

    return modelr_plot( model, colourmap, args )
Exemplo n.º 5
0
def run_script(args):
    from modelr.constants import dt, duration

    matplotlib.interactive(False)

    args.ntraces = 300
    args.pad = 150
    args.reflectivity_method = zoeppritz
    args.title = 'Wedge Model - RGB Blended Cross Section'
    #args.colourmap = 'Greys' #might need to overwrite this (it may not be used)
    args.wavelet = ricker
    args.wiggle_skips = 10
    args.aspect_ratio = 1
    args.margin = 1
    args.left = (0, 0)
    args.right = (0, 50)
    args.slice = 'spatial'
    args.trace = 0

    left = (args.left[0], args.left[1])
    right = (args.right[0], args.right[1])

    l1 = (150, 110, 110)
    l2 = (110, 150, 110)
    l3 = (110, 110, 150)
    layers = [l1, l2]

    # This is a hack to conserve colors
    colourmap = {
        rgb(l1[0], l1[1], l1[2]): args.Rock0,
        rgb(l2[0], l2[1], l2[2]): args.Rock1
    }

    if not isinstance(args.Rock2, str):
        colourmap[rgb(l3[0], l3[1], l3[2])] = args.Rock2
        layers.append(l3)

    model = mb.body(traces=args.ntraces,
                    pad=args.pad,
                    margin=args.margin,
                    left=left,
                    right=right,
                    layers=layers)

    return modelr_plot(model, colourmap, args)
Exemplo n.º 6
0
def run_script(args):

    matplotlib.interactive(False)
    """if args.transparent == 'False' or args.transparent == 'No':
        transparent = False
    else:
        transparent = True"""

    args.ntraces = 300
    args.pad = 150
    args.reflectivity_method = zoeppritz
    args.title = 'Wedge - AVA gather'
    args.theta = (0.0, 50, .5)
    args.wavelet = ricker
    args.wiggle_skips = 10
    args.aspect_ratio = 1
    args.left = (0, 0)
    args.right = (0, 50)
    args.margin = 1
    args.slice = 'angle'

    # This is a hack to conserve colors
    l1 = (150, 110, 110)
    l2 = (110, 150, 110)
    l3 = (110, 110, 150)
    layers = [l1, l2]
    colourmap = {
        rgb(150, 110, 110): args.Rock0,
        rgb(110, 150, 110): args.Rock1
    }

    if not isinstance(args.Rock2, str):
        colourmap[rgb(110, 110, 150)] = args.Rock2
        layers.append(l3)

    # Get the physical model (an array of rocks)
    model = mb.body(traces=args.ntraces,
                    pad=args.pad,
                    margin=args.margin,
                    left=args.left,
                    right=args.right,
                    layers=layers)

    return modelr_plot(model, colourmap, args)