gszbdlg.cpp

大地测量高斯投影坐标计算。可以正算

// gszbDlg.cpp : 实现文件
//

#include "stdafx.h"
#include "gszb.h"
#include "gszbDlg.h"
#include ".\gszbdlg.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#endif


// 用于应用程序“关于”菜单项的 CAboutDlg 对话框

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

// 对话框数据
	enum { IDD = IDD_ABOUTBOX };

	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV 支持

// 实现
protected:
	DECLARE_MESSAGE_MAP()
};

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

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

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
END_MESSAGE_MAP()


// CgszbDlg 对话框



CgszbDlg::CgszbDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CgszbDlg::IDD, pParent)
	, io_L_Degree(0)
	, io_L_Minute(0)
	, io_L_Second(0)
	, io_B_Degree(0)
	, io_B_Minute(0)
	, io_B_Second(0)
	, io_x(0)
	, io_y(0)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CgszbDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT1, io_L_Degree);
	DDX_Text(pDX, IDC_EDIT2, io_L_Minute);
	DDX_Text(pDX, IDC_EDIT3, io_L_Second);
	DDX_Text(pDX, IDC_EDIT4, io_B_Degree);
	DDX_Text(pDX, IDC_EDIT5, io_B_Minute);
	DDX_Text(pDX, IDC_EDIT6, io_B_Second);
	DDX_Text(pDX, IDC_EDIT7, io_x);
	DDX_Text(pDX, IDC_EDIT8, io_y);
}

BEGIN_MESSAGE_MAP(CgszbDlg, CDialog)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	//}}AFX_MSG_MAP
	ON_BN_CLICKED(IDC_BUTTON_Reset, OnBnClickedButtonReset)
	ON_BN_CLICKED(IDC_RADIO_GtoE, OnBnClickedRadioGtoe)
	ON_BN_CLICKED(IDC_RADIO_EtoG, OnBnClickedRadioEtog)
	ON_BN_CLICKED(IDOK, OnBnClickedOk)
	ON_EN_SETFOCUS(IDC_EDIT7, OnEnSetfocusEdit7)
END_MESSAGE_MAP()


// CgszbDlg 消息处理程序

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

	// 将\“关于...\”菜单项添加到系统菜单中。

	// IDM_ABOUTBOX 必须在系统命令范围内。
	ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
	ASSERT(IDM_ABOUTBOX < 0xF000);

	CMenu* pSysMenu = GetSystemMenu(FALSE);
	if (pSysMenu != NULL)
	{
		CString strAboutMenu;
		strAboutMenu.LoadString(IDS_ABOUTBOX);
		if (!strAboutMenu.IsEmpty())
		{
			pSysMenu->AppendMenu(MF_SEPARATOR);
			pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
		}
	}

	// 设置此对话框的图标。当应用程序主窗口不是对话框时，框架将自动
	//  执行此操作
	SetIcon(m_hIcon, TRUE);			// 设置大图标
	SetIcon(m_hIcon, FALSE);		// 设置小图标

	// TODO: 在此添加额外的初始化代码
	io_L_Degree=0;
	io_L_Minute=0;
	io_L_Second=0;
	io_B_Degree=0;
	io_B_Minute=0;
	io_B_Second=0;
	io_x=0;
	io_y=0;


	return TRUE;  // 除非设置了控件的焦点，否则返回 TRUE
}

void CgszbDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialog::OnSysCommand(nID, lParam);
	}
}

// 如果向对话框添加最小化按钮，则需要下面的代码
//  来绘制该图标。对于使用文档/视图模型的 MFC 应用程序，
//  这将由框架自动完成。

void CgszbDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // 用于绘制的设备上下文

		SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);

		// 使图标在工作矩形中居中
		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;

		// 绘制图标
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

//当用户拖动最小化窗口时系统调用此函数取得光标显示。
HCURSOR CgszbDlg::OnQueryDragIcon()
{
	return static_cast<HCURSOR>(m_hIcon);
}

void CgszbDlg::OnBnClickedButtonReset()
{
	// TODO: 在此添加控件通知处理程序代码
	io_L_Degree=0;
	io_L_Minute=0;
	io_L_Second=0;
	io_B_Degree=0;
	io_B_Minute=0;
	io_B_Second=0;
	io_x=0;
	io_y=0;
	UpdateData(FALSE);

}

void CgszbDlg::OnBnClickedRadioGtoe()
{
	// TODO: 在此添加控件通知处理程序代码
	calculate_state=0;
	CEdit *m_Edit=(CEdit *)GetDlgItem(IDC_EDIT7);
	m_Edit->SetFocus();
	m_Edit->SetSel(0,-1);
}

void CgszbDlg::OnBnClickedRadioEtog()
{
	// TODO: 在此添加控件通知处理程序代码
	calculate_state=1;
	CEdit *m_Edit=(CEdit *)GetDlgItem(IDC_EDIT1);
	m_Edit->SetFocus();
	m_Edit->SetSel(0,-1);;
}

