core.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 MATHCORE_CORE_H
#define MATHCORE_CORE_H


#include "cimpl.h"
using namespace CIMPL;


#include <vector>
#include <cmath>

#define PI 3.141592653589793115997963468544185161590576171875

#define SQR(X)  ((X)*(X))
#define ROUND(X)  floor((X)+0.5)



namespace MathCore
{

	// Create class specializations


	// Eye() goes into the classes...
	//Eye(int N);


	
	
	// Core functions
	
	
	template< class T >
	void Plot(Vector<T>& v, int height = 30);

	template <class T>
	Vector<int> Find(Vector<T>& m);
	template <class T>
	Vector<int> Find(Matrix<T>& m);
	template <class T>
	void Find(Matrix<T>& m, Vector<int>& I, Vector<int>& J);
	
	void Ind2Sub(int rows, int cols, Vector<int>& ind, Vector<int>& I, Vector<int>& J);
	
	Vector<int> Sub2Ind(int rows, int cols, Vector<int>& I, Vector<int>& J);
	
	template <class T>
	Vector<T> Cross(Vector<T>& x, Vector<T>& y);
	template <class T>
	Matrix<T> Cross(Matrix<T>& x, Matrix<T>& y); 
	template <class T>
	Matrix<T> Cross(Matrix<T>& x, Matrix<T>& y, int dimension); 
	
	template <class T>
	Vector<T> CumSum(Vector<T>& m);
	template <class T>
	Matrix<T> CumSum(Matrix<T>& m);
	template <class T>
	Matrix<T> CumSum(Matrix<T>& m, int dimension);
	template <class T>
	Vector<T>& CumSumI(Vector<T>& m);
	template <class T>
	Matrix<T>& CumSumI(Matrix<T>& m);
	template <class T>
	Matrix<T>& CumSumI(Matrix<T>& m, int dimension);

	template <class T>
	Vector<T> CumProd(Vector<T>& m);
	template <class T>
	Matrix<T> CumProd(Matrix<T>& m);
	template <class T>
	Matrix<T> CumProd(Matrix<T>& m, int dimension);
	template <class T>
	Vector<T>& CumProdI(Vector<T>& m);
	template <class T>
	Matrix<T>& CumProdI(Matrix<T>& m);
	template <class T>
	Matrix<T>& CumProdI(Matrix<T>& m, int dimension);

	template <class T>
	Matrix<T> Diag(Vector<T>& m);
	template <class T>
	Matrix<T> Diag(Vector<T>& m, int k);
	template <class T>
	Vector<T> Diag(Matrix<T>& m);
	template <class T>
	Vector<T> Diag(Matrix<T>& m, int k);

	template <class T>
	Matrix<T> FlipLR(Matrix<T>& m);
	template <class T>
	Matrix<T>& FlipLRI(Matrix<T>& m);

	template <class T>
	Matrix<T> FlipUD(Matrix<T>& m);
	template <class T>
	Matrix<T>& FlipUDI(Matrix<T>& m);

	template <class T>
	Vector<T> Reverse(Vector<T>& m);
	template <class T>
	Vector<T>& ReverseI(Vector<T>& m);

	template <class T>
	Matrix<T> Rot90(Matrix<T>& m);

	template <class T>
	Matrix<T> Rot90(Matrix<T>& m, int k);

	//Kron();

	template <class T>
	Vector<T> Maximum(Matrix<T>& m);
	template <class T>
	T Maximum(Vector<T>& m);
	
	template <class T>
	Vector<T> Minimum(Matrix<T>& m);
	template <class T>
	T Minimum(Vector<T>& m);

	template <class T>
	Vector<double> Mean(Matrix<T>& m);
	template <class T>
	double Mean(Vector<T>& m);
	
	template <class T>
	inline void Swap(T& a, T& b)
	{
		T temp = a;
		a = b;
		b = temp;
	}

	template <class T>
	Vector<T> Median(Matrix<T>& m);
	template <class T>
	T Median(Vector<T>& m);

	template <class T>
	Vector<double> Std(Matrix<T>& m);
	template <class T>
	double Std(Vector<T>& m);

	template <class T>
	Vector<double> Var(Matrix<T>& m);
	template <class T>
	double Var(Vector<T>& m);


	template <class T>
	Vector<T> Sum(Matrix<T>& m);
	template <class T>
	Vector<T> Sum(Matrix<T>& m, int dimension);
	template <class T>
	T Sum(Vector<T>& m);

	
	template <class T>
	Vector<T> Prod(Matrix<T>& m);
	template <class T>
	Vector<T> Prod(Matrix<T>& m, int dimension);
	template <class T>
	T Prod(Vector<T>& m);

	template <class T>
	Vector<T> Diff(Vector<T>& m);
	template <class T>
	Matrix<T> Diff(Matrix<T>& m);

