albersdlg.cpp

albers投影转换代码,使用vc开发的com

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

#include "stdafx.h"
#include "albers.h"
#include "albersDlg.h"
#include "math.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()


// CalbersDlg 对话框




CalbersDlg::CalbersDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CalbersDlg::IDD, pParent)
	, PI(3.1415926535897)
	, dLat(0)
	, dPrecision(0)
	, dX(0)
	, dY(0)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
	DA=6378245;
	DE2=0.006693422;
}

void CalbersDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT_LATITUDE, dLat);
	DDX_Text(pDX, IDC_EDIT_PRECISION, dPrecision);
	DDX_Text(pDX, IDC_EDIT_X, dX);
	DDX_Text(pDX, IDC_EDIT_Y, dY);
}

BEGIN_MESSAGE_MAP(CalbersDlg, CDialog)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	//}}AFX_MSG_MAP
	ON_BN_CLICKED(IDC_BUTTON_TRANSForm, &CalbersDlg::OnBnClickedButtonTransform)
END_MESSAGE_MAP()


// CalbersDlg 消息处理程序

BOOL CalbersDlg::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: 在此添加额外的初始化代码

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

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

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

void CalbersDlg::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 CalbersDlg::OnQueryDragIcon()
{
	return static_cast<HCURSOR>(m_hIcon);
}


void CalbersDlg::OnBnClickedButtonTransform()
{
	// TODO: 在此添加控件通知处理程序代码
	UpdateData(true);
	U0=100;
	double a=0.0,C=0.0,Ps=0.0;//,P1,P2;
	double r1,r2,F1,F2,P,D;//r1第一纬线半径，r2第二纬线半径，F1第一纬线体形面面积，F2第二纬线体形面面积，P要计算经纬度的极半径，D要计算经纬度的极角
	r1=computeR(PI*25/180);
	r2=computeR(PI*47/180);
	F1=computeF(PI*25/180);
	F2=computeF(PI*47/180);
	a=(pow(r1,2)-pow(r2,2))/(2*(F2-F1));//0.577108986
	C=(pow(r1,2)*F2-pow(r2,2)*F1)/(pow(r1,2)-pow(r2,2));//4.607731902E+13
	Ps=100*sqrt(2/a)*sqrt(C-computeF(PI*0/180))/U0;//最底纬度0的极半径
	P=100*sqrt(2/a)*sqrt(C-computeF(PI*dLat/180))/U0;//极坐标的极半径
	D=a*(dPrecision-105)*PI/180;//极坐标的极角
	dX=Ps-P*cos(D);
	dY=P*sin(D);
	UpdateData(false);
}

double CalbersDlg::computeR(double lat)
{
	double r=0.0,n;
	n=sqrt(1-DE2*pow(sin(lat),2));
	r=DA*cos(lat)/n;
	return r;
}

double CalbersDlg::computeF(double lat)
{
	double item1,item2,f;
	item1=sin(lat)/(2*(1-DE2*pow(sin(lat),2)));
	item2=(1/(4*sqrt(DE2)))*log((1+sqrt(DE2)*sin(lat))/(1-sqrt(DE2)*sin(lat)));
	f=pow(DA,2)*(1-DE2)*(item1+item2);
	return f;
}