单像空间后方交会摄影测量

是用C#编写的摄影测量【单像空间后方交会】程序，压缩包内包含数据文件（按此格式），还有源程序。

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;


namespace 单像空间后方交会_20111003DataGridView_
{
    class BackIntersect
    {
        /****************************************************
         * 方法   提取但对像点误差方程的系数阵和常数阵
         * 形参   地面点和像点坐标、外方位元素、摄影主距
         * 返回值 误差方程的系数阵和常数阵
         ****************************************************/

        public double[] Want_AL（double[] ConPoint double[] Position double f）
        {

            /***要用到数据的申明***/

            //旋转矩阵R（采用awk转角系统）
            double[] Matrix_R = new double[3 3];
            double a1 a2 a3 b1 b2 b3 c1 c2 c3;

            //误差方程系数阵
            double[] Matrix_A = new double[2 6];
            double a11 a12 a13 a14 a15 a16;
            double a21 a22 a23 a24 a25 a26;

            //误差方程式的常数项
            double Lx;
            double Ly;

            //误差方程系数阵与常数阵的【组合阵】
            double[] Matrix_AL = new double[2 7];

            //控制点像点坐标的近似值
            double xo;
            double yo;

            //【控制点像点（xy-f）】在像空间坐标系中坐标
            //【控制点物点（XYZ）】在摄影测量测量坐标系中坐标
            //【六个外方位元素】假设值
            double x = ConPoint[0];
            double y = ConPoint[1];
            double X = ConPoint[2];
            double Y = ConPoint[3];
            double Z = ConPoint[4];

            double Xso = Position[0];
            double Yso = Position[1];
            double Zso = Position[2];
            double ao = Position[3];
            double wo = Position[4];
            double ko = Position[5];

            /***计算模块***/
            //计算旋转矩阵（采用awk转角系统）
            Matrix_R[0 0] = a1 = Math.Cos（ao） * Math.Cos（ko） - Math.Sin（ao） * Math.Sin（wo） * Math.Sin（ko）;
            Matrix_R[0 1] = a2 = -Math.Cos（ao） * Math.Sin（ko） - Math.Sin（ao） * Math.Sin（wo） * Math.Cos（ko）;
            Matrix_R[0 2] = a3 = -Math.Sin（ao） * Math.Cos（wo）;
            Matrix_R[1 0] = b1 = Math.Cos（wo） * Math.Sin（ko）;
            Matrix_R[1 1] = b2 = Math.Cos（wo） * Math.Cos（ko）;
            Matrix_R[1 2] = b3 = -Math.Sin（wo）;
            Matrix_R[2 0] = c1 = Math.Sin（ao） * Math.Cos（ko） + Math.Cos（ao） * Math.Sin（wo） * Math.Sin（ko）;
            Matrix_R[2 1] = c2 = -Math.Sin（ao） * Math.Sin（ko） + Math.Cos（ao） * Math.Sin（wo） * Math.Cos（ko）;
            Matrix_R[2 2] = c3 = Math.Cos（ao） * Math.Cos（wo）;

            //计算控制点像点坐标的近似值
            //共线条件方程
            double temp1 = a1 * （X - Xso） + b1 * （Y - Yso） + c1 * （Z - Zso）;
            double temp2 = a2 * （X - Xso） + b2 * （Y - Yso） + c2 * （Z - Zso）;
            double temp3 = a3 * （X - Xso） + b3 * （Y - Yso） + c3 * （Z - Zso）;
            xo = -f * temp1 / temp3;
            yo = -f * temp2 / temp3;

            //误差方程常数项计算
            Lx = x - xo;
            Ly = y - yo;

            //误差方程系数项计算
            //公式5-8
            Matrix_A[0 0] = a11 = （a1 * f + a3 * x） / temp3;
            Matrix_

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-07-01 16:22  后方交会\
     文件         733  2011-10-21 20:36  后方交会\dat.txt
     文件         257  2011-10-21 17:06  后方交会\单像空间后方交会数据 - 副本.txt
     目录           0  2012-07-01 16:21  后方交会\后方交会\
     目录           0  2012-07-01 16:21  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\
     文件        8620  2011-10-21 20:28  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\BackIntersect.cs
     文件       16126  2011-10-21 17:18  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Form1.Designer.cs
     文件        7852  2011-10-21 20:34  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Form1.cs
     文件        8819  2011-10-21 17:18  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Form1.resx
     文件        6584  2011-10-04 16:19  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Matrix.cs
     文件         527  2011-10-03 22:04  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Program.cs
     目录           0  2012-07-01 16:21  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Properties\
     文件        1432  2011-10-03 22:04  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Properties\AssemblyInfo.cs
     文件        2942  2011-10-03 22:04  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Properties\Resources.Designer.cs
     文件        5612  2011-10-03 22:04  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Properties\Resources.resx
     文件        1131  2011-10-03 22:04  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Properties\Settings.Designer.cs
     文件         249  2011-10-03 22:04  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\Properties\Settings.settings
     目录           0  2012-07-01 16:21  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\bin\
     目录           0  2012-07-01 16:21  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\bin\Debug\
     文件       20992  2011-10-21 20:28  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\bin\Debug\单像空间后方交会（20111003DataGridView）.exe
     文件       50688  2011-10-21 20:28  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\bin\Debug\单像空间后方交会（20111003DataGridView）.pdb
     文件       11600  2011-10-21 20:20  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\bin\Debug\单像空间后方交会（20111003DataGridView）.vshost.exe
     目录           0  2012-07-01 16:22  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\bin\Release\
     目录           0  2012-07-01 16:21  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\
     目录           0  2012-07-01 16:21  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\x86\
     目录           0  2012-07-01 16:21  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\x86\Debug\
     文件        4440  2011-10-03 23:40  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\x86\Debug\DesignTimeResolveAssemblyReferences.cache
     文件        6478  2011-10-21 20:28  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\x86\Debug\DesignTimeResolveAssemblyReferencesInput.cache
     文件         428  2011-10-21 19:57  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\x86\Debug\GenerateResource.read.1.tlog
     文件        1046  2011-10-21 19:57  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\x86\Debug\GenerateResource.write.1.tlog
     目录           0  2012-07-01 16:22  后方交会\后方交会\单像空间后方交会（20111003DataGridView）\obj\x86\Debug\TempPE\
............此处省略8个文件信息