Leica DNA03水准测量数据处理程序

Leica DNA03水准测量数据处理程序

1 数据格式------------------------------------------------------------------------------------ 3

2 程序设计--------------------------------------------------------------------------------------- 3

3 源代码----------------------------------------------------------------------------- 4

4 编程总结--------------------------------------------------------------------------- 8

1数据格式

Leica数据格式为GSI(Geo Serial Interface)，根据其数据位是8或16位分为GSI8和GSI16，其他格式完全相同。数据以行为单位，每行包含若干个16位或24位的字段，每个字段以开头两位数字表示数据意义，称为WI(Word Index)，如41表示测量模式，11表示点号，32表示视距读数，331表示后视第一次中丝读数，332表示前视第一次中丝读数，336表示前视第二次中丝读数，335表示后视第二次中丝读数。

本次作业采用Leica DNA03数字水准仪原始数据。如下图：

第一行：观测模式为BFFB；

第二行：测段起点为BM3，起始高程为480.1230m；

第三行：测站后视点号BM3，视距读数28.58072m，第一次后视中丝读数0.69018m，中丝读数为2次读数的平均值，标准偏差为0.04mm；

第四行：测站前视点号Z1，视距读数28.61654m，第一次前视中丝读数1.69615m，中丝读数为2次读数的平均值，标准偏差为0.02mm；

第五行：测站前视点号Z1，视距读数28.62065m，第二次前视中丝读数1.69610m，中丝读数为2次读数的平均值，标准偏差为0.11mm；

第六行：测站后视点号BM3，视距读数28.58569m，第二次后视中丝读数0.69017m，中丝读数为2次读数的平均值，标准偏差为0.04mm。

2程序设计

本次数据处理程序的窗口设计非常简单，仅一个窗口，一个按钮。“读取数据”按钮读入原始水准仪数据，程序自动转换，以表格形式输出。

表头设计为“测站数”、“后视点号”、“前视点号”、“第一次后视中丝读数”、“第二次后视中丝读数”、“第一次前视中丝读数”、“第二次前视中丝读数”“测站两次前后高差较差”、“高差之差累计值”、“累计视距差”、“累计视距”、“前视点高程”共12项，表头之前输出“起始点点号”和“起始点高程”。

3源代码

using System;

using System.Collections.Generic; using https://www.360docs.net/doc/e911117357.html,ponentModel;

using System.Data;

using System.Drawing;

using System.Text;

using System.Windows.Forms;

using System.IO;

namespace Leica_DNA03水准数据处理

