STL like template based coding with the MMX/SSE extension

#ifdef _CH_
#define WIN32
#error "The file needs cvaux, which is not wrapped yet. Sorry"
#endif



#include "vigra/bordertreatment.hxx"
#include "vigra/stdconvolution.hxx"
#include "vigra/recursiveconvolution.hxx"
#include "vigra/multi_array.hxx" 
#include "vigra/imageiterator.hxx"
#include "vigra/stdimage.hxx"
#include "vigra/convolution.hxx"
//#include "vigra/fftw.hxx"
#include "vigra/fftw3.hxx"
#include "vigra/basicimage.hxx"

#include "vigra/inspectimage.hxx" 
#include "vigra/transformimage.hxx"


#include "stllcv/iplvfftw.hxx"
#include "stllcv/vimageoperation.hxx"


#include "stllcv/iplvbasicimageoperation.hxx"

#include "stllcv/iplvbasicimage.hxx"


//#include "iplvbasicimageex.hxx"

//#include "stllcv/iplvhistogram.hxx"
               //iplvhistogram.hxx

#include "stllcv/iplvinspectimage.hxx"

#include "vigra/localminmax.hxx"
//#include <ipl.h>



#ifndef _EiC
#include "cv.h"
#include "highgui.h"

#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>
#include <math.h>
#include <float.h>
#include <limits.h>
#include <time.h>
#include <ctype.h>
#endif

#ifdef _EiC
#define WIN32
#endif

//#include "..\cookroom\macros.inc"
//#include "..\cookroom\cookroom.h"
//#include "..\cookroom\cookdraw.h"


//#include <ipl.h>
//#include "cv.h"


static const bool is_modal = true;  /// ipDraw parameter

//static char msg[256];
//static __int64 start, stop;


#define DebugFileRead 1
#define DebugShowImage 1
#define SelfResize 0
//#define PhaseOnly 0
#define SMOOTHBorder 1

#define SMOOTHPhase 1
#define SMOOTHPhaseAverage 0

//#define PEXFFT char
//#define PEXFFT float
#define PEXFFT float
//#define PEXORGIMG unsigned short
//#define PEXORGIMG int
//#define PEXORGIMG float
#define PEXORGIMG unsigned char
//#define PEXORGIMG double
//#define PEXORGIMG char

//#define TileLen 100
#define TileLen 50






#if !DebugFileRead
#define Filename argv[1]
#define Filename2 argv[2] 
#define OutImageName argv[3] 
#endif

using namespace stllcv; 
int main( int argc, char** argv )
{

	const char* input_name;
	

	input_name =0;

#if DebugFileRead
    const char* Filename = input_name ? input_name : (char*)"lena.jpg";
	const char* Filename2 = input_name ? input_name : (char*)"lena2.jpg";
#endif


	CiplvBasicImage<PEXORGIMG> image1org(Filename);
	CiplvBasicImage<PEXORGIMG> image1(Filename);
	CiplvBasicImage<PEXORGIMG> image2org(Filename2);
	CiplvBasicImage<PEXORGIMG> image2(Filename2);

	int width1=image1org.width();
	int height1=image1org.height();
	CiplvBasicImage<PEXORGIMG> imageCorrelation(width1,height1);
	int xmax, ymax;

#if DebugShowImage
		bool is_modal = true;  /// ipDraw parameter
		char winName[]="srcImg";
		cvNamedWindow( winName, 1 );
		showIplvBasicImag<PEXORGIMG>(&image1,winName);	
		showIplvBasicImag<PEXORGIMG>(&image2,winName);	
#endif	
	phaseOnlyCorrelation(image1, image2,imageCorrelation , xmax, ymax , true, 6.0 );




	std::cout << "xmax " << xmax << " ymax " << ymax << std::endl;


	CiplvBasicImage<PEXORGIMG> image3(width1,height1);
	iplvTileImageNoCycle<PEXORGIMG>(&image1org,&image2org,				   
				   &image3,
				   TileLen,xmax,ymax
				   );

#if DebugShowImage
		showIplvBasicImag<PEXORGIMG>(&image1,winName);	
		showIplvBasicImag<PEXORGIMG>(&image2,winName);	
		showIplvBasicImag<PEXORGIMG>(&imageCorrelation,winName);
		showIplvBasicImag<PEXORGIMG>(&image3,winName);	
		cvDestroyWindow("srcImg");
#endif

}