Normalization.hpp

Go to the documentation of this file.

/* ImLib3D
 * Copyright (c) 2001, ULP-IPB Strasbourg.
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or (at
 * your option) any later version.
 * 
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 */
#ifndef _Normalization_hpp
#define _Normalization_hpp

#include<ImLib3D/ImageProcessor.hpp>
#include<ImLib3D/ImageStats.hpp>
#include<ImLib3D/Signal1D.hpp>

namespace IP3D
{
    void NormalizeAverage(const Image3Df& src,const Image3Df& ref,Image3Df& res);
};

namespace IP3D
{
        void NormalizeAverageAndVariance(const Image3Df& src,const Image3Df& ref,Image3Df& res);
};



class JointHistogram
{
    //FIXME should not be public
public:
    int nbSrcBins; // on x-axis
    int nbRefBins; // on y-axis
    Image3Df histogramImage;
    double srcMin;
    double srcMax;
    double refMin;
    double refMax;


    inline float  SrcBinF(float v)  const   {return nbSrcBins  *(v-srcMin)/(srcMax-srcMin);}
    inline float  RefBinF(float v)  const   {return nbRefBins  *(v-refMin)/(refMax-refMin);}
    inline int    SrcBin(float v)   const   {return (int)rint(SrcBinF(v));}
    inline int    RefBin(float v)   const   {return (int)rint(RefBinF(v));}
    inline double SrcVal(int x)     const   {return (srcMin+x*(srcMax-srcMin)/(double)nbSrcBins );}
    inline double SrcVal(double x)  const   {return (srcMin+x*(srcMax-srcMin)/(double)nbSrcBins );}
    inline double RefVal(int y)     const   {return (refMin+y*(refMax-refMin)/(double)nbRefBins);}
    inline double RefVal(double y)  const   {return (refMin+y*(refMax-refMin)/(double)nbRefBins);}
    inline int    SSrcBin(float v)  const   {return min(nbSrcBins -1,max(0,SrcBin(v)));}
    inline int    SRefBin(float v)  const   {return min(nbRefBins-1 ,max(0,RefBin(v)));}
    float SrcLimit(float v) const {return min((float)srcMax,max(v,(float)srcMin));}
    float RefLimit(float v) const {return min((float)refMax,max(v,(float)refMin));}
protected:
    void ComputeFromIndependentSamples(const Image3Df &src,const Image3Df &ref);
    void ComputeFromDirectionalLinearInterpolation(const Image3Df &src,const Image3Df &ref,int dir);
    virtual void FindRanges(const Image3Df& src, const Image3Df &ref);
    JointHistogram(int _nbSrcBins=100,int _nbRefBins=400) :
        nbSrcBins(_nbSrcBins),nbRefBins(_nbRefBins),histogramImage(_nbSrcBins,_nbRefBins,1) {;}
public:
    virtual ~JointHistogram(){;}
public:
    void AddDisplayGrid();
    float GetSrcBinTotalProba(int x) const;
    Signal1Df MedianTransferFct() const ;
    Signal1Df MaximumProbabilityTransferFct() const;

    Signal1Df BestFitTransferFct(double smoothingCoef=1,bool allowDecreasing=false) const;

    static Signal1Df SmoothAndExtend(const Signal1Df &src,const Signal1Df &weights,double smoothSigma,double outlierSigma);
    Signal1Df ConfidenceWeights() const;


    JointHistogram(const string &filename);
    JointHistogram(const Image3Df &src,const Image3Df &ref,int _nbSrcBins=100,int _nbRefBins=400);
};

namespace IP3D
{
    void ComputeJointHistogramImage(const Image3Df& src,const Image3Df& ref,Image3Df& res,int drawTransferFunction=0,int nbSrcBins=100,int nbRefBins=200);
};


namespace IP3D
{
    void NormalizeJointHistogram(const Image3Df& src,const Image3Df& ref,Image3Df& res,int nbSrcBins=100,int nbRefBins=200);
};

#endif// _Normalization_hpp

---