Hi,
good... but not yet very good !
I suspect a problem of numerical accuracy with such high values of X and
Y coordinates.
You can test this hypothesis by subtracting the mean of X and Y
coordinates to have 0 centered data...
If it's true, it should be better to correct the algorithm, I know...
D.
Le 07/04/2023 à 08:41, JUAN LORENZO DIAZ NAVEAS a écrit :
Dear Glenn and Dominique,
Thanks a lot for your help. Indeed the keyword scalco=-10000 was
guilty for scaling sx,sy,gx and gy coordinates. So now I am able to
produce an su file with cdp headers.
But now unfortunately I found another problem. The su file with cdp
headers has about 66% of cdp=0. The cdp=0 traces are at the
"beginning" of the su file.
And the intriguing aspect is that I have a crooked line with 10
points. See the attached map. There the black dots are the source
coordinates, the red dots are the receiver coordinates and the green
circles are the points of the crooked line.
As you can guess I smoothed both the source and receiver coordinates
using the GMT trend1d command.
The command line for my binning looks like this:
*sucdpbin
xline=723193.2,744776.011111,766358.822222,787941.633333,809524.444444,831107.255556,852690.066667,874272.877778,895855.688889,920934
\
yline=7242472.4,7242099.3532,7241693.29184,7241254.1793,7240782.11164,7240277.08676,7239739.07882,7239168.07407,7238564.07513,7237820
\
< SO104-17_geom.su verbose=2 dcdp=25 cdpmin=1 distmax=100 >
SO104-17_geom-cmp.su
*
See below the ranges of the file SO104-17_geom-cmp.su and the
attached log file with the verbosity of the sucdpbin command.
I tried also with distmax less than 100m, the smallest value being
15m, but the result was similar.
Any idea on what is going on?
Best wishes
Juan Díaz Naveas
-------------------------------------------------------------------
Juan Diaz Naveas (Dr. rer. nat.)
Escuela de Ciencias del Mar
www.cienciasdelmar.pucv.cl
<http://www.cienciasdelmar.pucv.cl>
Pontificia Universidad Catolica de Valparaiso Tel.: +56-32-2274269
Avenida Universidad 330, Curauma Cel.: +56-9-91630936
Valparaiso, CHILE E-Mail: juan.diaz(a)pucv.cl
-------------------------------------------------------------------
*466560 traces:
tracl 1 513204 (1 - 513204)
tracr 1 513204 (1 - 513204)
fldr 1 3888 (1 - 3888)
tracf 1 120 (1 - 120)
ep 1 3888 (1 - 3888)
cdp 0 2572 (0 - 14)
trid 1
nhs 1
offset 413 3388 (3388 - 413)
scalel -10000
scalco 1
sx 723302 917547 (917547 - 723302)
sy 7237924 7242471 (7237924 - 7242471)
gx 723715 920934 (920934 - 723715)
gy 7237823 7242464 (7237823 - 7242464)
counit 1
ns 3585
dt 4000
gain 9
afilf 102
afils 209
lcf 3
hcf 102
lcs 6
hcs 209
year 1995
Shot coordinate limits:
North(723302,7.24247e+06) South(917547,7.23792e+06)
East(917547,7.23792e+06) West(723302,7.24247e+06)
Receiver coordinate limits:
North(723740,7.24246e+06) South(920934,7.23782e+06)
East(920934,7.23782e+06) West(723715,7.24246e+06)
Midpoint coordinate limits:
North(723521,7.24247e+06) South(919240,7.23787e+06)
East(919240,7.23787e+06) West(723508,7.24247e+06)*
El jue, 6 abr 2023 a las 3:10, Dominique ROUSSET
(<dominique.rousset(a)univ-pau.fr>) escribió:
Hi Juan,
after a quick look at sucdpbin.c, i note that this tools honors
the scalco keyword. So, since your data bear a -10000 scalco, you
have to divide xline and yline parameters by 10000 to work with
actual units, not with sx and sy units.
Let me know it it works
Good luck !
D.
Le 06/04/2023 à 02:47, JUAN LORENZO DIAZ NAVEAS a écrit :
Dear SU users,
I am trying to make the CMP binning of a seismic line (see the
headers ranges below).
The command line is:
*sucdpbin < so104-17_geom.su <http://so104-17_geom.su>
xline=723715,920883 yline=7242464,7237825 verbose=2 dcdp=25
distmax=500 > so104-17_cdp.su <http://so104-17_cdp.su>*
As you see I have used a distmax of 500m and still I am not able
to get the binning.
Just in case I also tried again the same command, but sorting the
date so that the first traces in the su file correspond to the
minimum x coordinate and the last traces to the maximum
coordinate (in the field the line was shot from east to west),
but it failed as well.
Any help would be very much appreciated.
Best wishes
Juan
*466560 traces:
tracl 1 513204 (1 - 513204)
tracr 1 513204 (1 - 513204)
fldr 1 3888 (1 - 3888)
tracf 1 120 (1 - 120)
ep 1 3888 (1 - 3888)
trid 1
nhs 1
offset 413 3388 (3388 - 413)
scalel -10000
scalco -10000
sx 723302 917547 (917547 - 723302)
sy 7237924 7242471 (7237924 - 7242471)
gx 723715 920934 (920934 - 723715)
gy 7237823 7242464 (7237823 - 7242464)
counit 1
ns 3585
dt 4000
gain 9
afilf 102
afils 209
lcf 3
hcf 102
lcs 6
hcs 209
year 1995
Shot coordinate limits:
North(-91.7547,-723.792) South(-72.3302,-724.247)
East(-72.3302,-724.247) West(-91.7547,-723.792)
Receiver coordinate limits:
North(-92.0934,-723.782) South(-72.374,-724.246)
East(-72.3715,-724.246) West(-92.0934,-723.782)
Midpoint coordinate limits:
North(-91.924,-723.787) South(-72.3521,-724.247)
East(-72.3509,-724.247) West(-91.924,-723.787)*
-------------------------------------------------------------------
Juan Diaz Naveas (Dr. rer. nat.)
Escuela de Ciencias del Mar
www.cienciasdelmar.pucv.cl
<http://www.cienciasdelmar.pucv.cl>
Pontificia Universidad Catolica de Valparaiso Tel.: +56-32-2274269
Avenida Universidad 330, Curauma Cel.: +56-9-91630936
Valparaiso, CHILE E-Mail: juan.diaz(a)pucv.cl
-------------------------------------------------------------------
_______________________________________________
Seisunix mailing list --seisunix(a)mailman.seismic-unix.org
To unsubscribe send an email toseisunix-leave(a)mailman.seismic-unix.org
--
|Dominique Rousset LFCR UMR5150 / E2S UPPA / Dpt des géosciences|
|Univ. Pau et des Pays de l'Adour BP 1155 F-64013 Pau Cedex - France|
|mailto:dominique.rousset@univ-pau.fr <mailto:dominique.rousset@univ-pau.fr>
http://lfc.univ-pau.fr|
| Off: +33 559 407 423 Cell: +33 660 861 391 |
_______________________________________________
Seisunix mailing list -- seisunix(a)mailman.seismic-unix.org
To unsubscribe send an email to
seisunix-leave(a)mailman.seismic-unix.org
_______________________________________________
Seisunix mailing list --seisunix(a)mailman.seismic-unix.org
To unsubscribe send an email toseisunix-leave(a)mailman.seismic-unix.org
--
|Dominique Rousset LFCR UMR5150 / E2S UPPA / Dpt des géosciences|
|Univ. Pau et des Pays de l'Adour BP 1155 F-64013 Pau Cedex - France|
|mailto:dominique.rousset@univ-pau.fr
| Off: +33 559 407 423 Cell: +33 660 861 391 |