本发明专利技术涉及一种基于总体最小二乘的空间数据精度提高方法，1)计算迭代初值；2)迭代计算和的值；3)重复第二步直到满足条件其中ε是给定的阈值；4)满足迭代条件后，计算参数平差值及观测值改正数；5)通过单位权方差与检核点点位中误差来检验进度。与现有技术相比，本发明专利技术具有同时考虑实测点和图上点的误差，并能进一步提高空间数据的精度等优点。
【技术实现步骤摘要】
本专利技术涉及一种精度提高方法，尤其是涉及。
技术介绍
空间数据几何纠正中，需要首先利用控制点来计算转换模型参数，然后根据转换参数实现地籍图上所有界址点的坐标转换，实现地籍图上空间数据的更新和精度提高。其中控制点由GPS实测的地面控制点坐标(X，Y)和对应地籍图上的同名点图上坐标(x，y)构成。常用转换模型包括平移变换、平移加缩放变换、相似变换、仿射变换和多项式变换模型等。平移变换模型为权利要求1. ，该方法同时考虑实测点和图上点的误差，考虑系数阵误差的总体最小二乘平差模型为2.根据权利要求1所述的，其特征在于，所述的步骤5)中的单位权方差如下3.根据权利要求1所述的一种基于总体最/J· 于，所述的步骤5)中的单位权方差如下 由权利要求1推导过程可知4.根据权利要求1所述的，其特征在于，所述的步骤5)中的检核点点位中误差如下把地面控制点划分成两部分，一部分作为参与平差计算转换参数的控制点，一部分作为未参与平差的检核点；令图上点纠正后坐标为(χ/，y/ )，设纠正后坐标与检核点实测坐标偏差为(AXi， 由此可以得到X，Y方向的坐标中误差以及点位中误差为全文摘要本专利技术涉及，1)计算迭代初值；2)迭代计算和的值；3)重复第二步直到满足条件其中ε是给定的阈值；4)满足迭代条件后，计算参数平差值及观测值改正数；5)通过单位权方差与检核点点位中误差来检验进度。与现有技术相比，本专利技术具有同时考虑实测点和图上点的误差，并能进一步提高空间数据的精度等优点。文档编号G06F17/00GK102222058SQ201010148258公开日2011年10月19日 申请日期2010年4月15日 优先权日2010年4月15日专利技术者梁丹, 童小华, 金雁敏 申请人:同济大学本文档来自技高网...
【技术保护点】
１．一种基于总体最小二乘的空间数据精度提高方法，该方法同时考虑实测点和图上点的误差，考虑系数阵误差的总体最小二乘平差模型为：（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）ｌ（／ｍｉ）?（ｍｏ）＋（／ｍｏ）?（ｍｉ）ｖ（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）Ａ（／ｍｉ）?（ｍｏ）＋（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｅ（／ｍｉ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏｖｅｒ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ）＾（／ｍｏ）?（／ｍｏｖｅｒ）?（／ｍｒｏｗ）?（／ｍａｔｈ）其中，ｌ为ｎ×１的观测值向量，ｖ为ｎ×１观测值改正数向量，Ａ为ｎ×ｍ的平差模型系数阵，ＥＡ为Ａ中元素对应的改正数构成的ｎ×ｍ的矩阵，是ｍ×１的参数向量；根据总体最小二乘的思想，利用拉格朗日乘数法设立目标方程如下：（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｉ）＆ｐｈｉ；（／ｍｉ）?（ｍｏ）＝（／ｍｏ）?（ｍｓｕｐ）?（ｍｉ）ｖ（／ｍｉ）?（ｍｉ）Ｔ（／ｍｉ）?（／ｍｓｕｐ）?（ｍｉ）Ｐｖ（／ｍｉ）?（ｍｏ）＋（／ｍｏ）?（ｍｓｕｐ）?（ｍｉ）ｖｅｃ（／ｍｉ）?（ｍｉ）Ｔ（／ｍｉ）?（／ｍｓｕｐ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｅ（／ｍｉ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｓｕｂ）?（ｍｉ）Ｐ（／ｍｉ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）ｖｅｃ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｅ（／ｍｉ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＋（／ｍｏ）?（ｍｎ）２（／ｍｎ）?（ｍｓｕｐ）?（ｍｉ）＆ｌａｍｂｄａ；（／ｍｉ）?（ｍｉ）Ｔ（／ｍｉ）?（／ｍｓｕｐ）?（ｍｏ）［（／ｍｏ）?（ｍｉ）ｌ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｉ）Ａ（／ｍｉ）?（ｍｏｖｅｒ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ）＾（／ｍｏ）?（／ｍｏｖｅｒ）?（ｍｏ）＋（／ｍｏ）?（ｍｉ）ｖ（／ｍｉ）?（ｍｏ）－（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｐ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｉ）Ｔ（／ｍｉ）?（／ｍｓｕｐ）?（ｍｏ）＆ｔｉｍｅｓ；（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｉ（／ｍｉ）?（ｍｉ）ｎ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｉ）ｖｅｃ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｅ（／ｍｉ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）］（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｏ）－（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（／ｍａｔｈ）其中，ｖｅｃ代表矩阵拉直运算，即把矩阵按列方向拉直为一个列向量，Ｐ为ｎ×ｎ的观测值向量的权矩阵，ＰＡ是把系数阵Ａ拉直后对应的ｎｍ×ｎｍ的权矩阵，λ是拉格朗日乘系数，Ｉｎ是ｎ×ｎ的单位向量；对（１）式求偏导数，得到（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｆｅｎｃｅｄｏｐｅｎ＝＇＇ｃｌｏｓｅ＝＇｝＇）?（ｍｔａｂｌｅ）?（ｍｔｒ）?（ｍｔｄ）?（ｍｉ）Ｐ（／ｍｉ）?（ｍｏｖｅｒ）?（ｍｉ）ｖ（／ｍｉ）?（ｍｏ）～（／ｍｏ）?（／ｍｏｖｅｒ）?（ｍｏ）＋（／ｍｏ）?（ｍｏｖｅｒ）?（ｍｉ）＆ｌａｍｂｄａ；（／ｍｉ）?（ｍｏ）＾（／ｍｏ）?（／ｍｏｖｅｒ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（／ｍｔｄ）?（／ｍｔｒ）?（ｍｔｒ）?（ｍｔｄ）?（ｍｓｕｂ）?（ｍｉ）Ｐ（／ｍｉ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｉ）ｖｅｃ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｓｕｂ）?（ｍｏｖｅｒ）?（ｍｉ）Ｅ（／ｍｉ）?（ｍｏ）～（／ｍｏ）?（／ｍｏｖｅｒ）?（ｍｉ）Ａ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｏｖｅｒ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ）＾（／ｍｏ）?（／ｍｏｖｅｒ）?（ｍｏ）＆ｔｉｍｅｓ；（／ｍｏ）?（ｍｓｕｂ）?（ｍｉ）Ｉ（／ｍｉ）?（ｍｉ）ｎ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏｖｅｒ）?（ｍｉ）＆ｌａｍｂｄａ；（／ｍｉ）?（ｍｏ）＾（／ｍｏ）?（／ｍｏｖｅｒ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（／ｍｔｄ）?（／ｍｔｒ）?（ｍｔｒｏｗ）?（ｍｏ）］（／ｍｏ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｏ）－（／ｍｏ）?（ｍｎ）１（／ｍｎ）?（／ｍｒｏｗ）?（／ｍｓｕｐ）?（ｍｉ）ｌ（／ｍｉ）?...
【技术特征摘要】
【专利技术属性】
技术研发人员：童小华，梁丹，金雁敏，
申请(专利权)人：同济大学，
类型：发明
国别省市：31
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。