relativitydlg.cpp

对不同图像之间的相关性可以进行较好的比较.对后续的图像分析很有帮助.

// RelativityDlg.cpp : implementation file
//

#include "stdafx.h"
#include "Relativity.h"
#include "RelativityDlg.h"
#include "math.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CRelativityDlg dialog

CRelativityDlg::CRelativityDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CRelativityDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CRelativityDlg)
	m_Msg = _T("");
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_bCanAnalysis=false;
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CRelativityDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CRelativityDlg)
	DDX_Text(pDX, IDC_EDIT1, m_Msg);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CRelativityDlg, CDialog)
	//{{AFX_MSG_MAP(CRelativityDlg)
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_DATA1, OnData1)
	ON_BN_CLICKED(IDC_DATA2, OnData2)
	ON_BN_CLICKED(IDC_ANALYSIS, OnAnalysis)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CRelativityDlg message handlers

BOOL CRelativityDlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	// Set the icon for this dialog.  The framework does this automatically
	//  when the application's main window is not a dialog
	SetIcon(m_hIcon, TRUE);			// Set big icon
	SetIcon(m_hIcon, FALSE);		// Set small icon
	
	// TODO: Add extra initialization here
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.

void CRelativityDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

		// Center icon in client rectangle
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon + 1) / 2;
		int y = (rect.Height() - cyIcon + 1) / 2;

		// Draw the icon
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

// The system calls this to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CRelativityDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CRelativityDlg::OnData1() 
{
	CString FileName="";
	CFileDialog dlg(TRUE,"dat","*.dat",OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,"信号数据文件(*.dat)|*.dat|所有文件(*.*)|*.*||",NULL);
	if(dlg.DoModal()==IDOK)
	{
		m_bCanAnalysis=true;
		FileName=dlg.GetPathName();
		CFile file;
		file.Open(FileName,CFile::modeReadWrite);
		buf1=new char [file.GetLength()];
		file.Read(buf1,file.GetLength());
		m_nData1Len=file.GetLength();
		file.Close();
	}
	else
	{
		m_bCanAnalysis=false;
	}
}

void CRelativityDlg::OnData2() 
{
	if(m_bCanAnalysis==false)
	{
		AfxMessageBox("请先输入第一序列信号数据！");
		return;
	}
	CString FileName="";
	CFileDialog dlg(TRUE,"dat","*.dat",OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,"信号数据文件(*.dat)|*.dat|所有文件(*.*)|*.*||",NULL);
	if(dlg.DoModal()==IDOK)
	{
		m_bCanAnalysis=true;
		FileName=dlg.GetPathName();
		CFile file;
		file.Open(FileName,CFile::modeReadWrite);
		buf2=new char [file.GetLength()];
		file.Read(buf2,file.GetLength());
		m_nData2Len=file.GetLength();
		file.Close();
	}
	else
	{
		m_bCanAnalysis=false;
	}
}

void CRelativityDlg::OnAnalysis() 
{
	if(m_bCanAnalysis==false)
	{
		AfxMessageBox("请先选择信号数据！");
		return;
	}
	int N=m_nData1Len>m_nData2Len?m_nData2Len:m_nData1Len;
	
	float A,B,C,Pxy;
	A=B=C=Pxy=0.0f;
	
	for(int i=0;i<N;i++)
	{
		A+=buf1[i]*buf2[i];
		B+=buf1[i]*buf1[i];
		C+=buf2[i]*buf2[i];
	}

	Pxy=A/(sqrt(B*C));
	
	m_Msg.Format("%f",Pxy);
	UpdateData(FALSE);
	m_bCanAnalysis=false;

	delete[] buf1;
	delete[] buf2;	
}