|
|
|
|
Marmousi Model |
| n1=726 | d1= 0.004 | o1=0.000 | label1=Depth | unit1=s | |
| n2=96 | d2=-0.025 | o2=2.575 | label2=Offset | unit2=km | |
| n3=240 | d3= 0.025 | o3=3.000 | label3=Shot-coord | unit3=km |
The file marmousi/shots/SConstruct gathers shot data, appends the header as necessary and produces several plots of the data. This file is reproduced here in Figure 5
from rsf.proj import *
# Fetch Files from repository
Fetch("marmrefl.hh","marm")
# Convert Files to RSF and update header
Flow('marmrefl','marmrefl.hh', '''dd form=native | put
label1=Depth Z label2=Offset unit2=km unit1=sec
label3=Position unit3=km d2=-.025 d3=.025 o2=2.575 o3=3''')
# Plotting Section
#Result('marmrefl','window j3=10 | grey')
Result('nearOffset','marmrefl','''window $SOURCE min2=.200 max2=.200 size2=1 |
grey color=I gainpanel=a title=Near Offset Data label1=Time label2=Position
screenratio=1.3 labelsz=4 titlesz=6 wheretitle=t''')
Result('shot20','marmrefl','''window $SOURCE min3=3.500 max3=3.500 size3=1 |
grey color=I gainpanel=a title=Shot # 20 wantframenum=y label2=Offset
label1=Time screenratio=1.3 labelsz=4 titlesz=6 wheretitle=t''')
End()
|
To run the script type Command 2 within the marmousishots directory.
|
|---|
|
shot20
Figure 2. Shot number 20 of Marmousi synthetic survey dataset. Shot position in km is shown in the lower left hand corner. |
|
|
|
|---|
|
nearOffset
Figure 3. Near offset data for Marmousi model |
|
|
|
|
|
|
Marmousi Model |