PhaseOnlyCorrelation.cpp

Go to the documentation of this file.

#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"





static const bool is_modal = true;  

//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


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;  
                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

}

---