基于VC++GPS高程曲面拟合方法实现.docVIP

基于VC++GPS高程曲面拟合方法实现摘要：GPS定位测量技术可以得到高精度的平面位置和大地高差,但在实际应用中,地面点的高程常采用的是正常高高程系统,因此,需将GPS大地高转换成正常高。利用VC++6.0实现GPS高程多项式曲面拟合程序的实现。 关键词：GPS高程，高程拟合，VC++实现 一 前言 GPS以其全天候、高精度、自动化、高效率、多功能的显著特点已经广泛的应用于大地测量、工程测量、航空摄影测量等,给测绘领域带来一场深刻的革命。大量的实践数据表明，GPS测量的平面坐标精度是可靠的，能达到工程测量的要求，其平面相对定位的精度能够达到(0.1~1)10-6D mm甚至更高(D是基线长度)。对于G P S高程,为了满足工程的需要,须将GPS中的大地高转换为工程中直接应用的正常高,才能真正实现GPS的三维测量的目的。目前,多项式曲面拟合GPS水准高程,以其计算简单、容易编程实现,已经广泛应用于工程实践。 二GPS高程多项式曲面拟合原理 当测区GPS点布设成一定区域面时,用几何水准联测部分GPS点的正常高,用数值拟合的方法求出测区似大地水准面,计算出未联测几何水准GPS点的高程异常,即:根据测区已知点的平面坐标x,y(或大地坐标B,L)和ζ值,用数值拟合原理拟合出测区似大地水准面,再内插出待求点的ζ值,从而求出待求点的正常高。 则其高程异常一般表示为ξ=f（x，y），根据测区的实际情况，f（x，y）对应有平面拟合、二次曲面拟合、三次曲面拟合3种常用数学模型。 一般来说，根据测区实际情况，确定适宜的高程异常拟合数学模型，并根据已知联测控制点的高程异常，用最小二乘拟合法确定拟合数学模型的系数。 假设某GPS测区内有i个高程联测点，异常拟合数学模型采用多项式曲面拟合： f(x,y)= …….（1） 式中-------多项式系数，i=0,1,…,5；（x,y）-------GPS高程控制点的平面坐标。 根据n个GPS点的可以确定多项式(1)中的系数。当n≥3且n KNOWDATE; 定义待求数据类CDATEGPS，并定义数组储存数据。 class CDATEGPS : public CObject { public: BOOL ReaddategpsCASS(CString str);//读取数据 BOOL ReadDateGps(CString StrLine); double x,y,h,v,Hv; CString pro; CString Pointname ; CDATEGPS(); virtual ~CDATEGPS(); }; typedef CTypedPtrArray DATEGPS; 定义矩阵类CMatrix，用于矩阵计算。 class CMatrix { CMatrix(); CMatrix(int nRows, int nCols); // 指定行列构造函数 BOOLInit(int nRows, int nCols); // 初始化矩阵 BOOL SetElement(int nRow, int nCol, double value); // 设置指定元素的值 double GetElement(int nRow, int nCol) const; // 获取指定元素的值 voidSetData(double value[]); // 设置矩阵的值 int GetNumColumns() const; // 获取矩阵的列数 int GetNumRows() const; // 获取矩阵的行数 // 数学操作 CMatrix operator=(const CMatrix other); BOOL operator==(const CMatrix other) const; BOOL operator!=(const CMatrix other) const; CMatrix operator+(const CMatrix other) const; CMatrix operator-(const CMatrix other) const; CMatrix operator*(double value) const; CMatrix operator*(const CMatrix other) const; // 复矩阵乘法 BOOL CMul(const CMatrix AR, const CMatrix AI, const CMatrix BR