correlation.cpp

基于小波的图像配准

/*This file is for computing correlation of two image data,also for 
  Mean of one image data and sum of square of image data*/

#include "correlation.h"

double correlation(double Sum_AxB,double Mean_A,double Mean_B,
				double Sum_SquareA,double Sum_SquareB,long dataSize)
{
	double corre;
	
	corre=(Sum_AxB-Mean_A*Mean_B*(double)dataSize)/
		   sqrt((Sum_SquareA-Mean_A*Mean_A*(double)dataSize)*
		        (Sum_SquareB-Mean_B*Mean_B*(double)dataSize));
	return corre;
}

void c_Mean_SumSquare(unsigned char* dataArray,double* Mean, 
				   double* Sum_Square,long dataSize)
{
	double mean=0, sumSquare=0;
	unsigned char* tmp;
	long i;

	tmp=dataArray;
	for(i=0;i<dataSize;i++){
		mean=mean+(double)*tmp;
		sumSquare=sumSquare+(double)(*tmp)*(double)(*tmp);
		tmp++;
	}
	mean=mean/(double)dataSize;
	*Mean=mean;
	*Sum_Square=sumSquare;

	return;
}

void d_Mean_SumSquare(double* dataArray,double* Mean, 
				   double* Sum_Square,long dataSize)
{
	double mean=0, sumSquare=0;
	double* tmp;
	long i;

	tmp=dataArray;
	for(i=0;i<dataSize;i++){
		mean = mean + *tmp;
		sumSquare = sumSquare + (*tmp) * (*tmp);
		tmp++;
	}
	mean=mean/(double)dataSize;
	*Mean=mean;
	*Sum_Square=sumSquare;

	return;
}