文章目錄

• 一、前言
• 二、代碼界面展示
• 三、代碼運算結(jié)果展示
• • 3.1 閉合導(dǎo)線
  • • 3.1.1 通過txt文件導(dǎo)入原始數(shù)據(jù)
    • 3.1.2 設(shè)置
    • 3.1.3 計算
    • 3.1.4 繪圖
    • 3.1.5 數(shù)據(jù)導(dǎo)出
    • 3.1.6 整體展示
  • 3. 2 附和導(dǎo)線
• 四、源代碼
• 五、結(jié)語

一、前言

本程序是使用C# 編寫的一個附和閉合導(dǎo)線平差程序
閉合導(dǎo)線：兩個已知點，一個已知點+方位角；
附和導(dǎo)線：雙定向附和導(dǎo)線；
前面都不是重點，有需要的直接看代碼

二、代碼界面展示

Form1（大窗體）中控件有：Menu，tabControl，dataGridView。
Form2（小窗體）中控件有：radioButton，comboBox，label，textbox button。

三、代碼運算結(jié)果展示

3.1 閉合導(dǎo)線

3.1.1 通過txt文件導(dǎo)入原始數(shù)據(jù)

沒有的數(shù)據(jù)我用的0代替，沿用的水準測量程序的讀取代碼。
我這里的數(shù)據(jù)是直接給出了水平角了。
有些同學實習給的原始數(shù)據(jù)是盤左盤右數(shù)據(jù)。如果想要程序計算，需要多設(shè)置個幾變量，如果想要數(shù)據(jù)顯示出來的話，表格，數(shù)據(jù)導(dǎo)入那里也需要修改一下，但是運算邏輯是沒有變的，自己改一下應(yīng)該不難。
實在不會的話，建議自己手動先把水平角算出來，再使用程序計算。

數(shù)據(jù)一：（2個已知點-右角）

數(shù)據(jù)二：（2個已知點-左角）
數(shù)據(jù)三：（方位角(48.4319)-右角）
數(shù)據(jù)四：（方位角（100）-右角）

3.1.2 設(shè)置

數(shù)據(jù)一：（2個已知點-右角）

數(shù)據(jù)二：（2個已知點-左角）

數(shù)據(jù)三：（方位角(48.4319)-右角）

數(shù)據(jù)四：（方位角（100）-右角）

3.1.3 計算

數(shù)據(jù)一：（2個已知點-右角）
計算結(jié)果和工程測量大師的是差不多的，mm位上有一些數(shù)據(jù)有差別，應(yīng)該主要是我整個計算的過程中都沒有去避免double的精度誤差（不明白精度誤差的可以去搜搜看看）。如果你要參考我的代碼可以自己添加使用round來保留位數(shù)。
這個界面上面是datagridview，下面是richtextbox，表格設(shè)計成這樣其實不符合導(dǎo)線平差計算表的樣式的。主要是自己沒去學習表格的設(shè)計，能力不足，所以按照工程測量大師的表格設(shè)計的表。限差這些我隨便找的一個標準。但是這個數(shù)據(jù)本身計算出來就是不合格的。
數(shù)據(jù)二：（2個已知點-左角）
與他自己計算的數(shù)據(jù)可以進行對比
數(shù)據(jù)差別的來源
1.角度改正數(shù)的分配，我程序分配的邏輯是先均分，多出來的分配給邊長較長的。
2.位數(shù)保留，和數(shù)據(jù)一 一樣
不過最后計算的B點坐標，程序和原始數(shù)據(jù)一樣。

數(shù)據(jù)三：（方位角(48.4319)-右角）
數(shù)據(jù)四：（方位角（100）-右角）

3.1.4 繪圖

這個功能只能大概畫出一個圖形，電腦畫圖的坐標系和測量坐標系不一樣，所以需要轉(zhuǎn)換。程序畫的時候已經(jīng)轉(zhuǎn)換了
說明一下，這個程序能畫的坐標范圍為x∈[-224,224],y∈[-496,496]，如果坐標不屬于這里面的就畫不出來。坐標太大的話需要你自己在繪圖那段代碼里面修改n，讓坐標同時縮小n倍。
但是如果你的坐標是（3465128.993,552965.750）這種的，x坐標和y坐標位數(shù)不一樣，設(shè)置n就沒有用了，就畫不出來了。
當然你可以自己發(fā)力，畫出圖形來。
數(shù)據(jù)一：（2個已知點+右角）

n=5

數(shù)據(jù)二：（2個已知點-左角）
n=5
畫出來確實像一個廠房的結(jié)構(gòu)
數(shù)據(jù)三：（方位角(48.4319)-右角）
n=5

數(shù)據(jù)四：（方位角（100）-右角）
n=5

