Skip to content

Commit

Permalink
Oz Yilmaz gathers
Browse files Browse the repository at this point in the history
  • Loading branch information
sfomel committed Apr 9, 2024
1 parent 6be505c commit 856b671
Showing 1 changed file with 15 additions and 11 deletions.
26 changes: 15 additions & 11 deletions book/data/oz/gathers/SConstruct
Original file line number Diff line number Diff line change
@@ -1,22 +1,26 @@
from rsf.proj import *
import os

# Download a collection of shot gathers and convert them from SU to RSF
# See https://seismicrocks.com/oz40.html for description
# Download a collection of shot gathers from Oz Yilmaz
# and convert them from SEGY to RSF
# Courtesy of David Forel
# See https://seismicrocks.com/oz40.html

tgz = 'ozdata._SU.tgz'
tgz = 'oz.forty.sgy.tgz'

Fetch(tgz, 'download',
server='https://nextcloud.seismic-unix.org',
top='s/wfbgjReLLjRjLwk')
Fetch(tgz, 'seismicunix',
server='https://seismicrocks.com',
top='')

files = [f'ozdata.{i:>02}.su' for i in range(1, 41)]
files = [f'oz.{i:>02}.sgy' for i in range(1, 41)]
Flow(files,tgz,
'gunzip -c $SOURCE | tar -xvf -',stdin=0,stdout=-1)

for su in files:
rsf = os.path.splitext(su)[0]+'.rsf'
Flow([rsf,'t-'+rsf],su,'suread tfile=${TARGETS[1]} suxdr=y')
Result(rsf,f'grey title={rsf}')
for sgy in files:
name = os.path.splitext(sgy)[0]
rsf = name+'.rsf'
Flow([rsf,'t-'+rsf,name+'.asc',name+'.bin'],sgy,
'segyread tfile=${TARGETS[1]} hfile=${TARGETS[2]} bfile=${TARGETS[3]}')
Result(rsf,f'grey title={name}')

End()

0 comments on commit 856b671

Please sign in to comment.