/** @file psImageStats.h
*  \brief Routines for calculating statistics on images.
*  @ingroup ImageStats
*
*  This file will hold the prototypes for procedures which calculate
*  statistic on images, histograms on images, and fit/evaluate Chebyshev
*  polynomials to images.
*
*  @author GLG, MHPCC
*
*  @version $Revision: 1.19 $ $Name: not supported by cvs2svn $
*  @date $Date: 2004-12-02 21:12:52 $
*
*  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
*/
#if !defined(PS_IMAGE_STATS_H)
#define PS_IMAGE_STATS_H

#include "psType.h"
#include "psVector.h"
#include "psImage.h"
#include "psStats.h"
#include "psFunctions.h"

/// @addtogroup ImageStats
/// @{

/** This routine must determine the various statistics for the image.
 *
 *  Determine statistics for image (or subimage). The statistics to be 
 *  determined are specified by stats. The mask allows pixels to be excluded 
 *  if their corresponding mask pixel value matches the value of maskVal. 
 *  This function must be defined for the following types: psS8, psU16, psF32, 
 *  psF64.
 *
 *  @return psStats*    the resulting statistics result(s)
 */
psStats* psImageStats(
    psStats* stats,                    ///< defines statistics to be calculated
    psImage* in,                       ///< image (or subimage) to calculate stats
    psImage* mask,                     ///< mask data for image (NULL ok)
    psS32 maskVal                        ///< mask Mask for mask
);

/** Construct a histogram from an image (or subimage).
 *
 *  The histogram to generate is specified by psHistogram hist (see section 
 *  4.3.2 in SDRS). This function must be defined for the following types: 
 *  psS8, psU16, psF32, psF64.
 *
 *  @return psHistogram*     the resulting histogram
 */
psHistogram* psImageHistogram(
    psHistogram* out,                  ///< input histogram description & target
    psImage* in,                       ///< Image data to be histogramed.
    psImage* mask,                     ///< mask data for image (NULL ok)
    psU32 maskVal               ///< mask Mask for mask
);

/** Fit a 2-D polynomial surface to an image.
 *
 *  The input structure coeffs contains the desired order and terms of 
 *  interest. This function must be defined for the following types: psS8, 
 *  psU16, psF32, psF64.
 *
 *  @return psPolynomial2D*     fitted polynomial result
 *
 */
psPolynomial2D* psImageFitPolynomial(
    psPolynomial2D* coeffs,            ///< coefficient structure carries in desired terms & target
    const psImage* input
);

/** Evaluate a 2-D polynomial surface for the image pixels.
 *
 *  Given the input polynomial coefficients, set the image pixel values on the 
 *  basis of the polynomial function. This function must be defined for the 
 *  following types: psS8, psU16, psF32, psF64.
 *
 *  @return psImage*    the resulting image
 */
psImage* psImageEvalPolynomial(
    psImage* input,                    ///< input image
    const psPolynomial2D* coeffs       ///< coefficient structure carries in desired terms
);

/// @}

#endif
