derivatemi3d.h

3D reconstruction, medical image processing from colons, using intel image processing for based clas

// DerivateMI3D.h: interface for the RxDerivateMI3D class.
//
//////////////////////////////////////////////////////////////////////
//
//	Title: Header file for Calculating of the Derivative of the MI
//
//////////////////////////////////////////////////////////////////////
//
//	Author: Helen Hong, 3DMed co. LTD
//	138-dong 417-ho Seoul National Univ.
//	San 56-1 Shinlim-dong Kwanak-gu, Seoul Korea
//	Email. [e-mail]
//
//	Date	: 2002. 9. 10.
//	Update	: 2002. 10. 8.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_DERIVATEMI3D_H__8AE66687_577C_4031_A059_0AC48A97C882__INCLUDED_)
#define AFX_DERIVATEMI3D_H__8AE66687_577C_4031_A059_0AC48A97C882__INCLUDED_

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000

#include "Transform3DInfo.h"

class RxVolumeInfo;
class RxMatrix3D;
class RxAutoSample3DInfo;

class RxDerivateMI3D  
{
public:
	RxDerivateMI3D();
	virtual ~RxDerivateMI3D();

protected:
	unsigned short		*m_pnVolume;			// reference volume data
	int					m_SampleNum;			// sample number
	double				m_varFloat,				// variance of float volume
						m_covarRefer,			// covariance of reference volume
						m_covarFloat;			// covariance of float volume
	RxMatrix3D			*m_pMatrix3D;			// matrix calculation
	RxVolumeInfo		*m_pResol;				// reference volume information
	RxTransform3DInfo	m_Trans;				// transformation
	RxAutoSample3DInfo	*sampleA;				// automatic sampling

public:
	BOOL				SetVolumeData(unsigned short *Buf);
	BOOL				SetReferVolInfo(RxVolumeInfo ReferInfo);
	BOOL				SetSampleNumber(int number);
	BOOL				SetSampleA(RxAutoSample3DInfo *sampleA);	
	BOOL				SetTransform(RxTransform3DInfo *Trans);
	BOOL				SetMIParameters(double varFloat,double covarRefer,double covarFloat);
	BOOL				Calculate(RxAutoSample3DInfo *sampleB,double (*totalfactor)[6]);

protected:
	BOOL				DerivateTrans3D(RxAutoSample3DInfo vi,RxAutoSample3DInfo vj,double (*result)[6]);
	double				Gauss_point(unsigned short v1,unsigned short v2,double variance);
	double				Gauss_pair(unsigned short u1,unsigned short v1,unsigned short u2, unsigned short v2,double cov1,double cov2);
	unsigned short		GetIntensity3D(int x,int y,int z);
};

#endif // !defined(AFX_DERIVATEMI3D_H__8AE66687_577C_4031_A059_0AC48A97C882__INCLUDED_)