Re: Use of pssegy

In setting up the -R the "X" coordinates are distance along the line,
in most basic case the trace numbers, in your case the range of offsets
you want to plot. "Y" corresponds
to time in seconds unless you do something weird with the sample interval.

Using -So as you have done then the trace headers for offset will be read
and interpreted as the X coordinate.

The reason for using -J is that is what GMT always uses for conversion between
data coordinates and distance on a page - for pssegy only -JX or -Jx makes
much sense.

I suggest you start with something simple like
pssegy a.SEGY -R0/1000/0/10 -JX15/-15 -W -D0.3 -V -N >> segy.ps
with -V it will print to the screen the X coordinate it is trying to plot
each trace at and the RMS amplitude of each trace. My guess is that the failure is because of
your Y coordinates, although it could also fail if your seismic data aren't
stored as native floating point.

Regards,

Tim
------------------------------------------------------
Tim Henstock, Reader in Geophysics
National Oceanography Centre Southampton
University of Southampton
European Way
Southampton

[email protected]
+44 (0)23 8059 6491
________________________________________
From: World-Wide GMT Usage and Help Mailing List [[email protected]] on behalf of Andreas B [[email protected]]
Sent: 23 April 2014 20:21
To: [email protected]
Subject: Re: [GMT-HELP] Use of pssegy

Hi Gery,

Thank you so much, it helps! I'm a bit confused though; does the values in -R refer to the Xmin, Xmax, Ymin and Ymax of the spatial extent of the line (like an ordinary -R would)? However, I ran into another problem; psegy crashes when doing a "pssegy a.SEGY -R17052/27803/731058/732041 -Jx1:1000000 -D0.3 -W -So -Fgray > segy.ps<http://segy.ps>". Tried to change all parameters, crashes everytime.

Andreas


2014-04-23 20:55 GMT+02:00 Gery . <[email protected]<mailto:[email protected]>>:
Hello Andreas,

Would this help? http://www.geo.utep.edu/gmt/seismic.html

Cheers,

Sent from my iRon

On Apr 23, 2014, at 20:43, "Andreas B" <[email protected]<mailto:[email protected]>> wrote:

Hey all,

Does anyone have a working code snippet of how to use pssegy(/psegyz)?

I'm trying a very basic pssegy a.SEGY -W -R17052/27804/732041/731059 -D1  > test.ps<http://test.ps>, but GMT complains that the -R is not correctly given (I've also tried with z-range in all flavours). Also, why is -J a required parameter? Why would you use a map projection when you want to plot a SEGY?

Andreas
To unsubscribe, send the message "signoff gmt-help" to [email protected]<mailto:[email protected]>

To unsubscribe, send the message "signoff gmt-help" to [email protected]
--
This message has been scanned for viruses and
dangerous content by MailScanner<http://www.mailscanner.info/>, and is
believed to be clean.

To unsubscribe, send the message "signoff gmt-help" to [email protected]