Пример #1
0
        private double delta;    // 仪器常数
        //private double oriPointHorCoor;   // 定向点横坐标

        /// <summary>
        /// 构造器
        /// </summary>
        public ObserveData(string[] rowText)
        {
            //坐标边方位角
            coorAzi = AngleTransUtil.DMSToDegree(Convert.ToDouble(rowText[1]));
            //陀螺边方位角
            gyroAzi = AngleTransUtil.DMSToDegree(Convert.ToDouble(rowText[2]));
            // 定向点横坐标
            //OriPointHorCoor = Convert.ToDouble(rowText[3]);
        }
Пример #2
0
        /// <summary>
        /// 打开并读取数据文件
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void OpenDataFileToolStripMenuItem_Click(object sender, EventArgs e)
        {
            MyOpenFileDialog.InitialDirectory = Application.StartupPath;
            MyOpenFileDialog.Filter           = "txt file(*.txt)|*.txt"; // 文件类型过滤器设置为txt文件
            bool estiThroErrData = false;

            if (MyOpenFileDialog.ShowDialog() == DialogResult.OK)
            {
                string[] allText = File.ReadAllLines(MyOpenFileDialog.FileName, Encoding.UTF8); // 读取txt文件中全部的数据
                for (int i = 0; i < allText.Length - 1; i++)
                {
                    string[] rowText = allText[i].Replace(" ", "").Replace(",,", ",").Split(',');   // 每一行的数据
                    if (rowText[0] == "")
                    {
                        estiThroErrData = true;
                        i       = i + 1;
                        rowText = allText[i].Replace(" ", "").Replace(",,", ",").Split(',');        // 每一行的数据
                        ml      = Convert.ToDouble(rowText[0]);
                        mb      = AngleTransUtil.DegreeToRand(Convert.ToDouble(rowText[1]) / 3600); //将mb转成弧度
                        //ma = DegreeToRand(DMSToDegree(Convert.ToDouble(rowText[2])));
                        continue;
                    }
                    if (estiThroErrData)
                    {   // 贯通数据
                        bool          isDirEdge  = Convert.ToInt16(rowText[3]) == 1;
                        int           surveyNum  = Convert.ToInt16(rowText[4]);
                        TraversePoint startPoint = new TraversePoint(rowText);
                        startPoint.Pointindex = PointNum;
                        rowText = allText[i + 1].Replace(" ", "").Replace(",,", ",").Split(',');   // 每一行的数据
                        TraversePoint endPoint = new TraversePoint(rowText);
                        endPoint.Pointindex = PointNum + 1;
                        TraverseEdge traverEdge = new TraverseEdge(startPoint, endPoint, isDirEdge, surveyNum);
                        traverseEdge.Add(traverEdge);
                        PointNum++;
                    }
                    else
                    {
                        // 陀螺经纬仪观测数据
                        ObserveData obData = new ObserveData(rowText); // 观测数据
                        if (rowText[0] == "DM")
                        {                                              //地面
                            groundData.Add(obData);
                        }
                        if (rowText[0] == "JX")
                        {
                            //井下
                            downholeData.Add(obData);
                        }
                    }
                }
            }
            isReadData = true;
            InitState();
            // 更新表格数据
            InsertObDataIntoTable();
        }
Пример #3
0
        /// <summary>
        /// 陀螺经纬仪定向精度评定
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void CalDownholeCoorAziToolStripMenuItem_Click(object sender, EventArgs e)
        {
            if (isCalma)
            {
                return;
            }

            //计算仪器常数中误差
            mv = 0;  //仪器常数中误差清零
            double vv       = 0;
            double deltaSum = 0;

            for (int i = 0; i < groundData.Count; i++)
            {
                groundData[i].Delta = groundData[i].CoorAzi - groundData[i].GyroAzi;
                DataTableGridView.Rows[i].Cells[3].Value = AngleTransUtil.DMSToStr(AngleTransUtil.DegreeToDMS(groundData[i].Delta));
                deltaSum = deltaSum + groundData[i].Delta;
            }
            double deltaAverage = deltaSum / groundData.Count;  //仪器常数平均值

            for (int i = 0; i < groundData.Count; i++)
            {
                DataTableGridView.Rows[i].Cells[4].Value = Math.Round(-(groundData[i].CoorAzi - groundData[i].GyroAzi - deltaAverage) * 3600, 2) + "″";
                vv = vv + (groundData[i].CoorAzi - groundData[i].GyroAzi - deltaAverage) * (groundData[i].CoorAzi - groundData[i].GyroAzi - deltaAverage) * 3600 * 3600; //单位:秒
            }
            mv = Math.Sqrt(vv / (groundData.Count - 1));                                                                                                                 //仪器常数一次测定中误差

            ma = Math.Sqrt(mv * mv / groundData.Count + mv * mv / 2);                                                                                                    //井下一次测定中误差

            //计算各个陀螺边的误差
            for (int i = 0; i < gyroEdge.Count; i++)
            {
                gyroEdge[i].Ma = Math.Sqrt(mv * mv / groundData.Count + mv * mv / gyroEdge[i].SurveyNum);
            }

            MessageBox.Show("陀螺经纬仪定向精度评定完成!", "提示",
                            MessageBoxButtons.OK, MessageBoxIcon.Information, MessageBoxDefaultButton.Button1);
            isCalma = true;
            OutMaReport();
        }
Пример #4
0
        /// <summary>
        /// 将观测的数据插入到表格中
        /// </summary>
        private void InsertObDataIntoTable()
        {
            int dataCount = groundData.Count;

            for (int i = 0; i < dataCount; i++)
            {
                DataTableGridView.Rows.Add(1);
                DataTableGridView.Rows[i].Cells[0].Value = i + 1;
                DataTableGridView.Rows[i].Cells[1].Value = AngleTransUtil.DMSToStr(groundData[i].CoorAzi);
                DataTableGridView.Rows[i].Cells[2].Value = AngleTransUtil.DMSToStr(groundData[i].GyroAzi);
            }
            int traverseEdgeCount = traverseEdge.Count;

            for (int i = 0; i < traverseEdgeCount; i++)
            {
                dataGridView1.Rows.Add(1);
                dataGridView1.Rows[i].Cells[0].Value = i + 1;
                dataGridView1.Rows[i].Cells[1].Value = traverseEdge[i].StartPoint.PointName;
                dataGridView1.Rows[i].Cells[2].Value = traverseEdge[i].EndPoint.PointName;
                dataGridView1.Rows[i].Cells[3].Value = traverseEdge[i].IsGyroDirEdge;
                dataGridView1.Rows[i].Cells[4].Value = traverseEdge[i].SurveyNum;
            }
        }