registration.h

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

// Registration.h: interface for the RxxRegistration class.
//
//////////////////////////////////////////////////////////////////////

#if !defined(AFX_REGISTRATION_H__BE3A01E0_0423_4C4C_A781_857758270AF0__INCLUDED_)
#define AFX_REGISTRATION_H__BE3A01E0_0423_4C4C_A781_857758270AF0__INCLUDED_

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

#include "RxxPowellMethod.h"

class RxTransform3DInfo;

#include "GlobalRegistPar.h"
#include "JJDistanceMap.h"
#include "JJSkeleton.h"
#include "RxxOpticalFlow.h"
#include "Rxx3DTPS.h"

class JJBinaryVolume;

class RxxRegistration
{
private:
	RxxOpticalFlow	m_OpticalFlow;

	JJBinaryVolume	m_RefBinaryVolume;
	JJBinaryVolume	m_FloatBinaryVolume;
	JJDistanceMap	m_RefDistanceMap;
	JJSkeleton		m_Skeleton;

	double			m_lfBoundingBoxX;
	double			m_lfBoundingBoxY;

	BOOL			toggle;

	int				m_iDilationLow;
	int				m_iDilationHigh;

	RxMatrix4D		m_mxRegistration;
	RxMatrix4D		m_mxInverse;

	Rxx3DTPS		m_3DTPS_R;
	Rxx3DTPS		m_3DTPS_L;

	double			m_lfCPX;
	double			m_lfCPY;
	double			m_lfCPZ;

	BOOL			m_bREFSegR_RealR;
	BOOL			m_bFLOATSegR_RealR;

	unsigned short* pRefData;

public:
    BOOL			m_isBOUNDING_BOX_Z_REGISTRATION;
	BOOL			m_isBOUNDING_BOX_XY_REGISTRATION;
	BOOL			m_isREGISTRATION;

	float			m_fDistanceMapTime[3];
	float			m_fRegistrationTime[3];
	float			m_fIRTime[3];
	float			m_fNRTime[3];
	float			m_fSurfaceAccuracy[4];
	float			m_fSurfaceAccuracyAfterNR[3];
	float			m_fNoduleError[3];
	float			m_fNoduleAccuracyBeforeNR[3];
	float			m_fNoduleAccuracy[4];

	float			m_fCPDistFirst[3][2];
	float			m_fCPDistIR[3][2];
	float			m_fCPDistRIG[3][2];
	float			m_fCPDistTPS[3][2];

	float			m_fEachCPDistFirst[3][20];
	float			m_fEachCPDistIR[3][20];
	float			m_fEachCPDistRIG[3][20];
	float			m_fEachCPDistTPS[3][20];

	short			m_nTPS_CP_R[15];
	short			m_nTPS_CP_L[15];
	short			m_nEval_CP_R[5];
	short			m_nEval_CP_L[5];

	double			POST[20];
	double			PRE[20];

	double			RCP[120];
	double			LCP[120];
	double			RtApex[6];
	double			LtApex[6];
	double			RtApexAuto[6];
	double			LtApexAuto[6];
	double			RH[6];
	double			LH[6];
	double			RCen[6];
	double			LCen[6];

	double			m_dbTime[3][4][2];

private:
	void InitializeParameter(void);
	void InitializeData(short nMode);
	void PrintRegResult(RxTransform3DInfo* pTransform);
	void Dilation(int m_nDilationMode);
	void SetRegistrationMatrix(RxTransform3DInfo* pTransform);
	void LoadSegmentVolume(void);

	void Warp3DAfterIR_Right(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL);
	void Warp3DAfterIR_Left(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL);
	void AverageCPDistAfterTPS(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL, double* dbRight, double* dbLeft);
	void VesselCenterline(void);

	RxMatrix4D MakeRotationMatrix(_FPoint Rot);

	void Smoothing(short nLevel);
	double AverageCPDistance(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL, short nStep,  double* dbRight, double* dbLeft);

protected:
	RxxPowellMethod *m_pPM;

public:
	_MODAL_INFO* m_pFloat;
	_MODAL_INFO* m_pRefer;

	BOOL InitialRegistration(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL);
	BOOL RigidReg3DofDistanceMap(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL);

	BOOL RigidReg3DofCC(RxTransform3DInfo* pTransform);
	RxxRegistration();
	virtual ~RxxRegistration();	

	void SetPatientData(void);
	void SetControlPoint(void);

	void MakeRefBinaryVolumeFromSeg(short nMode);
	void MakeFloatBinaryVolumeFromSeg(short nMode);

	void SubtractionWithOutReg(void);
	void RigidRegSubtraction(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL);
	void OpticalFlow(RxTransform3DInfo* pTransform);
	void Subtraction(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL, BOOL bReg);
	void SubtractionForTPS(void);

	void CalTPSAfterIR(RxTransform3DInfo* pTransformR, RxTransform3DInfo* pTransformL, BOOL bWarp);

	void PrintPerfusionResult(void);

	short	*	m_pFloatEdgeListX;
	short	*	m_pFloatEdgeListY;
	short	*	m_pFloatEdgeListZ;
	unsigned long	m_nNumofFloatEdge;

	void MakeColorTable(int r1, int r2, int r3, int r4);

	void PrintHistogram(short nMode);
	void EdgeDilation(int nLevel);
	void SaveVolumeDicom(void);
};

#endif // !defined(AFX_REGISTRATION_H__BE3A01E0_0423_4C4C_A781_857758270AF0__INCLUDED_)