DEM内插程序

摄影测量学 数字高程模型 DEM内插，运用C#编写的，界面整洁，使用方便

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

namespace DEM内插程序
{
    public partial class DEMForm : Form
    {
        public DEMForm（）
        {
            InitializeComponent（）;
        }

        #region 定义变量
        double[] XYZ;//存储原始数据
        double[] X Y Z;//存储已知坐标

        //定义需要的矩阵及变量
        double Xp Yp;//待求点的XYZ坐标
        double[] _x = new double[10] _y = new double[10];//用于存放已知点与待定点的差值
        double[] A = new double[10 6];//用于存放系数
        double[] P = new double[10 10];//用于存放权值
        double[] AT = new double[6 10];//A系数阵转置
        double[] ATP = new double[6 10];//A的转置与P的乘积
        double[] ATPA = new double[6 6];//ATP与A的乘积
        double[] temp = new double[6 12];//临时矩阵
        double[] ATPAR = new double[6 6];//ATP与A的乘积的逆矩阵
        double[] ATPZ = new double[6];//ATP与Z的乘积
        double[] xx = new double[6];//X向量

        #endregion

        #region 定义函数
        public void QinNi（double[] c）
        {
            int i j h k n = 6;
            double p;
            double[] q = new double[6 12];
            for （i = 0; i < n; i++）//构造高斯矩阵
                for （j = 0; j < n; j++）
                    q[i j] = c[i j];
            for （i = 0; i < n; i++）
                for （j = n; j < 12; j++）
                {
                    if （i + n == j）
                        q[i j] = 1;
                    else
                        q[i j] = 0;
                }
            for （h = k = 0; k < n - 1; k++ h++）//消去对角线以下的数据
                for （i = k + 1; i < n; i++）
                {
                    if （q[i h] == 0）
                        continue;
                    p = q[k h] / q[i h];
                    for （j = 0; j < 12; j++）
                    {
                        q[i j] *= p;
                        q[i j] -= q[k j];
                    }
                }
            for （h = k = n - 1; k > 0; k-- h--） // 消去对角线以上的数据
                for （i = k - 1; i >= 0; i--）
                {
                    if （q[i h] == 0）
                        continue;
                    p = q[k h] / q[i h];
                    for （j = 0; j < 12; j++）
                    {
                        q[i j] *= p;
                        q[i j] -= q[k j];
                    }
                }
            for （i = 0; i < n; i++）//将对角线上数据化为1
            {
                p = 1.0 / q[i i];
                for （j = 0; j < 12; j++）
                    q[i j] *= p;
            }
            for （i = 0; i < n; i++） //提取逆矩阵
                for （j = 0; j < n; j++）
                    c[i j] = q[i j + n];
        }
        #endregion

        #region 导入数据
        private void btnInput_Click（object sender EventArgs e）
        {

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

     文件     186368  2011-11-18 15:13  DEM内插程序\DEM内插程序\bin\Debug\DEM内插程序.exe

     文件      36352  2011-11-18 15:13  DEM内插程序\DEM内插程序\bin\Debug\DEM内插程序.pdb

     文件      11600  2011-11-18 15:13  DEM内插程序\DEM内插程序\bin\Debug\DEM内插程序.vshost.exe

     文件        490  2010-03-17 22:39  DEM内插程序\DEM内插程序\bin\Debug\DEM内插程序.vshost.exe.manifest

     文件      10104  2011-11-15 10:27  DEM内插程序\DEM内插程序\DEMForm.cs

     文件       9596  2011-11-15 10:27  DEM内插程序\DEM内插程序\DEMForm.Designer.cs

     文件     260612  2011-11-15 10:27  DEM内插程序\DEM内插程序\DEMForm.resx

     文件       3699  2011-11-14 09:51  DEM内插程序\DEM内插程序\DEM内插程序.csproj

     文件       2088  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\DEM内插程序.csproj.FileListAbsolute.txt

     文件     170093  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\DEM内插程序.DEMForm.resources

     文件     186368  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\DEM内插程序.exe

     文件      36352  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\DEM内插程序.pdb

     文件        180  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\DEM内插程序.Properties.Resources.resources

     文件       4440  2011-11-14 14:03  DEM内插程序\DEM内插程序\obj\x86\Debug\DesignTimeResolveAssemblyReferences.cache

     文件       6289  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\DesignTimeResolveAssemblyReferencesInput.cache

     文件       1282  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\GenerateResource.read.1.tlog

     文件       1426  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug\GenerateResource.write.1.tlog

     文件       7897  2011-11-14 14:08  DEM内插程序\DEM内插程序\obj\x86\Debug\ResolveAssemblyReference.cache

     文件        498  2011-11-14 09:38  DEM内插程序\DEM内插程序\Program.cs

     文件       1362  2011-11-14 09:33  DEM内插程序\DEM内插程序\Properties\AssemblyInfo.cs

     文件       2880  2011-11-14 09:33  DEM内插程序\DEM内插程序\Properties\Resources.Designer.cs

     文件       5612  2011-11-14 09:33  DEM内插程序\DEM内插程序\Properties\Resources.resx

     文件       1100  2011-11-14 09:33  DEM内插程序\DEM内插程序\Properties\Settings.Designer.cs

     文件        249  2011-11-14 09:33  DEM内插程序\DEM内插程序\Properties\Settings.settings

     文件        887  2011-11-14 09:33  DEM内插程序\DEM内插程序.sln

    ..A..H.     18944  2011-11-18 15:14  DEM内插程序\DEM内插程序.suo

     文件      10240  2011-11-14 15:20  DEM内插程序\DEM数据.xls

     目录          0  2011-11-14 09:33  DEM内插程序\DEM内插程序\obj\x86\Debug\TempPE

     目录          0  2011-11-18 15:13  DEM内插程序\DEM内插程序\obj\x86\Debug

     目录          0  2011-11-18 15:05  DEM内插程序\DEM内插程序\bin\Debug

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