In [2]:

    

%%file SConstruct
from rsf.proj import *
from rsf.recipes.helderiv import Helderiv

# Download data

horizon = 'Penobscot_HorB.txt'

Fetch(horizon,'data',
      server='https://raw.github.com',
      top='seg/tutorials/master')

# Convert from xyz triples to horizon

Flow('xyz',horizon,
     '''
     echo n1=3 n2=254674 data_format=ascii_float in=$SOURCE |
     dd form=native
     ''')

Flow('xy','xyz','window n1=2 | dd type=int')
Flow('horizon','xyz xy',
     '''
     window f1=2 squeeze=n | 
     intbin head=${SOURCES[1]} xkey=0 ykey=1 |
     window | put label=Time unit=ms label1=Inline label2=Crossline
     ''')

Result('horizon','grey color=j bias=750 scalebar=y title=Horizon yreverse=n transp=n')

# Interpolate missing values

Flow('filled','horizon','lapfill grad=y verb=y niter=500')

# Display with the cube helix colormap and the correct aspect ratio

Plot('cubeh','filled',
     '''
     grey color=x bias=950 clip=133 scalebar=y title=Horizon 
     yreverse=n transp=n screenwd=11.9 screenht=9.62
     ''')

Result('cubeh','Overlay')

# Take a window

Flow('window','horizon','window f1=200 n1=125 f2=200 n2=100')

Result('window',
       'grey color=x bias=950 clip=133 scalebar=y title=Horizon yreverse=n transp=n screenratio=0.8')

Flow('wind.asc',None,
     '''
     echo 
     200 200 
     200 300
     325 300
     325 200
     200 200
     n1=2 n2=5 data_format=ascii_int in=$TARGET
     ''')
Plot('wind','wind.asc',
     '''
     dd type=complex form=native | window |
     graph wanttitle=n wantaxis=n plotcol=7 plotfat=10 
     min1=1 max1=595 min2=1 max2=463 scalebar=y screenwd=11.9 screenht=9.62
     ''')

Result('wind','cubeh wind','Overlay')

# Add noise
Flow('noisy','filled','noise type=n seed=2014 range=15')

# Regular smoothing
Flow('smoothed','noisy','smooth rect1=4 rect2=4')

# Anisotropic diffusion
Flow('diffused','noisy','impl2 rect1=10 rect2=10 tau=1')

# Median filter
Flow('median','noisy','despike2 wide1=10 wide2=10')

# Bilateral filter
Flow('bilateral','noisy','bilat2 a1=0.1 a3=0.001 r1=10 r2=10')

for case in Split('filled noisy smoothed diffused median bilateral'):
    Result(case,
           '''
           grey color=j bias=950 clip=133 scalebar=y 
           title="%s Horizon" yreverse=n transp=n
           minval=800 maxval=1100
           ''' % case.capitalize())
    if case != 'filled':
        Result(case+'-slice',['filled',case],
               '''
               cat axis=3 ${SOURCES[1]} |
               window n1=1 f1=300 n2=450 |
               graph plotfat=7,1 wanttitle=n 
               yreverse=y label2=Time unit2=ms
               ''')

# Shaded relief
Flow('dx','filled','deriv')
Flow('dy','filled','transp | deriv | transp')
Flow('shaded','dx dy','add ${SOURCES[1]}')

# Helix derivative
Helderiv('hel',na=10)

Flow('helix','filled hel','helicon filt=${SOURCES[1]}')

for case in Split('shaded helix'):
      Result(case,
           '''
           grey title="%s Horizon" yreverse=n transp=n mean=y
           ''' % case.capitalize())

End()


    

    




Writing SConstruct

In [3]:

    

!scons -Q


    

    




