sfpow multiplies the input data by a gain function of the form

$$G\left({x}_{1},{x}_{2},\dots ,{x}_{n}\right)={x}_{1}^{{p}_{1}}\phantom{\rule{0.167em}{0ex}}{x}_{2}^{{p}_{2}}\phantom{\rule{0.167em}{0ex}}\cdots \phantom{\rule{0.167em}{0ex}}{x}_{n}^{{p}_{n}}$$ The powers

$${p}_{1},{p}_{2},\dots ,{p}_{n}$$ are given by

**pow1=**,

**pow2=**, etc. parameters.

For backward compatibility,

**sftpow tpow=** is an alias for

**sfpow pow1=**.

The following example from

geo391/hw1/tpow shows an application of

**sfpow pow1=2 **to a shot gather. The gain of

**pow1=2** for seismic data was

advocated by Jon Claerbout.

To estimate an appropriate power gain from the data, you can try

sffpow.

### 10 previous programs of the month:

sfnmo implements normal moveout (NMO) correction, one of the most fundamental operations in seismic reflection data processing. The following example from jsg/avo/avo shows synthetic data before and after NMO correction. NMO transforms presta

Weblog:Madagascar development blogTracked:Apr 08, 06:02sfvscan implements seismic velocity analysis by scanning stacking velocities. This transformation is also known as the velocity transform or the hyperbolic Radon transform. The following example from bei/vela/vscan shows an example for transforming a C

Weblog:Madagascar development blogTracked:May 11, 12:25sfwiggle plots data using the traditional seismic method of wiggly traces. The following example from rsf/rsf/rsftour shows a typical output: Similarly to other plotting programs, there are multiple parameters that control the output. For examp

Weblog:Madagascar development blogTracked:Jun 12, 10:51sfwiggle plots data using the traditional seismic method of wiggly traces. The following example from rsf/rsf/rsftour shows a typical output: Similarly to other plotting programs, there are multiple parameters that control the output. For examp

Weblog:Madagascar development blogTracked:Jun 12, 10:59sftime2depth converts the input from vertical time to depth coordinates. The following example from rsf/su/rsflab9 shows a seismic image converted from time to depth by this transformation: The example is borrowed from John Stockwell's lecture not

Weblog:Madagascar development blogTracked:Jul 01, 22:28sftime2depth converts the input from vertical time to depth coordinates. The following example from rsf/su/rsflab9 shows a seismic image converted from time to depth by this transformation: The example is borrowed from John Stockwell's lecture not

Weblog:Madagascar development blogTracked:Jul 01, 22:29sfai2refl converts acoustic impedance to normal-incidence PP reflectivity using the simple equation R(t)=I(t+Δt)-I(t)I(t+Δt)+I(t) The program is useful for convolution modeling. The following example from rsf/rsf/wedge shows a classic e

Weblog:Madagascar development blogTracked:Aug 02, 15:22sfpatch breaks the input data into local windows or "patches", possibly with overlap. The patching technique is explained by Jon Claerbout in Nonstationarity: patching chapter from Image Estimation by Example. Suppose you have a 1-D signal with

Weblog:Madagascar development blogTracked:Sep 14, 12:33sfremap1 interpolates the input data to a different grid along the first axis. Here is an elementary example: making some data and interpolating it to a denser grid bash$ sfmath n1=5 o1=0 d1=1 output=x1 | sfdisfil 0: 0 1

Weblog:Madagascar development blogTracked:Nov 03, 06:22sfcausint implements an operation of causal numerical integration. This is a simple operation, which mathematically amounts to recursion yn=yn-1+xn or to inversion of a simple bidiagonal matrix. See Geophysical Image Estimation by Example for more expl

Weblog:Madagascar development blogTracked:Dec 01, 19:50sfinttest1 performs forward interpolation from a regular grid to irregular locations (in 1-D). The following example from sep/forwd/chirp shows regularly sampled values of a variable-frequency signal and the error of its interpolation using linear and

Weblog:Madagascar development blogTracked:Jan 12, 10:56