测绘类C程序代码

常用测量程序设计代码（1）用全站仪在A点观测了B点斜边和垂直角，求A到B的高差。

（提示：22sin(1)cos2ABDh D a k a i vR=+-+-，D--斜边，a--垂直角，i--仪器高，v--反光镜高，k--大气折光系数）using System;using ;namespace ConsoleApplication1{class Application{static void Main(string[] args){("请输入斜边=");double D = ());("请输入垂直角[]=");double a = DEG()));("请输入仪器高=");double i = ());("请输入反光镜高=");double v = ());double h = D * (a) + (1 - * D / * D / * (a) * (a)/ + i - v;("高差为{0}",h);}MainMainaMaindd(DEG(str)));}else{break;}} while (true);i = 1;foreach (double a in导线转角集合) {a0 += a + ;if (a0 > 2 *{a0 -= 2 * ;}else if (a0 <{a0 += 2 * ;}("第{0}条边的方位角为{1}",i++,DMS(a0)); }}Maindd(DEG(str)));}else{break;}("请输入第{0}条边长值=", i++);导线边长集合.Add( ()));} while (true);ount; j++){a0 += 导线转角集合[j] + ;if (a0 > 2 * a0 -= 2 * ;x0 = x0 + 导线边长集合[j] * (a0);y0 = y0 + 导线边长集合[j] * (a0);("P{0}点的坐标是:{1},{2}",j+2,x0,y0); }}Maindd(str));}else{break;}("请输入第{0}段高差=", i++);路线段高差集合.Add()));} while (true);ount; j++){w += 路线段高差集合[j];总长度 += 路线段长度集合[j];}w = w - (Hb - Ha); ount; j++){H += 路线段高差集合[j] + v * 路线段长度集合[j]; ("P{0}点的高程是:{1}",j+1,H);}("最后一点就是B点！");}}}（7）在如图所示的前方交会中，ABJ三点按逆时针方向排列，已知AB两点的坐标和两个交会角度，求J点的坐标。

