分段累加法大地主题解算及高斯投影.docx

下载文档 降价啦

27
0
约8.6千字
约 11页
2016-12-30 发布于重庆
举报
版权申诉
保障服务

分段累加法大地主题解算及高斯投影.docx

1、本文档共11页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

大地测量编程实习报告班号：XXXX 学号：XXXXXXXXXX 姓名：XXX大地主题解算（分段累加法正算）结果截图：主要代码：privatevoid button1_Click(object sender, EventArgs e) {double bx, by, bz, B1, lx, ly, lz, L1, ax, ay, az, A1, S;double dB, dL;double e2 = 0.006693421622966, a = 6378245.0000000000;double M, N, W, C;double B2 = 0, L2 = 0, A2 = 0;int Bx, By, Lx, Ly, Ax, Ay;double Bz, Lz, Az; bx = Convert.ToDouble(du1.Text); by = Convert.ToDouble(fen1.Text); bz = Convert.ToDouble(miao1.Text); lx = Convert.ToDouble(du2.Text); ly = Convert.ToDouble(fen2.Text); lz = Convert.ToDouble(miao2.Text); ax = Convert.ToDouble(du3.Text); ay = Convert.ToDouble(fen3.Text); az = Convert.ToDouble(miao3.Text); S = Convert.ToDouble(m.Text);double PI = Math.PI; B1 = (bx + by / 60 + bz / 3600) * PI / 180; L1 = (lx + ly / 60 + lz / 3600) * PI / 180; A1 = (ax + ay / 60 + az / 3600) * PI / 180; W = Math.Sqrt(1 - e2 * (Math.Sin(B1) * Math.Sin(B1))); M = a * (1 - e2) / (W * W * W); N = a / W; C = N * Math.Cos(B1) * Math.Sin(A1);int n ,i; double s; n = (int)(S / 4000 + 1); s = S / n;for (i = 1; i = n; i++) { dB = Math.Cos(A1) * s / M; dL = Math.Sin(A1) * s / N / Math.Cos(B1); B2 = B1 + dB; L2 = L1 + dL; W = Math.Sqrt(1 - e2 * (Math.Sin(B2) * Math.Sin(B2))); M = a * (1 - e2) / (W * W * W); N = a / W; A2 = Math.Asin(C / N / Math.Cos(B2)); B1 = B2; L1 = L2; A1 = A2; } Bx =(int)( B2 / PI * 180); By = (int)((B2 / PI * 180 - Bx) * 60); Bz = Math.Abs(((B2 / PI * 180 - Bx) * 60 - By) * 60); By = Math.Abs(By); Lx = (int)(L2 / PI * 180); Ly = (int)((L2 / PI * 180 - Lx) * 60); Lz = Math.Abs(((L2 / PI * 180 - Lx) * 60 - Ly) * 60); Ly = Math.Abs(Ly); Ax = (int)(A2 / P