SURFLIB_STATS_GRID - Calculate statistics of binned data

Description:

Calculate the statistics of the binned data and return the median and upper and lower limits of acceptable data.

Invocation:

CALL SURFLIB_STATS_GRID(SMODE, NX, NY, NMAX, NSIGMA, IPOS, JPOS, BINS, PNTS, STATS, STATUS)

Arguments:

SMODE = CHARACTER (Given)

Smoothing mode. The statistics can be smoothed by adjacent pixels in order to remove spikes from the clipping envelope. Options are: NONE - No smoothing (ie the statistics of each bin) HANN - Hanning smoothing. Triangular envolope across 3 points. Note that smoothing depends on the way the grid was unwrapped. since that determines which pixels are adjacent (except for SMODE=NONE.

NX = INTEGER (Given)

Size of X dimension

NY = INTEGER (Given)

Size of Y dimension

NMAX = INTEGER (Given)

Maximum value allowed for third dimension of BINS

NSIGMA = REAL (Given)

Standard deviation limit.

IPOS(NX $*$ NY) = INTEGER (Given)

I coordinate for each pixel

JPOS(NX $*$ NY) = INTEGER (Given)

J coordinate for each pixel

BINS(NX, NY, NMAX) = REAL (Given)

The data stored in relation to its position

PNTS(NMAX) = REAL (Given)

Scratch space for copying in the data from each I,J

STATS(NX, NY, 3) = REAL (Returned)

Statistics of each bin. 3 components are: Median, Upper limit lower limit.

STATUS = INTEGER (Given & Returned)

Global status.

Copyright

Copyright ©1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.