C#单像空间后方交会

摄影测量学的单像空间后方交会程序，用c＃实现，望共同学习

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Text;
using System.Windows.Forms;
using System.IO;
using System.Data.OleDb;


namespace WindowsApplication1
{
    public partial class Form1 : Form
    {

        #region//全局变量

        int i;
        int K;
        double ψ ω κ m Xs0 Ys0 Zs0 f S1 = 0.0 S2 = 0.0 S3 = 0.0;
        double[] x = new double[4];//{ -86.15 -53.40 -14.78 10.46 };
        double[] y = new double[4];// { -68.99 82.21 -76.63 64.43 };
        double[] X = new double[4];//{ 36589.41 37631.08 39100.97 40426.54 };
        double[] Y = new double[4];// { 25273.32 31324.51 24934.98 30319.81 };
        double[] Z = new double[4];// { 2195.17 728.69 2385.50 757.31 };


        double[] H = new double[4];
        double[] I = new double[4];
      
        double[] R = new double[3 3];  

        double[] x_0 = new double[4];  //xy 的近似值
        double[] y_0 = new double[4];

        double[] L = new double[6 1];  //系数矩阵A ※ L
        double[] A = new double[6 6];
         
           //矩阵计算中间矩阵
        double[] AT = new double[6 6];
        double[] AT_A = new double[6 6];
        double[] ATA_ni_AT = new double[6 6];
        double[] E = new double[6 6];
        double[] det = new double[6 1];

        bool Iscontiue = true; //当次数超过500时 用于判断是否继续的bool值

        #endregion

        public Form1（）
        {
            InitializeComponent（）;
        }

        private void button1_Click（object sender EventArgs e）
        {
            try
            {
                for （int j = 0; j < 5; j++）
            {
                for （int i = 0; i < 3; i++）
                {
                    X[i] = Convert.ToDouble（dataGridView1[0 i].Value）;
                    Y[i] = Convert.ToDouble（dataGridView1[1 i].Value）;
                    Z[i] = Convert.ToDouble（dataGridView1[2 i].Value）;
                    x[i] = （Convert.ToDouble（dataGridView1[3 i].Value））;
                    y[i] = （Convert.ToDouble（dataGridView1[4 i].Value））;

                }
            }

            #region// 初始化

            

                for （i = 0; i < 3; i++）
                {
                    S1 += X[i];
                    S2 += Y[i];
                    S3 += Z[i];
                }
                    //确定外方位元素的初始值
                ψ = 0;
                ω = 0;
                κ = 0;
                if （tb_m.Text == ““ || tb_f.Text == ““）
                {
                    MessageBox.Show（“情输入m和f“）;
                    return;
                }
                m = Convert.ToDouble（tb_m.Text）;
               // m = 50000;
                f = Convert.ToDouble（tb_f.Text）;
               // f = 0.15324;
                Xs0 = S1 / 3.0;
                Ys0 = S2 / 3.0;
                Zs0 = m * f + S3 / 3.0;

                    //像片坐标的

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     406528  2008-12-05 22:36  sheying\后方交会.msi

     文件     118784  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\bin\Debug\WindowsApplication1.exe

     文件      52736  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\bin\Debug\WindowsApplication1.pdb

     文件       5632  2005-12-08 14:51  sheying\后方交会源码\WindowsApplication1\bin\Debug\WindowsApplication1.vshost.exe

     文件      18337  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\Form1.cs

     文件      32150  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\Form1.Designer.cs

     文件      19155  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\Form1.resx

     文件        342  2008-12-05 17:40  sheying\后方交会源码\WindowsApplication1\help.cs

     文件       3142  2008-12-05 17:40  sheying\后方交会源码\WindowsApplication1\help.Designer.cs

     文件     100665  2008-12-05 17:40  sheying\后方交会源码\WindowsApplication1\help.resx

     文件       5930  2008-12-04 22:44  sheying\后方交会源码\WindowsApplication1\Matrix.cs

     文件        899  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\obj\Debug\WindowsApplication1.csproj.GenerateResource.Cache

     文件     118784  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\obj\Debug\WindowsApplication1.exe

     文件       9384  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\obj\Debug\WindowsApplication1.Form1.resources

     文件      63519  2008-12-05 17:40  sheying\后方交会源码\WindowsApplication1\obj\Debug\WindowsApplication1.help.resources

     文件      52736  2008-12-05 22:16  sheying\后方交会源码\WindowsApplication1\obj\Debug\WindowsApplication1.pdb

     文件        180  2008-12-04 22:24  sheying\后方交会源码\WindowsApplication1\obj\Debug\WindowsApplication1.Properties.Resources.resources

     文件       1226  2008-12-08 13:14  sheying\后方交会源码\WindowsApplication1\obj\WindowsApplication1.csproj.FileListAbsolute.txt

     文件        478  2008-12-04 22:23  sheying\后方交会源码\WindowsApplication1\Program.cs

     文件       1192  2008-12-04 22:23  sheying\后方交会源码\WindowsApplication1\Properties\AssemblyInfo.cs

     文件       2894  2008-12-04 22:23  sheying\后方交会源码\WindowsApplication1\Properties\Resources.Designer.cs

     文件       5612  2008-12-04 22:23  sheying\后方交会源码\WindowsApplication1\Properties\Resources.resx

     文件       1104  2008-12-04 22:23  sheying\后方交会源码\WindowsApplication1\Properties\Settings.Designer.cs

     文件        249  2008-12-04 22:23  sheying\后方交会源码\WindowsApplication1\Properties\Settings.settings

     文件       3616  2008-12-05 08:43  sheying\后方交会源码\WindowsApplication1\WindowsApplication1.csproj

     文件        946  2008-12-08 13:14  sheying\后方交会源码\WindowsApplication1.sln

    ..A..H.     19456  2008-12-08 13:15  sheying\后方交会源码\WindowsApplication1.suo

     文件       7680  2008-12-05 18:18  sheying\测试数据.xls

     目录          0  2008-12-05 08:10  sheying\后方交会源码\WindowsApplication1\obj\Debug\Refactor

     目录          0  2008-12-04 22:23  sheying\后方交会源码\WindowsApplication1\obj\Debug\TempPE

............此处省略11个文件信息