analysis.h

图像分割算法

//Copyright (c) 2004-2005, Baris Sumengen
//All rights reserved.
//
// CIMPL Matrix Performance Library
//
//Redistribution and use in source and binary
//forms, with or without modification, are
//permitted provided that the following
//conditions are met:
//
//    * No commercial use is allowed. 
//    This software can only be used
//    for non-commercial purposes. This 
//    distribution is mainly intended for
//    academic research and teaching.
//    * Redistributions of source code must
//    retain the above copyright notice, this
//    list of conditions and the following
//    disclaimer.
//    * Redistributions of binary form must
//    mention the above copyright notice, this
//    list of conditions and the following
//    disclaimer in a clearly visible part 
//    in associated product manual, 
//    readme, and web site of the redistributed 
//    software.
//    * Redistributions in binary form must
//    reproduce the above copyright notice,
//    this list of conditions and the
//    following disclaimer in the
//    documentation and/or other materials
//    provided with the distribution.
//    * The name of Baris Sumengen may not be
//    used to endorse or promote products
//    derived from this software without
//    specific prior written permission.
//
//THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT
//HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
//EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT
//NOT LIMITED TO, THE IMPLIED WARRANTIES OF
//MERCHANTABILITY AND FITNESS FOR A PARTICULAR
//PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
//CONTRIBUTORS BE LIABLE FOR ANY
//DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
//EXEMPLARY, OR CONSEQUENTIAL DAMAGES
//(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT
//OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
//DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
//HOWEVER CAUSED AND ON ANY THEORY OF
//LIABILITY, WHETHER IN CONTRACT, STRICT
//LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
//OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
//OF THIS SOFTWARE, EVEN IF ADVISED OF THE
//POSSIBILITY OF SUCH DAMAGE.



#pragma once
#ifndef ANALYSIS_H
#define ANALYSIS_H

#include "cimpl.h"
using namespace CIMPL;

#include "./IO.h"


#include <time.h>

using namespace MathCore;
using namespace IO;


/// \brief Analysis Toolbox
namespace Analysis
{
	void StatusDisplay(long status);

	
	enum Border
	{
		ZeroPad,
		Symmetric,
		Replicate,
		Circular
	};

	enum FilterArea
	{
		Full,
		Same,
		Valid
	};

	Vector<ComplexFloat> FFT(Vector<ComplexFloat>& x);
	Vector<ComplexFloat>& FFTI(Vector<ComplexFloat>& x);
	Vector<ComplexDouble> FFT(Vector<ComplexDouble>& x);
	Vector<ComplexDouble>& FFTI(Vector<ComplexDouble>& x);
	Vector<ComplexFloat> FFT(Vector<float>& x);
	Vector<ComplexDouble> FFT(Vector<double>& x);

	Vector<ComplexFloat> IFFT(Vector<ComplexFloat>& x);
	Vector<ComplexFloat>& IFFTI(Vector<ComplexFloat>& x);
	Vector<ComplexDouble> IFFT(Vector<ComplexDouble>& x);
	Vector<ComplexDouble>& IFFTI(Vector<ComplexDouble>& x);

	Matrix<ComplexFloat> FFT2(Matrix<ComplexFloat>& x);
	Matrix<ComplexFloat>& FFT2I(Matrix<ComplexFloat>& x);
	Matrix<ComplexDouble> FFT2(Matrix<ComplexDouble>& x);
	Matrix<ComplexDouble>& FFT2I(Matrix<ComplexDouble>& x);

	Matrix<ComplexFloat> IFFT2(Matrix<ComplexFloat>& x);
	Matrix<ComplexFloat>& IFFT2I(Matrix<ComplexFloat>& x);
	Matrix<ComplexDouble> IFFT2(Matrix<ComplexDouble>& x);
	Matrix<ComplexDouble>& IFFT2I(Matrix<ComplexDouble>& x);


	Vector<float> Conv(Vector<float>& input, Vector<float>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);
	Vector<double> Conv(Vector<double>& input, Vector<double>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);
	Vector<ComplexFloat> Conv(Vector<ComplexFloat>& input, Vector<ComplexFloat>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);
	Vector<ComplexDouble> Conv(Vector<ComplexDouble>& input, Vector<ComplexDouble>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);


	Matrix<float> Conv2(Matrix<float>& input, Matrix<float>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
	Matrix<float> Conv2(Vector<float>& filter1, Vector<float>& filter2, Matrix<float>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
	Matrix<double> Conv2(Matrix<double>& input, Matrix<double>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
	Matrix<double> Conv2(Vector<double>& filter1, Vector<double>& filter2, Matrix<double>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
	Matrix<ComplexFloat> Conv2(Matrix<ComplexFloat>& input, Matrix<ComplexFloat>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
	Matrix<ComplexFloat> Conv2(Vector<ComplexFloat>& filter1, Vector<ComplexFloat>& filter2, Matrix<ComplexFloat>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
	Matrix<ComplexDouble> Conv2(Matrix<ComplexDouble>& input, Matrix<ComplexDouble>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
	Matrix<ComplexDouble> Conv2(Vector<ComplexDouble>& filter1, Vector<ComplexDouble>& filter2, Matrix<ComplexDouble>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);


	//FFTN();
	//IFFTN();
	//
	//ConvN();
	//Deconv();
	//Filter();
	//Filter2();
	//
	//Cov();
	//Corr2();
	//CorrCoef();
	//XCorr();
	//XCorr2();
	

};







#endif