《常用测绘C程序设计(源代码).doc》由会员分享，可在线阅读，更多相关《常用测绘C程序设计(源代码).doc（32页珍藏版）》请在三一文库上搜索。
1、常用测量程序设计（1）用全站仪在A点观测了B点斜边和垂直角，求A到B的高差。（提示：，D-斜边，-垂直角，-仪器高，-反光镜高，-大气折光系数）using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1 class Application static void Main(string args) Console.Write(请输入斜边=); double D = double.Parse(Console.ReadLine(); Console.Write(请输入垂直
2、角ddd.mmss=); double a = DEG(double.Parse(Console.ReadLine(); Console.Write(请输入仪器高=); double i = double.Parse(Console.ReadLine(); Console.Write(请输入反光镜高=); double v = double.Parse(Console.ReadLine(); double h = D * Math.Sin(a) + (1 - 0.13) * D / 6371000.0 * D / 6371000.0 * Math.Cos(a) * Math.Cos(a) /
3、2.0 + i - v; Console.WriteLine(高差为0,h); /将ddd.mmss转为弧度 static public double DEG(double ang) int fuhao = (int)(ang / Math.Abs(ang); ang = Math.Abs(ang); int d = (int)ang; int m = (int)(ang * 100) - d * 100; double s = ang * 10000 - m * 100 - d * 10000; return (d + m / 60.0 + s / 3600.0) * fuhao) / 18
4、0.0 * Math.PI; （2）如图所示，已知A点的坐标及A点到B点的边长及方位角，计算B点的坐标。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1 class Application static void Main(string args) Console.Write(请输入A点的X坐标=); double XA = double.Parse(Console.ReadLine(); Console.Write(请输入A点的Y坐标=); double
5、YA = double.Parse(Console.ReadLine(); Console.Write(请输入A到B的方位角ddd.mmss=); double a = DEG(double.Parse(Console.ReadLine(); Console.Write(请输入A到B的水平距离=); double S = double.Parse(Console.ReadLine(); double XB = XA + S*Math.Cos(a); double YB = YA + S*Math.Sin(a); Console.WriteLine(B点的坐标（0,1）,XB,YB); /将dd
6、d.mmss转为弧度 static public double DEG(double ang) int fuhao = (int)(ang / Math.Abs(ang); ang = Math.Abs(ang); int d = (int)ang; int m = (int)(ang * 100) - d * 100; double s = ang * 10000 - m * 100 - d * 10000; return (d + m / 60.0 + s / 3600.0) * fuhao) / 180.0 * Math.PI; （3）如图所示，已知A点和B点的坐标，计算A点到的边长及方
7、位角。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1 class Application static void Main(string args) Console.Write(请输入A点的X坐标=); double XA = double.Parse(Console.ReadLine(); Console.Write(请输入A点的Y坐标=); double YA = double.Parse(Console.ReadLine(); Console.Wri
8、te(请输入B点的X坐标=); double XB = double.Parse(Console.ReadLine(); Console.Write(请输入B点的Y坐标=); double YB = double.Parse(Console.ReadLine(); double S = 距离(XA, YA, XB, YB); double a = 方位角(XA, YA, XB, YB); Console.WriteLine(AB间的距离=0，从A到B的方位角=1,S,DMS(a); /将弧度转为ddd.mmss static public double DMS(double ang) ang
9、+= 1.0E-15;/加上一个小量，以保证进位 int fuhao = (int)(ang / Math.Abs(ang); ang = Math.Abs(ang) * 180.0 / Math.PI; int d = (int)ang; ang = (ang - d) * 60.0; int m = (int)ang; double s = (ang - m) * 60.0; return (d + m / 100.0 + s / 10000.0) * fuhao; /计算方位角,返回弧度值 public static double 方位角(double x1, double y1, do
10、uble x2, double y2) double deltaX = x2 - x1; double deltaY = y2 - y1; double angle = Math.PI * 0.5; if (Math.Abs(deltaX) 0.000000001) angle = Math.Atan2(deltaY, deltaX); if (angle 0) angle += Math.PI; if (deltaY 2 * Math.PI) a0 -= 2 * Math.PI; List 导线转角集合 = new List(); int i = 1; do Console.Write(请输
11、入第0个转角的水平角左角为正，右角为负 =,i+); string str = Console.ReadLine(); if (str != ) 导线转角集合.Add(DEG(double.Parse(str); else break; while (true); i = 1; foreach (double a in 导线转角集合) a0 += a + Math.PI; if (a0 2 * Math.PI) a0 -= 2 * Math.PI; else if (a0 0.0) a0 += 2 * Math.PI; Console.WriteLine(第0条边的方位角为1,i+,DMS(a
12、0); /将弧度转为ddd.mmss static public double DMS(double ang) ang += 1.0E-15;/加上一个小量，以保证进位 int fuhao = (int)(ang / Math.Abs(ang); ang = Math.Abs(ang) * 180.0 / Math.PI; int d = (int)ang; ang = (ang - d) * 60.0; int m = (int)ang; double s = (ang - m) * 60.0; return (d + m / 100.0 + s / 10000.0) * fuhao; /将
13、ddd.mmss转为弧度 static public double DEG(double ang) int fuhao = (int)(ang / Math.Abs(ang); ang = Math.Abs(ang); int d = (int)ang; int m = (int)(ang * 100) - d * 100; double s = ang * 10000 - m * 100 - d * 10000; return (d + m / 60.0 + s / 3600.0) * fuhao) / 180.0 * Math.PI; （5）在如图所示的支中导线，已知A点和M点的坐标及每个
14、左角和每条边长，求每个点的坐标。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1 class Application static void Main(string args) Console.Write(请输入M点的x坐标=); double Mx = double.Parse(Console.ReadLine(); Console.Write(请输入M点的y坐标=); double My = double.Parse(Console.ReadLine()
15、; Console.Write(请输入A点的x坐标=); double Ax = double.Parse(Console.ReadLine(); Console.Write(请输入A点的y坐标=); double Ay = double.Parse(Console.ReadLine(); List 导线转角集合 = new List(); List 导线边长集合 = new List(); int i = 1; do Console.Write(请输入第0个转角的水平角左角为正，右角为负=,i); string str = Console.ReadLine(); if (str != ) 导
16、线转角集合.Add(DEG(double.Parse(str); else break; Console.Write(请输入第0条边长值=, i+); 导线边长集合.Add( double.Parse(Console.ReadLine(); while (true); /计算M到A的坐标方位角 double a0 = 方位角(Mx, My, Ax, Ay); /计算每个点的坐标 double x0 = Ax; double y0 = Ay; for (int j = 0; j 2 * Math.PI) a0 -= 2 * Math.PI; x0 = x0 + 导线边长集合j * Math.Co
17、s(a0); y0 = y0 + 导线边长集合j * Math.Sin(a0); Console.WriteLine(P0点的坐标是:1,2,j+2,x0,y0); /将ddd.mmss转为弧度 static public double DEG(double ang) int fuhao = (int)(ang / Math.Abs(ang); ang = Math.Abs(ang); int d = (int)ang; int m = (int)(ang * 100) - d * 100; double s = ang * 10000 - m * 100 - d * 10000; retur
18、n (d + m / 60.0 + s / 3600.0) * fuhao) / 180.0 * Math.PI; /计算方位角,返回弧度值 public static double 方位角(double x1, double y1, double x2, double y2) double deltaX = x2 - x1; double deltaY = y2 - y1; double angle = Math.PI * 0.5; if (Math.Abs(deltaX) 0.000000001) angle = Math.Atan2(deltaY, deltaX); if (angle
19、0) angle += Math.PI; if (deltaY 0.0) angle += Math.PI; return angle; （6）在如图所示的单一附合水准路线中，已知A点和B点的高程及每段的长度和高差，求每个点的高程。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1 class Application static void Main(string args) Console.Write(请输入A点的高程=); double Ha = doub
20、le.Parse(Console.ReadLine(); Console.Write(请输入B点的高程=); double Hb = double.Parse(Console.ReadLine(); List 路线段长度集合 = new List(); List 路线段高差集合 = new List(); int i = 1; do Console.Write(请输入第0段长度公里为单位=,i); string str = Console.ReadLine(); if (str != ) 路线段长度集合.Add(double.Parse(str); else break; Console.Wr
21、ite(请输入第0段高差=, i+); 路线段高差集合.Add(double.Parse(Console.ReadLine(); while (true); /计算水准路线闭合差 double w = 0; double 总长度 = 0; for (int j = 0; j 路线段长度集合.Count; j+) w += 路线段高差集合j; 总长度 += 路线段长度集合j; w = w - (Hb - Ha); /闭合差 double v = -w / 总长度; /每公里改正数 /计算每个点的高程 double H = Ha; for (int j = 0; j 0.000000001) an
22、gle = Math.Atan2(deltaY, deltaX); if (angle 0) angle += Math.PI; if (deltaY 0.0) angle += Math.PI; return angle; /计算距离 public static double 距离(double x1, double y1, double x2, double y2) return Math.Sqrt(x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1); （8）在如图所示的单个多边形中，已知每个点的坐标，求多边形的面积。提示面积的计算公式为：using
23、System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1 class Application static void Main(string args) List 多边形X坐标集 = new List(); List 多边形Y坐标集 = new List(); int i = 1; do Console.Write(请输入第0点的x坐标=,i); string str = Console.ReadLine(); if (str != ) 多边形X坐标集.Add(double.Parse(str); else break; Console.Write(请输入第0点的y坐标=,i+); 多边形Y坐标集.Add( double.Parse(Console.ReadLine(); while (true); /计算面积 double area = 0; /点的个数 int Count = 多边形X坐标集.Count; for (int j = 0; j Count-1; j+) area += 多边形X坐标集j * 多边形Y坐标集j + 1 - 多边形X坐标集j + 1 * 多边形Y坐标集j; /计算最后一点