retrieve(["Penobscot_HorB.txt"], [])
< Penobscot_HorB.txt echo n1=3 n2=254674 data_format=ascii_float in=Penobscot_HorB.txt | /home/sergey/RSFROOT/bin/sfdd form=native > xyz.rsf
< xyz.rsf /home/sergey/RSFROOT/bin/sfwindow n1=2 | /home/sergey/RSFROOT/bin/sfdd type=int > xy.rsf
< xyz.rsf /home/sergey/RSFROOT/bin/sfwindow f1=2 squeeze=n | /home/sergey/RSFROOT/bin/sfintbin head=xy.rsf xkey=0 ykey=1 | /home/sergey/RSFROOT/bin/sfwindow | /home/sergey/RSFROOT/bin/sfput label=Time unit=ms label1=Inline label2=Crossline > horizon.rsf
sfintbin: 463 of 463;
< horizon.rsf /home/sergey/RSFROOT/bin/sfgrey color=j bias=750 scalebar=y title=Horizon yreverse=n transp=n > Fig/horizon.vpl
< horizon.rsf /home/sergey/RSFROOT/bin/sflapfill grad=y verb=y niter=500 > filled.rsf
sflapfill: iteration 1 res 1.000000 mod 479436.968750 grad 1.000000
sflapfill: iteration 2 res 0.624057 mod 480367.031250 grad 0.554027
sflapfill: iteration 3 res 0.440545 mod 481365.687500 grad 0.352207
sflapfill: iteration 4 res 0.334763 mod 482159.656250 grad 0.239700
sflapfill: iteration 5 res 0.267765 mod 482771.312500 grad 0.173057
sflapfill: iteration 6 res 0.224722 mod 483225.500000 grad 0.130411
sflapfill: iteration 7 res 0.194830 mod 483583.031250 grad 0.100369
sflapfill: iteration 8 res 0.174097 mod 483861.937500 grad 0.080679
sflapfill: iteration 9 res 0.159915 mod 484075.218750 grad 0.063897
sflapfill: iteration 10 res 0.149630 mod 484248.000000 grad 0.053979
sflapfill: iteration 11 res 0.141611 mod 484399.062500 grad 0.047501
sflapfill: iteration 12 res 0.135059 mod 484537.812500 grad 0.041849
sflapfill: iteration 13 res 0.129297 mod 484675.281250 grad 0.039744
sflapfill: iteration 14 res 0.124048 mod 484816.562500 grad 0.036404
sflapfill: iteration 15 res 0.119261 mod 484961.500000 grad 0.034843
sflapfill: iteration 16 res 0.114710 mod 485116.000000 grad 0.033544
sflapfill: iteration 17 res 0.110409 mod 485278.843750 grad 0.031172
sflapfill: iteration 18 res 0.106556 mod 485440.062500 grad 0.028819
sflapfill: iteration 19 res 0.103127 mod 485597.062500 grad 0.027489
sflapfill: iteration 20 res 0.099868 mod 485759.343750 grad 0.026215
sflapfill: iteration 21 res 0.096664 mod 485932.468750 grad 0.025557
sflapfill: iteration 22 res 0.093778 mod 486100.781250 grad 0.023943
sflapfill: iteration 23 res 0.091089 mod 486268.906250 grad 0.022714
sflapfill: iteration 24 res 0.088578 mod 486436.500000 grad 0.021249
sflapfill: iteration 25 res 0.086274 mod 486599.937500 grad 0.020493
sflapfill: iteration 26 res 0.084198 mod 486755.750000 grad 0.018734
sflapfill: iteration 27 res 0.082366 mod 486900.562500 grad 0.018119
sflapfill: iteration 28 res 0.080568 mod 487049.937500 grad 0.017541
sflapfill: iteration 29 res 0.078784 mod 487206.000000 grad 0.017894
sflapfill: iteration 30 res 0.076916 mod 487378.093750 grad 0.017437
sflapfill: iteration 31 res 0.075157 mod 487548.718750 grad 0.016762
sflapfill: iteration 32 res 0.073549 mod 487712.218750 grad 0.015385
sflapfill: iteration 33 res 0.072123 mod 487863.531250 grad 0.014581
sflapfill: iteration 34 res 0.070845 mod 488004.562500 grad 0.013793
sflapfill: iteration 35 res 0.069603 mod 488146.906250 grad 0.013938
sflapfill: iteration 36 res 0.068344 mod 488296.875000 grad 0.013452
sflapfill: iteration 37 res 0.067131 mod 488446.968750 grad 0.013548
sflapfill: iteration 38 res 0.065887 mod 488607.000000 grad 0.013198
sflapfill: iteration 39 res 0.064647 mod 488772.906250 grad 0.013364
sflapfill: iteration 40 res 0.063393 mod 488947.187500 grad 0.012952
sflapfill: iteration 41 res 0.062178 mod 489122.625000 grad 0.012979
sflapfill: iteration 42 res 0.060934 mod 489308.875000 grad 0.012829
sflapfill: iteration 43 res 0.059740 mod 489494.000000 grad 0.012036
sflapfill: iteration 44 res 0.058703 mod 489660.031250 grad 0.011235
sflapfill: iteration 45 res 0.057706 mod 489824.125000 grad 0.011290
sflapfill: iteration 46 res 0.056722 mod 489990.875000 grad 0.010977
sflapfill: iteration 47 res 0.055727 mod 490163.968750 grad 0.010989
sflapfill: iteration 48 res 0.054776 mod 490334.031250 grad 0.010386
sflapfill: iteration 49 res 0.053867 mod 490500.750000 grad 0.010443
sflapfill: iteration 50 res 0.052939 mod 490675.406250 grad 0.010280
sflapfill: iteration 51 res 0.051992 mod 490858.093750 grad 0.010518
sflapfill: iteration 52 res 0.051000 mod 491054.312500 grad 0.010579
sflapfill: iteration 53 res 0.049969 mod 491263.656250 grad 0.010685
sflapfill: iteration 54 res 0.048930 mod 491480.031250 grad 0.010385
sflapfill: iteration 55 res 0.047930 mod 491693.062500 grad 0.010012
sflapfill: iteration 56 res 0.047015 mod 491892.093750 grad 0.009563
sflapfill: iteration 57 res 0.046134 mod 492087.406250 grad 0.009295
sflapfill: iteration 58 res 0.045305 mod 492274.406250 grad 0.009024
sflapfill: iteration 59 res 0.044465 mod 492466.875000 grad 0.009050
sflapfill: iteration 60 res 0.043618 mod 492664.375000 grad 0.008913
sflapfill: iteration 61 res 0.042766 mod 492865.937500 grad 0.008923
sflapfill: iteration 62 res 0.041943 mod 493063.562500 grad 0.008628
sflapfill: iteration 63 res 0.041104 mod 493267.843750 grad 0.008793
sflapfill: iteration 64 res 0.040258 mod 493476.500000 grad 0.008486
sflapfill: iteration 65 res 0.039447 mod 493679.218750 grad 0.008274
sflapfill: iteration 66 res 0.038658 mod 493878.500000 grad 0.008117
sflapfill: iteration 67 res 0.037898 mod 494072.375000 grad 0.007801
sflapfill: iteration 68 res 0.037184 mod 494256.375000 grad 0.007445
sflapfill: iteration 69 res 0.036524 mod 494427.718750 grad 0.007212
sflapfill: iteration 70 res 0.035870 mod 494598.781250 grad 0.007138
sflapfill: iteration 71 res 0.035218 mod 494770.625000 grad 0.007075
sflapfill: iteration 72 res 0.034583 mod 494939.125000 grad 0.006746
sflapfill: iteration 73 res 0.033999 mod 495095.187500 grad 0.006403
sflapfill: iteration 74 res 0.033481 mod 495234.218750 grad 0.005906
sflapfill: iteration 75 res 0.033013 mod 495360.656250 grad 0.005823
sflapfill: iteration 76 res 0.032553 mod 495485.468750 grad 0.005620
sflapfill: iteration 77 res 0.032114 mod 495605.093750 grad 0.005501
sflapfill: iteration 78 res 0.031699 mod 495718.875000 grad 0.005253
sflapfill: iteration 79 res 0.031307 mod 495826.937500 grad 0.005133
sflapfill: iteration 80 res 0.030919 mod 495934.437500 grad 0.005043
sflapfill: iteration 81 res 0.030547 mod 496038.093750 grad 0.005053
sflapfill: iteration 82 res 0.030167 mod 496144.718750 grad 0.004910
sflapfill: iteration 83 res 0.029806 mod 496246.406250 grad 0.004845
sflapfill: iteration 84 res 0.029460 mod 496344.218750 grad 0.004609
sflapfill: iteration 85 res 0.029120 mod 496441.281250 grad 0.004639
sflapfill: iteration 86 res 0.028805 mod 496531.531250 grad 0.004278
sflapfill: iteration 87 res 0.028523 mod 496612.812500 grad 0.004155
sflapfill: iteration 88 res 0.028258 mod 496689.875000 grad 0.003971
sflapfill: iteration 89 res 0.027993 mod 496767.093750 grad 0.004050
sflapfill: iteration 90 res 0.027719 mod 496847.375000 grad 0.004040
sflapfill: iteration 91 res 0.027457 mod 496924.968750 grad 0.003968
sflapfill: iteration 92 res 0.027208 mod 496998.750000 grad 0.003683
sflapfill: iteration 93 res 0.026981 mod 497066.687500 grad 0.003636
sflapfill: iteration 94 res 0.026760 mod 497133.375000 grad 0.003584
sflapfill: iteration 95 res 0.026541 mod 497199.625000 grad 0.003520
sflapfill: iteration 96 res 0.026337 mod 497261.968750 grad 0.003325
sflapfill: iteration 97 res 0.026142 mod 497321.843750 grad 0.003301
sflapfill: iteration 98 res 0.025954 mod 497380.218750 grad 0.003309
sflapfill: iteration 99 res 0.025754 mod 497442.500000 grad 0.003401
sflapfill: iteration 100 res 0.025550 mod 497506.718750 grad 0.003338
sflapfill: iteration 101 res 0.025352 mod 497569.406250 grad 0.003311
sflapfill: iteration 102 res 0.025162 mod 497629.843750 grad 0.003150
sflapfill: iteration 103 res 0.024986 mod 497686.593750 grad 0.003064
sflapfill: iteration 104 res 0.024814 mod 497742.593750 grad 0.003045
sflapfill: iteration 105 res 0.024646 mod 497797.406250 grad 0.002941
sflapfill: iteration 106 res 0.024488 mod 497849.468750 grad 0.002885
sflapfill: iteration 107 res 0.024329 mod 497902.406250 grad 0.002909
sflapfill: iteration 108 res 0.024176 mod 497953.718750 grad 0.002822
sflapfill: iteration 109 res 0.024023 mod 498005.250000 grad 0.002805
sflapfill: iteration 110 res 0.023875 mod 498055.812500 grad 0.002751
sflapfill: iteration 111 res 0.023727 mod 498106.562500 grad 0.002762
sflapfill: iteration 112 res 0.023583 mod 498156.531250 grad 0.002707
sflapfill: iteration 113 res 0.023437 mod 498207.531250 grad 0.002852
sflapfill: iteration 114 res 0.023277 mod 498264.187500 grad 0.002772
sflapfill: iteration 115 res 0.023127 mod 498317.781250 grad 0.002720
sflapfill: iteration 116 res 0.022984 mod 498369.062500 grad 0.002564
sflapfill: iteration 117 res 0.022859 mod 498414.625000 grad 0.002459
sflapfill: iteration 118 res 0.022743 mod 498457.281250 grad 0.002373
sflapfill: iteration 119 res 0.022631 mod 498498.343750 grad 0.002317
sflapfill: iteration 120 res 0.022526 mod 498537.468750 grad 0.002275
sflapfill: iteration 121 res 0.022419 mod 498577.593750 grad 0.002289
sflapfill: iteration 122 res 0.022316 mod 498616.937500 grad 0.002176
sflapfill: iteration 123 res 0.022221 mod 498653.187500 grad 0.002109
sflapfill: iteration 124 res 0.022133 mod 498686.750000 grad 0.002015
sflapfill: iteration 125 res 0.022052 mod 498718.187500 grad 0.001931
sflapfill: iteration 126 res 0.021980 mod 498746.093750 grad 0.001762
sflapfill: iteration 127 res 0.021918 mod 498770.625000 grad 0.001731
sflapfill: iteration 128 res 0.021857 mod 498794.687500 grad 0.001680
sflapfill: iteration 129 res 0.021799 mod 498817.906250 grad 0.001657
sflapfill: iteration 130 res 0.021742 mod 498840.750000 grad 0.001620
sflapfill: iteration 131 res 0.021689 mod 498861.937500 grad 0.001534
sflapfill: iteration 132 res 0.021642 mod 498881.156250 grad 0.001459
sflapfill: iteration 133 res 0.021597 mod 498899.593750 grad 0.001438
sflapfill: iteration 134 res 0.021554 mod 498916.968750 grad 0.001396
sflapfill: iteration 135 res 0.021516 mod 498932.781250 grad 0.001280
sflapfill: iteration 136 res 0.021484 mod 498946.062500 grad 0.001162
sflapfill: iteration 137 res 0.021456 mod 498957.843750 grad 0.001153
sflapfill: iteration 138 res 0.021428 mod 498969.718750 grad 0.001144
sflapfill: iteration 139 res 0.021399 mod 498981.593750 grad 0.001129
sflapfill: iteration 140 res 0.021373 mod 498992.875000 grad 0.001104
sflapfill: iteration 141 res 0.021347 mod 499004.062500 grad 0.001095
sflapfill: iteration 142 res 0.021323 mod 499014.375000 grad 0.001002
sflapfill: iteration 143 res 0.021302 mod 499023.187500 grad 0.000956
sflapfill: iteration 144 res 0.021284 mod 499031.187500 grad 0.000889
sflapfill: iteration 145 res 0.021267 mod 499038.437500 grad 0.000853
sflapfill: iteration 146 res 0.021252 mod 499045.062500 grad 0.000811
sflapfill: iteration 147 res 0.021239 mod 499051.031250 grad 0.000783
sflapfill: iteration 148 res 0.021226 mod 499056.593750 grad 0.000738
sflapfill: iteration 149 res 0.021215 mod 499061.656250 grad 0.000712
sflapfill: iteration 150 res 0.021204 mod 499066.468750 grad 0.000704
sflapfill: iteration 151 res 0.021193 mod 499071.406250 grad 0.000710
sflapfill: iteration 152 res 0.021183 mod 499076.125000 grad 0.000681
sflapfill: iteration 153 res 0.021173 mod 499080.625000 grad 0.000667
sflapfill: iteration 154 res 0.021163 mod 499085.031250 grad 0.000655
sflapfill: iteration 155 res 0.021154 mod 499089.156250 grad 0.000637
sflapfill: iteration 156 res 0.021146 mod 499093.156250 grad 0.000618
sflapfill: iteration 157 res 0.021138 mod 499096.906250 grad 0.000607
sflapfill: iteration 158 res 0.021130 mod 499100.593750 grad 0.000599
sflapfill: iteration 159 res 0.021122 mod 499104.250000 grad 0.000587
sflapfill: iteration 160 res 0.021115 mod 499107.687500 grad 0.000559
sflapfill: iteration 161 res 0.021108 mod 499110.812500 grad 0.000537
sflapfill: iteration 162 res 0.021103 mod 499113.562500 grad 0.000499
sflapfill: iteration 163 res 0.021097 mod 499116.156250 grad 0.000494
sflapfill: iteration 164 res 0.021092 mod 499118.687500 grad 0.000488
sflapfill: iteration 165 res 0.021087 mod 499121.187500 grad 0.000474
sflapfill: iteration 166 res 0.021082 mod 499123.468750 grad 0.000443
sflapfill: iteration 167 res 0.021078 mod 499125.500000 grad 0.000422
sflapfill: iteration 168 res 0.021075 mod 499127.375000 grad 0.000419
sflapfill: iteration 169 res 0.021071 mod 499129.250000 grad 0.000415
sflapfill: iteration 170 res 0.021067 mod 499131.156250 grad 0.000411
sflapfill: iteration 171 res 0.021063 mod 499133.031250 grad 0.000416
sflapfill: iteration 172 res 0.021060 mod 499134.906250 grad 0.000411
sflapfill: iteration 173 res 0.021056 mod 499136.812500 grad 0.000403
sflapfill: iteration 174 res 0.021053 mod 499138.531250 grad 0.000375
sflapfill: iteration 175 res 0.021050 mod 499140.000000 grad 0.000352
sflapfill: iteration 176 res 0.021047 mod 499141.343750 grad 0.000338
sflapfill: iteration 177 res 0.021045 mod 499142.656250 grad 0.000343
sflapfill: iteration 178 res 0.021042 mod 499144.031250 grad 0.000348
sflapfill: iteration 179 res 0.021040 mod 499145.375000 grad 0.000334
sflapfill: iteration 180 res 0.021037 mod 499146.656250 grad 0.000318
sflapfill: iteration 181 res 0.021035 mod 499147.781250 grad 0.000307
sflapfill: iteration 182 res 0.021033 mod 499148.843750 grad 0.000297
sflapfill: iteration 183 res 0.021031 mod 499149.875000 grad 0.000280
sflapfill: iteration 184 res 0.021030 mod 499150.781250 grad 0.000271
sflapfill: iteration 185 res 0.021028 mod 499151.656250 grad 0.000275
sflapfill: iteration 186 res 0.021026 mod 499152.562500 grad 0.000276
sflapfill: iteration 187 res 0.021025 mod 499153.500000 grad 0.000275
sflapfill: iteration 188 res 0.021023 mod 499154.406250 grad 0.000264
sflapfill: iteration 189 res 0.021022 mod 499155.218750 grad 0.000253
sflapfill: iteration 190 res 0.021020 mod 499155.968750 grad 0.000240
sflapfill: iteration 191 res 0.021019 mod 499156.687500 grad 0.000235
sflapfill: iteration 192 res 0.021018 mod 499157.343750 grad 0.000229
sflapfill: iteration 193 res 0.021017 mod 499158.031250 grad 0.000232
sflapfill: iteration 194 res 0.021016 mod 499158.718750 grad 0.000230
sflapfill: iteration 195 res 0.021014 mod 499159.375000 grad 0.000229
sflapfill: iteration 196 res 0.021013 mod 499160.031250 grad 0.000227
sflapfill: iteration 197 res 0.021012 mod 499160.718750 grad 0.000226
sflapfill: iteration 198 res 0.021011 mod 499161.343750 grad 0.000216
sflapfill: iteration 199 res 0.021010 mod 499161.968750 grad 0.000213
sflapfill: iteration 200 res 0.021009 mod 499162.531250 grad 0.000200
sflapfill: iteration 201 res 0.021008 mod 499163.031250 grad 0.000189
sflapfill: iteration 202 res 0.021008 mod 499163.468750 grad 0.000180
sflapfill: iteration 203 res 0.021007 mod 499163.906250 grad 0.000184
sflapfill: iteration 204 res 0.021006 mod 499164.406250 grad 0.000185
sflapfill: iteration 205 res 0.021005 mod 499164.843750 grad 0.000175
sflapfill: iteration 206 res 0.021005 mod 499165.218750 grad 0.000166
sflapfill: iteration 207 res 0.021004 mod 499165.593750 grad 0.000164
sflapfill: iteration 208 res 0.021004 mod 499165.968750 grad 0.000162
sflapfill: iteration 209 res 0.021003 mod 499166.343750 grad 0.000159
sflapfill: iteration 210 res 0.021002 mod 499166.687500 grad 0.000156
sflapfill: iteration 211 res 0.021002 mod 499167.031250 grad 0.000164
sflapfill: iteration 212 res 0.021001 mod 499167.437500 grad 0.000167
sflapfill: iteration 213 res 0.021001 mod 499167.843750 grad 0.000173
sflapfill: iteration 214 res 0.021000 mod 499168.281250 grad 0.000169
sflapfill: iteration 215 res 0.020999 mod 499168.687500 grad 0.000165
sflapfill: iteration 216 res 0.020999 mod 499169.062500 grad 0.000159
sflapfill: iteration 217 res 0.020998 mod 499169.437500 grad 0.000158
sflapfill: iteration 218 res 0.020998 mod 499169.781250 grad 0.000153
sflapfill: iteration 219 res 0.020997 mod 499170.125000 grad 0.000147
sflapfill: iteration 220 res 0.020997 mod 499170.437500 grad 0.000140
sflapfill: iteration 221 res 0.020996 mod 499170.718750 grad 0.000137
sflapfill: iteration 222 res 0.020996 mod 499171.000000 grad 0.000133
sflapfill: iteration 223 res 0.020996 mod 499171.281250 grad 0.000130
sflapfill: iteration 224 res 0.020995 mod 499171.531250 grad 0.000127
sflapfill: iteration 225 res 0.020995 mod 499171.781250 grad 0.000122
sflapfill: iteration 226 res 0.020995 mod 499172.000000 grad 0.000118
sflapfill: iteration 227 res 0.020994 mod 499172.218750 grad 0.000114
sflapfill: iteration 228 res 0.020994 mod 499172.406250 grad 0.000110
sflapfill: iteration 229 res 0.020994 mod 499172.593750 grad 0.000108
sflapfill: iteration 230 res 0.020993 mod 499172.781250 grad 0.000107
sflapfill: iteration 231 res 0.020993 mod 499172.968750 grad 0.000106
sflapfill: iteration 232 res 0.020993 mod 499173.125000 grad 0.000103
sflapfill: iteration 233 res 0.020993 mod 499173.312500 grad 0.000101
sflapfill: iteration 234 res 0.020993 mod 499173.468750 grad 0.000099
sflapfill: iteration 235 res 0.020992 mod 499173.625000 grad 0.000094
sflapfill: iteration 236 res 0.020992 mod 499173.750000 grad 0.000088
sflapfill: iteration 237 res 0.020992 mod 499173.875000 grad 0.000088
sflapfill: iteration 238 res 0.020992 mod 499174.031250 grad 0.000087
sflapfill: iteration 239 res 0.020992 mod 499174.156250 grad 0.000084
sflapfill: iteration 240 res 0.020992 mod 499174.250000 grad 0.000076
sflapfill: iteration 241 res 0.020991 mod 499174.343750 grad 0.000067
sflapfill: iteration 242 res 0.020991 mod 499174.406250 grad 0.000062
sflapfill: iteration 243 res 0.020991 mod 499174.468750 grad 0.000061
sflapfill: iteration 244 res 0.020991 mod 499174.531250 grad 0.000061
sflapfill: iteration 245 res 0.020991 mod 499174.593750 grad 0.000058
sflapfill: iteration 246 res 0.020991 mod 499174.656250 grad 0.000055
sflapfill: iteration 247 res 0.020991 mod 499174.687500 grad 0.000054
sflapfill: iteration 248 res 0.020991 mod 499174.750000 grad 0.000052
sflapfill: iteration 249 res 0.020991 mod 499174.812500 grad 0.000051
sflapfill: iteration 250 res 0.020991 mod 499174.843750 grad 0.000048
sflapfill: iteration 251 res 0.020991 mod 499174.875000 grad 0.000047
sflapfill: iteration 252 res 0.020991 mod 499174.906250 grad 0.000047
sflapfill: iteration 253 res 0.020991 mod 499174.968750 grad 0.000046
sflapfill: iteration 254 res 0.020991 mod 499175.000000 grad 0.000046
sflapfill: iteration 255 res 0.020991 mod 499175.031250 grad 0.000049
sflapfill: iteration 256 res 0.020990 mod 499175.062500 grad 0.000050
sflapfill: iteration 257 res 0.020990 mod 499175.125000 grad 0.000049
sflapfill: iteration 258 res 0.020990 mod 499175.156250 grad 0.000047
sflapfill: iteration 259 res 0.020990 mod 499175.187500 grad 0.000045
sflapfill: iteration 260 res 0.020990 mod 499175.250000 grad 0.000045
sflapfill: iteration 261 res 0.020990 mod 499175.281250 grad 0.000043
sflapfill: iteration 262 res 0.020990 mod 499175.312500 grad 0.000042
sflapfill: iteration 263 res 0.020990 mod 499175.343750 grad 0.000041
sflapfill: iteration 264 res 0.020990 mod 499175.375000 grad 0.000041
sflapfill: iteration 265 res 0.020990 mod 499175.406250 grad 0.000041
sflapfill: iteration 266 res 0.020990 mod 499175.437500 grad 0.000042
sflapfill: iteration 267 res 0.020990 mod 499175.468750 grad 0.000043
sflapfill: iteration 268 res 0.020990 mod 499175.500000 grad 0.000042
sflapfill: iteration 269 res 0.020990 mod 499175.531250 grad 0.000040
sflapfill: iteration 270 res 0.020990 mod 499175.562500 grad 0.000039
sflapfill: iteration 271 res 0.020990 mod 499175.593750 grad 0.000039
sflapfill: iteration 272 res 0.020990 mod 499175.625000 grad 0.000037
sflapfill: iteration 273 res 0.020990 mod 499175.656250 grad 0.000036
sflapfill: iteration 274 res 0.020990 mod 499175.656250 grad 0.000037
sflapfill: iteration 275 res 0.020990 mod 499175.687500 grad 0.000037
sflapfill: iteration 276 res 0.020990 mod 499175.718750 grad 0.000036
sflapfill: iteration 277 res 0.020990 mod 499175.750000 grad 0.000035
sflapfill: iteration 278 res 0.020990 mod 499175.781250 grad 0.000033
sflapfill: iteration 279 res 0.020990 mod 499175.781250 grad 0.000032
sflapfill: iteration 280 res 0.020990 mod 499175.812500 grad 0.000031
sflapfill: iteration 281 res 0.020990 mod 499175.843750 grad 0.000029
sflapfill: iteration 282 res 0.020990 mod 499175.843750 grad 0.000027
sflapfill: iteration 283 res 0.020990 mod 499175.875000 grad 0.000026
sflapfill: iteration 284 res 0.020990 mod 499175.875000 grad 0.000025
sflapfill: iteration 285 res 0.020990 mod 499175.875000 grad 0.000024
sflapfill: iteration 286 res 0.020990 mod 499175.906250 grad 0.000023
sflapfill: iteration 287 res 0.020990 mod 499175.906250 grad 0.000022
sflapfill: iteration 288 res 0.020990 mod 499175.906250 grad 0.000020
sflapfill: iteration 289 res 0.020990 mod 499175.937500 grad 0.000020
sflapfill: iteration 290 res 0.020989 mod 499175.937500 grad 0.000019
sflapfill: iteration 291 res 0.020989 mod 499175.937500 grad 0.000018
sflapfill: iteration 292 res 0.020989 mod 499175.937500 grad 0.000018
sflapfill: iteration 293 res 0.020989 mod 499175.968750 grad 0.000017
sflapfill: iteration 294 res 0.020989 mod 499175.968750 grad 0.000016
sflapfill: iteration 295 res 0.020989 mod 499175.968750 grad 0.000016
sflapfill: iteration 296 res 0.020989 mod 499175.968750 grad 0.000014
sflapfill: iteration 297 res 0.020989 mod 499175.968750 grad 0.000014
sflapfill: iteration 298 res 0.020989 mod 499175.968750 grad 0.000013
sflapfill: iteration 299 res 0.020989 mod 499176.000000 grad 0.000014
sflapfill: iteration 300 res 0.020989 mod 499176.000000 grad 0.000013
sflapfill: iteration 301 res 0.020989 mod 499176.000000 grad 0.000013
sflapfill: iteration 302 res 0.020989 mod 499176.000000 grad 0.000012
sflapfill: iteration 303 res 0.020989 mod 499176.000000 grad 0.000012
sflapfill: iteration 304 res 0.020989 mod 499176.000000 grad 0.000011
sflapfill: iteration 305 res 0.020989 mod 499176.000000 grad 0.000011
sflapfill: iteration 306 res 0.020989 mod 499176.000000 grad 0.000010
sflapfill: iteration 307 res 0.020989 mod 499176.000000 grad 0.000010
sflapfill: iteration 308 res 0.020989 mod 499176.000000 grad 0.000009
sflapfill: iteration 309 res 0.020989 mod 499176.000000 grad 0.000009
sflapfill: iteration 310 res 0.020989 mod 499176.031250 grad 0.000009
sflapfill: iteration 311 res 0.020989 mod 499176.031250 grad 0.000009
sflapfill: iteration 312 res 0.020989 mod 499176.031250 grad 0.000009
sflapfill: iteration 313 res 0.020989 mod 499176.031250 grad 0.000009
sflapfill: iteration 314 res 0.020989 mod 499176.031250 grad 0.000008
sflapfill: iteration 315 res 0.020989 mod 499176.031250 grad 0.000007
sflapfill: iteration 316 res 0.020989 mod 499176.031250 grad 0.000007
sflapfill: iteration 317 res 0.020989 mod 499176.031250 grad 0.000007
sflapfill: iteration 318 res 0.020989 mod 499176.031250 grad 0.000007
sflapfill: iteration 319 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 320 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 321 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 322 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 323 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 324 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 325 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 326 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 327 res 0.020989 mod 499176.031250 grad 0.000006
sflapfill: iteration 328 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 329 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 330 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 331 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 332 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 333 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 334 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 335 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 336 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 337 res 0.020989 mod 499176.031250 grad 0.000005
sflapfill: iteration 338 res 0.020989 mod 499176.031250 grad 0.000004
sflapfill: iteration 339 res 0.020989 mod 499176.031250 grad 0.000004
sflapfill: iteration 340 res 0.020989 mod 499176.031250 grad 0.000004
sflapfill: iteration 341 res 0.020989 mod 499176.031250 grad 0.000004
sflapfill: iteration 342 res 0.020989 mod 499176.031250 grad 0.000004
sflapfill: iteration 343 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 344 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 345 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 346 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 347 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 348 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 349 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 350 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 351 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 352 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 353 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 354 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 355 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 356 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 357 res 0.020989 mod 499176.031250 grad 0.000003
sflapfill: iteration 358 res 0.020989 mod 499176.031250 grad 0.000002
sflapfill: iteration 359 res 0.020989 mod 499176.031250 grad 0.000002
sflapfill: iteration 360 res 0.020989 mod 499176.031250 grad 0.000002
sflapfill: iteration 361 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 362 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 363 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 364 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 365 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 366 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 367 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 368 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 369 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 370 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 371 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 372 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 373 res 0.020989 mod 499176.062500 grad 0.000002
sflapfill: iteration 374 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 375 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 376 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 377 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 378 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 379 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 380 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 381 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 382 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 383 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 384 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 385 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 386 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 387 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 388 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 389 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 390 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 391 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 392 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 393 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 394 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 395 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 396 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 397 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 398 res 0.020989 mod 499176.062500 grad 0.000001
sflapfill: iteration 399 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 400 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 401 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 402 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 403 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 404 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 405 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 406 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 407 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 408 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 409 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 410 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 411 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 412 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 413 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 414 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 415 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 416 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 417 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 418 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 419 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 420 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 421 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 422 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 423 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 424 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 425 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 426 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 427 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 428 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 429 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 430 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 431 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 432 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 433 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 434 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 435 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 436 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 437 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 438 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 439 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 440 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 441 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 442 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 443 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 444 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 445 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 446 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 447 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 448 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 449 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 450 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 451 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 452 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 453 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 454 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 455 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 456 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 457 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 458 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 459 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 460 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 461 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 462 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 463 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 464 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 465 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 466 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 467 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 468 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 469 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 470 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 471 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 472 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 473 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 474 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 475 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 476 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 477 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 478 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 479 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 480 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 481 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 482 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 483 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 484 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 485 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 486 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 487 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 488 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 489 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 490 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 491 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 492 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 493 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 494 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 495 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 496 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 497 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 498 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 499 res 0.020989 mod 499176.062500 grad 0.000000
sflapfill: iteration 500 res 0.020989 mod 499176.062500 grad 0.000000
< filled.rsf /home/sergey/RSFROOT/bin/sfgrey color=x bias=950 clip=133 scalebar=y title=Horizon yreverse=n transp=n screenwd=11.9 screenht=9.62 > cubeh.vpl
/home/sergey/RSFROOT/bin/vppen erase=o vpstyle=n cubeh.vpl > Fig/cubeh.vpl
< horizon.rsf /home/sergey/RSFROOT/bin/sfwindow f1=200 n1=125 f2=200 n2=100 > window.rsf
< window.rsf /home/sergey/RSFROOT/bin/sfgrey color=x bias=950 clip=133 scalebar=y title=Horizon yreverse=n transp=n screenratio=0.8 > Fig/window.vpl
echo 200 200 200 300 325 300 325 200 200 200 n1=2 n2=5 data_format=ascii_int in=wind.asc > wind.asc
< wind.asc /home/sergey/RSFROOT/bin/sfdd type=complex form=native | /home/sergey/RSFROOT/bin/sfwindow | /home/sergey/RSFROOT/bin/sfgraph wanttitle=n wantaxis=n plotcol=7 plotfat=10 min1=1 max1=595 min2=1 max2=463 scalebar=y screenwd=11.9 screenht=9.62 > wind.vpl
/home/sergey/RSFROOT/bin/vppen erase=o vpstyle=n cubeh.vpl wind.vpl > Fig/wind.vpl
< filled.rsf /home/sergey/RSFROOT/bin/sfnoise type=n seed=2014 range=15 > noisy.rsf
< noisy.rsf /home/sergey/RSFROOT/bin/sfsmooth rect1=4 rect2=4 > smoothed.rsf
< noisy.rsf /home/sergey/RSFROOT/bin/sfimpl2 rect1=10 rect2=10 tau=1 > diffused.rsf
sfimpl2: build/api/c/impl2.c: nstep=8 tau=(1.03125,1.03125) nclip=137741
sfimpl2: slice 0 of 1;
< noisy.rsf /home/sergey/RSFROOT/bin/sfdespike2 wide1=10 wide2=10 > median.rsf
< noisy.rsf /home/sergey/RSFROOT/bin/sfbilat2 a1=0.1 a3=0.001 r1=10 r2=10 > bilateral.rsf
< filled.rsf /home/sergey/RSFROOT/bin/sfgrey color=j bias=950 clip=133 scalebar=y title="Filled Horizon" yreverse=n transp=n minval=800 maxval=1100 > Fig/filled.vpl
< noisy.rsf /home/sergey/RSFROOT/bin/sfgrey color=j bias=950 clip=133 scalebar=y title="Noisy Horizon" yreverse=n transp=n minval=800 maxval=1100 > Fig/noisy.vpl
< filled.rsf /home/sergey/RSFROOT/bin/sfcat axis=3 noisy.rsf | /home/sergey/RSFROOT/bin/sfwindow n1=1 f1=300 n2=450 | /home/sergey/RSFROOT/bin/sfgraph plotfat=7,1 wanttitle=n yreverse=y label2=Time unit2=ms > Fig/noisy-slice.vpl
< smoothed.rsf /home/sergey/RSFROOT/bin/sfgrey color=j bias=950 clip=133 scalebar=y title="Smoothed Horizon" yreverse=n transp=n minval=800 maxval=1100 > Fig/smoothed.vpl
< filled.rsf /home/sergey/RSFROOT/bin/sfcat axis=3 smoothed.rsf | /home/sergey/RSFROOT/bin/sfwindow n1=1 f1=300 n2=450 | /home/sergey/RSFROOT/bin/sfgraph plotfat=7,1 wanttitle=n yreverse=y label2=Time unit2=ms > Fig/smoothed-slice.vpl
< diffused.rsf /home/sergey/RSFROOT/bin/sfgrey color=j bias=950 clip=133 scalebar=y title="Diffused Horizon" yreverse=n transp=n minval=800 maxval=1100 > Fig/diffused.vpl
< filled.rsf /home/sergey/RSFROOT/bin/sfcat axis=3 diffused.rsf | /home/sergey/RSFROOT/bin/sfwindow n1=1 f1=300 n2=450 | /home/sergey/RSFROOT/bin/sfgraph plotfat=7,1 wanttitle=n yreverse=y label2=Time unit2=ms > Fig/diffused-slice.vpl
< median.rsf /home/sergey/RSFROOT/bin/sfgrey color=j bias=950 clip=133 scalebar=y title="Median Horizon" yreverse=n transp=n minval=800 maxval=1100 > Fig/median.vpl
< filled.rsf /home/sergey/RSFROOT/bin/sfcat axis=3 median.rsf | /home/sergey/RSFROOT/bin/sfwindow n1=1 f1=300 n2=450 | /home/sergey/RSFROOT/bin/sfgraph plotfat=7,1 wanttitle=n yreverse=y label2=Time unit2=ms > Fig/median-slice.vpl
< bilateral.rsf /home/sergey/RSFROOT/bin/sfgrey color=j bias=950 clip=133 scalebar=y title="Bilateral Horizon" yreverse=n transp=n minval=800 maxval=1100 > Fig/bilateral.vpl
< filled.rsf /home/sergey/RSFROOT/bin/sfcat axis=3 bilateral.rsf | /home/sergey/RSFROOT/bin/sfwindow n1=1 f1=300 n2=450 | /home/sergey/RSFROOT/bin/sfgraph plotfat=7,1 wanttitle=n yreverse=y label2=Time unit2=ms > Fig/bilateral-slice.vpl
< filled.rsf /home/sergey/RSFROOT/bin/sfderiv > dx.rsf
< filled.rsf /home/sergey/RSFROOT/bin/sftransp | /home/sergey/RSFROOT/bin/sfderiv | /home/sergey/RSFROOT/bin/sftransp > dy.rsf
< dx.rsf /home/sergey/RSFROOT/bin/sfadd dy.rsf > shaded.rsf
echo 1 1000 n=1000,1000 n1=2 in=hel-slag0.asc data_format=ascii_int > hel-slag0.asc
< hel-slag0.asc /home/sergey/RSFROOT/bin/sfdd form=native > hel-slag.rsf
< hel-slag.rsf echo -1 -1 a0=2.0005 n1=2 lag=hel-slag.rsf in=hel-ss0.asc data_format=ascii_float > hel-ss0.asc
< hel-ss0.asc /home/sergey/RSFROOT/bin/sfdd form=native > hel-ss.rsf
echo 1 2 3 4 5 6 7 8 9 10 991 992 993 994 995 996 997 998 999 1000 n=1000,1000 n1=20 in=hel-alag0.asc data_format=ascii_int > hel-alag0.asc
< hel-alag0.asc /home/sergey/RSFROOT/bin/sfdd form=native > hel-alag.rsf
< hel-ss.rsf /home/sergey/RSFROOT/bin/sfwilson lagin=hel-alag0.asc lagout=hel-lag.rsf > hel.rsf
sfwilson: wilson 0 0.249938
sfwilson: wilson 1 0.098000
sfwilson: wilson 2 0.030887
sfwilson: wilson 3 0.008770
sfwilson: wilson 4 0.002203
sfwilson: wilson 5 0.002526
sfwilson: wilson 6 0.002592
sfwilson: wilson 7 0.002604
sfwilson: wilson 8 0.002607
sfwilson: wilson 9 0.002607
sfwilson: wilson 10 0.002607
sfwilson: wilson 11 0.002607
sfwilson: wilson 12 0.002607
sfwilson: wilson 13 0.002607
sfwilson: wilson 14 0.002607
sfwilson: wilson 15 0.002607
sfwilson: wilson 16 0.002607
sfwilson: wilson 17 0.002607
sfwilson: wilson 18 0.002607
sfwilson: wilson 19 0.002607
< filled.rsf /home/sergey/RSFROOT/bin/sfhelicon filt=hel.rsf > helix.rsf
< shaded.rsf /home/sergey/RSFROOT/bin/sfgrey title="Shaded Horizon" yreverse=n transp=n mean=y > Fig/shaded.vpl
< helix.rsf /home/sergey/RSFROOT/bin/sfgrey title="Helix Horizon" yreverse=n transp=n mean=y > Fig/helix.vpl

In [8]:

    

import m8r


    

In [7]:

    

m8r.Vplot('Fig/horizon.vpl')


    

    
Out[7]:

In [11]:

    

m8r.Vplot('Fig/wind.vpl')


    

    
Out[11]:

In [22]:

    

m8r.Vplot('Fig/bilateral-slice.vpl')


    

    
Out[22]:

In [24]:

    

m8r.Vplot('Fig/helix.vpl')


    

    
Out[24]:

In [ ]: