測繪類C程序代碼

1、 常用測量程序設(shè)計代碼（1）用全站儀在A點觀測了B點斜邊和垂直角，求A到B的高差。（提示：，D-斜邊，-垂直角，-儀器高，-反光鏡高，-大氣折光系數(shù)）using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(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 / 6

3、371000.0 * Math.Cos(a) * Math.Cos(a) / 2.0 + i - v;Console.WriteLine("高差為0",h); /將ddd.mmss轉(zhuǎn)為弧度staticpublicdouble 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

4、n (d + m / 60.0 + s / 3600.0) * fuhao) / 180.0 * Math.PI; （2）如圖所示，已知A點的坐標(biāo)及A點到B點的邊長及方位角，計算B點的坐標(biāo)。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(string args) Console.Write("請輸入A點的X坐標(biāo)=");double XA = double.Parse(Cons

5、ole.ReadLine();Console.Write("請輸入A點的Y坐標(biāo)=");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.Co

6、s(a);double YB = YA + S*Math.Sin(a);Console.WriteLine("B點的坐標(biāo)（0,1）",XB,YB); /將ddd.mmss轉(zhuǎn)為弧度staticpublicdouble 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

7、+ m / 60.0 + s / 3600.0) * fuhao) / 180.0 * Math.PI; （3）如圖所示，已知A點和B點的坐標(biāo)，計算A點到的邊長及方位角。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(string args) Console.Write("請輸入A點的X坐標(biāo)=");double XA = double.Parse(Console.ReadLi

8、ne();Console.Write("請輸入A點的Y坐標(biāo)=");double YA = double.Parse(Console.ReadLine();Console.Write("請輸入B點的X坐標(biāo)=");double XB = double.Parse(Console.ReadLine();Console.Write("請輸入B點的Y坐標(biāo)=");double YB = double.Parse(Console.ReadLine();double S = 距離(XA, YA, XB, YB);double a = 方位角(XA,

9、YA, XB, YB);Console.WriteLine("AB間的距離=0，從A到B的方位角=1",S,DMS(a); /將弧度轉(zhuǎn)為ddd.mmssstaticpublicdouble DMS(double ang) ang += 1.0E-15;/加上一個小量，以保證進(jìn)位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

10、- m) * 60.0;return (d + m / 100.0 + s / 10000.0) * fuhao; /計算方位角,返回弧度值publicstaticdouble 方位角(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 &

11、lt; 0) angle += Math.PI; if (deltaY < 0.0) angle += Math.PI; return angle; /計算距離publicstaticdouble 距離(double x1, double y1, double x2, double y2) returnMath.Sqrt(x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1); （4）在如圖所示的支中導(dǎo)線，已知A點到M點的坐標(biāo)方位角及每個左角，求每條邊的坐標(biāo)方位角。using System;using System.Collections.Generic

12、;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(string args) Console.Write("請輸入A點到M點的坐標(biāo)方位角=");/將A點到M點的坐標(biāo)方位角換算為M點到A點的坐標(biāo)方位角，以便利用公式double a0 = DEG(double.Parse(Console.ReadLine() +180.0);if (a0 > 2 * Math.PI) a0 -= 2 * Math.PI; List<double> 導(dǎo)線轉(zhuǎn)角集合 = ne

13、wList<double>();int i = 1;do Console.Write("請輸入第0個轉(zhuǎn)角的水平角左角為正，右角為負(fù)<直接回車結(jié)束輸入>=",i+);string str = Console.ReadLine();if (str != "") 導(dǎo)線轉(zhuǎn)角集合.Add(DEG(double.Parse(str); else break; while (true); i = 1;foreach (double a in 導(dǎo)線轉(zhuǎn)角集合) a0 += a + Math.PI;if (a0 > 2 * Math.PI) a

14、0 -= 2 * Math.PI; elseif (a0 < 0.0) a0 += 2 * Math.PI; Console.WriteLine("第0條邊的方位角為1",i+,DMS(a0); /將弧度轉(zhuǎn)為ddd.mmssstaticpublicdouble DMS(double ang) ang += 1.0E-15;/加上一個小量，以保證進(jìn)位int fuhao = (int)(ang / Math.Abs(ang); ang = Math.Abs(ang) * 180.0 / Math.PI;int d = (int)ang; ang = (ang - d) *

15、 60.0;int m = (int)ang;double s = (ang - m) * 60.0;return (d + m / 100.0 + s / 10000.0) * fuhao; /將ddd.mmss轉(zhuǎn)為弧度staticpublicdouble 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 * 1000

16、0;return (d + m / 60.0 + s / 3600.0) * fuhao) / 180.0 * Math.PI; （5）在如圖所示的支中導(dǎo)線，已知A點和M點的坐標(biāo)及每個左角和每條邊長，求每個點的坐標(biāo)。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(string args) Console.Write("請輸入M點的x坐標(biāo)=");double Mx = doub

17、le.Parse(Console.ReadLine();Console.Write("請輸入M點的y坐標(biāo)=");double My = double.Parse(Console.ReadLine();Console.Write("請輸入A點的x坐標(biāo)=");double Ax = double.Parse(Console.ReadLine();Console.Write("請輸入A點的y坐標(biāo)=");double Ay = double.Parse(Console.ReadLine();List<double> 導(dǎo)線轉(zhuǎn)角集合

18、= newList<double>();List<double> 導(dǎo)線邊長集合 = newList<double>();int i = 1;do Console.Write("請輸入第0個轉(zhuǎn)角的水平角左角為正，右角為負(fù)<直接回車結(jié)束輸入>=",i);string str = Console.ReadLine();if (str != "") 導(dǎo)線轉(zhuǎn)角集合.Add(DEG(double.Parse(str); else break; Console.Write("請輸入第0條邊長值=",

19、i+); 導(dǎo)線邊長集合.Add( double.Parse(Console.ReadLine(); while (true);/計算M到A的坐標(biāo)方位角double a0 = 方位角(Mx, My, Ax, Ay);/計算每個點的坐標(biāo)double x0 = Ax;double y0 = Ay;for (int j = 0; j < 導(dǎo)線轉(zhuǎn)角集合.Count; j+) a0 += 導(dǎo)線轉(zhuǎn)角集合j + Math.PI;if (a0 > 2 * Math.PI) a0 -= 2 * Math.PI; x0 = x0 + 導(dǎo)線邊長集合j * Math.Cos(a0); y0 = y0 + 導(dǎo)

20、線邊長集合j * Math.Sin(a0);Console.WriteLine("P0點的坐標(biāo)是:1,2",j+2,x0,y0); /將ddd.mmss轉(zhuǎn)為弧度staticpublicdouble 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.

21、0 + s / 3600.0) * fuhao) / 180.0 * Math.PI; /計算方位角,返回弧度值publicstaticdouble 方位角(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 < 0) angle +

22、= Math.PI; if (deltaY < 0.0) angle += Math.PI; return angle; （6）在如圖所示的單一附合水準(zhǔn)路線中，已知A點和B點的高程及每段的長度和高差，求每個點的高程。using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(string args) Console.Write("請輸入A點的高程=");double Ha = do

23、uble.Parse(Console.ReadLine();Console.Write("請輸入B點的高程=");double Hb = double.Parse(Console.ReadLine();List<double> 路線段長度集合 = newList<double>();List<double> 路線段高差集合 =newList<double>();int i = 1;do Console.Write("請輸入第0段長度公里為單位<直接回車結(jié)束輸入>=",i);string str

24、= Console.ReadLine();if (str != "") 路線段長度集合.Add(double.Parse(str); else break; Console.Write("請輸入第0段高差=", i+); 路線段高差集合.Add(double.Parse(Console.ReadLine(); while (true);/計算水準(zhǔn)路線閉合差double w = 0;double 總長度 = 0;for (int j = 0; j < 路線段長度集合.Count; j+) w += 路線段高差集合j; 總長度 += 路線段長度集合j;

25、 w = w - (Hb - Ha); /閉合差double v = -w / 總長度; /每公里改正數(shù)/計算每個點的高程double H = Ha;for (int j = 0; j < 路線段長度集合.Count; j+) H += 路線段高差集合j + v * 路線段長度集合j;Console.WriteLine("P0點的高程是:1",j+1,H); Console.WriteLine("最后一點就是B點！"); （7）在如圖所示的前方交會中，ABJ三點按逆時針方向排列，已知AB兩點的坐標(biāo)和兩個交會角度，求J點的坐標(biāo)。using Syste

26、m;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(string args) Console.Write("請輸入A點的x坐標(biāo)=");double Ax = double.Parse(Console.ReadLine();Console.Write("請輸入A點的y坐標(biāo)=");doubleAy = double.Parse(Console.ReadLine();Console.Wr

27、ite("請輸入B點的x坐標(biāo)=");double Bx = double.Parse(Console.ReadLine();Console.Write("請輸入B點的y坐標(biāo)=");double By = double.Parse(Console.ReadLine();Console.Write("請輸入的角度值=");double = DEG(double.Parse(Console.ReadLine();Console.Write("請輸入的角度值=");double = DEG(double.Parse(Con

28、sole.ReadLine();/計算A到B的坐標(biāo)方位角double a0 = 方位角(Ax, Ay, Bx, By);/計算A到J的坐標(biāo)方位角 a0 -= ;/計算A到B的邊長double s = 距離(Ax, Ay, Bx, By);/計算A到J的邊長 s = s / Math.Sin(Math.PI - - ) * Math.Sin();/計算J的坐標(biāo)double Jx = Ax + s * Math.Cos(a0);double Jy = Ay + s * Math.Sin(a0);Console.WriteLine("J點的坐標(biāo)=0,1",Jx,Jy); /將dd

29、d.mmss轉(zhuǎn)為弧度staticpublicdouble 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; /計算方位角,返回弧度值publicstaticdouble 方位角(d

30、ouble 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 < 0) angle += Math.PI; if (deltaY < 0.0) angle += Math.PI; return angle; /計算距離publicstati

31、cdouble 距離(double x1, double y1, double x2, double y2) returnMath.Sqrt(x2 - x1) * (x2 - x1) + (y2 - y1) * (y2 - y1); （8）在如圖所示的單個多邊形中，已知每個點的坐標(biāo)，求多邊形的面積。提示面積的計算公式為：using System;using System.Collections.Generic;using System.Text;namespace ConsoleApplication1classApplication staticvoidMain(string args) L

32、ist<double> 多邊形X坐標(biāo)集 = newList<double>();List<double> 多邊形Y坐標(biāo)集 = newList<double>();int i = 1;do Console.Write("請輸入第0點的x坐標(biāo)<直接回車就結(jié)束輸入>=",i);string str = Console.ReadLine();if (str != "") 多邊形X坐標(biāo)集.Add(double.Parse(str); else break; Console.Write("請輸入第

33、0點的y坐標(biāo)=",i+); 多邊形Y坐標(biāo)集.Add( double.Parse(Console.ReadLine(); while (true);/計算面積double area = 0;/點的個數(shù)int Count = 多邊形X坐標(biāo)集.Count;for (int j = 0; j < Count-1; j+) area += 多邊形X坐標(biāo)集j * 多邊形Y坐標(biāo)集j + 1 - 多邊形X坐標(biāo)集j + 1 * 多邊形Y坐標(biāo)集j; /計算最后一點到起點的項 area += 多邊形X坐標(biāo)集Count - 1 * 多邊形Y坐標(biāo)集0 - 多邊形X坐標(biāo)集0 * 多邊形Y坐標(biāo)集Count

34、- 1;Console.WriteLine("多邊形的面積=0",Math.Abs(area * 0.5); （9）在如圖所示的僅有一個連接角的附合導(dǎo)線中，已知A點和B點的高程及每條邊的斜邊、垂直角、儀器高和反光鏡高，求每個點的高程。（提示：，D-斜邊，-垂直角，-儀器高，-反光鏡高，-大氣折光系數(shù)）using System;usingSystem.Collections.Generic;using System.Linq;using System.Text;namespace導(dǎo)線classProgram staticvoidMain(string args) Consol

35、e.Write("請輸入A點的高程=");double HA = double.Parse(Console.ReadLine();Console.Write("請輸入B點的高程=");double HB = double.Parse(Console.ReadLine();List<double> 高差集合 = newList<double>();List<double> 水平距離平方集合 = newList<double>();int i = 1;double高差和 = 0;double水平距離平方和 =

36、0;do Console.Write("請輸入0點到1點的斜距D<直接回車結(jié)束輸入>=",i,i+1);string str = Console.ReadLine();if (str = "") break;double D = double.Parse(str);Console.Write("請輸入0點到1點的垂直角=", i, i + 1);double = 導(dǎo)線.DEG(double.Parse(Console.ReadLine();/并化為弧度Console.Write("請輸入0的儀器高=",

37、 i);double儀器高 = double.Parse(Console.ReadLine();Console.Write("請輸入0的反光鏡高=", i+1);double反光鏡高 = double.Parse(Console.ReadLine(); i+;double水平距離;double h = 導(dǎo)線.高差(D, , 儀器高, 反光鏡高, out水平距離);水平距離平方集合.Add(水平距離 * 水平距離);高差集合.Add(h);水平距離平方和 += 水平距離 * 水平距離;高差和 += h; while (true);/計算高差閉合差double高差閉合差 = H

38、A + 高差和 - HB;/計算單位長度改正數(shù),三角高程的權(quán)與距離的平方成反比double單位長度改正數(shù) = -高差閉合差 / 水平距離平方和;/計算每點高程double H = HA;for (int j = 0; j < 高差集合.Count; j+) H += 高差集合j + 單位長度改正數(shù) * 水平距離平方集合j;Console.WriteLine("第0點的高程H=1",j+2,H); Console.WriteLine("最后一個點的高程應(yīng)該與B點的高程相同!"); class導(dǎo)線 /將ddd.mmss轉(zhuǎn)為弧度staticpublicd

39、ouble 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; /高差計算,返回高差publicstaticdouble高差(double斜距,double垂直角,double儀器高

40、,double反光鏡高,outdouble水平距離) double h = 斜距 * Math.Sin(垂直角) + 斜距 * Math.Cos(垂直角) / 6371000.0 * 斜距 * Math.Cos(垂直角) * (1 - 0.13) + 儀器高 - 反光鏡高;水平距離 = 斜距*Math.Cos(垂直角);return h; （10）在如圖所示的僅有一個連接角的附合導(dǎo)線中，已知M點、A點和B點的坐標(biāo)，每個點的左角及每條邊的水平距離，求每個點的坐標(biāo)。using System;using System.Collections.Generic;using System.Text;nam

41、espace導(dǎo)線classProgram staticvoidMain(string args) Console.Write("請輸入M點的X坐標(biāo)=");double MX = double.Parse(Console.ReadLine();Console.Write("請輸入M點的Y坐標(biāo)=");double MY = double.Parse(Console.ReadLine();Console.Write("請輸入A點的X坐標(biāo)=");double AX = double.Parse(Console.ReadLine();Conso

42、le.Write("請輸入A點的Y坐標(biāo)=");double AY = double.Parse(Console.ReadLine();Console.Write("請輸入B點的X坐標(biāo)=");double BX = double.Parse(Console.ReadLine();Console.Write("請輸入B點的Y坐標(biāo)=");double BY = double.Parse(Console.ReadLine();/計算M到A的方位角double0 = 導(dǎo)線.方位角(MX, MY, AX, AY);List<double&g

43、t; 坐標(biāo)增量X集合 = newList<double>();List<double> 坐標(biāo)增量Y集合 = newList<double>();int i = 1;double坐標(biāo)增量X和 = 0;double坐標(biāo)增量Y和 = 0;double坐標(biāo)增量X絕對值和 = 0;double坐標(biāo)增量Y絕對值和 = 0;double水平距離和 = 0;do Console.Write("請輸入0點到1點的水平距離S<直接回車結(jié)束輸入>=",i,i+1);string str = Console.ReadLine();if (str =

44、"") break;double S = double.Parse(str);Console.Write("請輸入0點的水平角<左角為正,右角為負(fù)>=", i);double = 導(dǎo)線.DEG(double.Parse(Console.ReadLine(); i+;/計算坐標(biāo)增量double dx, dy;0 = 導(dǎo)線.坐標(biāo)(0, 0, , S, 0, out dx, out dy);坐標(biāo)增量X集合.Add(dx);坐標(biāo)增量Y集合.Add(dy);坐標(biāo)增量X和 += dx; ;坐標(biāo)增量Y和 += dy; ;坐標(biāo)增量X絕對值和 += Math.

45、Abs(dx);坐標(biāo)增量Y絕對值和 += Math.Abs(dy);水平距離和 += S; while (true);/計算閉合差double閉合差X = AX + 坐標(biāo)增量X和 - BX;double閉合差Y = AY + 坐標(biāo)增量Y和 - BY;/計算單位長度改正數(shù)double X坐標(biāo)單位長度改正數(shù) = -閉合差X / 坐標(biāo)增量X絕對值和;double Y坐標(biāo)單位長度改正數(shù) = -閉合差Y / 坐標(biāo)增量Y絕對值和;/計算導(dǎo)線相對精度int f = (int)(水平距離和 / Math.Sqrt(閉合差X * 閉合差X + 閉合差Y * 閉合差Y);Console.WriteLine(&qu

46、ot;導(dǎo)線的相對精度為1/0",f);/計算每點坐標(biāo)double X = AX;double Y = AY;for (int j = 0; j < 坐標(biāo)增量X集合.Count; j+) X += 坐標(biāo)增量X集合j + X坐標(biāo)單位長度改正數(shù) * Math.Abs(坐標(biāo)增量X集合j); Y += 坐標(biāo)增量Y集合j + Y坐標(biāo)單位長度改正數(shù) * Math.Abs(坐標(biāo)增量Y集合j);Console.WriteLine("第0點的X坐標(biāo)=1,Y坐標(biāo)=2",j+2,X,Y); Console.WriteLine("最后一個點的坐標(biāo)應(yīng)該與B點的坐標(biāo)相同!&qu

47、ot;); class導(dǎo)線 /將ddd.mmss轉(zhuǎn)為弧度staticpublicdouble 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; /計算方位角,返回弧度值public

48、staticdouble方位角(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 < 0) angle += Math.PI; if (deltaY < 0.0) angle += Math.PI; return angle;

49、 /計算坐標(biāo),返回已知點到計算點的方位角publicstaticdouble坐標(biāo)(double x0, double y0, double左角, double水平距離,double已知方位角,outdouble x,outdouble y) double方位角 = 已知方位角 + 左角 + Math.PI; /將方位角調(diào)整到0到2之間if (方位角 >= Math.PI * 2) 方位角 -= Math.PI * 2;if (方位角 < 0.0) 方位角 += Math.PI * 2; x = x0 + 水平距離 * Math.Cos(方位角); y = y0 + 水平距離 * M

50、ath.Sin(方位角);return方位角; （11）在如圖所示的僅有一個連接角的附合導(dǎo)線中，已知M點、A點和B點的坐標(biāo)，A點和B點的高程，每個點的左角及每條邊的斜邊、垂直角、儀器高和反光鏡高，求每個點的坐標(biāo)和高程。（提示：，D-斜邊，-垂直角，-儀器高，-反光鏡高，-大氣折光系數(shù)）using System;using System.Collections.Generic;using System.Text;namespace導(dǎo)線classProgram staticvoidMain(string args) Console.Write("請輸入M點的X坐標(biāo)=");double MX = double.Parse(Console.ReadLine();Console.Write("請輸入M點的Y坐標(biāo)=");double MY = double.Parse(Console.ReadLine();Console.Write("請輸入A點的X坐標(biāo)=");double AX = double.Pars