// guass coordinateDlg.cpp : implementation file
//
#include "stdafx.h"
#include "guass coordinate.h"
#include "guass coordinateDlg.h"
#include "math.h"
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif
///////////////////////////////////////////////////////////////////////////// // 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)
//}}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()
///////////////////////////////////////////////////////////////////////////// // CGuasscoordinateDlg dialog
CGuasscoordinateDlg::CGuasscoordinateDlg(CWnd* pParent /*=NULL*/) : CDialog(CGuasscoordinateDlg::IDD, pParent)
{
//{{AFX_DATA_INIT(CGuasscoordinateDlg)
m_bdu = 0;
m_bfen = 0;
m_bmiao = 0.0;
m_x = 0.0;
m_y = 0.0;
m_ldu = 0;
m_lfen = 0;
m_lmiao = 0.0;
m_ZoneWidth = -1;
//}}AFX_DATA_INIT
// Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}
void CGuasscoordinateDlg::DoDataExchange(CDataExchange* pDX)
{
CDialog::DoDataExchange(pDX);
//{{AFX_DATA_MAP(CGuasscoordinateDlg)
DDX_Text(pDX, IDC_B_DU, m_bdu);
DDV_MinMaxInt(pDX, m_bdu, 0, 90);
DDX_Text(pDX, IDC_B_FEN, m_bfen);
DDV_MinMaxInt(pDX, m_bfen, 0, 60);
DDX_Text(pDX, IDC_B_MIAO, m_bmiao);
DDV_MinMaxDouble(pDX, m_bmiao, 0., 60.);
DDX_Text(pDX, IDC_x, m_x);
DDX_Text(pDX, IDC_y, m_y);
DDX_Text(pDX, IDC_L_DU, m_ldu);
DDV_MinMaxInt(pDX, m_ldu, 0, 180);
DDX_Text(pDX, IDC_L_FEN, m_lfen);
DDV_MinMaxInt(pDX, m_lfen, 0, 60);
DDX_Text(pDX, IDC_L_MIAO, m_lmiao);
DDV_MinMaxDouble(pDX, m_lmiao, 0., 60.);
DDX_Radio(pDX, IDC_RADIO1, m_ZoneWidth);
//}}AFX_DATA_MAP
}
BEGIN_MESSAGE_MAP(CGuasscoordinateDlg, CDialog)
//{{AFX_MSG_MAP(CGuasscoordinateDlg)
ON_WM_SYSCOMMAND()
ON_WM_PAINT()
ON_WM_QUERYDRAGICON()
ON_COMMAND(ID_ABOUT, OnAbout)
ON_COMMAND(ID_ZHENG, OnZheng)
ON_COMMAND(ID_FAN, OnFan)
ON_BN_CLICKED(ID_CAL, OnCal)
ON_BN_CLICKED(IDC_RADIO1, OnRadio1)
ON_BN_CLICKED(IDC_RADIO2, OnRadio2)
//}}AFX_MSG_MAP
END_MESSAGE_MAP()
///////////////////////////////////////////////////////////////////////////// // CGuasscoordinateDlg message handlers
BOOL CGuasscoordinateDlg::OnInitDialog()
{
CDialog::OnInitDialog();
// Add "About..." menu item to system menu.
// IDM_ABOUTBOX must be in the system command range.
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);
}
}
// 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
CMenu menu;
menu.LoadMenu(IDR_MENU); //IDR_MENU在资源中创建的菜单
SetMenu(&menu);
return TRUE; // return TRUE unless you set the focus to a control
}
void CGuasscoordinateDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
if ((nID & 0xFFF0) == IDM_ABOUTBOX)
{
CAboutDlg dlgAbout;
dlgAbout.DoModal();
}
else
{
CDialog::OnSysCommand(nID, lParam);
}
}
// 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 CGuasscoordinateDlg::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 CGuasscoordinateDlg::OnQueryDragIcon()
{
return (HCURSOR) m_hIcon;
}
void CGuasscoordinateDlg::OnAbout()
{
// TODO: Add your command handler code here
CAboutDlg Dlg;
Dlg.DoModal();
}
void CGuasscoordinateDlg::OnZheng()
{
// TODO: Add your command handler code here
GetDlgItem(IDC_B_DU)->EnableWindow(TRUE);
GetDlgItem(IDC_B_FEN)->EnableWindow(TRUE);
GetDlgItem(IDC_B_MIAO)->EnableWindow(TRUE);
GetDlgItem(IDC_L_DU)->EnableWindow(TRUE);
GetDlgItem(IDC_L_FEN)->EnableWindow(TRUE);
GetDlgItem(IDC_L_MIAO)->EnableWindow(TRUE);
GetDlgItem(IDC_x)->EnableWindow(FALSE);
GetDlgItem(IDC_y)->EnableWindow(FALSE);
SetDlgItemText(IDC_STATIC_INPUT,"请输入B,L坐标值：");
m_x=m_y=0;
//m_bdu=m_bfen=m_ldu=m_lfen=0;
// m_bmiao=m_lmiao=0;
UpdateData(FALSE);
}
void CGuasscoordinateDlg::OnFan()
{
// TODO: Add your command handler code here
GetDlgItem(IDC_x)->EnableWindow(TRUE);
GetDlgItem(IDC_y)->EnableWindow(TRUE);
GetDlgItem(IDC_B_DU)->EnableWindow(FALSE);
GetDlgItem(IDC_B_FEN)->EnableWindow(FALSE);
GetDlgItem(IDC_B_MIAO)->EnableWindow(FALSE);
GetDlgItem(IDC_L_DU)->EnableWindow(FALSE);
GetDlgItem(IDC_L_FEN)->EnableWindow(FALSE);
GetDlgItem(IDC_L_MIAO)->EnableWindow(FALSE);
SetDlgItemText(IDC_STATIC_INPUT,"请输入x,y坐标值：");
m_bdu=m_bfen=m_ldu=m_lfen=0;
m_bmiao=m_lmiao=0;
// m_x=m_y=0;
UpdateData(FALSE);
}
void CGuasscoordinateDlg::OnCal()
{
// TODO: Add your control notification handler code here
UpdateData(TRUE);
double B,L,l,N,a0,a4,a6,a3,a5;
int n,Lo;
double b,Bf,Nf,b2,b3,b4,b5,Z;
if(!GetDlgItem(IDC_x)->IsWindowEnabled()&&(!GetDlgItem(IDC_B_DU)->IsWindow Enabled()))
{
MessageBox("请选择计算方式！","提示",MB_OK|MB_ICONWARNING);
return;
}
if (m_ZoneWidth==-1)
{
MessageBox("请选择投影分带方式！","提示",MB_OK|MB_ICONWARNING);
return;
}
if(GetDlgItem(IDC_x)->IsWindowEnabled())
{
b=m_x/6367558.4969;
n=int(m_y/1000000.0);
Bf=b+(50221746+(293622+(2350+22*cos(b)*cos(b))*cos(b)*cos(b))*cos(b)*cos(b))*p ow(10,-10)*sin(b)*cos(b);
Nf=6399698.902-(21562.267-(108.973-0.612*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf)) *cos(Bf)*cos(Bf);
Z=(m_y-500000-n*1000000.0)/(Nf*cos(Bf));
b2=(0.5+0.00336975*cos(Bf)*cos(Bf))*sin(Bf)*cos(Bf);
b3=0.333333-(0.1666667-0.001123*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
b4=0.25+(0.161612+0.005617*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
b5=0.2-(0.16667-0.00878*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
B=(Bf-(1-(b4-0.147*Z*Z)*Z*Z)*Z*Z*b2)*p;
l=(1-(b3-b5*Z*Z)*Z*Z)*Z*p;
if (m_ZoneWidth==0)
Lo=(n*6-3)*60*60;
else if(m_ZoneWidth=1)
Lo=n*3*60*60;
L=Lo+l;
m_bdu=int(B/60/60);
m_bfen=int(B/60-m_bdu*60);
m_bmiao=B-m_bdu*60*60-m_bfen*60;
m_ldu=int(L/60/60);
m_lfen=int(L/60-m_ldu*60);
m_lmiao=L-m_ldu*60*60-m_lfen*60;
}
if(GetDlgItem(IDC_B_DU)->IsWindowEnabled())
{
B=(m_bdu*60*60+m_bfen*60+m_bmiao)/p;
L=m_ldu*60*60+m_lfen*60+m_lmiao;
if (m_ZoneWidth==0)
{
n=int((m_ldu+m_lfen/60.0+m_lmiao/60.0/60.0)/6.0)+1;
Lo=(6*n-3)*60*60;
}
else if (m_ZoneWidth==1)
{
n=int((m_ldu+m_lfen/60.0+m_lmiao/60.0/60.0)/3.0);
Lo=3*n*60*60;
}
l=(L-Lo)/p;
N=6399698.902-(21562.267-(108.973-0.612*cos(B)*cos(B))*cos(B)*cos(B))*cos( B)*cos(B);
a0=32140.404-(135.3302-(0.7092-0.0040*cos(B)*cos(B))*cos(B)*cos(B))*cos(B) *cos(B);
a4=(0.25+0.00252*cos(B)*cos(B))*cos(B)*cos(B)-0.04166;
a6=(0.166*cos(B)*cos(B)-0.084)*cos(B)*cos(B);
a3=(0.3333333+0.001123*cos(B)*cos(B))*cos(B)*cos(B)-0.1666667;
a5=0.0083-(0.1667-(0.1968+0.0040*cos(B)*cos(B))*cos(B)*cos(B))*cos(B)*cos( B);
m_x=6367558.4969*B-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(B)*cos(B);
m_y=(1+(a3+a5*l*l)*l*l)*l*N*cos(B)+500000.0+n*1000000.0;
}
UpdateData(FALSE);
}
void CGuasscoordinateDlg::OnRadio1()
{
// TODO: Add your control notification handler code here
m_ZoneWidth=0;
}
void CGuasscoordinateDlg::OnRadio2()
{
// TODO: Add your control notification handler code here
m_ZoneWidth=1;
}
高斯投影坐标正反算 V B程序 文件编码（GHTU-UITID-GGBKT-POIU-WUUI-8968）
高斯投影坐标正反算 学院： 班级： 学号： 姓名： 课程名称: 指导老师：
实验目的： 1.了解高斯投影坐标正反算的基本思想； 2.学会编写高斯正反算程序，加深了解。 实验原理： 高斯投影正算公式中应满足的三个条件： 1. 中央子午线投影后为直线； 2. 中央子午线投影后长度不变； 3. 投影具有正形性质，即正形投影条件。 高斯投影反算公式中应满足的三个条件： 1. x坐标轴投影成中央子午线，是投影的对称轴； 2. x轴上的长度投影保持不变； 3. 正形投影条件，即高斯面上的角度投影到椭球面上后角度没有 变形，仍然相等。 操作工具： 计算机中的 代码： Dim a As Double, b As Double, x As Double, y As Double, y_#
Dim l_ As Double, b_ As Double, a0#, a2#, a4#, a6#, a8#, m2#, m4#, m6#, m8#, m0#, l0#, e#, e1# Dim deg1 As Double, min1 As Double, sec1 As Double, deg2 As Double, min2 As Double, sec2 As Double Private Sub Command1_Click() Dim x_ As Double, t#, eta#, N#, W#, k1#, k2#, ik1%, ik2%, dh% deg1 = Val min1 = Val sec1 = Val deg2 = Val min2 = Val sec2 = Val l_ = (deg1 * 3600 + min1 * 60 + sec1) / 206265 b_ = (deg2 * 3600 + min2 * 60 + sec2) / 206265 dh = Val k1 = ((l_ * 180 / + 3) / 6) k2 = (l_ * 180 / / 3) ik1 = Round(k1, 0) ik2 = Round(k2, 0) If dh = 6 Then l0 = 6 * ik1 - 3 Else
GPS ３度、６度带高斯投影如何区分 择投影的目的在于使所选投影的性质、特点适合于地图的用途，同时考虑地图在图廓范围内变形较小而且变形分布均匀。海域使用的地图多采用保角投影，因其能保持方位角度的正确。 我国的基本比例尺地形图(1:5千，1:1万，1:2.5万，1:5万，1:10万，1:25万，1:50万，1:100万)中，大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger)，这是一个等角横切椭圆柱投影，又叫横轴墨卡托投影(Transverse Mercator)；小于50万的地形图采用等角正轴割园锥投影，又叫兰勃特投影(Lambert Conformal Conic)；海上小于50万的地形图多用等角正轴圆柱投影，又叫墨卡托投影(Mercator)。一般应该采用与我国基本比例尺地形图系列一致的地图投影系统。 地图坐标系由大地基准面和地图投影确定，大地基准面是利用特定椭球体对特定地区地球表面的逼近，因此每个国家或地区均有各自的大地基准面，我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系，1978年采用国际大地测量协会推荐的IAG 75地球椭球体建立了我
国新的大地坐标系--西安80坐标系，目前GPS定位所得出的结果都属于WGS84坐标系统，WGS84基准面采用WGS84椭球体，它是一地心坐标系，即以地心作为椭球体中心的坐标系。因此相对同一地理位置，不同的大地基准面，它们的经纬度坐标是有差异的。 采用的3个椭球体参数如下（源自“全球定位系统测量规范GB/T 18314-2001”）： 椭球体与大地基准面之间的关系是一对多的关系，也就是基准面是在椭球体基础上建立的，但椭球体不能代表基准面，同样的椭球体能定义不同的基准面，如前苏联的Pulkovo 1942、非洲索马里的Afgooye基准面都采用了Krassovsky椭球体，但它们的大地基准面显然是不同的。在目前的GIS商用软件中，大地基准面都通过当地基准面向WGS84的转换7参数来定义，即三个平移参数ΔX、ΔY、ΔZ表示两坐标原点的平移值；三个旋转参数εx、εy、εz表示当地坐标系旋转至与地心坐标系平行时，分别绕Xt、Yt、Zt的旋转角；最后是比例校正因子，用于调整椭球大小。北京54、西安80相对WGS84的转换参数至今没有公开，实际工作中可利用工作区内已知的北京54或西安80坐标控制点进行与WGS84坐标值的转换，在
３度６度带高斯投影 选择投影的目的在于使所选投影的性质、特点适合于地图的用途，同时考虑地图在图廓范围内变形较小而且变形分布均匀。海域使用的地图多采用保角投影，因其能保持方位角度的正确。 我国的基本比例尺地形图(1:5千，1:1万，1:万，1:5万，1:10万，1:25万，1:50万，1:100万)中，大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger)，这是一个等角横切椭圆柱投影，又叫横轴墨卡托投影(Transverse Mercator)；小于50万的地形图采用等角正轴割园锥投影，又叫兰勃特投影(Lambert Conformal Conic)；海上小于50万的地形图多用等角正轴圆柱投影，又叫墨卡托投影(Mercator)。一般应该采用与我国基本比例尺地形图系列一致的地图投影系统。 地图坐标系由大地基准面和地图投影确定，大地基准面是利用特定椭球体对特定地区地球表面的逼近，因此每个国家或地区均有各自的大地基准面，我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系，1978年采用国际大地测量协会推荐的IAG 75地球椭球体建立了我国新的大地坐标系--西安80坐标系，目前GPS定位所得出的结果都属于WGS84坐标系统，WGS84基准面采用WGS84椭球体，它是一地心坐标系，即以地心作为椭球体中心的坐标系。因此相对同一地理位置，不同的大地基准面，它们的经纬度坐标是有差异的。 采用的3个椭球体参数如下（源自“全球定位系统测量规范 GB/T 8314-2001”）：
高斯投影坐标正反算编程报告 1. 编程思想 进行高斯投影坐标正反算的编程需要牵涉到大量的公式,为了使程序条理更清楚,各块的数据复用性更强,这里采取了结构化的编程思想。 程序由四大块组成。 GeodesyHomework 、cpp 文件用于存放main()函数,就是整个程序的入口。通过结构化的编程尽力使main()函数变得简单。 MyFunction 、h 与MyFunction 、cpp 用于存放计算过程中进行角度弧度换算时所要用到的一些自定的转换函数。 Zhengsuan 、h 与Zhengsuan 、cpp 用于存放Zhengsuan 类,在Zhengsuan 类中声明了高斯投影坐标正算所要用到的所有变量,在类的构造函数中进行成员变量的初始化及正算计算。通过get 函数获得相应的正算结果。 Fansuan 、h 与Fansuan 、cpp 用于存放Fansuan 类,类似于Zhengsuan 类,Fansuan 类中声明了高斯投影坐标反算所要用到的所有变量,在类的构造函数中进行成员变量的初始化及反算计算。通过get 函数获得相应的反算结果。 2. 计算模型 高斯投影正算公式 6 4256 4 42234 22)5861(cos sin 720)495(cos 24cos sin 2l t t B B N l t B simB N l B B N X x ''+-''+ ''++-''+''?''+=ρηηρρ 5 2224255 3 2233 )5814185(cos 120)1(cos 6cos l t t t B N l t B N l B N y ''-++-''+ ''+-''+''?''=ηηρηρρ 高斯投影反算公式 () () ()( ) 2 22425 52 23 36 4254 222232 8624285cos 12021cos 6cos 459061720935242f f f f f f f f f f f f f f f f f f f f f f f f f f f f f t t t B N y t B N y B N y l y t t y N M t y t t N M t y N M t B B ηηηηη+++++++-=++- -+++ -= 3. 程序框图
class Gauss { #region 高斯投影正反算 ///
高斯投影正反算编程一．高斯投影正反算基本公式（1）高斯正算基本公式 （2）高斯反算基本公式
以上主要通过大地测量学基础课程得到，这不进行详细的推导，只是列出基本公式指导编程的进行。 二．编程的基本方法和流程图 （1）编程的基本方法 高斯投影正反算基本上运用了所有的编程基本语句，本文中是利用C++语言进行基本的设计。高斯正算中对椭球参数和带宽的选择主要运用了选择语句。而高斯反算中除了选择语句的应用，在利用迭代算法求底点纬度还应用了循环语句。编程中还应特别注意相关的度分秒和弧度之间的相互转换，这是极其重要的。 （2）相关流程图 1）正算
选择带宽 3/6 度带 计算带号 输入大地坐标 B ，L 和经差 L0 6 度带 3 度带 选择椭球参数 计算带号 计算弧长 计算平面坐标 x,y 打印 x,y 开始 计算平面坐标 x,y 计算弧长 打印 x,y
开始 输入自然值坐标x,y 和经差L0 选择椭球参数 利用迭代算法 求解底点纬度 利用公式计算B 和L 打印B 和L 2）反算
三．编程的相关代码（1）正算 # include "stdio.h" # include "stdlib.h" # include "math.h" # include "assert.h" #define pi (4*atan(1.0)) int i; struct jin { double B; double L; double L0; }; struct jin g[100]; main(int argc, double *argv[]) { FILE *r=fopen("a.txt","r"); assert(r!=NULL); FILE *w=fopen("b.txt","w"); assert(r!=NULL); int i=0;
几种地图投影的特点及分带方法 做空间分析之前数据准备的时候，将多源数据（如DEM，遥感影像，土地利用图，土壤图，行政区划图等等）转换到统一的坐标系下，让它们能叠在一起，这个工作繁琐，经常让俺头疼，每次得摸索一阵子，好不容易搞明白了，下次做的时候，又因为好久不做，忘得一干二净，为了防止下次做的时候又重新再摸索，就在博客里记一下笔记，供以后用到的时候参考。 在ARCGIS下经纬度坐标的数据和公里格网数据是能自动叠加在一起的——在公里格网数据的投影设臵正确的情况下。而且，6度带的数据与3度带的数据也能自动叠加在一起。只要投影设臵正确了，所有图层都能在ArcGIS里面叠加在一起，整个Data Frame的坐标系统以第一个添加的图层为准，全部自动统一到这个坐标系统下。拿到数据第一件事情，先看X，Y坐标的整数位数，有以下两种情况：(东阳何生的经验总结) 1、X坐标6位，Y坐标7位(东阳何生的经验总结) 没有加带号的坐标，坐标单位是米，假偏东500公里。(东阳何生的经验总结) 2、X坐标8位，Y坐标7位(东阳何生的经验总结) 加了带号的坐标，坐标单位是米。X坐标最前面两位就是添加的带号，这时就要判断是3度带还是6度带，我国幅员辽阔，经度从东经72度到135度，有经验的人一看带号就能大致知道是6度分带还是3度分带；没有经验的，就随便假设一个，然后根据下面的公式算出其中央经线，再与研究区域所在的经度对照，看是否相符，从而判断出是3度分带还是6度分带。带号与中央经线一一对应，知道两者中的任何一个，都能推算出另外一个的值，计算公式如下：(东阳何生的经验总结)
6度带中央经线经度的计算：当地中央经线经度＝6°×当地带号－3°（适用于1∶2．5万和1∶5万地形图） 3度带中央经线经度的计算：中央经线经度＝3°×当地带号（适用于1∶1万地形图） 搞清楚数据坐标的投影之后，就可以在ARCGIS里面定义，此方法可以解决大部分数据叠加问题，采用地方坐标系的特例另当别论，这里只讨论通常情况。(东阳何生的经验总结) 附1：6度分带和3 度分带是怎么回事_百度 知道 1．我国采用6度分 带和3度分带： 1∶2.5万及1∶5万 的地形图采用6度分带 投影，即经差为6度，从 零度子午线开始，自西向 东每个经差6度为一投 影带，全球共分60个带， 用1，2，3，4，5，…… 表示．即东经0～6度为 第一带，其中央经线的经 度为东经3度，东经6～ 12度为第二带，其中央经线的经度为9度。 1∶1万的地形图采用3度分带，从东经1.5度的经线开始，每隔3度为一带，用1，2，3，……表示，全球共划分120个投影带，即东经1.5～ 4.5度为第1带，其中央经线的经度为东经3度，东经4.5～7.5度为第2带，其中央经线的经度为东经6度．我省位于东经113度－东经120度之间，跨第38、39、40共计3个带，其中东经115.5度以西为第38带，其中央经线为东经114度；东经115.5～118.5度为39带，其中央经线为东经117度；东经118.5度以东到山海关为40带，其中央经线为东经120度。 地形图上公里网横坐标前2位就是带号，例如：1∶5万地形图上的横坐标为20345486，其中20即为带号，345486为横坐标值。 2．当地中央经线经度的计算 六度带中央经线经度的计算：当地中央经线经度＝6°×当地带号－3°，例如：地形图上的横坐标为20345，其所处的六度带的中央经线经度为：6°×20－3°＝117°（适用于1∶2．5万和1∶5万地形图）。 三度带中央经线经度的计算：中央经线经度＝3°×当地带号（适用于1∶1万地形图）。 附２：几种投影的特点及分带方法 一、只谈比较常用的几种：“墨卡托投影”、“高斯-克吕格投影”、“UTM 投影”、“兰勃特等角投影” 1．墨卡托(Mercator)投影 1.1 墨卡托投影简介
高斯投影坐标正反算编程报告 10021班 张鑫 学号：2010302590040 1. 编程思想 进行高斯投影坐标正反算的编程需要牵涉到大量的公式，为了使程序条理更清楚，各块的数据复用性更强，这里采取了结构化的编程思想。 程序由四大块组成。 GeodesyHomework.cpp 文件用于存放main()函数，是整个程序的入口。通过结构化的编程尽力使main()函数变得简单。 MyFunction.h 和MyFunction.cpp 用于存放计算过程中进行角度弧度换算时所要用到的一些自定的转换函数。 Zhengsuan.h 和Zhengsuan.cpp 用于存放Zhengsuan 类，在Zhengsuan 类中声明了高斯投影坐标正算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及正算计算。通过get 函数获得相应的正算结果。 Fansuan.h 和Fansuan.cpp 用于存放Fansuan 类，类似于Zhengsuan 类，Fansuan 类中声明了高斯投影坐标反算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及反算计算。通过get 函数获得相应的反算结果。 2. 计算模型 高斯投影正算公式 6 4256 4 42234 22)5861(cos sin 720)495(cos 24cos sin 2l t t B B N l t B simB N l B B N X x ''+-''+ ''++-''+''?''+=ρηηρρ 5 2224255 3 2233)5814185(cos 120)1(cos 6cos l t t t B N l t B N l B N y ''-++-''+''+-''+''?''=ηηρηρρ 高斯投影反算公式 () () ()( ) 222425 52 2336 4254 2222 32 8624285cos 12021cos 6cos 459061720935242f f f f f f f f f f f f f f f f f f f f f f f f f f f f f t t t B N y t B N y B N y l y t t y N M t y t t N M t y N M t B B ηηηηη+++++++-=++- -+++ -=
高斯投影坐标正反算 学院： 班级： 学号： 姓名： 课程名称: 指导老师：
实验目的： 1.了解高斯投影坐标正反算的基本思想； 2.学会编写高斯正反算程序，加深了解。 实验原理： 高斯投影正算公式中应满足的三个条件： 1. 中央子午线投影后为直线； 2. 中央子午线投影后长度不变； 3. 投影具有正形性质，即正形投影条件。 高斯投影反算公式中应满足的三个条件： 1. x坐标轴投影成中央子午线，是投影的对称轴； 2. x轴上的长度投影保持不变； 3. 正形投影条件，即高斯面上的角度投影到椭球面上后角度没 有变形，仍然相等。 操作工具： 计算机中的VB6.0 代码： Dim a As Double, b As Double, x As Double, y As Double, y_# Dim l_ As Double, b_ As Double, a0#, a2#, a4#, a6#, a8#, m2#, m4#,
m6#, m8#, m0#, l0#, e#, e1# Dim deg1 As Double, min1 As Double, sec1 As Double, deg2 As Double, min2 As Double, sec2 As Double Private Sub Command1_Click() Dim x_ As Double, t#, eta#, N#, W#, k1#, k2#, ik1%, ik2%, dh% deg1 = Val(Text1.Text) min1 = Val(Text2.Text) sec1 = Val(Text3.Text) deg2 = Val(Text4.Text) min2 = Val(Text5.Text) sec2 = Val(Text6.Text) l_ = (deg1 * 3600 + min1 * 60 + sec1) / 206265 b_ = (deg2 * 3600 + min2 * 60 + sec2) / 206265 dh = Val(Text9.Text) k1 = ((l_ * 180 / 3.14159 + 3) / 6) k2 = (l_ * 180 / 3.14159 / 3) ik1 = Round(k1, 0) ik2 = Round(k2, 0) If dh = 6 Then l0 = 6 * ik1 - 3 Else If dh = 3 Then
高斯投影坐标正反算 一、相关概念 大地坐标系由大地基准面和地图投影确定，由地图投影到特定椭圆柱面后在南北两极剪开展开而成，是对地球表面的逼近，各国或地区有各自的大地基准面，我国目前主要采用的基准面为：基准面，为GPS基准面，17届国际大地测量协会上推荐，椭圆柱长半轴a=6378137m，短半轴b=； 2.西安80坐标系，1975年国际大地测量协会上推荐，椭圆柱长半轴a=6378140m，短半轴b=; 3.北京54坐标系，参照前苏联克拉索夫斯基椭球体建立,椭圆柱长半轴a=6378245m, 短半轴b=; 通常所说的高斯投影有三种，即投影后： a)角度不变（正角投影），投影后经线和纬线仍然垂直； b)长度不变； c)面积不变； 大地坐标一般采用高斯正角投影，即在地球球心放一点光源，地图投影到过与中央经线相切的椭圆柱面上而成；可分带投影，按中央经线经度值分带，有每6度一带或每3度一带两种(起始带中央经线经度为均为3度，即：6度带1带位置0-6度，3度带1带位置度)，即所谓的高斯-克吕格投影。
图表11高斯投影和分带 地球某点经度（L）为过该点和地球自转轴的半圆与子午线所在半圆夹角，东半球为东经，西半球为西经；地球某点纬度（B）为所在水平面法线与赤道圆面的线面角。 正算是已知大地坐标（L，B），求解高斯平面坐标（X，Y）,为确保Y值为正，Y增加500公里；反算则是由高斯平面坐标（X，Y）求解大地坐标（L，B）。 二、计算模型： 地球椭球面由椭圆绕地球自转轴旋转180度而成。 图表 1 椭圆 椭圆长半轴a，椭圆短半轴b, 椭圆方程：
(1) 图表2椭球面 椭球面方程： y2 a2+ x2 b2 + z2 a2 =1 /*************************************** 与网上充斥的将函数关系先展开为泰勒级数，再依据投影规则确定各参数不同，本文直接依据空间立体三角函数关系得出结果。 *****/ （一）正算 由图表1，
一、高斯-克吕格投影 1、高斯-克吕格简介 高斯-克吕格(Gauss-Kruger)投影简称“高斯投影”，又名"等角横切椭圆柱投影”，地球椭球面和平面间正形投影的一种。德国数学家、物理学家、天文学家高斯（Carl FriedrichＧauss，1777一1855）于十九世纪二十年代拟定，后经德国大地测量学家克吕格（Johannes Kruger，1857～1928）于1912年对投影公式加以补充，故名。该投影按照投影带中央子午线投影为直线且长度不变和赤道投影为直线的条件，确定函数的形式，从而得到高斯一克吕格投影公式。投影后，除中央子午线和赤道为直线外，其他子午线均为对称于中央子午线的曲线。设想用一个椭圆柱横切于椭球面上投影带的中央子午线，按上述投影条件，将中央子午线两侧一定经差范围内的椭球面正形投影于椭圆柱面。将椭圆柱面沿过南北极的母线剪开展平，即为高斯投影平面。取中央子午线与赤道交点的投影为原点，中央子午线的投影为纵坐标x 轴，赤道的投影为横坐标y轴，构成高斯克吕格平面直角坐标系。 2、高斯-克吕格特性 （1）等角投影——投影前后的角度相等，但长度和面积有变形； （2）等距投影——投影前后的长度相等，但角度和面积有变形； （3）等积投影——投影前后的面积相等，但角度和长度有变形。 3、投影的基本概念 它是一种横轴等角切圆柱投影。它把地球视为球体,假想一个平面卷成一个横圆柱面并把它套在球体外面,使横轴圆柱的轴心通过球的中心,球面上一根子午线与横轴圆柱面相切。这样,该子午线在圆柱面上的投影为一直线,赤道面与圆柱面的交线是一条与该子午线投影垂直的直线。将横圆柱面展开成平面,由这两条正交直线就构成高斯-克吕格平面直角坐标系。为减少投影变形,高斯-克吕格投影分为3o带和6o带投影。
三度带与六度带的概念源于高斯平面直角坐标，高斯平面直角坐标适用于：测区范围较大，不能将测区曲面当作平面看待。 当测区范围较大，若将曲面当作平面来看待，则把地球椭球面上的图形展绘到平面上来，必然产生变形，为减小变形，必须采用适当的方法来解决。测量上常采用的方法是高斯投影方法。 高斯投影方法是将地球划分成若干带，然后将每带投影到平面上。 1．我国采用6度带和3度带： 1∶2.5万及1∶5万的地形图采用6度分带投影，即经差为6度，从零度子午线开始，自西向东每个经差6度为一投影带，全球共分60个带，用1，2，3，4，5，……表示．即东经0～6度为第一带，其中央经线的经度为东经3度，东经6～12度为第二带，其中央经线的经度为9度。 1∶1万的地形图采用3度分带，从东经1.5度的经线开始，每隔3度为一带，用1，2，3，……表示，全球共划分120个投影带，即东经1.5～ 4.5度为第1带，其中央经线的经度为东经3度，东经4.5～7.5度为第2带，其中央经线的经度为东经6度．我省位于东经113度－东经120度之间，跨第38、39、40共计3个带，其中东经115.5度以西为第38带，其中央经线为东经114度；东经115.5～118.5度为39带，其中央经线为东经117度；东经118.5度以东到山海关为40带，其中央经线为东经120度。 地形图上公里网横坐标前2位就是带号，例如：1∶5万地形图上的横坐标为20345486，其中20即为带号，345486为横坐标值。 2．当地中央经线经度的计算 六度带中央经线经度的计算：当地中央经线经度＝6°×当地带号－3°，例如：地形图上的横坐标为20345，其所处的六度带的中央经线经度为：6°×20－3°＝117°（适用于1∶2.5万和1∶5万地形图）。 三度带中央经线经度的计算：中央经线经度＝3°×当地带号（适用于1∶1万地形图）。 例如在河北境内大部分地区： 三度带中央经线经度的计算：114°=3°×38 （横坐标为38******） 六度带中央经线经度的计算：117°=6°×20-3°（横坐标为20******）
高斯投影坐标正算公式 高斯投影坐标正反算公式 2.2.2. 1高斯投影坐标正算公式： B, x,y 高斯投影必须满足以下三个条件： ⑴中央子午线投影后为直线；⑵中央子午线投影后长度不变；⑶投影具有正形性质，即正形投影条件。 由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即 式中，x为的偶函数，y为的奇函数；，即， 如展开为的级数，收敛。 （2-10） 式中是待定系数，它们都是纬度B的函数。 由第三个条件知： 分别对和q求偏导数并代入上式 (2-11) 上两式两边相等，其必要充分条件是同次幂前的系数应相等，即
(2-12) (2-12)是一种递推公式，只要确定了就可依次确定其余各系数。 由第二条件知:位于中央子午线上的点，投影后的纵坐标x应等于投影前从赤道量至该点的子午线弧长X，即(2-10)式第一式中，当时有： (2-13) 顾及(对于中央子午线) 得： (2-14,15) (2-16) 依次求得并代入(2-10)式，得到高斯投影正算公式
(2-17) 2.2.2. 2高斯投影坐标反算公式 x,y B, 投影方程： (2-18) 高斯投影坐标反算公式推导要复杂些。 ⑴由x求底点纬度(垂足纬度),对应的有底点处的等量纬度，求x,y与 的关系式，仿照式有， 由于y和椭球半径相比较小(1/16.37)，可将展开为y的幂级数；又由于是对称投影，q必是y的偶函数，必是y的奇函数。 (2-19) 是待定系数，它们都是x的函数. 由第三条件知： ，
， (2-20) (2-19)式分别对x和y求偏导数并代入上式 上式相等必要充分条件，是同次幂y前的系数相等， 第二条件，当y=0时，点在中央子午线上，即x=X，对应的点称为底点，其纬度为底点纬度，也就是x=X时的子午线弧长所对应的纬度，设所对应的等量纬度为。也就是在底点展开为y的幂级数。 由(2-19)1式 依次求得其它各系数 (2-21) (2-21)1 ………… 将代入(2-19)1式得
如何判断投影坐标是3度带坐标还是6度带坐标标题1．我国采用6度分带和3度分带： 1∶2.5万及1∶5万的地形图采用6度分带投影，即经差为6度，从零度子午线开始，自西向东每个经差6度为一投影带，全球共分60个带，用1，2，3，4，5，……表示．即东经0～6度为第一带，其中央经线的经度为东经3度，东经6～12度为第二带，其中央经线的经度为9度。 1∶1万的地形图采用3度分带，从东经1.5度的经线开始，每隔3度为一带，用1，2，3，……表示，全球共划分120个投影带，即东经1.5～4.5度为第1带，其中央经线的经度为东经3度，东经4.5～7.5度为第2带，其中央经线的经度为东经6度．我省位于东经113度－东经120度之间，跨第38、39、40共计3个带，其中东经115.5度以西为第38带，其中央经线为东经114度；东经115.5～118.5度为39带，其中央经线为东经117度；东经118.5度以东到山海关为40带，其中央经线为东经120度。 地形图上公里网横坐标前2位就是带号，例如：1∶5万地形图上的横坐标为，其中20即为带号，345486为横坐标值。 2．当地中央经线经度的计算 六度带中央经线经度的计算：当地中央经线经度＝6°×当地带号－3°，例如：地形图上的横坐标为20345，其所处的六度带的中央经线经度为：6°×20－3°＝117°（适用于1∶2．5万和1∶5万地形图）。 三度带中央经线经度的计算：中央经线经度＝3°×当地带号（适用于1∶1万地形图）。 3、如何计算当地的中央子午线？ 当地中央子午线决定于当地的直角坐标系统，首先确定您的直角坐标系统是3度带还是6度带投影公式推算: 6度带中央子午线计算公式：当地经度/6=N；中央子午线L=6 * N （带号）
高斯投影坐标正反算V B 程序 Jenny was compiled in January 2021
高斯投影坐标正反算 学院： 班级： 学号： 姓名： 课程名称: 指导老师：
实验目的： 1.了解高斯投影坐标正反算的基本思想； 2.学会编写高斯正反算程序，加深了解。 实验原理： 高斯投影正算公式中应满足的三个条件： 1. 中央子午线投影后为直线； 2. 中央子午线投影后长度不变； 3. 投影具有正形性质，即正形投影条件。 高斯投影反算公式中应满足的三个条件： 1. x坐标轴投影成中央子午线，是投影的对称轴； 2. x轴上的长度投影保持不变； 3. 正形投影条件，即高斯面上的角度投影到椭球面上后角度没有 变形，仍然相等。 操作工具： 计算机中的 代码： Dim a As Double, b As Double, x As Double, y As Double, y_#
Dim l_ As Double, b_ As Double, a0#, a2#, a4#, a6#, a8#, m2#, m4#, m6#, m8#, m0#, l0#, e#, e1# Dim deg1 As Double, min1 As Double, sec1 As Double, deg2 As Double, min2 As Double, sec2 As Double Private Sub Command1_Click() Dim x_ As Double, t#, eta#, N#, W#, k1#, k2#, ik1%, ik2%, dh% deg1 = Val min1 = Val sec1 = Val deg2 = Val min2 = Val sec2 = Val l_ = (deg1 * 3600 + min1 * 60 + sec1) / 206265 b_ = (deg2 * 3600 + min2 * 60 + sec2) / 206265 dh = Val k1 = ((l_ * 180 / + 3) / 6) k2 = (l_ * 180 / / 3) ik1 = Round(k1, 0) ik2 = Round(k2, 0) If dh = 6 Then l0 = 6 * ik1 - 3 Else
高斯投影正、反算代码 //高斯投影正、反算 //////6度带宽 54年北京坐标系 //高斯投影由经纬度(Unit:DD)反算大地坐标(含带号，Unit:Metres) void GaussProjCal(double longitude, double latitude, double *X, double *Y) { int ProjNo=0; int ZoneWide; ////带宽 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; double a,f, e2,ee, NN, T,C,A, M, iPI; iPI = 0.0174532925199433; ////3.1415926535898/180.0; ZoneWide = 6; ////6度带宽 a=6378245.0; f=1.0/298.3; //54年北京坐标系参数 ////a=6378140.0; f=1/298.257; //80年西安坐标系参数 ProjNo = (int)(longitude / ZoneWide) ; longitude0 = ProjNo * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ; latitude0=0; longitude1 = longitude * iPI ; //经度转换为弧度 latitude1 = latitude * iPI ; //纬度转换为弧度 e2=2*f-f*f; ee=e2*(1.0-e2); NN=a/sqrt(1.0-e2*sin(latitude1)*sin(latitude1)); T=tan(latitude1)*tan(latitude1); C=ee*cos(latitude1)*cos(latitude1); A=(longitude1-longitude0)*cos(latitude1); M=a*((1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)*latitude1-(3*e2/8+3*e2*e2 /32+45*e2*e2 *e2/1024)*sin(2*latitude1)
在采用分带的投影坐标系统中，我们最常用的是高斯-克吕格投影，它是由德国数学家、物理学家、天文学家高斯（Carl Friedrich Ｇauss，1777—1855）于十九世纪二十年代拟定，后经德国大地测量学家克吕格（Johannes Kruger，1857～1928）于1912年对投影公式加以补充，所以因此而得名。它是横轴墨卡托投影的一个变种，高斯-克吕格只是它通俗的名称，比较专业的名称叫做横轴等角切椭圆柱投影。 设想用一个圆柱横切于球面上投影带的中央经线，按照投影带中央经线投影为直线且长度不变和赤道投影为直线的条件，将中央经线两侧一定经差范围内的球面正形投影于圆柱面。然后将圆柱面沿过南北极的母线剪开展平，即获高斯—克吕格投影平面。高斯—克吕格投影后，除中央经线和赤道为直线外，其他经线均为对称于中央经线的曲线。高斯-克吕格投影没有角度变形，在长度和面积上变形也很小，中央经线无变形，自中央经线向投影带边缘，变形逐渐增加，变形最大处在投影带内赤道的两端。下图是高斯—克吕格投影方式示意图。 图一高斯克吕格投影的投影方式 高斯—克吕格投影按一定经差将地球椭球面划分成若干投影带，这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差，又要使带数不致过多以减少换带计算工作，据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带，以便分带投影。通常按经差6度或3度分为六度带或三度带。六度带自0度子午线起每隔经差6度自西向东分带，带号依次编为第1、2…60带。三度带是在六度带的基础上分成的，它的中央子午线与六度带的中央子午线和分带子午线重合，即自1.5度子午线起每隔经差3度自西
§8.3高斯投影坐标正反算公式 任何一种投影①坐标对应关系是最主要的；②如果是正形投影，除了满足正形投影的条件外（C-R 偏微分方程），还有它本身的特殊条件。 8.3.1高斯投影坐标正算公式： B,l ? x,y 高斯投影必须满足以下三个条件： ①中央子午线投影后为直线；②中央子午线投影后长度不变；③投影具有正形性质，即正形投影条件。 由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即 (8-10)式中，x 为l 的偶函数，y 为l 的奇函数；0330 '≤l ，即20/1/≈''''ρl ， 如展开为l 的级数，收敛。 +++=++++=553316644220l m l m l m y l m l m l m m x （8-33） 式中 ,,10m m 是待定系数，它们都是纬度B 的函数。 由第三个条件知： q y l x l y q x ??-=????=??, (8-33)式分别对l 和q 求偏导数并代入上式 ----=++++++=+++553315 63424 42204 52 3164253l dq dm l dq dm l dq dm l m l m l m l dq dm l dq dm dq dm l m l m m (8-34) 上两式两边相等，其必要充分条件是同次幂l 前的系数应相等，即
dq dm m dq dm m dq dm m 231 20 13121? =? -== (8-35) (8-35)是一种递推公式，只要确定了 0m 就可依次确定其余各系数。 由第二条件知:位于中央子午线上的点，投影后的纵坐标x 应等于投影前从赤道量至该点的子午线弧长X ，即(8-33)式第一式中，当0=l 时有： 0m X x == (8-36) 顾及(对于中央子午线) B V M r M B N dq dB M dB dX cos cos 2 ==== 得： B V c B N r dq dB dB dX dq dX dq dm m cos cos 01===?===(8-37,38) B B N dq dB dB dm dq dm m cos sin 2 2121112=?-=?-= (8-39) 依次求得6543,,,m m m m 并代入(8-33)式，得到高斯投影正 算公式
高斯投影正反算c代码 Coca-cola standardization office【ZZ5AB-ZZSYT-ZZ2C-ZZ682T-ZZT18】
高斯投影正反算程序设计 一．程序设计流程 本程序的设计思路如下： （1），程序采用VS08版本作为开发平台，并采用C#语言作为开发语言，设计为WindowsForm窗体程序形式。 （2），程序主要的算法来自于教材。但是本程序为了更加实用，添加了更多的解算基准，包括：WGS-84，国际椭球1975，克氏椭球，和2000国家大地坐标系。 （3），程序为了更方便的读取数据和输出数据，故需要自己定义了固定的数据输入格式和数据输出格式或形式，请老师注意查看。 二．代码 using System; using ; using ; using ; using ; using Gauss { public partial class Form1 : Form { double b = (a * a * (1 - ee * ee)); double c = a * a / b; double epp = ((a * a - b * b) / b / b); CRDGEODETIC pcrdGeo; CRDCARTESIAN pcrdCar; double midlong = 0;
//求X，Y和带号 = ; ytext = ; string temp = (0, 2); num = (temp); ytext = (0, 2); = (ytext) - 500000; try { tt = } catch { ("Choose 3/6 error!"); return; } if ("3度带") == 0) { midlong = num * 3 * pai / 180; } if ("6度带") == 0) { midlong = (6 * num - 3) * pai / 180; } b = (a * a * (1 - ee * ee)); c = a * a / b; epp = (a * a - b * b) / b; double m0, m2, m4, m6, m8; double a0, a2, a4, a6, a8; m0 = a * (1 - ee * ee);