sumean.c

su 的源代码库

/* Copyright (c) Colorado School of Mines, 2006.*/
/* All rights reserved.                       */

/* SUMEAN: $Revision: 1.3 $ ; $Date: 2003/08/20 18:32:49 $	*/

#include "su.h"
#include "segy.h"

/*************************** self documentation ******************************/
char *sdoc[] = {
"									",
" SUMEAN - get the mean values of data traces				",	
"									",
" sumean < stdin > stdout [optional parameters] 			",
"									",
" Required parameters:							",
"   power = 2.0		mean to the power				",
"			(e.g. = 1.0 mean amplitude, = 2.0 mean energy)	",
"									",
" Optional parameters: 							",
"   verbose = 0		writes mean value of section to outpar	   	",
"			= 1 writes mean value of each trace / section to",
"				outpar					",
"   outpar=/dev/tty   output parameter file				",
"   abs = 1             average absolute value ",
"                       = 0 preserve sign if power=1.0",
"			 						",
" Notes:			 					",
" Each sample is raised to the requested power, and the sum of all those",
" values is averaged for each trace (verbose=1) and the section.	",
" The values power=1.0 and power=2.0 are physical, however other powers	",
" represent other mathematical L-p norms and may be of use, as well.	",
"			 						",
NULL};

/* Credits:
 *  Bjoern E. Rommel, IKU, Petroleumsforskning / October 1997
 *		    [e-mail]
 */

/**************** end self doc ***********************************************/


/* Globals */
segy tr;

int
main (int argc, char **argv)
{
  char *outpar ="";	   /* name of output file			*/
  FILE *outparfp = NULL;   /* file pointer of output file		*/
  int verbose = 0;	   /* flag for printing extra information	*/
  int itr = 0;   	   /* trace number				*/
  int it = 0;		   /* sample number				*/
  float power = 0.0;       /* mean to the power of			*/
  float tmean = 0.0;	   /* average mean of trace			*/
  float gmean = 0.0;	   /* average mean of section			*/
  int abs = 1;             /* absolute value flag */


  /* Initialize */
  initargs (argc, argv);
  requestdoc (1);

  /* Get optional parameters */
  if (!getparint ("verbose", &verbose))   verbose = 0;
  if ((verbose < 0) || (verbose > 1))   err ("unknown mode for verbose!");
  if (!getparstring ("outpar", &outpar))   outpar = "/dev/tty" ;
  outparfp = efopen (outpar, "w");
  if (!getparint ("abs", &abs))   abs = 0;
  if (!getparfloat ("power", &power))   power = 2.0;

  /* Get info from first trace */
  if (!gettr(&tr))  err ("can't get first trace");

  /* Initialize mean value of section */
  gmean = 0.0;

  /* Loop through traces */
  itr = 0;
  do {
    /* Reset */
    itr++ ;	/* count traces */
    tmean = 0.0;   /* reset mean value of trace */

    /* Loop through samples */
    for (it = 0; it < tr.ns; ++it) {
    
	/* Raise sample to the requested power, add to mean value of trace */
        if( abs ){

	   tmean += pow (fabs (tr.data[it]), power);

        }else{

	   tmean += pow ((tr.data[it]), power);

        }
    
    }
    
    /* Average mean value of trace */
    tmean = pow (tmean / tr.ns, 1.0 / power);

    /* Add to the mean value of section */
    gmean += tmean;

    /* Print mean value of trace */
    if (verbose == 1)   
	fprintf (outparfp, "trace: %i   mean value: %e\n", itr, tmean);

  } while (gettr(&tr)); 

  /* Average mean value of section */
  gmean = gmean / itr;

  /* Print mean value of section */
  fprintf (outparfp, "global mean: %e\n", gmean);
  
  return(CWP_Exit());
}