MisrRegression.h

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                              MisrRegression                               =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

            Copyright 2005, California Institute of Technology.
                           ALL RIGHTS RESERVED.
                 U.S. Government Sponsorship acknowledged.

============================================================================*/

#ifndef MISRREGRESSION_H
#define MISRREGRESSION_H

#include "MisrUtil.h"
#include "MisrError.h"
#include "MisrProjParam.h"
#include "MisrMapQuery.h"

typedef struct {
       MTKt_DataBuffer valid_mask;
       MTKt_DataBuffer slope;
       MTKt_DataBuffer intercept;
       MTKt_DataBuffer correlation;
} MTKt_RegressionCoeff;

#define MTKT_REGRESSION_COEFF_INIT { \
    MTKT_DATABUFFER_INIT, \
    MTKT_DATABUFFER_INIT, \
    MTKT_DATABUFFER_INIT, \
}

MTKt_status MtkRegressionCoeffAllocate(
  int nline,   
  int nsample, 
  MTKt_RegressionCoeff *regressbuf 
);

MTKt_status MtkRegressionCoeffFree(
  MTKt_RegressionCoeff *regressbuf 
);

MTKt_status MtkLinearRegressionCalc(
  int Size,   
  const double *X, 
  const double *Y, 
  const double *YSigma, 
  double *A, 
  double *B, 
  double *Correlation 
);

MTKt_status MtkSmoothData(
  const MTKt_DataBuffer *Data, 
  const MTKt_DataBuffer *Valid_mask, 
  int Width_line, 
  int Width_sample, 
  MTKt_DataBuffer *Data_smoothed 
);

MTKt_status MtkRegressionCoeffCalc(
  const MTKt_DataBuffer *Data1, 
  const MTKt_DataBuffer *Valid_mask1, 
  const MTKt_DataBuffer *Data2, 
  const MTKt_DataBuffer *Data2_sigma, 
  const MTKt_DataBuffer *Valid_mask2, 
  const MTKt_MapInfo *Map_info, 
  int Regression_size_factor, 
  MTKt_RegressionCoeff *Regression_coeff, 
  MTKt_MapInfo *Regression_map_info 
) ;

MTKt_status MtkResampleRegressionCoeff(
  const MTKt_RegressionCoeff *Regression_coeff, 
  const MTKt_MapInfo *Regression_coeff_map_info, 
  const MTKt_MapInfo *Target_map_info, 
  MTKt_RegressionCoeff *Regression_coeff_out 
) ;

MTKt_status MtkApplyRegression(
  const MTKt_DataBuffer *Source, 
  const MTKt_DataBuffer *Source_mask, 
  const MTKt_MapInfo    *Source_map_info, 
  const MTKt_RegressionCoeff *Regression_coeff, 
  const MTKt_MapInfo *Regression_coeff_map_info, 
  MTKt_DataBuffer *Regressed,  
  MTKt_DataBuffer *Regressed_mask 
);

MTKt_status MtkDownsample(
  const MTKt_DataBuffer *Source, 
  const MTKt_DataBuffer *Source_mask, 
  int Size_factor, 
  MTKt_DataBuffer *Result, 
  MTKt_DataBuffer *Result_mask 
              );

MTKt_status MtkUpsampleMask(
  const MTKt_DataBuffer *Source_mask, 
  int Size_factor, 
  MTKt_DataBuffer *Result_mask 
            );

#endif /* MISRREGRESSION_H */