geoinmodel.cpp

利用多面函数拟合法来拟合似大地水准面

// GeoinModel.cpp : Defines the class behaviors for the application.
//

#include "stdafx.h"
#include "GeoinModel.h"

#include "MainFrm.h"
#include "GeoinModelDoc.h"
#include "GeoinModelView.h"

#if ( defined(__BORLANDC__) || _MSC_VER <= 1000 ) && !defined( __GNUG__ )
#  include <stdio.h>
#  include <stdlib.h>
#  include <math.h>
#  include <iostream.h>
#  include <string.h>
#else
#  include <cmath>
#  include <cstdio>
#  include <cstdlib>
#  include <string>
#  include <iostream>
#endif


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

#define PI (4.*atan(1.0));

//中心化经度
double CenterL(double dl)
{
	return (dl-2.0748967728807)*100.0;
}
//中心化纬度
double CenterB(double db)
{
	return (db-0.444297153027155)*100.0;
}

//还原中心化经度
double ReturnCenterL(double dl)
{
	return (dl+2.0748967728807*100.0)/100.0;
}
//还原中心化纬度
double ReturnCenterB(double db)
{
	return (db+0.444297153027155*100.0)/100.0;
}

//弧度——dd.mmsssssss
double Arc2DMS(double dd)
{
	int iDegree,iMinute;
	double dSecond;

	dd=dd*180.0/PI;  //度
	iDegree=int(dd);
	iMinute=int((dd-iDegree)*60.0);
	dSecond=(dd-iDegree-iMinute/60.0)*3600.0;
	return iDegree + iMinute/100.0+ dSecond/10000.0;

}

// 十进制度 -- 弧度 
double int_radian(double intdeg)
{
	 return intdeg * 0.0174532925;
}

// 度分秒 -- 十进制度
double deg_int(double gms)
{
	 short g, m;	 
	 double m_gms, s, intdeg;

	 if(gms > -.000001 && gms < .000001) return 0;

	 g = (short)gms;
	 m_gms = gms - g;
	 m = short(100 * m_gms);
	 s = 100 * (100 * m_gms - m);

	 intdeg = g + m / 60.0 + s / 3600.0 + .000005;
	 return intdeg;
}


double myu(double a, double b) // Uniform Distribution
{
  double y;
  if(a>b) {
	printf("\nThe first parameter should be less than the second!");
	exit(1);
  }
  y = (double)rand()/(RAND_MAX);
  return (a+(b-a)*y); 
}

/////////////////////////////////////////////////////////////////////////////
// CGeoinModelApp

BEGIN_MESSAGE_MAP(CGeoinModelApp, CWinApp)
	//{{AFX_MSG_MAP(CGeoinModelApp)
	ON_COMMAND(ID_APP_ABOUT, OnAppAbout)
	//}}AFX_MSG_MAP
	// Standard file based document commands
	ON_COMMAND(ID_FILE_NEW, CWinApp::OnFileNew)
	ON_COMMAND(ID_FILE_OPEN, CWinApp::OnFileOpen)
	// Standard print setup command
	ON_COMMAND(ID_FILE_PRINT_SETUP, CWinApp::OnFilePrintSetup)
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CGeoinModelApp construction

CGeoinModelApp::CGeoinModelApp()
{
	// TODO: add construction code here,
	// Place all significant initialization in InitInstance
}

/////////////////////////////////////////////////////////////////////////////
// The one and only CGeoinModelApp object

CGeoinModelApp theApp;

/////////////////////////////////////////////////////////////////////////////
// CGeoinModelApp initialization

BOOL CGeoinModelApp::InitInstance()
{
	AfxEnableControlContainer();

	// Standard initialization
	// If you are not using these features and wish to reduce the size
	//  of your final executable, you should remove from the following
	//  the specific initialization routines you do not need.

#ifdef _AFXDLL
	Enable3dControls();			// Call this when using MFC in a shared DLL
#else
	Enable3dControlsStatic();	// Call this when linking to MFC statically
#endif

	// Change the registry key under which our settings are stored.
	// TODO: You should modify this string to be something appropriate
	// such as the name of your company or organization.
	SetRegistryKey(_T("Local AppWizard-Generated Applications"));

	LoadStdProfileSettings();  // Load standard INI file options (including MRU)

	// Register the application's document templates.  Document templates
	//  serve as the connection between documents, frame windows and views.

	CSingleDocTemplate* pDocTemplate;
	pDocTemplate = new CSingleDocTemplate(
		IDR_MAINFRAME,
		RUNTIME_CLASS(CGeoinModelDoc),
		RUNTIME_CLASS(CMainFrame),       // main SDI frame window
		RUNTIME_CLASS(CGeoinModelView));
	AddDocTemplate(pDocTemplate);

	// Parse command line for standard shell commands, DDE, file open
	CCommandLineInfo cmdInfo;
	ParseCommandLine(cmdInfo);

	// Dispatch commands specified on the command line
	if (!ProcessShellCommand(cmdInfo))
		return FALSE;

	// The one and only window has been initialized, so show and update it.
	m_pMainWnd->ShowWindow(SW_SHOW);
	m_pMainWnd->UpdateWindow();

	return TRUE;
}


/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
		// No message handlers
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAboutDlg)
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

// App command to run the dialog
void CGeoinModelApp::OnAppAbout()
{
	CAboutDlg aboutDlg;
	aboutDlg.DoModal();
}

/////////////////////////////////////////////////////////////////////////////
// CGeoinModelApp message handlers