/*-------------------------------------------------------------------- * The MB-system: mb_io.h 1/19/93 * $Id: mb_segy.h 1942 2012-03-16 19:20:44Z caress $ * * Copyright (c) 2004-2012 by * David W. Caress (caress@mbari.org) * Monterey Bay Aquarium Research Institute * Moss Landing, CA 95039 * and Dale N. Chayes (dale@ldeo.columbia.edu) * Lamont-Doherty Earth Observatory * Palisades, NY 10964 * * See README file for copying and redistribution conditions. *--------------------------------------------------------------------*/ /* * mb_segy.h defines the SEG-Y trace header used by MB-System when * extracting seismic reflection or subbottom profiler data from * swath mapping data files. The SIOSEIS implementation of the * trace header is used because it incorporates a deep water delay * value not found in most other implementations.See the SIOSEIS * web pages for information:. * http://sioseis.ucsd.edu/ * http://sioseis.ucsd.edu/segy.header.html * * Author: D. W. Caress * Date: April 13, 2004 * * $Log: mb_segy.h,v $ * Revision 5.7 2008/09/13 06:08:09 caress * Updates to apply suggested patches to segy handling. Also fixes to remove compiler warnings. * * Revision 5.6 2006/11/10 22:36:04 caress * Working towards release 5.1.0 * * Revision 5.5 2005/05/05 23:53:01 caress * Just added a space in a line of code... * * Revision 5.4 2004/09/16 19:02:33 caress * Changes to better support segy data. * * Revision 5.3 2004/07/27 19:44:38 caress * Working on handling subbottom data. * * Revision 5.2 2004/06/18 05:22:32 caress * Working on adding support for segy i/o and for Reson 7k format 88. * * Revision 5.1 2004/05/21 23:44:49 caress * Progress supporting Reson 7k data, including support for extracing subbottom profiler data. * * Revision 5.0 2004/04/27 01:50:15 caress * Adding support for Reson 7k sonar data, including segy extensions. * * */ /* Standard SEGY format sizes */ #define MB_SEGY_ASCIIHEADER_LENGTH 3200 #define MB_SEGY_FILEHEADER_LENGTH 400 #define MB_SEGY_TRACEHEADER_LENGTH 240 /* Flags used to specify desired data type in mb_extract_segy() calls */ #define MB_SEGY_SAMPLEFORMAT_NONE 1 #define MB_SEGY_SAMPLEFORMAT_TRACE 2 #define MB_SEGY_SAMPLEFORMAT_ENVELOPE 3 #define MB_SEGY_SAMPLEFORMAT_ANALYTIC 4 /* SEGY structures */ struct mb_segyasciiheader_struct { char line[40][80]; }; struct mb_segyfileheader_struct { int jobid; int line; int reel; unsigned short channels; unsigned short aux_channels; unsigned short sample_interval; unsigned short sample_interval_org; unsigned short number_samples; unsigned short number_samples_org; unsigned short format; /* 1 IBM 32 bit floating point 2 32 bit integer 3 16 bit integer 5 IEEE 32 bit floating point 6 IEEE 32 bit floating point 8 8 bit integer 11 Little-endian IEEE 32 bit floating point */ unsigned short cdp_fold; unsigned short trace_sort; unsigned short vertical_sum; unsigned short sweep_start; unsigned short sweep_end; unsigned short sweep_length; unsigned short sweep_type; unsigned short sweep_trace; unsigned short sweep_taper_start; unsigned short sweep_taper_end; unsigned short sweep_taper; unsigned short correlated; unsigned short binary_gain; unsigned short amplitude; unsigned short units; unsigned short impulse_polarity; unsigned short vibrate_polarity; short domain; short rev; short fixed_length; short num_ext_headers; char extra[338]; }; struct mb_segytraceheader_struct { int seq_num; /* bytes 0-3, trace sequence number in the line */ int seq_reel; /* bytes 4-7, trace sequence number in the reel */ int shot_num; /* bytes 8-11, shot number or stacked trace number "Original field record number" */ int shot_tr; /* bytes 12-15, trace number within the shot */ int espn; /* bytes 16-19, "Energy source point number -- used when more than one record occurs at the same effective surface location". */ int rp_num; /* bytes 20-23, rp or cdp number */ int rp_tr; /* bytes 24-27, trace number within the cdp */ unsigned short trc_id; /* bytes 28-29, trace id: 1= live, 2=dead */ unsigned short num_vstk; /* bytes 30-31, Number of traces vertically stacked */ unsigned short cdp_fold; /* bytes 32-33, cdp fold (coverage) */ unsigned short use; /* bytes 34-35, Data use: 1=production, 2=test */ int range; /* bytes 36-39, source to receiver distance (range) */ int grp_elev; /* bytes 40-43, Receiver group elevation w.r.t. sea level (depth is negative) */ int src_elev; /* bytes 44-47, Source elevation */ int src_depth; /* bytes 48-51, Source depth below surface. (depth is a positive number!) */ int grp_datum; /* bytes 51-55, Datum elevation at receiver group. */ int src_datum; /* bytes 56-59, Datum elevation at source */ int src_wbd; /* bytes 60-63, water depth at the source */ int grp_wbd; /* bytes 64-67, water depth at the receiver group. */ short elev_scalar; /* bytes 68-69, Scalar to be applied to elevations and depths in bytes 41-68. If positive use as a multiplier. If negative, use as a divisor. */ short coord_scalar; /* bytes 70-71, Scalar to be applied to coordinates in bytes 72-87. If positive use as a multiplier. If negative, use as a divisor. */ int src_long; /* bytes 72-75, longitude in seconds of arc. Source X coordinate */ int src_lat; /* bytes 76-79, latitude in seconds of arc. Source Y coordinate */ int grp_long; /* bytes 80-83, Receiver longitude or X coordinate */ int grp_lat; /* bytes 84-87, Receiver latitude or Y coordinate */ short coord_units; /* bytes 88-89, = 2, coordinate units = seconds of arc */ short wvel; /* bytes 90-91, weathering or water velocity */ short sbvel; /* bytes 92-93, subweathering velocity */ short src_up_vel; /* bytes 94-95, uphole time at source */ short grp_up_vel; /* bytes 96-97, uphole time at group */ short src_static; /* bytes 98-99, Source static correction */ short grp_static; /* bytes 100-101, Group static correction */ short tot_static; /* bytes 102-103, Total static applied */ short laga; /* bytes 104-105, Lag time A in ms. before time 0 */ short lagb; /* bytes 106-107, Lag time B in ms. before time 0 */ int delay_mils; /* bytes 106-109, deep water delay in ms. (or meters) */ short smute_mils; /* bytes 110-111, start mute time in ms. */ short emute_mils; /* bytes 112-113, end mutes time in ms. */ unsigned short nsamps; /* bytes 114-115, "Number of data samples in this trace" - excludes header */ unsigned short si_micros; /* bytes 116-117, Sample interval in us for this trace */ short other_1[19]; /* bytes 118-155, Other short integer stuff */ short year; /* bytes 156-157, year data was recorded. */ short day_of_yr; /* bytes 158-159, recording day of year */ short hour; /* bytes 160-161, recording hour of day */ short min; /* bytes 162-163, recording minute of hour */ short sec; /* bytes 164-165, recording second of minute */ short mils; /* bytes 166-167, recording millisecond */ /* OFFICIAL SEGY says: "time basis code" */ short tr_weight; /* bytes 168-169, Trace weighting factor */ short other_2[5]; /* bytes 170-179, Other short integer stuff */ float delay; /* bytes 180-183, deep water delay in seconds (or meters) */ float smute_sec; /* bytes 184-187, start mute time in seconds */ float emute_sec; /* bytes 188-191, end mute time in seconds */ float si_secs; /* bytes 192-195, sample interval in seconds */ float wbt_secs; /* bytes 196-199, water bottom time in seconds */ int end_of_rp; /* bytes 200-203, <0, indicates the end of a gather >0, the number of traces stacked Also EdgeTech's Trace Scalar. */ float dummy1; /* bytes 204-207 */ float dummy2; /* bytes 208-211 */ float dummy3; /* bytes 212-215 */ float dummy4; /* bytes 216-219 */ float soundspeed; /* bytes 220-223, water sound speed in m/sec (MB-System only) */ float distance; /* bytes 224-227, distance from previous shot along track in meters (MB-System only) */ float roll; /* bytes 228-231, roll in degrees (MB-System only) */ float pitch; /* bytes 232-235, pitch in degrees (MB-System only) */ float heading; /* bytes 236-239, heading in degrees (MB-System only) */ }; struct mb_segyio_struct { FILE *fp; char segyfile[MB_PATH_MAXLINE]; int bufferalloc; char *buffer; int asciiheader_set; int fileheader_set; struct mb_segyasciiheader_struct asciiheader; struct mb_segyfileheader_struct fileheader; struct mb_segytraceheader_struct traceheader; int tracealloc; float *trace; }; /* function prototypes */ int mb_segy_read_init(int verbose, char *segyfile, void **mbsegyio_ptr, struct mb_segyasciiheader_struct *segyasciiheader, struct mb_segyfileheader_struct *segyfileheader, int *error); int mb_segy_write_init(int verbose, char *segyfile, struct mb_segyasciiheader_struct *asciiheader, struct mb_segyfileheader_struct *fileheader, void **mbsegyio_ptr, int *error); int mb_segy_close(int verbose,void **mbsegyio_ptr, int *error); int mb_segy_read_trace(int verbose, void *mbsegyio_ptr, struct mb_segytraceheader_struct *traceheaderptr, float **traceptr, int *error); int mb_segy_write_trace(int verbose, void *mbsegyio_ptr, struct mb_segytraceheader_struct *traceheader, float *trace, int *error); void hilbert(int n, double delta[], double kappa[]); void hilbert2(int n, double data[]);