using System;using ;namespace ConsoleApplication1{class Application{static void Main(string[] args) {("请输入A点的x坐标=");double Ax = ());("请输入A点的y坐标=");double Ay = ());("请输入B点的x坐标=");double Bx = ());("请输入B点的y坐标=");double By = ());("请输入α的角度值=");doubleα = DEG()));("请输入β的角度值=");doubleβ = DEG()));1111()2ni i i i i P x y x y ++==-∑Maindd(str));}else{break ;}("请输入第{0}点的y 坐标=",i++);多边形Y 坐标集.Add( ())); } while (true );ount;for (int j = 0; j < Count-1; j++){area += 多边形X 坐标集[j] * 多边形Y 坐标集[j + 1] - 多边形X 坐标集[j + 1] * 多边形Y 坐标集[j];}22sin (1)cos 2AB D h D a k a i v R=+-+-a i v kMainEG()));差(D, α, 仪器高, 反光镜高, out 水平距离);水平距离平方集合.Add(水平距离 * 水平距离); 高差集合.Add(h);水平距离平方和 += 水平距离 * 水平距离;高差和 += h;} while (true );ount; j++){H += 高差集合[j] + 单位长度改正数 * 水平距离平方集合[j];("第{0}点的高程H={1}",j+2,H);}("最后一个点的高程应该与B 点的高程相同!");}}class 导线{iMain位角(MX, MY, AX, AY);List<double> 坐标增量X集合 = new List<double>();List<double> 坐标增量Y集合 = new List<double>();int i = 1;double坐标增量X和 = 0;double坐标增量Y和 = 0;double坐标增量X绝对值和 = 0;double坐标增量Y绝对值和 = 0;double水平距离和 = 0;do{("请输入{0}点到{1}点的水平距离S<直接回车结束输入>=",i,i+1);string str = ();if (str == "") break;double S = (str);("请输入{0}点的水平角<左角为正,右角为负>=", i);doubleβ = 导线.DEG()));i++;标(0, 0, β, S, α0, out dx, out dy);坐标增量X集合.Add(dx);坐标增量Y集合.Add(dy);坐标增量X和 += dx; ;坐标增量Y和 += dy; ;坐标增量X绝对值和 += (dx);坐标增量Y绝对值和 += (dy);水平距离和 += S;} while (true);ount; j++){X += 坐标增量X集合[j] + X坐标单位长度改正数 * (坐标增量X集合[j]);Y += 坐标增量Y集合[j] + Y坐标单位长度改正数 * (坐标增量Y集合[j]);("第{0}点的X坐标={1},Y坐标={2}",j+2,X,Y);}("最后一个点的坐标应该与B点的坐标相同!");}}class 导线{i β22sin (1)cos 2AB D h D a k a i v R=+-+-a i v kMain 位角(MX, MY, AX, AY);List <double > 坐标增量X 集合 = new List <double >(); List <double > 坐标增量Y 集合 = new List <double >(); List <double > 高差集合 = new List <double >();List <double > 水平距离平方集合 = new List <double >(); int i = 1;double 坐标增量X 和 = 0;double 坐标增量Y 和 = 0;double 坐标增量X 绝对值和 = 0;double 坐标增量Y 绝对值和 = 0;double 水平距离和 = 0;double 高差和 = 0;double 水平距离平方和 = 0;do{("请输入{0}点到{1}点的斜距D<直接回车结束输入>=", i, i + 1);string str = ();if (str == "") break;double D = (str);("请输入{0}点的水平角<左角为正,右角为负>=", i);doubleβ = 导线.DEG()));("请输入{0}点到{1}点的垂直角=", i, i + 1);doubleα = 导线.DEG()));差(D, α, 仪器高, 反光镜高, out水平距离);水平距离平方集合.Add(水平距离 * 水平距离);高差集合.Add(h);水平距离平方和 += 水平距离 * 水平距离;高差和 += h;标(0, 0, β, 水平距离, α0, out dx, out dy);坐标增量X集合.Add(dx);坐标增量Y集合.Add(dy);坐标增量X和 += dx; ;坐标增量Y和 += dy; ;坐标增量X绝对值和 += (dx);坐标增量Y绝对值和 += (dy);水平距离和 += 水平距离;} while (true );ount; j++){H += 高差集合[j] + 单位长度改正数 * 水平距离平方集合[j];X += 坐标增量X 集合[j] + X 坐标单位长度改正数 * (坐标增量X 集合[j]);Y += 坐标增量Y 集合[j] + Y 坐标单位长度改正数 * (坐标增量Y 集合[j]);("第{0}点的X 坐标={1},Y 坐标={2},高程H={3}",j+2,X,Y,H);}("最后一个点的坐标和高程应该与B 点的坐标和高程相同!");}}class 导线{i β22sin (1)cos 2AB D h D a k a i v R=+-+-a i v kMain位角(MX, MY, AX, AY);位角(BX, BY, NX, NY);List<double> 坐标增量X集合 = new List<double>();List<double> 坐标增量Y集合 = new List<double>();List<double> 高差集合 = new List<double>();List<double> 水平距离平方集合 = new List<double>();List<double> 方位角集合 = new List<double>();int i = 1;double坐标增量X和 = 0;double坐标增量Y和 = 0;double坐标增量X绝对值和 = 0;double坐标增量Y绝对值和 = 0;double水平距离和 = 0;double高差和 = 0;double水平距离平方和 = 0;doubleα = α0;do{("请输入{0}点到{1}点的斜距D<直接回车结束输入>=",i, i + 1);string str = ();if (str == "") break;double D = (str);("请输入{0}点的水平角<左角为正,右角为负>=", i);doubleβ = 导线.DEG()));("请输入{0}点到{1}点的垂直角=", i, i + 1);double垂直角 = 导线.DEG()));差(D, 垂直角, 仪器高, 反光镜高, out水平距离);水平距离平方集合.Add(水平距离 * 水平距离);高差集合.Add(h);水平距离平方和 += 水平距离 * 水平距离;高差和 += h;标(0, 0, β, 水平距离, α, out dx, out dy);方位角集合.Add(α);坐标增量X集合.Add(dx);坐标增量Y集合.Add(dy);坐标增量X和 += dx; ;坐标增量Y和 += dy; ;坐标增量X绝对值和 += (dx);坐标增量Y绝对值和 += (dy);水平距离和 += 水平距离;} while (true);("请输入最后一个连接角<左角为正,右角为负>=");doubleβn = 导线.DEG()));位角(方位角集合[方位角集合.Count-1],βn);ount + 1);lear();坐标增量Y集合.Clear();坐标增量X和 = 0;坐标增量Y和 = 0;坐标增量X绝对值和 = 0;坐标增量Y绝对值和 = 0;for (int j = 0; j < 方位角集合.Count; j++){方位角集合[j] += 方位角改正数 * (j + 1); dd(dx);坐标增量Y集合.Add(dy);坐标增量X和 += dx; ;坐标增量Y和 += dy; ;坐标增量X绝对值和 += (dx);坐标增量Y绝对值和 += (dy);}ount; j++){H += 高差集合[j] + 单位长度改正数 * 水平距离平方集合[j];X += 坐标增量X集合[j] + X坐标单位长度改正数 * (坐标增量X集合[j]);Y += 坐标增量Y集合[j] + Y坐标单位长度改正数 * (坐标增量Y集合[j]);("第{0}点的X坐标={1},Y坐标={2},高程H={3}",j+2,X,Y,H);}("最后一个点的坐标和高程应该与B点的坐标和高程相同!");}}class导线{iMainEG()));("请输入B点到N点的坐标方位角=");doubleα1 = 导线.DEG()));List<double> 方位角集合 = new List<double>();int i = 1;doubleα = α0;do{("请输入{0}点的水平角[直接回车结束输入]<左角为正,右角为负>=", i);string str = ();if (str == "") break;doubleβ = 导线.DEG(str));i++;α = 导线.方位角(α,β);方位角集合.Add(α);} while (true);ount-1] - α1;double方位角改正数 = -方位角闭合差 / 方位角集合.Count;ount; j++){方位角集合[j] += 方位角改正数 * (j + 1);("{0}点到{1}点的坐标方位角={2}", j, j + 2, 导线.DMS(方位角集合[j]));}("最后一个坐标方位角应该与B点到N的坐标方位角相同!");}}class导线{iMain位角(MX, MY, AX, AY);位角(BX, BY, NX, NY);List<double> 坐标增量X集合 = new List<double>();List<double> 坐标增量Y集合 = new List<double>();List<double> 水平距离平方集合 = new List<double>();List<double> 方位角集合 = new List<double>();int i = 1;double坐标增量X和 = 0;double坐标增量Y和 = 0;double坐标增量X绝对值和 = 0;double坐标增量Y绝对值和 = 0;double水平距离和 = 0;double水平距离平方和 = 0;doubleα = α0;do{("请输入{0}点到{1}点的平距S<直接回车结束输入>=", i, i + 1);string str = ();if (str == "") break;double S = (str);("请输入{0}点的水平角<左角为正,右角为负>=", i);doubleβ = 导线.DEG()));i++;水平距离平方集合.Add(S * S);水平距离平方和 += S * S;标(0, 0, β, S, α, out dx, out dy);方位角集合.Add(α);坐标增量X集合.Add(dx);坐标增量Y集合.Add(dy);坐标增量X和 += dx; ;坐标增量Y和 += dy; ;坐标增量X绝对值和 += (dx);坐标增量Y绝对值和 += (dy);} while (true);("请输入最后一个连接角<左角为正,右角为负>=");doubleβn = 导线.DEG()));位角(方位角集合[方位角集合.Count-1],βn);ount + 1);lear();坐标增量Y集合.Clear();坐标增量X和 = 0;坐标增量Y和 = 0;坐标增量X绝对值和 = 0;坐标增量Y绝对值和 = 0;for (int j = 0; j < 方位角集合.Count; j++){方位角集合[j] += 方位角改正数 * (j + 1); dd(dx);坐标增量Y集合.Add(dy);坐标增量X和 += dx; ;坐标增量Y和 += dy; ;坐标增量X绝对值和 += (dx);坐标增量Y绝对值和 += (dy);}ount; j++){X += 坐标增量X集合[j] + X坐标单位长度改正数 * (坐标增量X集合[j]);Y += 坐标增量Y集合[j] + Y坐标单位长度改正数 * (坐标增量Y集合[j]);("第{0}点的X坐标={1},Y坐标={2}",j+2,X,Y);}("最后一个点的坐标应该与B点的坐标相同!");}}class导线{//将转为弧度static public double DEG(double ang){int fuhao = (int)(ang / (ang));ang = (ang);int d = (int)ang;int m = ((int)(ang * 100)) - d * 100;double s = ang * 10000 - m * 100 - d * 10000;return ((d + m / + s / * fuhao) / * ;}//计算方位角,返回弧度值public static double方位角(double x1, double y1, double x2, double y2){double deltaX = x2 - x1;double deltaY = y2 - y1;double angle = * ;if (deltaX) >{angle = (deltaY, deltaX);}if (angle < 0){angle += ;}if (deltaY <{angle += ;}return angle;}//计算坐标,返回已知点到计算点的方位角public static double坐标(double x0, double y0, double左角, double水平距离,double已知方位角,out double x,out double y) {double方位角 = 已知方位角 + 左角 + ;//将方位角调整到0到2π之间if (方位角 >= * 2) 方位角 -= * 2;if (方位角 < 方位角 += * 2;x = x0 + 水平距离 * (方位角);y = y0 + 水平距离 * (方位角);return方位角;}//根据后视边的方位角与左角,计算前进边的方位角public static double方位角(double后视边方位角, double 左角){double方位角 = 后视边方位角 + 左角 + ;//将方位角调整到0到2π之间if (方位角 >= * 2) 方位角 -= * 2;if (方位角 < 方位角 += * 2;return方位角;}}}。