3.1.5 數(shù)據(jù)導(dǎo)出

這個沒什么好聊的，就是把計算的數(shù)據(jù)導(dǎo)出成excel
數(shù)據(jù)一：（2個已知點+右角）

3.1.6 整體展示

3. 2 附和導(dǎo)線

這里我就不過多描述了，就是計算雙定向附和導(dǎo)線
實驗數(shù)據(jù)都是搜索的，經(jīng)過了驗證的。
數(shù)據(jù)一


數(shù)據(jù)二

四、源代碼

namespace 附和閉合導(dǎo)線平差程序 {public partial class Form1 : Form{public Form1(){InitializeComponent();}//畫圖設(shè)置Graphics g;int midx, midy;List<Point> curvelist;Form2 form2 = new Form2();List<RowData> rdatas = new List<RowData>();List<ResultData> rtdatas = new List<ResultData>();ResultData rtdata = new ResultData();RowData rdata = new RowData();Tolerance tolerance = new Tolerance();double D = 0, M = 0, S = 0;double d1 = 0, m1 = 0, s1 = 0;double d2 = 0, m2 = 0;private void 文件ToolStripMenuItem1_Click(object sender, EventArgs e){tabControl1.SelectedTab = tabControl1.TabPages[0];dataGridView1.Rows.Clear();rdatas.Clear();string[] temp;OpenFileDialog ofd = new OpenFileDialog();ofd.Filter = "文本文件(*.txt)|*.txt|所有文件(*.*)|*.*";ofd.Title = "打開數(shù)據(jù)文件";if (ofd.ShowDialog(this) == DialogResult.OK){StreamReader sr = new StreamReader(ofd.FileName);while (!sr.EndOfStream){temp = sr.ReadLine().Split(new[] { ',', ' ' }, StringSplitOptions.RemoveEmptyEntries);RowData rdata = new RowData();rdata.DH = temp[0];rdata.DMS = Convert.ToDouble(temp[1]);rdata.L = Convert.ToDouble(temp[2]);rdata.X = Convert.ToDouble(temp[3]);rdata.Y = Convert.ToDouble(temp[4]);rdatas.Add(rdata);}sr.Close();GridBuild1();}else{MessageBox.Show("輸入失敗");}}public void GridBuild1(){for (int i = 0; i < rdatas.Count; i++){dataGridView1.Rows.Add();dataGridView1.Rows[i].Cells[0].Value = rdatas[i].DH;dataGridView1.Rows[i].Cells[1].Value = rdatas[i].DMS;dataGridView1.Rows[i].Cells[2].Value = rdatas[i].L;dataGridView1.Rows[i].Cells[3].Value = rdatas[i].X;dataGridView1.Rows[i].Cells[4].Value = rdatas[i].Y;}}private void 設(shè)置ToolStripMenuItem_Click(object sender, EventArgs e){form2.ShowDialog();}private void 計算ToolStripMenuItem_Click(object sender, EventArgs e){tabControl1.SelectedTab = tabControl1.TabPages[1];ClearData();ShowHAandVandCHA();ShowL();ShowCA();ShowΔxandΔy();ShowCΔxandCΔy();ShowXandY();ShowTolerance();Tolerance tolerance1 = new Tolerance();}public void ShowHAandVandCHA(){//閉合if (form2.radioButton1.Checked){int n = 0, i = 0;if (form2.comboBox1.Text == "兩個已知點"){n = rdatas.Count - 2;i = 2;}else if (form2.comboBox1.Text == "一個已知點+方位角"){n = rdatas.Count - 1;i = 1;}double sumli = 180 * (n - 2);while (i < rdatas.Count){D += ReturnDMS(rdatas[i].DMS).d;M += ReturnDMS(rdatas[i].DMS).m;S += ReturnDMS(rdatas[i].DMS).s;i++;}tolerance.fβ = (D + M / 60.0 + S / 3600.0 - sumli) * 3600;for (i = 0; i < rdatas.Count; i++){dataGridView2.Rows.Add();dataGridView2.Rows[i].Cells[0].Value = rtdatas[i].DH;if (i == 0){rtdatas[i].X = rdatas[i].X;rtdatas[i].Y = rdatas[i].Y;dataGridView2.Rows[i].Cells[10].Value = rtdatas[i].X;dataGridView2.Rows[i].Cells[11].Value = rtdatas[i].Y;}else if (i == 1 && form2.comboBox1.Text == "兩個已知點"){rtdatas[i].X = rdatas[i].X;rtdatas[i].Y = rdatas[i].Y;dataGridView2.Rows[i].Cells[10].Value = rtdatas[i].X;dataGridView2.Rows[i].Cells[11].Value = rtdatas[i].Y;d1 = ReturnDMS(rdatas[i].DMS).d;m1 = ReturnDMS(rdatas[i].DMS).m;s1 = ReturnDMS(rdatas[i].DMS).s;double HAb = toDu(d1, m1, s1);rtdatas[i].HA = HAb;dataGridView2.Rows[i].Cells[1].Value = d1.ToString() + " " + m1.ToString() + " " + s1.ToString();}else{d1 = ReturnDMS(rdatas[i].DMS).d;m1 = ReturnDMS(rdatas[i].DMS).m;s1 = ReturnDMS(rdatas[i].DMS).s;double HAx = toDu(d1, m1, s1);rtdatas[i].HA = HAx;double S = VDistribution(s1, i);d2 = d1;m2 = m1;if (S < 0){S += 60;m2 = m1 - 1;}else if (S >= 60){S -= 60;m2 = m1 + 1;}else if (m2 < 0){m2 += 60;d2 = d1 - 1;}else if (m2 >= 60){m2 -= 60;d2 = d1 + 1;}double CHAx = toDu(d2, m2, S);rtdatas[i].CHA = CHAx;dataGridView2.Rows[i].Cells[1].Value = d1.ToString() + " " + m1.ToString() + " " + s1.ToString();dataGridView2.Rows[i].Cells[2].Value = Math.Round(rtdatas[i].V, 0).ToString();dataGridView2.Rows[i].Cells[3].Value = d2.ToString() + " " + m2.ToString() + " " + Math.Round(S, 0).ToString();}}}//附和else if (form2.radioButton2.Checked){double Xa = rdatas[0].X;double Ya = rdatas[0].Y;double Xb = rdatas[1].X;double Yb = rdatas[1].Y;double Theta1 = Math.Atan(Math.Abs((Yb - Ya) / (Xb - Xa)));double radCA = JudgementQuadrant(Theta1, 1, 0);double CAa = radCA * 180 / Math.PI;rtdatas[0].CA = CAa;double Xc = rdatas[rtdatas.Count - 2].X;double Yc = rdatas[rtdatas.Count - 2].Y;double Xd = rdatas[rtdatas.Count - 1].X;double Yd = rdatas[rtdatas.Count - 1].Y;double Theta2 = Math.Atan(Math.Abs((Yd - Yc) / (Xd - Xc)));double radCc = JudgementQuadrant(Theta2, rtdatas.Count - 1, rtdatas.Count - 2);double CAc = radCc * 180 / Math.PI;rtdatas[rtdatas.Count - 2].CA = CAc;for (int i = 1; i < rtdatas.Count - 1; i++){D += ReturnDMS(rdatas[i].DMS).d;M += ReturnDMS(rdatas[i].DMS).m;S += ReturnDMS(rdatas[i].DMS).s;}if (form2.comboBox2.Text=="右角"){double βli = rtdatas[0].CA - rtdatas[rtdatas.Count - 2].CA + ((rtdatas.Count - 2) * 180);double βce = D + M / 60.0 + S / 3600.0;tolerance.fβ = βce - βli;}else if (form2.comboBox2.Text == "左角"){double βli = rtdatas[rtdatas.Count - 2].CA - rtdatas[0].CA+360 + ((rtdatas.Count - 2) * 180);double βce = D + M / 60.0 + S / 3600.0;tolerance.fβ = βce - βli; }while (tolerance.fβ > 360){tolerance.fβ -= 360;}while (tolerance.fβ < -360){tolerance.fβ += 360;}tolerance.fβ =Math.Round( tolerance.fβ * 3600,0);for (int i = 0; i < rdatas.Count; i++){dataGridView2.Rows.Add();dataGridView2.Rows[i].Cells[0].Value = rtdatas[i].DH;if (i == 0 || i == rdatas.Count - 1){rtdatas[i].X = rdatas[i].X;rtdatas[i].Y = rdatas[i].Y;dataGridView2.Rows[i].Cells[10].Value = rtdatas[i].X;dataGridView2.Rows[i].Cells[11].Value = rtdatas[i].Y;}else{if (i == 1 || i == rtdatas.Count - 2){rtdatas[i].X = rdatas[i].X;rtdatas[i].Y = rdatas[i].Y;dataGridView2.Rows[i].Cells[10].Value = rtdatas[i].X;dataGridView2.Rows[i].Cells[11].Value = rtdatas[i].Y;}d1 = ReturnDMS(rdatas[i].DMS).d;m1 = ReturnDMS(rdatas[i].DMS).m;s1 = ReturnDMS(rdatas[i].DMS).s;double HAx = toDu(d1, m1, s1);rtdatas[i].HA = HAx;double S = VDistribution(s1, i);d2 = d1;m2 = m1;if (S < 0){S += 60;m2 =m1- 1;}else if (S>=60){S -= 60;m2=m1+ 1;}else if (m2 < 0){m2 += 60;d2 =d1- 1;}else if (m2>=60){m2 -= 60;d2 =d1+ 1;}double CHAx = toDu(d2, m2, S);rtdatas[i].CHA = CHAx;dataGridView2.Rows[i].Cells[1].Value = d1.ToString() + " " + m1.ToString() + " " + s1.ToString();dataGridView2.Rows[i].Cells[2].Value = Math.Round(rtdatas[i].V, 0).ToString();dataGridView2.Rows[i].Cells[3].Value = d2.ToString() + " " + m2.ToString() + " " + Math.Round(S, 0).ToString();}}}}//點號和水平角和改正后的水平角public void ShowL(){if (form2.radioButton1.Checked){int n = 0;if (form2.comboBox1.Text == "兩個已知點"){n = rdatas.Count - 1;}else if (form2.comboBox1.Text == "一個已知點+方位角"){n = rdatas.Count;}for (int i = 1; i < n; i++){rtdatas[i].L = rdatas[i].L;tolerance.fl += rtdatas[i].L;dataGridView2.Rows[i].Cells[5].Value = rtdatas[i].L;}}else if (form2.radioButton2.Checked){for (int i = 1; i < rdatas.Count - 2; i++){rtdatas[i].L = rdatas[i].L;tolerance.fl += rtdatas[i].L;dataGridView2.Rows[i].Cells[5].Value = rtdatas[i].L;}}}//展示距離public void ShowCA(){//閉合if (form2.radioButton1.Checked){int n = 0;if (form2.comboBox1.Text == "兩個已知點"){n = rdatas.Count - 1;}else if (form2.comboBox1.Text == "一個已知點+方位角"){n = rdatas.Count;}//已知點坐標（X1，Y1）,(X2，Y2)double CAa = 0;double CAb = 0;for (int i = 0; i < n; i++){if (i == 0 && form2.comboBox1.Text == "兩個已知點"){double X1 = rdatas[0].X;double Y1 = rdatas[0].Y;double X2 = rdatas[1].X;double Y2 = rdatas[1].Y;double Theta = Math.Atan(Math.Abs((Y2 - Y1) / (X2 - X1)));double radCA = JudgementQuadrant(Theta, 1, 0);CAa = radCA * 180 / Math.PI;rtdatas[i].CA = CAa;int[] dms = toDMS(rtdatas[i].CA);dataGridView2.Rows[0].Cells[4].Value = dms[0].ToString() + " " + dms[1].ToString() + " " + dms[2].ToString();}else if (i == 0 && form2.comboBox1.Text == "一個已知點+方位角"){CAa = Convert.ToDouble(form2.textBox1.Text) + Convert.ToDouble(form2.textBox2.Text) / 60 + Convert.ToDouble(form2.textBox3.Text) / 3600;rtdatas[i].CA = CAa;int[] dms = toDMS(rtdatas[i].CA);dataGridView2.Rows[0].Cells[4].Value = dms[0].ToString() + " " + dms[1].ToString() + " " + dms[2].ToString();}else if (i == 1 && form2.comboBox1.Text == "兩個已知點"){if (form2.comboBox2.Text == "右角"){CAb = CAa - rtdatas[i].HA + 180;}if (form2.comboBox2.Text == "左角"){CAb = CAa + rtdatas[i].HA + 180;}if (CAb > 360){CAb -= 360;}if (CAb<0){CAb += 360;}rtdatas[i].CA = CAb;int[] dms1 = toDMS(rtdatas[i].CA);dataGridView2.Rows[i].Cells[4].Value = dms1[0].ToString() + " " + dms1[1].ToString() + " " + dms1[2].ToString();}else if (i == 1 && form2.comboBox1.Text == "一個已知點+方位角"){CAb = Calculate_CA(CAa, i);rtdatas[i].CA = CAb;int[] dms1 = toDMS(rtdatas[i].CA);dataGridView2.Rows[i].Cells[4].Value = dms1[0].ToString() + " " + dms1[1].ToString() + " " + dms1[2].ToString();}else{CAb = Calculate_CA(CAb, i);rtdatas[i].CA = CAb;int[] dms1 = toDMS(rtdatas[i].CA);dataGridView2.Rows[i].Cells[4].Value = dms1[0].ToString() + " " + dms1[1].ToString() + " " + dms1[2].ToString();}}}//附和else if (form2.radioButton2.Checked){double CAx=0;int[] dms = toDMS(rtdatas[0].CA);dataGridView2.Rows[0].Cells[4].Value = dms[0].ToString() + " " + dms[1].ToString() + " " + dms[2].ToString();int[] dms1= toDMS(rtdatas[rtdatas.Count-2].CA);dataGridView2.Rows[rtdatas.Count - 2].Cells[4].Value = dms1[0].ToString() + " " + dms1[1].ToString() + " " + dms1[2].ToString();for (int i = 1; i < rtdatas.Count-1; i++){if (i==1){CAx = Calculate_CA(rtdatas[0].CA, i);rtdatas[i].CA = CAx;int[] dms2 = toDMS(rtdatas[i].CA);dataGridView2.Rows[i].Cells[4].Value = dms2[0].ToString() + " " + dms2[1].ToString() + " " + dms2[2].ToString();}else{CAx = Calculate_CA(CAx, i);rtdatas[i].CA = CAx;int[] dms2 = toDMS(rtdatas[i].CA);dataGridView2.Rows[i].Cells[4].Value = dms2[0].ToString() + " " + dms2[1].ToString() + " " + dms2[2].ToString();}}}}//計算展示坐標方位角//計算坐標方位角public double Calculate_CA(double CA, int i){double CA1 = 0;if (form2.comboBox2.Text == "右角"){CA1 =CA - rtdatas[i].CHA + 180;}if (form2.comboBox2.Text == "左角"){CA1 = CA + rtdatas[i].CHA + 180;}while (CA1 >= 360){CA1 -= 360;}while (CA1 <0){CA1 += 360;}return CA1;}public void ShowΔxandΔy(){if (form2.radioButton1.Checked){int i = 0;if (form2.comboBox1.Text == "兩個已知點"){i = 2;}else if (form2.comboBox1.Text == "一個已知點+方位角"){i = 1;}while (i < rdatas.Count){if (form2.comboBox1.Text == "兩個已知點"){rtdatas[i].Δx = Math.Cos(rtdatas[i - 1].CA * Math.PI / 180) * rtdatas[i - 1].L;rtdatas[i].Δy = Math.Sin(rtdatas[i - 1].CA * Math.PI / 180) * rtdatas[i - 1].L;}else if (form2.comboBox1.Text == "一個已知點+方位角"){rtdatas[i].Δx = Math.Cos(rtdatas[i - 1].CA * Math.PI / 180) * rtdatas[i].L;rtdatas[i].Δy = Math.Sin(rtdatas[i - 1].CA * Math.PI / 180) * rtdatas[i].L;}tolerance.fx += rtdatas[i].Δx;tolerance.fy += rtdatas[i].Δy;dataGridView2.Rows[i].Cells[6].Value = Math.Round(rtdatas[i].Δx, 4);dataGridView2.Rows[i].Cells[7].Value = Math.Round(rtdatas[i].Δy, 4);i++;}tolerance.fD = Math.Sqrt(tolerance.fx * tolerance.fx + tolerance.fy * tolerance.fy);tolerance.K = tolerance.fD / tolerance.fl;}else if (form2.radioButton2.Checked){double fxce=0, fyce=0,fxli=0,fyli=0;fxli = rdatas[rdatas.Count - 2].X - rdatas[1].X;fyli = rdatas[rdatas.Count - 2].Y - rdatas[1].Y;for (int i = 1; i < rdatas.Count-2; i++){rtdatas[i].Δx = Math.Cos(rtdatas[i].CA * Math.PI / 180) * rtdatas[i].L;rtdatas[i].Δy = Math.Sin(rtdatas[i].CA * Math.PI / 180) * rtdatas[i].L;fxce += rtdatas[i].Δx;fyce += rtdatas[i].Δy;dataGridView2.Rows[i].Cells[6].Value = Math.Round(rtdatas[i].Δx, 4);dataGridView2.Rows[i].Cells[7].Value = Math.Round(rtdatas[i].Δy, 4);}tolerance.fx = fxce - fxli;tolerance.fy = fyce - fyli;tolerance.fD = Math.Sqrt(tolerance.fx * tolerance.fx + tolerance.fy * tolerance.fy);tolerance.K = tolerance.fD / tolerance.fl;}}//計算展示Δx和Δy和fd，fl，fx，fx，Kpublic void ShowCΔxandCΔy(){if (form2.radioButton1.Checked){int i = 0;double temp = 0, temp1 = 0;if (form2.comboBox1.Text == "兩個已知點"){i = 2;}else if (form2.comboBox1.Text == "一個已知點+方位角"){i = 1;}while (i < rdatas.Count){if (form2.comboBox1.Text == "兩個已知點"){temp = tolerance.fx * rtdatas[i - 1].L / tolerance.fl;temp1 = tolerance.fy * rtdatas[i - 1].L / tolerance.fl;}else if (form2.comboBox1.Text == "一個已知點+方位角"){temp = tolerance.fx * rtdatas[i].L / tolerance.fl;temp1 = tolerance.fy * rtdatas[i].L / tolerance.fl;}rtdatas[i].CΔx = rtdatas[i].Δx - temp;rtdatas[i].CΔy = rtdatas[i].Δy - temp1;dataGridView2.Rows[i].Cells[8].Value = Math.Round(rtdatas[i].CΔx, 4);dataGridView2.Rows[i].Cells[9].Value = Math.Round(rtdatas[i].CΔy, 4);i++;}}else if (form2.radioButton2.Checked){for (int i = 1; i < rdatas.Count-2; i++){double temp = 0, temp1 = 0;temp = tolerance.fx * rtdatas[i].L / tolerance.fl;temp1 = tolerance.fy * rtdatas[i].L / tolerance.fl;rtdatas[i].CΔx = rtdatas[i].Δx - temp;rtdatas[i].CΔy = rtdatas[i].Δy - temp1;dataGridView2.Rows[i].Cells[8].Value = Math.Round(rtdatas[i].CΔx, 4);dataGridView2.Rows[i].Cells[9].Value = Math.Round(rtdatas[i].CΔy, 4);}}}//計算展示CΔ和CΔypublic void ShowXandY(){if (form2.radioButton1.Checked){int i = 0;if (form2.comboBox1.Text == "兩個已知點"){i = 2;}else if (form2.comboBox1.Text == "一個已知點+方位角"){i = 1;}while (i < rdatas.Count){rtdatas[i].X = rtdatas[i - 1].X + rtdatas[i].CΔx;rtdatas[i].Y = rtdatas[i - 1].Y + rtdatas[i].CΔy;dataGridView2.Rows[i].Cells[10].Value = Math.Round(rtdatas[i].X, 4);dataGridView2.Rows[i].Cells[11].Value = Math.Round(rtdatas[i].Y, 4);i++;}}else if (form2.radioButton2.Checked){for (int i = 2; i < rtdatas.Count-1; i++){rtdatas[i].X = rtdatas[i - 1].X + rtdatas[i-1].CΔx;rtdatas[i].Y = rtdatas[i - 1].Y + rtdatas[i-1].CΔy;dataGridView2.Rows[i].Cells[10].Value = Math.Round(rtdatas[i].X, 4);dataGridView2.Rows[i].Cells[11].Value = Math.Round(rtdatas[i].Y, 4);}}}//計算展示X和Ypublic void ShowTolerance(){if (form2.radioButton1.Checked){int n = 0;if (form2.comboBox1.Text == "兩個已知點"){n = rdatas.Count - 2;}else if (form2.comboBox1.Text == "一個已知點+方位角"){n = rdatas.Count - 1;}richTextBox1.AppendText(Environment.NewLine);if (Math.Abs( Math.Round(tolerance.fβ, 0))> Math.Round(24 * Math.Sqrt(n), 4)){richTextBox1.Text = "角度閉合差為 fβ=" + Math.Round(tolerance.fβ, 0).ToString() + "″" + " fβ允=±24√n=±" + Math.Round(24 * Math.Sqrt(n), 4) + " 角度閉合差超限"+"\n";}else if (Math.Abs(Math.Round(tolerance.fβ, 0)) < Math.Round(24 * Math.Sqrt(n), 4)){richTextBox1.Text = "角度閉合差為 fβ=" + Math.Round(tolerance.fβ, 0).ToString() + "″" + " fβ允=±24√n=±" + Math.Round(24 * Math.Sqrt(n), 4) + "\n";}richTextBox1.Text += "fx=" + Math.Round(tolerance.fx, 4).ToString() + "m" + " fy=" + Math.Round(tolerance.fy, 4).ToString() + "m" + " fl=" + Math.Round(tolerance.fl, 4).ToString() + "m" + "\n";richTextBox1.Text += "導(dǎo)線全長閉合差fD =" + Math.Round(tolerance.fD, 4) + "m" + "\n";richTextBox1.Text += "導(dǎo)線全長相對閉合差容許值=" + (1.0 / 2000) + "\n";if (tolerance.K > 1.0 / 2000){richTextBox1.Text += "導(dǎo)線全長相對閉合差K =" + Math.Round(tolerance.K, 4) + ">" + (1.0 / 2000) + " " + "結(jié)果超限" + "\n";}else if (tolerance.K < 1.0 / 2000){richTextBox1.Text += "導(dǎo)線全長相對閉合差K =" + Math.Round(tolerance.K, 4) + "<" + (1.0 / 2000) + " " + "結(jié)果不超限" + "\n";}}else if (form2.radioButton2.Checked){int n = rtdatas.Count - 2;richTextBox1.AppendText(Environment.NewLine);richTextBox1.Text = "角度閉合差為 fβ=" + Math.Round(tolerance.fβ, 0).ToString() + "″" + " fβ允=±24√n=±" + Math.Round(24 * Math.Sqrt(n), 4) + "\n";richTextBox1.Text += "fx=" + Math.Round(tolerance.fx, 4).ToString() + "m" + " fy=" + Math.Round(tolerance.fy, 4).ToString() + "m" + " fl=" + Math.Round(tolerance.fl, 4).ToString() + "m" + "\n";richTextBox1.Text += "導(dǎo)線全長閉合差fD =" + Math.Round(tolerance.fD, 4) + "m" + "\n";richTextBox1.Text += "導(dǎo)線全長相對閉合差容許值=" + (1.0 / 2000) + "\n";if (tolerance.K > 1.0 / 2000){richTextBox1.Text += "導(dǎo)線全長相對閉合差K =" + Math.Round(tolerance.K, 4) + ">" + (1.0 / 2000) + " " + "結(jié)果超限" + "\n";}else if (tolerance.K < 1.0 / 2000){richTextBox1.Text += "導(dǎo)線全長相對閉合差K =" + Math.Round(tolerance.K, 4) + "<" + (1.0 / 2000) + " " + "結(jié)果不超限" + "\n";}}}//展示限差public void ClearData(){dataGridView2.Rows.Clear();rtdatas.Clear();D = 0; M = 0; S = 0;d1 = 0; m1 = 0; s1 = 0;tolerance.fx = 0;tolerance.fy = 0;tolerance.fl = 0;for (int i = 0; i < rdatas.Count; i++){ResultData rtdata = new ResultData();rtdata.DH = rdatas[i].DH;rtdatas.Add(rtdata);}}//計算前清理數(shù)據(jù)，避免數(shù)據(jù)重復(fù)計算public (double d, double m, double s) ReturnDMS(double DMS){double d, m, s;d = Math.Floor(DMS);m = Math.Floor(Math.Round(((DMS - d) * 100),6));s = Math.Round(((Math.Round(DMS * 100, 8) - Math.Floor(Math.Round(DMS * 100,8))) * 100), 8);return (d, m, s);}//格式轉(zhuǎn)換public double JudgementQuadrant(double Theta, int i, int j) //判斷方位角象限{double CA = 0;double m = rdatas[i].X - rdatas[j].X;double n = rdatas[i].Y - rdatas[j].Y;if (m > 0 && n >= 0){CA = Theta;}else if (m <= 0 && n > 0){CA = Math.PI - Theta;}else if (m < 0 && n <= 0){CA = Theta + Math.PI;}else if (m > 0 && n < 0){CA = 2 * Math.PI - Theta;}else if (m == 0 && n > 0){CA = Math.PI / 2;}else if (m == 0 && n < 0){CA = Math.PI * 1.5;}return CA;}public int[] AngularCarry(int d, int m, int s)//角度進位{if (s >= 60){s -= 60;m += 1;}else if (m >= 60){m -= 60;d += 1;}else if (s < 0){s += 60;m -= 1;}else if (m < 0){m += 60;d -= 1;}return new int[] { d, m, s };}public double VDistribution(double s, int i)//改正數(shù)分配{double temp2 = 0;//閉合if (form2.radioButton1.Checked){ArrayList arrayList = new ArrayList();int n = 0;if (form2.comboBox1.Text == "兩個已知點"){n = rdatas.Count - 2;for (int j = 1; j < rdatas.Count - 1; j++){arrayList.Add(rdatas[j].L);arrayList.Sort();}}else if (form2.comboBox1.Text == "一個已知點+方位角"){n = rdatas.Count - 1;for (int j = 1; j < rdatas.Count ; j++){arrayList.Add(rdatas[j].L);arrayList.Sort();}}int temp = Convert.ToInt32(Math.Round(tolerance.fβ, 5) % n);//取余int temp1 = Math.Abs(temp);double avernum = (tolerance.fβ - temp) / n;if (temp != 0){if (rdatas[i].L >= Convert.ToDouble(arrayList[n - temp1])){rtdatas[i].V = -(temp / temp1 + avernum);temp2 = s + rtdatas[i].V;}else{rtdatas[i].V = -(avernum);temp2 = s + rtdatas[i].V;}}else{rtdatas[i].V = -(tolerance.fβ / n);temp2 = s + rtdatas[i].V;}}//附和 else if (form2.radioButton2.Checked){int n=rtdatas.Count-2;int temp = Convert.ToInt32(Math.Round(tolerance.fβ, 5) % n);//取余int temp1 = Math.Abs(temp);double avernum = (tolerance.fβ - temp) / n;ArrayList arrayList = new ArrayList();for (int j = 1; j < rdatas.Count - 1; j++){arrayList.Add(rdatas[j].L);arrayList.Sort();}if (temp != 0){if (rdatas[i].L >= Convert.ToDouble(arrayList[n - temp1])){rtdatas[i].V = -(temp / temp1 + avernum);temp2 = s + rtdatas[i].V;}else{rtdatas[i].V = -(avernum);temp2 = s + rtdatas[i].V;}}else{rtdatas[i].V = -(tolerance.fβ / n);temp2 = s + rtdatas[i].V;}}return temp2;}public int[] toDMS(double CA){int d = (int)CA;int m = (int)((CA - d) * 60);int s = (int)Math.Round(((CA - d) * 3600 - m * 60), 0);int[] temp = AngularCarry(d, m, s);return new int[] { temp[0], temp[1], temp[2] };}//十進制度數(shù)轉(zhuǎn)度分秒public double toDu(double d1, double m1, double s1){double du = d1 + m1 / 60 + s1 / 3600;return du;}//度分秒轉(zhuǎn)十進制度數(shù)public void drawPoint(double X, double Y){g.FillRectangle(Brushes.Red, (float)Y, (float)X, 4, 4);}private void 繪圖ToolStripMenuItem_Click(object sender, EventArgs e){tabControl1.SelectedTab = tabControl1.TabPages[2];tabPage3.Show();g = tabPage3.CreateGraphics();midx = tabPage3.ClientRectangle.Height / 2;midy = tabPage3.ClientRectangle.Width / 2;g.DrawLine(Pens.Black, 0, midx, this.ClientRectangle.Width, midx);g.DrawLine(Pens.Black, midy, 0, midy, this.ClientRectangle.Height);curvelist = new List<Point>();int n = 20;for (int i = 0; i < rtdatas.Count; i++){if (i == 0){drawPoint((midx - rtdatas[i].X / n), (midy + rtdatas[i].Y / n));}else{drawPoint((midx - rtdatas[i].X / n), (midy + rtdatas[i].Y / n));g.DrawLine(Pens.Black, (float)(midy + rtdatas[i].Y / n), (float)(midx - rtdatas[i].X /n), (float)(midy + rtdatas[i - 1].Y / n), (float)(midx - rtdatas[i - 1].X / n));}}}private void 保存ToolStripMenuItem_Click(object sender, EventArgs e){//這里的代碼自己搜索一下，就是將datagridview 轉(zhuǎn)換成txt或者excel}} } //這里是創(chuàng)建的變量class RowData{public string DH;//點號public double DMS;//度分秒public double L;//距離public double X;public double Y;}class ResultData {public string DH;public double HA;//水平角public double V;//改正數(shù)public double CHA;//改正后水平角public double CA;//方位角public double L;public double Δx;public double Δy;public double CΔx;public double CΔy;public double X;public double Y;}class Tolerance //限差{public double fβ;public double fl;public double fx;public double fy;public double fD;public double K;}

Form2的代碼我就不復(fù)制了，如上圖一些設(shè)置的按鈕。

五、結(jié)語

不管是水準測量平差或者是導(dǎo)線測量，主要面臨的不是公式不會用（高斯正反算就是純粹的將公式套用編程），而是各種編程的基礎(chǔ)問題（數(shù)據(jù)導(dǎo)入，顯示，分隔，循環(huán)處理，debug，輸出等等），計算公式倒是好理解。其實最簡單的就是寫一個最死的程序，只處理一組數(shù)據(jù)或者單個數(shù)據(jù)，過后慢慢百度學習，久而久之就相對熟練了。
如果你有基礎(chǔ)，看源碼就可以了。如果你才學C#窗體，我覺得你下載這一個程序或者水準那個基本上可以滿足測繪程序編寫的要求了，不用每個問題都去搜索，這兩個程序基本都包括了。
我才學C#的時候，有很多自己也不會，也要看別人的，比如如何打開文件夾，讀取txt，用逗號或者空格分隔，換行讀取過后再顯示，跨窗體調(diào)用數(shù)據(jù)，保存數(shù)據(jù)，畫圖，方法的調(diào)用等等，遇到了各種報錯了就開始搜索再修改，過后再把它們?nèi)诤系阶约簩懙某绦颉?br /> 如果有一個現(xiàn)成的程序，會對你的程序?qū)W習有很大幫助。
我的其他文章：

利用Python編寫一個高斯正反算程序：https://blog.csdn.net/Zj1638/article/details/125740379

利用C#編寫一個高斯正反算程序：https://blog.csdn.net/Zj1638/article/details/125380593

利用C#編寫一個水準測量近似平差程序：https://blog.csdn.net/Zj1638/article/details/119303957

利用C#編寫一個GPS高程擬合（二次曲面擬合模型）程序：https://blog.csdn.net/Zj1638/article/details/125752243

總結(jié)

以上是生活随笔為你收集整理的利用C#编写一个附和闭合导线平差程序的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。