{

public partial class Form1 : Form {

public Form1()

{

InitializeComponent();

}

double[,] Rb = null;

double[,] Rf = null;

string[] sb = null; //后视点号

string[] sf = null; //前视点号

double[,] Rbs = null;//后视视距

double[,] Rfs = null;//前视视距

string st2,bmn1,bmn2;

double bm;

double n;

double[,] z1 = null;

string[] s=null;

public void Open_LeicaDNA03() //打开egm96文件的子函数

{

Rb = new double[100,2]; //后视中丝读数

Rf = new double[100,2];

sb = new string[100];

sf = new string[100];

Rbs = new double[100,2];

Rfs = new double[100,2];

z1 = new double[100,5];

string s1,s2,s3,s4,s5,s6,s7,s8,s9,s10,s11;

double bm;

s=new string[100];

openFileDialog1.Title = "选择文件";

openFileDialog1.InitialDirectory = "";

openFileDialog1.Filter = "文本文件(*txt)|*txt|数据文件(*dat)|*dat|数据文件

(*gsi)|*gsi";

if (openFileDialog1.ShowDialog() == DialogResult.OK)

{

System.IO.StreamReader sr = new

System.IO.StreamReader(openFileDialog1.FileName);

string st = sr.ReadLine();

string st2 = sr.ReadLine();

bmn1 = st2.Substring(12, 3); //起始点点号

bm = (double.Parse(st2.Substring(25, 6))) / 10000; //起始点高程

s1=bm.ToString("F4");

s[0] = "起始点点号：" + bmn1 + " " + "起始点高程:" + s1+" m";

int i=0;

textBox1.AppendText(s[0] + "\n");

s[1]="测站"+" "+"后视点" +" "+ "前视点"+" "+"第一次后视"+" "+"第二次后视"+" "+"第一次前视"+" "+"第二次前视"+" "+" "+"测站两次前后"+" "+"高差之差"+" "+" 累

计 "+" "+" 累计 "+" "+"前视点";

s[2]=" 数 "+" "+" 点号 " +" " +" 点号 "+" "+"中丝读数/m"+" "+"中丝读数

/m"+" "+"中丝读数/m"+" "+"中丝读数/m"+" "+" "+" 高差较差/m "+" "+"累计值/m"+" "+"视

距差/m"+" "+" 视距/m "+" "+"高程/m";

textBox1.AppendText("测站"+" "+"后视点" +" "+ "前视点"+" "+"第一次后视"+" "+"第二次后视"+" "+"第一次前视"+" "+"第二次前视"+" "+" "+"测站两次前后"+" "+"高差之差"+" "+" 累计 "+" "+" 累计 "+" "+"前视点"+"\n");

textBox1.AppendText(" 数 "+" "+" 点号 " +" " +" 点号 "+" "+"中丝读数/m"+" "+"中丝读数/m"+" "+"中丝读数/m"+" "+"中丝读数/m"+" "+" "+" 高差较差/m "+" "+"累计值

/m"+" "+"视距差/m"+" "+" 视距/m "+" "+"高程/m"+"\n");

st2 = sr.ReadLine();

do

{

for (int j = 0; j < 5; j++)

{

if (j == 0)

{

sb[i] = st2.Substring(12, 3); //后视点号

Rbs[i, 0] = (double.Parse(st2.Substring(25, 6))) / 10000.0; //后视视距

Rb[i, 0] = (double.Parse(st2.Substring(38, 9))) / 10000.0; //后视中丝读数

}

if (j == 1)

{

st2 = sr.ReadLine();

sf[i] = st2.Substring(12, 3);

bmn2 = st2.Substring(12, 3);

Rfs[i, 0] = (double.Parse(st2.Substring(25, 6))) / 10000;

Rf[i, 0] = (double.Parse(st2.Substring(38, 9))) / 10000;

}

if (j == 2)

{

st2 = sr.ReadLine();

Rfs[i, 1] = (double.Parse(st2.Substring(25, 6))) / 10000; //前视二次视距

Rf[i, 1] = (double.Parse(st2.Substring(38, 9))) / 10000; //前视二次中丝读数

}

if (j == 3)

{

st2 = sr.ReadLine();

Rbs[i, 1] = (double.Parse(st2.Substring(25, 6))) / 10000;

Rb[i, 1] = (double.Parse(st2.Substring(38, 9))) / 10000;

}

if (j == 4) //测站总结

{

st2 = sr.ReadLine();

z1[i, 0] = (double.Parse(st2.Substring(25, 6))) / 10000; //测站两次高差之差

z1[i, 1] = (double.Parse(st2.Substring(38, 9))) / 10000;//高差之差累计值

z1[i, 2] = (double.Parse(st2.Substring(54, 9))) / 10000;//测站累计视距差

z1[i, 3] = (double.Parse(st2.Substring(70, 9))) / 10000;//测站累计视距

z1[i, 4] = (double.Parse(st2.Substring(86, 9))) / 10000;//前视点高程

}

}

n = i + 1;

s2 = Rb[i, 0].ToString("F4");

s3 = Rb[i, 1].ToString("F4");

s4 = Rf[i, 0].ToString("F4");

s5 = Rf[i, 1].ToString("F4");

s6 = z1[i, 0].ToString("F4");

s7 = z1[i, 1].ToString("F4");//高差之差累计值

s8 = z1[i, 2].ToString("F4");

s9 = z1[i, 3].ToString("F4");

s10 = z1[i, 4].ToString("F4");

s11 = n.ToString("F0");

s[i+3] = s11.PadLeft(3, ' ') + sb[i].PadLeft(8, ' ') + sf[i].PadLeft(8, ' ') + s2.PadLeft(12, ' ') + s3.PadLeft(12, ' ') + s4.PadLeft(12, ' ') + s5.PadLeft(10, ' ') + s6.PadLeft(16, ' ')

+ s7.PadLeft(12, ' ') + s8.PadLeft(12, ' ') + s9.PadLeft(14, ' ') + s10.PadLeft(11, ' ');

textBox1.AppendText(s[i+3] + "\n");

i++;

st2 = sr.ReadLine();

}

while (st2!=null);

sr.Close();

}

}

private void button1_Click(object sender, EventArgs e)

{

Open_LeicaDNA03();

int i=0;

System.IO.StreamWriter sw = new StreamWriter("从"+bmn1+"到"+bmn2+"水准数据处

理.txt");

do

{

sw.Write(s[i]);

sw.WriteLine();

i++;

}

while (s[i] != null);

sw.Close();

}

private void Form1_Load(object sender, EventArgs e)

{

}

}

}

4编程总结

本次作业的原始数据中出现的WI如571、572、573、574一开始并不了解其意思，但在查阅Leica公司说明书GSI ONLINE for Leica TPS and DNA 后弄清楚了。

整个编程中，并无大量计算，主要在于数据的读入与重新输出。重点就是原始数据与输出结果每一项的对应要正确无误，代码编写时数据的读入、创建数组、调用数据等都是重点，由于数据繁多，并且格式特殊，在读入数据时每一行每个字符位置都要数清楚，应细致反复检查。