From 856b67139148861b60ed3afee37e0aa72632f964 Mon Sep 17 00:00:00 2001 From: sfomel Date: Tue, 9 Apr 2024 14:49:13 -0500 Subject: [PATCH] Oz Yilmaz gathers --- book/data/oz/gathers/SConstruct | 26 +++++++++++++++----------- 1 file changed, 15 insertions(+), 11 deletions(-) diff --git a/book/data/oz/gathers/SConstruct b/book/data/oz/gathers/SConstruct index 1e9fdd800..a4f15277e 100644 --- a/book/data/oz/gathers/SConstruct +++ b/book/data/oz/gathers/SConstruct @@ -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()