Category

Development Platform
Home > SourceCode/Document > Other systems > View Code
fgettra.c
Package: seimiclabtools.rar　[view]
Upload User: qiu126
Upload Date: 2015-02-06
Package Size: 3950k
Code Size: 3k
Category:

Other systems

Development Platform:

Visual C++

fgettra.c：Code Content
  1. /* FGETTRA: $Revision: 1.6 $; $Date: 90/12/23 17:03:09 $ */
  3. /*----------------------------------------------------------------------
  4.  * Copyright (c) Colorado School of Mines, 1989.
  5.  * All rights reserved.
  6.  *
  7.  * This code is part of SU.  SU stands for Seismic Unix, a processing line
  8.  * developed at the Colorado School of Mines, partially based on Stanford
  9.  * Exploration Project (SEP) software.  Inquiries should be addressed to:
  10.  *
  11.  *  Jack K. Cohen, Center for Wave Phenomena, Colorado School of Mines,
  12.  *  Golden, CO 80401  (jkc@dix.mines.colorado.edu)
  13.  *----------------------------------------------------------------------
  14.  */
  16. #include "su.h"
  17. #include "segy.h"
  18. #include "header.h"
  20. /* fgettra - get disk trace by trace number
  21.  *
  22.  * Returns:
  23.  * integer number of traces in the disk file
  24.  *
  25.  * Synopsis:
  26.  * #include "segy.h"
  27.  * fgettra(FILE *fp, segy *traceptr, int itr)
  28.  *
  29.  * gettra(tp, itr) is a macro for fgettra(stdin, tp, itr)
  30.  *
  31.  * Credits:
  32.  * SEP: Einar
  33.  * CWP: Jack
  34.  *
  35.  */
  38. int fgettra(FILE *fp, segy *tp, int itr)
  39. {
  40. static int nsfirst; /* samples on first trace */
  41. static int nsegy; /* size of trace in bytes */
  42. static int ntr; /* number of traces in file */
  43. static bool first=true; /* flag for first entry */
  44. int nread; /* bytes read */
  45. static char card[EBCBYTES];
  46. static char bhdr[BNYBYTES];
  47. static int hdbytes;
  49. if (first) { /* first entry; get number of traces */
  50. unsigned short bytesper;/* bytes per float (packed?) */
  51. long length; /* length of file in bytes */
  52. filetype ftype = filestat(fileno(fp));
  54. first = false;
  56. switch(ftype) {
  57. case TTY:
  58. warn("stdin not redirected");
  59. selfdoc();
  60. break;
  61. case DISK: /* correct */
  62. break;
  63. default:
  64. err("%s: input must be disk file", __FILE__);
  65. break;
  66. }
  68. efread(card, 1, 10, fp);
  69.                 if (strncmp(card, "C 1 CLIENT",10)==0 ) {
  70. efseek(fp, EBCBYTES+BNYBYTES, 0);
  71. hdbytes = EBCBYTES + BNYBYTES;
  72. }
  73. else {
  74. efseek(fp, 0, 0);
  75. hdbytes = 0;
  77. }
  78.    
  80. if (HDRBYTES != (nread = efread(tp, 1, HDRBYTES, fp))) {
  81. err("%s: read only %d bytes of header",
  82. __FILE__, nread);
  83. }
  85. nsfirst = tp->ns;
  86. if (nsfirst > SU_NFLTS) {
  87. err("%s: trace too long: nsfirst=%d>SU_NFLTS=%d",
  88. __FILE__, nsfirst, SU_NFLTS);
  89. }
  91. if      (tp->trid == CHARPACK)  bytesper = sizeof(char);
  92. else if (tp->trid == SHORTPACK) bytesper = sizeof(short);
  93. else                            bytesper = sizeof(float);
  95. nsegy = HDRBYTES + nsfirst * bytesper;
  96. efseek(fp, 0, SEEK_END);
  97. length = eftell(fp);
  98. ntr = (length-hdbytes) / nsegy;
  100. } /* end first entry */
  103. /* Check on requested trace number */
  104. if (itr >= ntr)  err("%s, trying to read off end of file", __FILE__);
  107. /* Position file pointer at start of requested trace */
  108. efseek(fp, itr*nsegy+hdbytes, SEEK_SET);
  110. nread = efread(tp, 1, nsegy, fp);
  111. if (nread != nsegy)
  112. err("%s: read only %d of %d bytes in trace",
  113. __FILE__, nread, nsegy);
  115. if (tp->ns != nsfirst)
  116.     warn("%s: header ns field = %d differs from first trace = %d",
  117. __FILE__, tp->ns, nsfirst);
  119. return ntr;
  120. }