	template <class T>
	Matrix<T> RepMat(Matrix<T>& m, int M, int N);
	
	template <class T>
	Matrix<T> Reshape(Matrix<T>& m, int M, int N);
	template <class T>
	Matrix<T> Reshape(Vector<T>& m, int M, int N);
	
	template <class T>
	Vector<T>& QuickSort(Vector<T>& m);
	
	template <class T>
	Vector<T> Sort(Vector<T>& m);
	template <class T>
	Matrix<T> Sort(Matrix<T>& m, int dimension);
	template <class T>
	Matrix<T> Sort(Matrix<T>& m);
	template <class T>
	Vector<T>& SortI(Vector<T>& m);
	template <class T>
	Matrix<T>& SortI(Matrix<T>& m, int dimension);
	template <class T>
	Matrix<T>& SortI(Matrix<T>& m);
	//SortRows(); needed...
	//
	//
	//FreqSpace();
	Vector<double> LinSpace(double x, double y, int N);
	Vector<double> LinSpace(double x, double y);

	Vector<double> LogSpace(double x, double y, int N);
	Vector<double> LogSpace(double x, double y);

	template <class T>
	void MeshGrid(Vector<T>& x, Vector<T>& y, Matrix<T>& X, Matrix<T>& Y);
	
	
	Matrix<float>& RandN(Matrix<float>& m);
	Vector<float>& RandN(Vector<float>& m);
	Matrix<double>& RandN(Matrix<double>& m);
	Vector<double>& RandN(Vector<double>& m);




	int NextPow2(int n);
	

	template <class T>
	Vector<T> Abs(Vector<T>& m);
	template <class T>
	Vector<T>& AbsI(Vector<T>& m);

	template <class T>
	Matrix<T> Abs(Matrix<T>& m);
	template <class T>
	Matrix<T>& AbsI(Matrix<T>& m);
	
	//Abs();
	//Angle();

	// OTHER
	
	template <class T>
	Vector<int> Floor(Vector<T>& m);
	template <class T>
	Matrix<int> Floor(Matrix<T>& m);

	template <class T>
	Vector<int> Ceil(Vector<T>& m);
	template <class T>
	Matrix<int> Ceil(Matrix<T>& m);
	
	template <class T>
	Vector<int> Fix(Vector<T>& m);
	template <class T>
	Matrix<int> Fix(Matrix<T>& m);

	template <class T>
	Vector<int> Round(Vector<T>& m);
	template <class T>
	Matrix<int> Round(Matrix<T>& m);

	//template <class T>
	//Vector<double> Mod(Vector<T>& m);
	//template <class T>
	//Matrix<double> Mod(Matrix<T>& m);

	//template <class T>
	//Vector<double> Rem(Vector<T>& m);
	//template <class T>
	//Matrix<double> Rem(Matrix<T>& m);


	template <class T>
	Vector<T> UnitScale(Vector<T>& m);
	template <class T>
	Matrix<T> UnitScale(Matrix<T>& m);
	
	template <class T>
	Vector<T>& UnitScaleI(Vector<T>& m);
	template <class T>
	Matrix<T>& UnitScaleI(Matrix<T>& m);

	template <class T>
	Vector<T> ImageScale(Vector<T>& m);
	template <class T>
	Matrix<T> ImageScale(Matrix<T>& m);

	template <class T>
	Vector<T>& ImageScaleI(Vector<T>& m);
	template <class T>
	Matrix<T>& ImageScaleI(Matrix<T>& m);

	


	template <class T>
	Vector<int> Sign(Vector<T>& m);
	template <class T>
	Matrix<int> Sign(Matrix<T>& m);
	
	template <class T>
	Vector<T> Gradient(Vector<T>& m);
	template <class T>
	void Gradient(Matrix<T>& m, Matrix<T>& X, Matrix<T>& Y);

	template <class T>
	Matrix<T> Divergence(Matrix<T>& u, Matrix<T>& v);

	template <class T>
	Vector<T> Del2(Vector<T>& m);
	template <class T>
	Matrix<T> Del2(Matrix<T>& m);

	
	template <class T>
	Vector<T> CircShift(Vector<T>& m, int shift);
	template <class T>
	Matrix<T> CircShift(Matrix<T>& m, int rshift, int cshift);

	template <class T>
	Vector<T> FFTShift(Vector<T>& m);
	template <class T>
	Matrix<T> FFTShift(Matrix<T>& m);
	

	template <class T>
	Vector<T> IFFTShift(Vector<T>& m);
	template <class T>
	Matrix<T> IFFTShift(Matrix<T>& m);


	
	
	
	
	
	
	
	
};	
	
	
	
	
#include "Core.inl"










#endif