void CgszbDlg::OnBnClickedOk()
{
	// TODO: 在此添加控件通知处理程序代码
	UpdateData(TRUE);
	sec_L=io_L_Degree*3600.0+io_L_Minute*60.0+io_L_Second;
	sec_B=io_B_Degree*3600.0+io_B_Minute*60.0+io_B_Second;
	if(calculate_state==1)
		EllipsoidToGauss(sec_B,sec_L);
	if(calculate_state==0)
		if(io_y<1000000)
		{
			MessageBox("y坐标有误","出错了！",MB_ICONERROR|MB_OK);
			CEdit *m_Edit=(CEdit *)GetDlgItem(IDC_EDIT8);
			m_Edit->SetFocus();
			m_Edit->SetSel(0,-1);;
		}
		else
			GaussToEllipsoid(io_x,io_y);

}
int CgszbDlg::EllipsoidToGauss(double sec_B,double sec_L)
{
	double x,y,rad_B,l,deg_L0,N,a0,a3,a4,a5,a6;
	int n;
	const double sec_rho=206264.8062;
	n=int(sec_L/(6*3600))+1;//计算6度带投影带号
	deg_L0=6*n-3.0;			//计算中央经线经度
	rad_B=sec_B/sec_rho;
	l=(sec_L-deg_L0*3600)/sec_rho;
	N=6399698.902-(21562.267-(108.973-0.612*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B);
	a0=32140.404-(135.3302-(0.7092-0.0040*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B);
	a3=(0.3333333+0.001123*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B)-0.1666667;
	a4=(0.25+0.00252*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B)-0.04166;
	a5=0.0083-(0.1667-(0.1968+0.0040*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B))*cos(rad_B)*cos(rad_B);
	a6=(0.166*cos(rad_B)*cos(rad_B)-0.084)*cos(rad_B)*cos(rad_B);
	x=6367558.4969*sec_B/sec_rho-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(rad_B)*cos(rad_B);
	y=(1+(a3+a5*l*l)*l*l)*l*N*cos(rad_B);
	io_x=x;
	io_y=y+500000+n*1000000;
	UpdateData(FALSE);
	return 0;
}
int CgszbDlg::GaussToEllipsoid(double x,double Y)
{
	double y,deg_B,deg_L,rad_Bf,rad_beta,Z,Nf,b2,b3,b4,b5,deg_L0,sec_l;
	int n;
	int Degree,Minute;
	double Second;
	n=int(Y/1000000);
	y=Y-n*1000000-500000;
	deg_L0=6*n-3.0;
	const double sec_rho=206264.8062;
	rad_beta=x/6367558.4969;
	rad_Bf=rad_beta+(50221746+(293622+(2350+22*cos(rad_beta)*cos(rad_beta))*cos(rad_beta)*cos(rad_beta))*cos(rad_beta)*cos(rad_beta))*sin(rad_beta)*cos(rad_beta)*1e-10;
	Nf=6399698.902-(21562.267-(108.973-0.612*cos(rad_Bf)*cos(rad_Bf))*cos(rad_Bf)*cos(rad_Bf))*cos(rad_Bf)*cos(rad_Bf);
	Z=y/(Nf*cos(rad_Bf));
	b2=(0.5+0.003369*cos(rad_Bf)*cos(rad_Bf))*sin(rad_Bf)*cos(rad_Bf);
	b3=0.333333-(0.166667-0.001123*cos(rad_Bf)*cos(rad_Bf))*cos(rad_Bf)*cos(rad_Bf);
	b4=0.25+(0.16161+0.00562*cos(rad_Bf)*cos(rad_Bf))*cos(rad_Bf)*cos(rad_Bf);
	b5=0.2-(0.1667-0.0088*cos(rad_Bf)*cos(rad_Bf))*cos(rad_Bf)*cos(rad_Bf);
	deg_B=(rad_Bf*sec_rho-(1-(b4-0.12*Z*Z)*Z*Z)*Z*Z*b2*sec_rho)/3600;
	sec_l=(1-(b3-b5*Z*Z)*Z*Z)*Z*sec_rho;
	deg_L=deg_L0+sec_l/3600;
	Degree=int(deg_B);
	Minute=int ((deg_B-Degree)*60);
	Second=((deg_B-Degree)*60-Minute)*60;
	io_B_Degree=Degree;
	io_B_Minute=abs(Minute);
	if(Minute<0&&Degree==0)
		io_B_Minute=Minute;
	io_B_Second=abs(Second);
	if(Second<0&&Minute==0&&Degree==0)
		io_B_Second=Second;
	Degree=int(deg_L);
	Minute=int ((deg_L-Degree)*60);
	Second=((deg_L-Degree)*60-Minute)*60;
	io_L_Degree=Degree;
	io_L_Minute=abs(Minute);
	if(Minute<0&&Degree==0)
		io_L_Minute=Minute;
	io_L_Second=abs(Second);
	if(Second<0&&Minute==0&&Degree==0)
		io_L_Second=Second;
	UpdateData(FALSE);
	return 0;
}

void CgszbDlg::OnEnSetfocusEdit7()
{
	// TODO: 在此添加控件通知处理程序代码
}