常用测绘C程序设计.docx

1、常用测绘C程序设计常用测量程序设计（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.Wr

2、ite(请输入垂直角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.

3、Cos(a) / 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) * fu

4、hao) / 180.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坐标=)

5、; double 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,

6、YB); /将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; （3）如图所示，已知A点和B点的坐标，计

7、算A点到的边长及方位角。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(); Co

8、nsole.Write(请输入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

9、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; /计算方位角,返回弧度值 public static double 方位角(double x1, doub

10、le 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 0) angle += Math.PI; if (deltaY 2 * Math.PI) a0 -= 2 * Math.PI; List 导线转角集合 = new List(); int i = 1; do Console

11、.Write(请输入第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

12、,i+,DMS(a0); /将弧度转为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) *

13、fuhao; /将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点

14、和M点的坐标及每个左角和每条边长，求每个点的坐标。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.R

15、eadLine(); 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 (s

16、tr != ) 导线转角集合.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

17、* Math.Cos(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 * 100

18、00; return (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); i

19、f (angle 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

20、Ha = double.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; C

21、onsole.Write(请输入第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) angle = 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 -