segy2grd
segy2grd - Converting SEGY data to a GMT grid
Synopsis
segy2grd segyfile -Ggrdfile
-Ixinc[unit][=|+][/yinc[unit][=|+]]
-R[unit]xmin/xmax/ymin/ymax[r]
[ -A[n|z] ]
[ -Dxname/yname/zname/scale/offset/title/remark ]
[ -M[flags] ] [ -Nnodata ] [ -S[zfile] ]
[ -V[level] ]
[ -bi[ncols][type][w][+L|+B] ]
[ -:[i|o] ]
Note: No space is allowed between the option flag and the associated arguments.
Description
segy2grd reads an IEEE SEGY file and creates a binary grid file.
Either a simple mapping (equivalent to xyz2grd -Z) or a more complicated
averaging where a particular grid cell includes values from more than
one sample in the SEGY file can be done. segy2grd will report if
some of the nodes are not filled in with data. Such unconstrained nodes
are set to a value specified by the user [Default is NaN]. Nodes with
more than one value will be set to the average value.
Required Arguments
segyfile is an IEEE floating point SEGY file. Traces are all assumed to start at 0 time/depth.
- -Ggrdfile
- grdfile is the name of the binary output grid file.
- -I
- x_inc [and optionally y_inc] is the grid spacing. Append m
to indicate minutes or c to indicate seconds.
- -Rwest/east/south/north[/zmin/zmax][r]
- west, east, south, and north specify the region of interest,
and you may specify them in decimal degrees or in
[+-]dd:mm[:ss.xxx][W|E|S|N] format. Append r if lower left
and upper right map coordinates are given instead of w/e/s/n. The
two shorthands -Rg and -Rd stand for global domain (0/360
and -180/+180 in longitude respectively, with -90/+90 in latitude).
Alternatively, specify the name of an existing grid file and the
-R settings (and grid spacing, if applicable) are copied from
the grid. Using -Runit expects projected (Cartesian)
coordinates compatible with chosen -J and we inversely project
to determine actual rectangular geographic region.
For perspective view (-p), optionally append /zmin/zmax.
In case of perspective view (-p), a z-range (zmin, zmax)
can be appended to indicate the third dimension. This needs to be
done only when using the -Jz option, not when using only the
-p option. In the latter case a perspective view of the plane is
plotted, with no third dimension.
Optional Arguments
- -A[n|z]
- Add up multiple values that belong to the same node (same as
-Az). Append n to simply count the number of data points
that were assigned to each node. [Default (no -A option) will
calculate mean value]. Not used for simple mapping.
- -Dxname/yname/zname/scale/offset/title/remark
- Give values for xname, yname, zname, scale, offset,
title, and remark. To leave some of these values untouched,
specify = as the value.
- -M[flags]
- Fix number of traces to read in. Default tries to read 10000 traces.
-M0 will read number in binary header, -Mn will
attempt to read only n traces.
- -Nnodata
- No data. Set nodes with no input sample to this value [Default is
NaN].
- -S[zfile]
- set variable spacing header is c for cdp, o for offset, b<number>
for 4-byte float starting at byte number If -S not set, assumes even
spacing of samples at the dx, dy supplied with -I
- -L
- Override number of samples in each trace
- -X
- applies scalar x-scale to coordinates in trace header to match the
coordinates specified in -R
- -Y
- Specifies sample interval as s_int if incorrect in the SEGY file
- -V[level] (more ...)
- Select verbosity level [c].
- -r (more ...)
- Set pixel node registration [gridline].
- -^
- Print a short message about the syntax of the command, then exits (NOTE: on Windows use just -).
- +
- Print an extensive usage (help) message, including the explanation of
any module-specific option (but not the GMT common options), then exits.
- -?
- Print a complete usage (help) message, including the explanation of
options, then exits.
- --version
- Print GMT version and exit.
- --show-datadir
- Print full path to GMT share directory and exit.
Examples
To create a grid file from an even spaced SEGY file test.segy, try
gmt segy2grd test.segy -I0.1/0.1 -Gtest.nc -R198/208/18/25 -V
Note that this will read in 18-25s (or km) on each trace, but the
first trace will be assumed to be at X=198
To create a grid file from the SEGY file test.segy, locating traces
according to the CDP number, where there are 10 CDPs per km and the
sample interval is 0.1, try
gmt segy2grd test.segy -Gtest.nc -R0/100/0/10 -I0.5/0.2 -V -X0.1 -Y0.1
Because the grid interval is larger than the SEGY file sampling, the
individual samples will be averaged in bins