from rsf.proj import *
# Fetch Files from repository
Fetch("obc_vz_1.segy","marm2")
Fetch("obc_vz_2.segy","marm2")
# Convert Files to RSF and update header
Flow('obc_vz_1','obc_vz_1.segy',
'''segyread read=d |
put n2=1381 n3=240 o1=0 o2=0 o3=3000
d2=12.32 d3=25''')
Flow('obc_vz_2','obc_vz_2.segy',
'''segyread read=d |
put n2=1381 n3=240 o1=0 o2=0 o3=11000
d2=12.32 d3=25
''')
# Concatinate Datasets
Flow('vz',['obc_vz_1','obc_vz_2'],'cat ${SOURCES[0:2]} axis=3',stdin=0)
# Plot Data
Result('zero','vz',
'''window $SOURCE min2=0 max2=0 size2=1 |
grey color=I gainpanel=a label2=Position X unit2=m
title=Zero Offset Data''')
Result('zero2','vz',
'''window $SOURCE min2=0 max2=0 size2=1 |
grey color=I gainpanel=a label2=Position X unit2=m
title=Zero Offset Data''')
End()
|