 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
一种新的增量式关联规则数据挖掘方法研究
下载积分：100
内容提示：
文档格式：PDF|
浏览次数：3|
上传日期： 15:47:04|
文档星级：
该用户还上传了这些文档
一种新的增量式关联规则数据挖掘方法研究.PDF
官方公共微信海尔BCD-100AJ
海尔BCD-100AJ(Haier海尔BCD-100AJ)
海尔BCD-100AJ
产品类别：
产品类别：
温控方式：
制冷方式：
额定耗电量：
能效等级：
海尔冰箱相关文章行情经典时尚样式的海尔冰箱值得选择，今天为大家介绍这款海尔BCD-192KTJX，单开门冰箱。1级能效等级，......行情海尔BCD-649WDCE冰箱促销中，原价7999元，目前这款海尔BCD-649WDCE冰箱在商家“石家庄鹏飞腾业”的......行情海尔BCD-521WDPW是一款性价比颇高的冰箱产品，适合家庭是使用。目前这款海尔BCD-521WDPW冰箱在商家......行情海尔BCD-626WABCB可以说是海尔的一款高端产品，功能不错，外观设计也是十分大气。目前这款海尔BCD-......
与本产品有相同属性的产品
温控方式：
额定耗电量：
您可能感兴趣的商品
相同品牌&,&&相同产品类别&
好评率83.33%
同价位&,&&相同品牌&
好评率100%
关注 海尔BCD-100AJ 的网友还关注
海尔BCD-100AJ 网友点评
4-500字，简明扼要的概括您对这个产品的感受
海尔最新冰箱
安装ZOL新版手机客户端，随时随地查报价
相关类别：
热门类别：
联系电话：400-699-9999
官方网站：
联系地址：青岛中国.青岛海尔路1号海尔工业园内 &每一次都在徘徊孤单中坚强，每一次就算很受伤也不闪泪光，我一直有双隐形的翅膀，带我飞飞过绝望&&&，《隐形的翅膀》这首歌相信时下的年轻一代都非常熟悉，它所呈现出的那种对于生活的顽强与拼搏态度总是能够激励我们不断向前，没有不可能的事情能够阻挡脚步，也许这一切的故事终将过去，而感恩生活的态度却会始终如一。近年来游戏市场迎来了井喷式的发展，而对于广大的电竞玩家来说，DIY的趋势也越来越强，究竟该如何选购称心的装备，以此来更好为实战提供借鉴呢?今天小编将要为您介绍的就是一款适应紧张刺激的对战，可以轻松收获一血的黑爵AJ100激光游戏，接下来请看详细介绍。
黑爵AJ 100 FirstBlood&一血&鼠标采用了魅力专业的设计，承载的价值已经远远的超过了游戏本身，绝对能够为游戏进行助力，享受胜利的氛围。在包装的搭配上，本款鼠标采用了坚固耐用的硬卡天地盖包装设计，进口黑牛卡材质提升了产品质量，专业塑封膜确保安全完整性，标配的绒布袋、驱动盘、脚贴及说明书等细节配件都呈现出了人性化的魅力，绝对能够满足不同用户的实际需求。
作为一款传递生活价值的专属产品，黑爵AJ100 FirstBlood鼠标在外观设计和整体构造上就让人眼前一亮，个性十足的&歪把子&曲面造型更为符合人体工学，极致舒适体验方便用户手掌紧密贴合，电竞游戏手感更加出众，用户指令都得到良好反馈，每一个指令都会得到完整贯彻，出击与战斗收获&一血&很容易。本款鼠标表面通过钢琴烤漆技术结合镭雕图案打造，&AJ&，& FirstBlood&等简美标识充满了时尚意味。
黑爵AJ100 FirstBlood凭借安华高A9800高精度光学引擎的搭载，能够轻松实现150ips的移动速度和30G极限加速度，鼠标移动及定位能力得到增强，桌面适应能力出色，不同的环境都是运行自如，绝对可以让你的战斗如虎添翼。黑爵AJ100 鼠标在超频模式下最高可实现16000DPI的巅峰效果，如此卓越的效果不在只存在于理论之中，将真真正正的呈现在您面前。
没魅力，再好的鼠标也不会受人追捧，黑爵AJ100 FirstBlood鼠标内置1600万色APM呼吸灯颜色及点亮模式可通过驱动自定义调节，魅力灯光随着游戏节奏而自由变换，让您的鼠标更加炫酷;鼠标搭载顶级欧姆龙白点微动，使用寿命长达1000万次，手感反馈清晰;鼠标MCU内置了可刷写闪存，预先写好配置在鼠标内，免去繁琐的驱动和安装程序;此外本鼠标搭载了内置上下盖双配重设计，合理调整重心，重量达到40g，适应出色魅力随行。
黑爵官方网站：http://www..cn/
黑爵天猫旗舰店：/?spm=a220o..0.XtLkjX
黑爵AJ100 FirstBlood一血激光游戏鼠标选购地址：
/item.htm?spm=a220o.03-.1.ZOWUAu&id=&mt=&scene=taobao_shop&refType=3c
【产品名称】黑爵【参考价格】199元【推荐商家】
【联 系 人】吴灿【联系电话】9【商家Q Q】(请注明来自中关村在线福建站的网友)【网店地址】【联系地址】福建省福州市大利嘉城三楼A3-30-31号商铺【友情提示】扫描以上二维码，添加中关村在线福建站官方微信【微信导购】微信搜索添加"fjzol520",即可掌握福建最新最全手机、笔记本等数码资讯。【购买提示】购买中提及中关村在线福建站(福建IT资讯网)的网友均可享受更低优惠！
适用类型 鼠标大小
硬件论坛精选
下载中关村在线Android 客户端
下载中关村在线 iPhone 客户端
下载中关村在线Windows8客户端
成为中关村在线微信好友
4￥2255￥2496￥597￥6208￥3209￥8510￥909
立即登录，签到赢金豆~
登录后即可查看通知，和作者互动中国矿业大学环境与测绘学院gps测量原理与利用 第九章GPS测量数据处理_简单学习网
当前位置： >>
中国矿业大学环境与测绘学院gps测量原理与利用 第九章GPS测量数据处理
第九章 GPS测量数据处理王 坚中国矿业大学环境与测绘学院 第九章GPS测量数据处理§9-1 GPS测量数据处理概述 §9-2 GPS基线向量的解算 §9-3 GPS定位成果的坐标转换 §9-4 GPS基线向量网平差 §9-5 GPS高程 9-1 GPS测量数据处理概述GPS接收机采集记录的是GPS接收机天线至卫星伪距、载波相位 和卫星星历等数据。如果采样间隔为20秒，则每20秒记录一组观 测值，一台接收机连续观测一小时将有180组观测值。观测值中有 对4颗以上卫星的观测数据以及地面气象观测数据。GPS数据处理 要从原始的观测值出发得到最终的测量定位成果，其中数据处理 过程大致分为GPS测量数据的基线向量解算、GPS基线向量网平 差以及GPS网平差或与地面网联合平差等几个阶段。数据处理的 基本流程如图9-1所示。 9.1.1 数据传输大多数的GPS接收机（如ASHTECH，TRIMBLE等型号），采集 的数据记录在接收机的内存模块上。数据传输是用专用电缆将接 收机与计算机连接，并在后处理软件的菜单中选择传输数据选项 后，便将观测数据传输至计算机。数据传输的同时进行数据分流， 生成四个数据文件：载波相位和伪距观测值文件、星历参数文件、 电离层参数和UTC参数文件、测站信息文件。 经数据分流后生成的四个数据文件中，除测站文件外，其余均 为二进制数据文件。为下一步预处理的方便，必须将它们解译成直 接识别的文件，将数据文件标准化。 9.1.2 预处理GPS数据预处理的目的是：①对数据进行平滑滤波检验，剔除粗 差；②统一数据文件格式并将各类数据文件加工成标准化文件 （如GPS卫星轨道方程的标准化，卫星时钟钟差标准化，观测值 文件标准化等）；③找出整周跳变点并修复观测值(整周跳变的修 复见5.3.3)；④对观测值进行各种模型改正。1.GPS卫星轨道方程的标准化 数据处理中要多次进行卫星位置的计算，而GPS广播星历每小时有一 组独立的星历参数，使得计算工作十分繁杂，需要将卫星轨道方程标准化 ，以便计算简便，节省内存空间。GPS卫星轨道方程标准化一般采用以时 间为变量的多项式进行拟合处理。 将已知的多组不同历元的星历参数所对应的卫星位置Pi(t)表达为时间t 的多项式形式： 利用拟合法求解多项式系数。解出的系数ain记入标准化星历文件， 用它们来计算任一时刻的卫星位Z。多项式的阶数n一般取8~10 就足以保证米级轨道拟合精度。 拟合计算时，时间t的单位需规格化，规格化时间T为：式中Ti为对应于ti的规格化时间；t1和tm分别为观测时段开始和 结束的时间。很显然，对应于t1和tm，T1和Tm分别为-1和+1。对任 意时刻|ti|≤1。 需指出的是，拟合时引进了规格化的时间，则在实际轨道计算时 也应使用规格化的时间 。 2.卫星钟差的标准化 来自广播星历的卫星钟差（即卫星钟钟面时间与GPS系统标准时 间之差Δts）是多组数值，需要通过多项式拟合求得唯一的，平滑 的钟差改正多项式。用于确定真正的信号发射时刻并计算该时刻 的卫星轨道位Z，同时也用于将各站对各卫星的时间基准统一起 来以估算它们之间的相对钟差。当多项式拟合精度优于±0.2ns时， 可精确探测整周跳变，估算整周未知数。 钟差的的多项式形式为：式中：a0,a1,a2为星钟参数，t0为星钟参数的参考历元。 由多个参考历元的卫星钟差，利用最小二乘法原理求定多项式系数 ai，再由(9-3)式计算任一时刻的钟差。因为GPS时间定义区间为一个星 期，即604800秒，故当t-t0&属于下一GPS周)时t应减去 604800，t-t0&-属于上一GPS周)时t应加上604800。 3.观测值文件的标准化 观测值文件是容量最大的文件。观测值记录中有对应的卫星号， 卫星高度角和方位角，C/A码伪距，L1、L2的相位观测值，观测 值对应的历元时间，积分多普勒记数，信噪比等。 不同的接收机提供的数据记录有不同的格式。例如观测时刻这个 记录，可能采用接收机参考历元，也可能是经过改正归算至GPS 标准时间。在进行平差（基线向量的解算）之前，观测值文件必 须规格化、标准化。具体项目包括：① 记录格式标准化。各种接收机输出的数据文件应在记录类型、 记录长度和存取方式方面采用同一记录格式。 ② 记录项目标准化。每一种记录应包含相同的数据项。如果某 些数据项缺项，则应以特定数据如“0”或空格填上。 ③ 采样密度标准化。各接收机的数据记录采样间隔可能不同， 如有的接收机每15秒钟记录一次，有的则20秒钟记录一次。标 准化后应将数据采样间隔统一成一个标准长度。标准长度应大 于或等于外业采样间隔的最长的标准值。采样密度标准化后， 数据量将成倍数地减少，所以这种标准化过程也称为数据压缩。 数据压缩应在周跳修复后进行。数据压缩常用多项式拟合法。 压缩后的数据应等价于被压缩区间内的全部数据，且保持各压 缩数据的误差独立。④ 数据单位的标准化。数据文件中，同一数据项的量纲和单 位应是统一的，例如，载波相位观测值统一以周为单位。 9.2GPS基线向量的解算在第五章中我们着重讨论了由双差观测值列出误差方程式，然后 利用最小二乘平差原理求解基线向量的方法。由于未知数个数和 误差方程个数很多，平差解算的工作量很大 . 本节重点将讨论9.2.2节法方程的组成及解算，双差观测值模型直 接从第五章引用，不作为必修内容。有兴趣的同学可自学，自学 方式为：课后在网上通过与教师讨论的形式进行。故9.2.1 双差观 测模型部分（省略）。9.2.2 法方程的组成及解算1、方程式的组成在第五章中我们着重讨论了由双差观测值列出误差方程 式，然后利用最小二乘平差原理求解基线向量的方法。由于 未知数个数和误差方程个数很多，平差解算的工作量很大。 本节将重点讨论9.2.2节法方程的组成及解算，双差观测值模 型直接从第五章引用即： kj j j DD (t ) ? ? (t ) ? ? (t ) ? ? k (t ) ? ? k (t ) 12 i 2 i 1 i 2 i 1 i ?? x ? ? 2? ? ?l k (t ) ?m k (t ) ?n k (t ) ? ?? y ? ? ??N k ? f k c[( ? k (t ) ? ? k (t ) ? ? j (t ) ? ? j (t )] ? ? f /c ? 2 2 2 ? ? 2? 20 1 20 1 ? ? ?? z ? ? 2?（9-4）令：??N k (t ) ? ?N k ? ?N j （9-5）j j jk ??Lk (t ) ? f / c( ? k (t ) ? ? k (t ) ? ? (t ) ? ? (t )) ? DD? (t ) 20 1 20 1 12（9-6）误差方程式的形式为：?? x ? ? 2? k (t ) ? ? f / c ? ?l k (t ) ?mk (t ) ?n k (t ) ? ?? y ? ? ??N k ? ??Lk (t ) V ? 2 2 2 ? ? 2? ? ? ?? z ? ? 2? （9-7）式中：T j ? 1 ? V (t ) ? ? V (t ) V 2 (t ) ? V n ? 1(t ) ? ? ? ? ? 当两站同步观测的卫星数为 n j 时，误差方程式如下：V (t ) ? A(t )? X 2 ? B(t )??N ? ??L(t )式中：? ?l1 (t ) 2 ? ? 2 ? ?l2 (t ) A(t ) ? ? f / c ? ? ? ? j ? ?l n ? 1 (t ) ? 2 0 ? 0? ?1 ?0 1 ? 0? ? B (t ) ? ? ?? ? ? ?? ? ? ?0 0 ? 1? ? 1 ??L(t ) ? ? ??L (t ) ? ? ? ??N ? ? ??N1 ? ? ?m1 (t ) 2 ?m 2 (t ) 2 ? ?m n 2 j ? 1 (t ) ? ? ? ? ? ? ? j n ? 1 (t ) ? ?n 2 ? ?n1 (t ) 2 ?n 2 (t ) 2 ?（9-8）?X? ?x 2 2??y2?z2?T??L2 (t ) ???N 2?T ? 1 (t ) ? ? ? ? T j ? ??N n ? 1 ? ? ? ??Ln j 如果在基线两端对同一组卫星观测的历元数为nt，相应的误差方程式组为：?? X ? V ? ? A B? ? 2 ? ? L ? ??N ? ? ? （9-9）式中：A(t ) ? A(t ) ? B ? ? B (t ) ? 2 nt 1 T V ? ?V (t ) V (t ) ? V (t ) ? 1 2 nt T L ? ? ??L(t ) ??L(t ) ? ??L(t ) ? 1 2 nt A ? A(t ) 1 T B (t ) ? B (t ) 2 nt?T相应的法方程式： 式中：N△X+U=0(9-10)?X ? ? XN ? ? AB ?T P ? AB ? U ? ? AB ?T PL?2??N?TP为双差观测量权矩阵 2、 权的确定 在上面的法方程式中权P应如何确定?各观测量是相互独立还是相 关？是我们必须关注的问题。下面我们从单差观测量的相关性讨 论出发，给出双差观测量的权的确定方法。1） 单差观测量的相关性由单差的定义可知：观测站T1、T2，与历元t同步观测卫星Sj的观测 量之差为：j j j ?? (t ) ? ? (t ) ? ? (t ) 2 1 （9-11）如果同一历元，还同步观测了另一颗卫星Sk，则同理可得：由矩 阵表示为：?? (t ) ? r (t ) ? ? (t )（9-12） 式中各量表示为：? (t ) ? ? (t )1r (t ) ??k? (t ) 21 0k? (t )1j? (t ) 20? 1? ?j?T??1 ?0 ?0 ?1观测量单差的方差和协方差阵：D (t ) ? r (t ) ? D (t ) ? rT (t ) ?? ? （9-13）由于： 由此得：D (t ) ? ? 2E(t )?E(t)：单位矩?1 0 ? D (t ) ? ? 2 2 ? ? ?? ?0 1 ? ?? k (t ) ? ? k (t ) ? ? k (t ) 2 1（9-14）从上面的协方差阵可知，两观测站同步观测两颗不同卫星的单差，其间是 不相关的。这一结论可推广到一般情况。 2） 双差观测量的相关性 设在观测站T1、T2，与历元t同步观测卫星Si、 Sj、 Sk,并取Si 作为参考星，则：??? j (t ) ? ?? j (t ) ? ?? i (t ) ??? k (t ) ? ?? k (t ) ? ?? i (t )（9-15）由矩阵表示为：??? (t ) ? r(t ) ??? (t )（9-16）?? j (t ) ?? k (t )? ??式中各量表示为：?? (t ) ? ???i (t ) ? ?Tr (t ) ? ???1 ??1 ?1 0? 1? ?0?观测量双差的方差和协方差阵:D (t ) ? r (t ) ? D (t ) ? rT (t ) ??? ??由于： D?? (t ) ? 2? 2 E (t )（9-17） 由此得：?2 D??? (t ) ? ? 2 2 ? ? ?1 ?1? ? ? 2? ?（9-18）从上面的协方差阵可知，两观测站同步观测两颗不同卫星的双差， 其间是相关的。由此可得到权阵:? 1? 2 P ? (t ) ? ? ?? 2 3 ? ?1 2? ?1?1? ? ? 2? ?（9-19）当同步观测nj颗卫星时，相应的权阵为：? n j ?1 ?1 ? n j ?1 1 1 ? P ? (t ) ? ? ?? 2? 2 n j ? ? ? ? ?1 ? ?1 ? ? ? ? ?1 ? ? ?1 ? ? ? ? ? n j ?1?（9-20） 如果同步观测历元数为nt时，则相应双差的权阵为：?P (t ) ? ??? 1 ? P ? ?? ?? ? ? ? 0 ?0 P ? (t2 ) ?? ? 0? ? ? ?? ? ? 0 ? ? ? P ? (tnt ) ? ?? ? 0（9-21） 9.2.3 精度评定1.单位权中误差估值m0 ? V T PV (n ? k ? 2) (9 ? 22)2.平差值的精度估计：未知数向量X中任一分量的精度估值为:mxi ? m0 1 p xi (9 ? 23)式中，Pxi由N-1中对角元素求得: p xi ? 1 Qxixi 基线长b???x0 12? ?X 12? ? ??y20 12? ?Y12? ? ??z20 12? ?Z12? ，在 ??X20 120 0 , ?Y12 , ?Z12? 处展开后:? b ? f T ?X(9 ? 24) 由协因数传播定律可得:Qbb ? f T Q?X f基线长度b中误差估值为:mb ? m0 Qbb (9 ? 25)基线长度相对中误差估值为:f b ? mb / b ?106 (9 ? 26)下面给出某一基线向量解算结果： 基线端点号：01～05 基线向量值：ΔX=-,ΔY=-,ΔZ=-,S= 向量标准差：MX=0.0010466 MY=0.0013094 MZ=0.0016443 MS=0.0009相关系数阵： dx dx 1.000000 dy -0.599817 dz -0.221641dy1..545551dz1.000000 9.2.4基线向量解算结果分析基线向量的解算是一个复杂的平差计算过程。实际处理时 要顾及 时段中信号间断引起的数据剔除、劣质观测数据的发现 及剔除、星座变化引起的整周未知参数的增加，进一步消除传 播延迟改正以及对接收机钟差重新评估等问题。基线处理完成 后应对其结果作以下分析和检核： 1 观测值残差分析。平差处理时假定观测值仅存在偶然误差。理论 上，载波相位观测精度为1%周，即对L1波段信号观测误差只有2mm。 因而当偶然误差达1cm时，应认为观测值质量存在系统误差或粗差。 当残差分布中出现突然的跳变时，表明周跳未处理成功。 2 基线向量环闭合差的计算及检核。由同时段的若干基线向量组 成的同步环和不同时段的若干基线向量组成的异步环，其闭合差应 能滞相应等级的精度要求。其闭合差值就小于相应等级的限差值。 基线向量检核合格后，便可进行基线向量网的平差计算（以解算的 基线向量作为观测值进行无约束平差）。平差后求得各GPS之间的相 对坐标差值，加上基准点的坐标值，求得各GPS点的坐标。 3.基线长度的精度。基线处理后基线长度中误差应在标称精度值内。 多数接收机的基线长度标称精度为5～10±1～2ppm?D（mm）. 对于20km以内的短基线，单频数据通过差分处理可有效地消除电离层 影响，从而确保相对定位结果的精度。当基线长度增长时，双频接收 机消除电离层的影响将明显优于单频接收机数据的处理结果。4.双差固定解与双差实数解。理论上整周未知数N是一整数，但平差解 算得的是一实数，称为双差实数解。将实数确定为整数在进一步平差时 不作为未知数求解时，这样的结果称为双差固定解。短基线情况下可以 精确确定整周未知数，因而其解算结果优于实数解，但两者之间的基线 向量坐标应符合良好（通常要求其差小于5cm）。当双差固定解与实数 解的向量坐标差达分米级时，则处理结果可能有疑，其中原因多为观测 值质量不佳。基线长度较长时，通常以双差实数解为佳。 9.3 GPS定位成果的坐标转换GPS坐标定位成果（包括单点定位的坐标以及相对定位中解 算的基线向量）属于WGS-84大地坐标系坐标（因为卫星星历是以 WGS-84坐标系为根据而建立的），而实用的测量成果往往是属于 某一国家坐标系或地方坐标系（或叫局部的，参考坐标系）。参 考坐标系与WGS-84坐标系之间一般存在着平移和旋转的关系。实 际应用中必须研究GPS成果与地面参考坐标系统的转换关系。 本节先介绍GPS定位结果的表示方法，然后介绍将GPS定位结 果转换为国家/地方独立坐标系的方法，最后讨论这几种转换方法的应用。 9.3.1GPS定位结果的表示方法WGS-84大地坐标系是GPS卫星定位系统采用的大地坐标系，因而， 所有利用GPS接收机进行测量计算的成果均属于WGS-84。我们知道，GPS定位有单点绝对定位和点间相对定位两种方法，定 位结果的表示形式也随结果的性质不同而不同，但都以WGS-84坐 标系作为参考体。 单点定位确定的是点在WGS-84坐标系中的位置。大地测量中点的 位置常用大地纬度B，大地经度L和大地高H表示，也常用三维直角 坐标X，Y，Z表示。 相对定位确定的是点之间的相对位置，因而可以用直角坐标差Δ X， Δ Y，Δ Z表示，也可以用大地坐标差Δ B、Δ L和Δ H表示。相对定 位时其中一个点是固定点。设为1号点，其坐标为X1、Y1、Z1或B1、 L1、H1，则另一点（2号点）的三维直角坐标和大地坐标可分别求 得如下： ?X ? ?X ? ? 2 ? ? 1 ? ? ?X ? ?Y ? ? ?Y ? ? ? ?Y ? ? 2 ? ?1 ? ? ? ? Z ? ? Z ? ? ?Z ? ? ? ? 2 ? ? 1?(9 ? 27)?B ? ?B ? ? 2 ? ? 1 ? ? ?B ? ? L ? ? ? L ? ? ? ?L ? ? ? 2 ? ? 1? ? ? ?H ? ?H ? ?H ? ? ? ? 2 ? ? 1?(9 ? 28)如果建立以固定点为原点的站心地平空间直角坐标系，参照（2-6） 式则2号点在该坐标系内的坐标X、Y、Z与基线向量Δ X，Δ Y，Δ Z关系 为：- sinB1sinL1 ? X ? ?sinB1 ? ?Y ? ? ? sinL cosL1 1 ? ? ? ? Z ? ?cosB1cosL1 cosB1sinL1 ? ? ? cosB1? ??X ? ? 0 ? ? ??Y ? ? ? sinB1 ? ??Z ? ? ? ? (9 ? 29)或??X ? ?- sinB1cosL1 ??Y ? ? ?- sinB sinL 1 1 ? ? ? ??Z ? ? cosB1 ? ? ?- sinB1 cosL1 0cosB1cosL1 ? ? X ? cosB1sinL1 ? ?Y ? ?? ? ? ?Z ? sinB1 ?? ?(9 ? 30) 如果以天顶距Z天，方位角A和水平距离D来表示2号点在该站心空间直角坐标系内的位置，则有:? X ? ? D ? cos A ? ?Y ? ? ? D ? sin A ? ? ? ? ? ? Z ? ? D ? ctgZ天 ? ? ? ? ?? D ? X 2 ?Y 2 ? ? A ? arctg(Y / X ) ? Z 天 ? arctg( Z / D )? ? ?(9 ? 31)或(9 ? 32)9.3.2GPS定位成果至国家/地方参考椭球的二维转换GPS定位技术用于测量，其所得结果多为三维的基线向量，而其基线向 量又构成一个GPS向量网，基线向量网是三维空间的。二维转换就是提取出 的三维基线向量中的二维信息，构成一个二维GPS向量网。 二维转换的目的是将三维的GPS基线向量网变换投影至国家大地 坐标系/地方独立坐标上去，或者说是将GPS基线网变换投影成与国 家大地测量网或与地方独立测量控制相匹配兼容。其要点是：使GPS 基线向量与常规地面测量控制网原点重合，起始方位一致，这样就 使两者在方位上具有可比性，而在坐标和边长上只存在两个系统间 尺度差影响。下面介绍二维转换的基本方法和步骤。1.GPS三维基线向量网的平移变换设常规地面测量控制网的原点在国家大地坐标系中的大地坐 标为B0、L0、H0（H0=h0+ξ 0），于是可求得该点在国家大地坐标 系中的直角坐标X0、Y0、Z0:? X ? N ? H cos B cos L ? 0 0 0 0 0 ? ? Y ? N ? H cos B sin L ? 0 0 0 0 0 ? ? N 1 ? e2 ? H ? sin B ? Z ?? 0 ? 0 0? 0? ? ??????(9 ? 33)? N0 ? a1 ? e2 sin 2 B 0(9 ? 34)其中，a、e2为国家大地坐标系参考椭球的长半径和第一偏心率。 再设GPS网在原点的三维直角坐标为 ，于是可求得GPS网平 移至地面测量控制网原点的平移参数为:?X ? X ? X 0 ? 0 ? ? ?Y ? Y ? Y 0 ? 0 ? ?Z ? Z ? Z 0 ? 0 ?(9 ? 35)于是，GPS网中各点坐标经下式变换就得到了在国家大地坐标系中 的三维直角坐标:X ? X 1 ? ?X ? 1 ? ? Y ? Y1 ? ?Y ? 1 ? Z ? Z1 ? ?Z ? 1 ? (9 ? 36)利用反算公式则得各点在国家大地坐标系中的大地坐标B1、L1、H1。 2.GPS网在国家大地坐标系内的二维投影变换由上述平移变换可知，GPS网各点已算得了在国家大地坐标系中的 大地坐标，为使GPS网与地面测量控制网在起始方位上一致，可利 用大地测量学中的赫里斯托夫第一微分公式，即使同一椭球面上 的网互相匹配。公式如下:dB1 ? P 1 dB0 ? P3 ? ds s ? ? P4 dA 0 ? ? ? dL1 ? Q1dB0 ? Q3 ? ds s ? ? Q4 dA0 ? dL0 ? ? (9 ? 37)其中，dB0、dL0为两网在原点上的纬、经度差。 Ds/s为两网在尺度上的差。 dA0为两网在起始方位上的差。 P1、P3、P4、Q1、Q3、Q4为微分公式的系数。 GPS网经平移变换后，已在原点上与地面网完全重合，因此有: dB ? 0 ? ? 0 ? dL ? 0? 0 ?(9 ? 38)在进行二维投影变换时，通常不确知两网在尺度上的差异（这一 问题1留待GPS网与地面网的联合平差或约束平差时论述），因而 可设:ds / s ? 0 (9 ? 39)两网在起始方位上的偏差需要计算。为此，须有地面网原点至起 始方位点的大地方位角A0和GPS网在相应方位上的大地方位角A0。A0和 A0可分别利用大地测量主题反解公式求得。于是:dA0 ? A0 ? A0 (9 ? 40)这样，赫里斯托夫第一类微分公式就简化成:dB1 ? P4 dA0 ? ? dL1 ? Q4 dA0 ? (9 ? 41) 最后得GPS网各点在国家大地坐标系内与此地面网点原点一致、起始方位一致的坐标为:B1 ? B1 ? dB1 ? ? L1 ? L1 ? dL1 ? (9 ? 42)利用高斯正算公式或其它平面投影变换公式可得GPS各点在国家平面坐标系 内的坐标X1和Y1。3.GPS网投影变换至地方独立坐标系将投影变换至国家大地坐标系上的GPS网再投影至地方独立坐标系， 方法如下：地方独立坐标系对应着一个地方参考椭球，该椭球与国 家参考椭球只存在长半径上的差异da, 因而，根据椭球变换的投影 公式有:dB ' ? ? a sin 2 B1 M 1 ? ? da dL ' ? 0 (9 ? 43) 其中a为椭球扁率。而M 1 ? a ?1 ? e 2 ? ? ? ??1 ? e2sin 2 B1 ? ? ? ?(9 ? 44)于是得GPS网点在地方参考椭球上的大地经纬度为:B L1 1? B 1? dB1 ? ? ? ? ? L1 ?(9 ? 45)4.地方独立坐标系内的平面投影公式我们知道，进行高斯投影变换或作其它平面变换时，需要知 道某一椭球的几何参数。由于地方参考椭球具有自己的几何参 数，它们不同于国家大地坐标系所参考椭球的几何参数，因而， 在进行平面投影计算时，必须利用该地方参考椭球的几何参数。 当然，大地经纬度也必须是属于地方参考椭球的。 9.3.3GPS定位成果转换为国家大地坐标系的三维坐标二维转换后GPS网与地面网在原点和起始大地方位上已完全重 合一致，而三维转换就是使GPS网与地面网在原点和空间起始方向上完全重合一致。首先也作GPS网的平移变换，使GPS网与地面网在原点上重合 一致，方法同上。然后，在原点建立站心空间直角坐标系，计算 出GPS网与地面网在起始方向上的方位角A和高度角β 。A和β 的 计算公式为:? A ? arctan ?Y X ? , ? ? arctan ? Z ? ? X2 ?Y2 ? ? (9 ? 46)于是，GPS网与地面网在起始方向上的方位角差和高度角差为:dA ? A地 ? a G ? ? d? ? ? 地 ? ? G ? (9 ? 47) 设GPS网各点相对于原点的三维直角坐标差为Δ X，Δ Y，Δ Z， 则各点经三维转换后相对于原点的三维直角坐标差Δ X1，Δ Y1， Δ Z1为:??X 1 ? ??X ? ? ? ? ? ??Y1 ? ? Q0 R地 ?dA, d? ?Q0 ??Y ? ??Z1 ? ??Z ? ? ? ? ?? 1 ? R地 ? dA, d ? ? ? ? dA ?-d /cosa n ? β dA 0 0 cosa n dβ ? ? sina n dβ ? ? 1 ?(9 ? 48)(9 ? 49)?- cosL0 sinB0 ? Q0 ? ?- sinL0 sinB0 ?cosB0 ?- sinL0 cosL0 cosL0 0? ? sinL0 cosB0 ? ? sinB0 ? cosB0(9 ? 50)最后得各点经三维变换后在国家大地坐标系内的三维直角坐标为:X 1 ? X ? ?X 1 ? ? Y1 ? Y ? ?Y1 ? Z 1 ? Z ? ?Z 1 ? ? (9 ? 5 1 ) 这样便可求得各点在国家大地坐标系内的大地坐标B1、L1、H1。若要再 将GPS网投影变换至地方坐标系内，可利用上述方法作类似转换。9.3.4二维和三维转换方法的应用GPS网测建完成后，为了进行与地面网的联合平差或约束平差，可首先利用上述转换方法将GPS网转换至国家大地坐标系内或地方独立坐标系内，再作平差计算。由于GPS网与地面网在原点和起始方位上已完全重合一致，因 此，可将转换后GPS网中各点的坐标，各基线向量的方位角和边长与地面网中的相应点的坐标，相应方向上的方位角和边长进行比较，可作为GPS外部检核的一种结果。 9.3.5利用不同空间直角坐标系的转换方法将GPS定位成果转换至国家大地坐标系以上所述GPS定位成果至国家大地坐标系的转换是通过三维大地坐标系统进行的。转换模型比较复杂，计算中要利用大地主题解算 公式求解大地方位角，利用大地测量微分公式进行大地坐标的变换，三维转换中还要通过建立站心坐标系计算有关转换量。下面提出用不同空产直角坐标系的转换方法实现GPS成果至国家大地坐 标系统的转换。关于不同空间直角坐标系的转换方法在第二章 2.3节已经详细论述，这里仅给出GPS定位成果转换至国家或地方坐标系的方法步骤。 1.用七参数法实现坐标转换应用七参数转换公式进行坐标转换时，GPS网与地面网应有三个以上 的重合点。当GPS网选定基准点的坐标后，便可由基准点的坐标值和基线向量的 平差值计算各GPS点的WGS-84坐标值（X Y Z）G，重合点在地面网中 的坐标由（B L H）D换算为（X Y Z）D，最后将重合点的两套坐标值 代入七参数公式解算转换参数（三个坐标平移参数，三个旋转参数， 一个尺度比参数）。重合点多于三个时，一般用平差的方法进行求 解转换参数。转换参数求出后，仍用七参数公式计算各GPS点在国家 坐标系中的坐标，便实现了GPS定位结果至国家坐标系的转换。应当指出的是，GPS定位结果中，随着基准点的坐标的不同，所求转 换参数会有很大差异。地面网重合点大地坐标中H值（大地高）往往不 能精确的给定，H=h+ζ 中高程异常最高精度为米级，所以会给转换后 的坐标带来一定误差。重合点的个数与几何图形结构也会影响转换精 度。所以求出的转换参数具有时间性和区域性。 当重合点较少时，如只有2个重合点，则只能求解部分转换参数， 如三个平移参数，三个旋转参数等。利用部分参数实现坐标转换， 检核少，精度不高。所以实际布测GPS网时，应尽量多联测地央面网 点。2.局部地区应用坐标差求解转换参数的方法因为GPS定位结果中，经基线向量网平差后获得高精度的基线向量 （Δ X Δ Y Δ Z）G，在重合点中选定一点为原点，分别求出各GPS 点对原点的坐标差，同时也求出地面网点对原点的坐标差，然后 利用式（2-24）求出尺度比与三个旋转角参数。求出4个转换参数 后，便可利用该式计算各GPS点转换至国家坐标系中的三维坐标值。 这种转换方法实践证明精度较高。3.在GPS网的约束平差中实现坐标转换 GPS基线向量网进行约束平差或GPS网与地面网联合平差时，将地 面网点的已知坐标、方位角和边长作为约束条件，坐标转换参数 也作为未知数，平差之后即得到各GPS点的地面网坐标系（国家 或地方坐标系）的坐标，平差的同时解算出坐标转换参数。具体 内容在下面几节中专述。9.3.6坐标转换中协因数阵的转换xQ 由协因数的定义可知，观测值的协因数阵 Q ? ? ? 0 2 ，其中Σ为协方差阵， σ02为单位权方差。由Q可求出观测值的权阵P=1/Q，当观测值由一个坐 标系转换为另一个坐标系时，其协因数阵也应进行转换。1.将空间直角坐标的协因数阵由式（2-3）可得:转化为大地坐标的协因数阵 QB(9 ? 52)?dB dL dH?T ? A?dX dY dZ?T 由协因数传播定律有:QB ? AQx AT (9 ? 53)2.将大地坐标的协因数阵 QB 转化为高斯平面直角坐标的协因数阵 QG 由高斯投影正算公式可得:(dX dY) T ? B(dB dL) T由协因数传播定律有:QG ? BQB B T(9 ? 54)(9 ? 55)3.直接由空间直角坐标的协因数阵 Qx 计算高斯平面坐标的协因数阵QGQG ? CQx CT(9 ? 56) 式中C=BA，简化后:?- sinB?cosL? sinL? l) C?? ? ?sin2BcosL l ? sinL - sinB(sinL cosL? l? ? sin2BsinLl ? cosL ? cosB ? ? sinBcosB l? ? (9 ? 57) 9.4基线向量网平差两观测站对GPS卫星的同步观测数据，经过平差后，解算出两 观测站间的基线向量及其方差与协方差。实际工作中，同时参加 作业的接收机可能多于两台，这样，在同一观测时间段中，便可 能在多个观测站上同步观测成网，称为GPS基线向量网。GPS基线 向量网的平差是以GPS基线向量为观测值，以其方差阵之逆阵为权， 进行平差计算，消除许多图形闭合条件不符值，求定各GPS网点的 坐标并进行精度评定。 在以两观测站之间的基线向量为观测量进行网的平差时，一般认为任 一基线向量的三个分量之间是相关的，其相关性的大小由基线向量各 自平差的结果确定；而不同的基线向量之间是相互独立的。 GPS基线向量网的平差分为三种类型：一是经典的自由网平差， 又叫无约束平差，平差时固定网中某一点的坐标，平差的主要目的 是检验网本身的内部符合精度以及基线向量之间有无明显的系统误 差和粗差，同时为用GPS大地高与公共点正高（或正常高）联合确定 GPS网点的正高（或正常高）提供平差处理后的大地高程数据；二是 非自由网平差，又叫约束平关，平差时以国家大地坐标系或地方坐 标系的某些点的坐标，边长和方位角为约束条件，顾及GPS网与地面 网之间的转换参数进行平差计算；三是GPS网与地面网联合平差，即 除了GPS基线向量观测值和约束数据以外，还有地面常规测量值如边 长、方向和高差等，将这些数据一并进行平差。非自由网平差与联 合平差一般是在国家坐标系或地方坐标系内进行，平差完成后网点 坐标已属于国家坐标系或地方坐标系或地方坐标系，因而这两种平 差方法是解决GPS成果转换的有效手段。 平差可以以三维模式进行，也可以以二维模式进行。当进行二 维平差时，应首先将三维GPS基线向量坐标及其方差阵转换至二维平 差计算面（椭球面或高斯投影平面等）。 为了解决GPS大地高的实际应用问题，需利用测区内若干公共点的正常高和高程异常将GPS大地高转换为实用的正常高。9.4.1GPS基线向量网的无约束平差GPS基线向量网的无约束平差常用的是三维无约束平差法。GPS基线向量提供的尺度和定向基准属于WGS-84坐标系，进行三维无约束平差时，需要引入位置基准，引入的位置基准不应引起观测值的变形和改正。引入位置基准的方法有三种，一种是网中有高级的GPS点时，将高级GPS点的坐标（属WGS-84坐标系） 作为网平差时的位置基准；第二种方法是网中无高级GPS点时， 取网中任一点的伪距定位坐标作为固定网点坐标的起算数据；第 三种方法是引入合适的近似坐标系统下的亏秩自由网基准。一般 采用前两种方法。 1.误差方程的列立 设网中固定点号为1，其坐标为: 基线向量观测值为: 其改正数为:?1? x1 , y1 , z1 ?TT?X ij ? ? ?xij , ?xij , ?xij ?TVij ? ?V ?xij ,V ?yij ,V ?zij ?基线向量平差值为: ? Xij? ?xij , ?xij , ?xij??T基线向量观测值的方差与协方差及其权阵分别为: D?x ,dX X 待定点近似坐标及其改正数分别为： ? ? dX , dX , dX ? ，T i i i i?1 P ? ? 2 D?xT0 i? ? xi0 , yi0 , zi0 ?待定点坐标平差值为:X i? ?x i , y i , z i ,?i ? 2,3,?, j ? 1,2,?n; i ? n 为网中点数。由?X ij ? X j ? X i、 ? X ij ? ?X ij?X ij? V?X ij以及X i ? X i0 ? dXi三式，不难得出基线向量观测值的误差方程为: 0 0 ?V ?Xij ? ?1 0 0 ? ? d xi ? ?1 0 0 ? ? d xj ? ? ?xij ? xi ? x j ? ?V ?Yij ? ? ? ? 0 1 0 ? ? d ? ? ?0 1 0 ? ? d ? ? ? ?y ? y 0 ? y 0 ? i j ? ? ? ? ? ? yi ? ? ? ? yj ? ? ij ?V ?Zij ? ? 0 0 1 ? ? d z ? ?0 0 1 ? ? d zj ? ? ?z ? z 0 ? z 0 ? ? ? ? ?? i ? ? ? ? ? ? ij i j ? ? ? ? ?(9 ? 58)写成矩阵形式为:Vij ? ? EdXi ? EdX j ? Lij (9 ? 59)¨ ,权Pij式中，E为单位阵，Lij 为（9-58）式的最后一项。 对于一端为固定点的基线向量，其误差方程式为 :0 ?V ?Xi1? ?1 0 0 ? ? d xi ? ? ?xi1 ? xi ? x1 ? ?V ?Yi1 ? ? ? ? 0 1 0 ? ? d ? ? ? ?y ? y 0 ? y ? i 1? ? ? ? ? ? yi ? ? i1 ? ? ? 0 ?V ?Zi1 ? ? 0 0 1 ? ? d z ? ?zi1 ? zi ? z1 ? ? ? ? ?? i ? ? ?(9 ? 60)写成矩阵形式为:Vi1 ? ? EdXi ? Li1 (9 ? 61i1 ，权为P) 2.法方程式的组成及解算由于各基线向量观测值之间互相独立，因而可分别对每个基线向量观测值的误差方程式组成法方程，将单个法方程的系数阵 及常数项加到总法方程的对应系数项的常数项上去。对应于(9-59)和（9-61）的法方程式分别为:? Pij dX i ? Pij dX j ? Pij Li j ? 0? Pi1 dX i ? Pi1 Li1 ? 0 (9 ? 63)(9 ? 62)总的法方程设为:NdX ?U ? 0 (9 ? 64)解算法方程后得到未知数dX为:dX ? N ?1U (9 ? 65) 各待定点坐标平差值Xi为:X i ? X i0 ? dX i (9 ? 66)3.精度评定 单位权方差估值为:? 02 ? V T PV /?3m ? 3(n ?1)?(9 ? 67)式中m为基线向量个数，n为网中点数。 平差未知数dX的方差估值为:2 Di ? ? 0 N ?1(9 ? 68)9.4.2 GPS基线向量网的约束平差 二维约束平差实际应用中以国家（或地方）坐标系的一个已知点和一个已知基 线的方向作为起算数据，平差时将GPS基线向量观测值及其方差阵 转换到国家（或地方）坐标系的二维平面（或球面）上，然后在 国家（或地方）坐标系中进行二维约束平差。转换后的GPS基线向 量网与地面网在一个起算点上位置重合，在一条空间基线方向上 重合。这种转换方法避免了三维基线网转换成二维向量时地面网 大地高不准确引起的尺度误差和变形，保证GPS网转换后整体及相 对几何关系的不变性。转换后，二维基线向量网与地面网之间只 存在尺度差和残余的定向差，因而进行二维约束平差时只要考虑 两网之间的尺度差参数和残余定向差参数。 1.GPS基线向量观测值的误差方程式V?xij ? ?dxi ? dx j ? ?xij dk ? ?yij / ? & d? ? ?xij ? xi0 ? x 0 ? j ? ? V?yij ? ?dyi ? dy j ? ?yij dk ? ?xij / ? & d? ? ?yij ? yi0 ? y 0 ? j? (9 ? 69)式中：Δx、Δy和dx、dy分别为转换后的二维基线向量观测值和待定点坐 标改正数，dκ和dα分别为尺度差和残余定向差数，当i点或j点为固定点时， 相应的改正数为0。 2.约束条件方程 (1)边长约束条件0 0 0 0 ? cos?ij ? sin ?ij cos?ij dx j ? sin ?ij dy j ? Wsij ? 0(9 ? 70)式中:0 ? ij ? arctg ??y 0 ? yi0 ?/ ?x 0 ? xi0 ?? j jW ???x0 j? xi0 ? y 0 ? y i0 j? ???1/ 2? S ij(2)坐标方位角约束条件:aij dxi ? bij dyi ? aij dx j ? b ij dy j ? Waij ? 0 (9 ? 71)式中:0 0 0 0 aij ? ?&sin aij / Sij , bij ? ?& cosaij / SijWaij ? arctg y 0 ? yi0 / x 0 ? xi0 ? aij j j?????? 9.4.3 GPS基线向量网与地面网联合平差当地面网除了已知数据（已知点坐标，已知边长和已知方位角 ）以外，还有常规观测值（如方向、边长等），则将GPS基线向量 观测值与地面网的已知数据和常规观测值一起进行的平差叫做GPS 基线向量网与地面网联合平差。 联合平差可以两网的原始观测量为根据，也可以两网单独平 差的结果为根据。平差时，引入坐标系统的转换参数，平差的同时 完成坐标系统的转换。1.三维联合平差GPS基线向量观测值的误差方程和条件方程同三维约束平 差。地面网观测值误差方程为：(1)空间弦长观测值的误差方程:VD ? ?Cij Ai dBi ? Cij A j dB j ? LDijij(9 ? 72) 式中:0 2 0 2 0 2 LDij ? Dij ? ?? ?xij ? ? ? ?yij ? ? ? ?yij ? ? ? ? ? ? 1/ 2Dij为空间弦长观测值，相应的权为PDij。 (2)方向观测值的误差方程:V?ij ? ?dzi ? Fij Ai dBi ? Fij Aj dB j ? L?ij (9 ? 73)式中:0 L?ij ? ? ij ? z i ? ? ijdzi为i点定向角未知数的改正数， zi 为i点定向角未知数近似值， 0 ? ij 为方向观测值，? ij 为大地方位角近似值。 法方程组成与解算以及精度评定与三维约束平差相同。求 单位权方差时自由度计算中应加上地面观测值个数。 三维联合平差也可以在三维直角坐标系中进行。0 由于地面网通常都有是在大地坐标系统或高斯平面坐标系统中进 行平差计算的，为计算网点的大地高程，必须以相应的精度确定点的 高程异常。但实际上高程异常的精度在东西沿海地区好于1m，而在西 北高山地区，只能保持数米的精度。这样，高程异常的误差直接影响 所求地面网点大地高的精度，从而影响据以计算的空间直角坐标的精 度，在这种情况下，大工业高的方差和协方差也难以比较可靠地确定， 这样一来便会对两网的联合平差造成不利影响。因此，通常应选择二 维联合平差的方案。 2.二维联合平差 二维联合平差时，GPS基线向量观测值的误差方程与约束条件方程 同二维约束平差。而地面网的方向和边长观测值的误差方程如下： 方向观测的误差方程为：Vl ij ? ?dzi ? aij dzi ? bij dyi ? aij dx j ? bij dy j ? Lij (9 ? 74)0 式中dzi为点i上定向角未知数改正数，其近似值为 zi ，而: 0 0 L ? zi ? lij ? aij ij边长观测值误差方程为:0 0 0 0 Vsij ? ? cos aij dxi ? sin aij dyi ? cos aij dx j ? sin aij dy j ? Lsij(9 ? 75)式中:0 Lsij ? Sij ? Sij 9.5GPS高程由GPS相对定位得到的三维基线向量，通过GPS网平差，可以 得到高精度的大地高差。如果网中有一点或多点具有精确的WGS84大地坐标系的大地高程，则在GPS网平差后，可求得各GPS点的 WGS-84大地高H84。 但在实际应用中，地面点 的高程采用正常高系统。 地面点的正常高Hr是地面 点沿铅垂线至似大地水准 面的距离。这种高程是通 过水准测量来确定的。这 就有必要找出GPS点的大 地高H84与正常高程Hr的 关系，并用一定的方法将 H84转换为Hr。 图9-3所示大地高与正常高之间的关系，其中，ζ 表示似大 地水准面至椭球面间的高差，叫做高程异常。显然，如果知道了 各GPS点的高程异常ζ 值，则不难由各GPS点的大地高H84求得各 GPS点的正常高Hr值。如果同时知道了各GPS点的大地高H84和正常 高Hr，则可以求得各点的高程异常ζ 。Hr ? H84 ? ? (9 ? 76)或? ? H84 ? Hr(9 ? 77)由此可见，研究GPS高程的意义有两个方面，一是精确求定 GPS点的正常高，一是求定高精度的似大地水准面的方法为GPS 水准。 实际上，很难获得高精度的高程异常ζ 值，而GPS单点定 位误差又较大，一般测区内缺少高精度的GPS基准点，GPS网平 差后，很难得到高精度的大地高H84。所以很难应用上式精确 地计算各GPS点的正常高。 精确计算各GPS点的正常高Hr，目前主要有GPS水准高程（简称 GPS水准）、GPS重力高程和GPS三角高程等方法。9.5.1 GPS水准高程 GPS水准高程，是目前GPS作业中最常用的一种方法。 目前，国内外有于GPS水准计算的各种方法主要有：绘等值 线图法；解析内插法（包括曲线内插法、样条函数法和Akima法） ；曲面拟合法（包括平面拟合法、多项式曲面拟合法、多面函数 拟合法，非参数回归曲面拟合法和移动曲面法）等。 下面介绍几种常用的GPS水准高程计算方法。1.绘等值线图法 这是最早的GPS水准方法。其原理是：设在某一测区，有m个GPS点， 用几何水准联测其中n个点的正常高（联测水准的点称为已知点，下 同），根据GPS观测获得的点的大地高，按（9-94）式求出n个已知点的 高程异常。然后，选定适合的比例尺，按n个已知点的平面坐标（平面 坐标经GPS网平差后获得），展绘在图纸上，并标注上相应的高程异常， 再用1～5cm的等高距，绘出测区的高程异常图。在图上内插出未联测几 何水准的(m～n)个点（未联测几何水准的GPS点称为待求点，下同）的 高程异常，从而求出这些待求点的正常高。 2.解析内插法 当GPS点布设成测线时，可应用以下曲线内插法，求定待求点 的正常高。其原理是：根据测线上已知点平面坐标和高程异常，用 数值拟合的方法，拟合出测线方向的似大地水准面曲线，再内插出 待求点的高程异常，从而求出点的正常高。(1)多项式曲线拟合法设点的ζ 与Xi（或Yi或拟合坐标）存在的函数关系（i=0,1,2,?n） 可以用下面（m≤n）次多项式(9-95) 来拟合。? ( x) ? a0 ? a1x ? a2 x2 ? ?am xmi m i i(9 ? 78)2 i在已知点处的高差 R ? ? ( x ) ? ? ，在 ? R 继而求出各点的ζ，从而获得点的Hr。(2)三次样条曲线拟合法? min条件下解各ai，当测线长，已知点多，ζ 变化大时，按 ? Ri2 ? min 解求的ai误差会增大，故通 常总采取分段计算。这样使曲线在分段点上不连续，也影响拟合精度。为此， 采用三次样条法来拟合。 设过n个已知点，ζ i和Xi（或Yi或拟合坐标）在区间[Xi,Xi+1](i=0,1,2,?,n-1)上有三次样条函数关系? ? x ? ? ? ? xi ? ? ( x ? xi )? ( xi , xi ?1 ) ? ( x ? xi )( x ? xi ?1 )? ( x, xi , xi ?1 )(9 ? 79)式中X为待点坐标，Xi、Xi+1为待求点两端已知点的坐标。其中 ? ?xi , xi?1 ? 为一阶差商。? ?xi , xi?1 ? ? (? i?1 ? ? i ) ?x i?1 ?xi ? ;? ( x, xi , xi ?1 ) ?? ( x, xi , xi ?1 )为二阶差商， ,而? &?xi ?(i ? 1,2,?, n ? 1)1 ?? & ?xi ? ? ? & ?x ? ? ? & ?x i ?1 ?? 6，满足系数矩阵为对称三角阵的线性方程组:( xi ? xi ?1 )? &( xi ?1 ) ? 2( xi ?1 ? xi ?1 )? &( xi ) ? ( xi ?1 ? xi )? &( x i ?1 ) ? 6 ?? &( x, xi ?1 ) ? ? ( xi ?1 , xi ) ? ? ? (9 ? 80)? ( x0 ) ? ? & ( x n ) ? 0 用追赶法解（9-80）方程组，可求出 ? & ( xi ) ? ( x , x ) ，而 和i i ?1? &( x) ? ? &( xi ) ? ( x ? xi )? &( x i , xi ?1 )(3)Akima法(9 ? 81)Akima法的原理是：在两个已知点间内插时，除用这两 个已知点外，还需用两已知点外二点，其目的是使曲线光滑，函 数连续。 设有六个已知点(i=1,2,3,4,5,6)，现需在3号和4号点 之间内插任一待求点，其计算公式为:? ( x) ? P ? P 1( x ? x3 ) ? P ( x ? x3 )2 ? P ( x ? x3 )3 0 2 3(9 ? 82)其中:P0 ? ? 3 ? ? P ? t3 1 ? ? P2 ? [3(? 4 ? ? 3 ) /( x4 ? x3 ) ? 2t3 ? t4 ]/( x 4 ? x3 ) ? P3 ? [t3 ? t4 ? 2(? 4 ? ? 3 ) /( x4 ? x3 )]/( x4 ? x3 ) 2 ? ? (9 ? 83) （9-83）式中的t3、t4为3号和4号点实测要素的斜率，t3用1、 2、3、4、5已知点计算，t4用2、3、4、5、6已知点计算，一般计 算公式为:ti ? (| mi?1 ? mi | ?mi?1 ? | mi?1 ? mi?2 | ?mi?1 /(| mi?1 ? mi ) | ? | mi ?1 ? mi ?2 |) (9 ? 84)式中i=3,4。mi ? ?? i ?1 ? ? i ? /( xi ?1 ? xi ) (9 ? 85)当（9-85）式分母为零时，t 3.曲面拟合法i? 1 / 2(mi ?1 ? mi )或t i ? mi。当GPS点布设成一定区域面状时，可以应用数学曲面拟合 法求待定点的正常高。其原理是，根据测区中已知点的平面坐 标x、y（或大地坐标B、L）和 ? 值，用数值拟合法，拟合出 测区似大地水准面，再内插出待求点的 ? ，从而求出待求点的 正常高。下面介绍几种常用的拟合方法。 (1)多项式曲面拟合法设点的 ? 与平面坐标x,y有以下关系:? ? f ( x, y) ? ?(9 ? 86)? 其中，f(x,y)为 ? 中趋势值， 为误差。设:f ( x, y) ? a0 ? a1x ? a2 y ? a3 x2 ? a4 y 2 ? a5 xy ? ? (9 ? 87)写成矩阵形式有：? ? XB ? ?(9 ? 88)x2 ? ? 1 ? ? x2 ? ? 2 ? ? ? ? ? x2 ? ? n ?其中：?1 x y 1 1 ? ? 1 x y 2 2 X ?? ? ?? ? ? ? ? 1 x n yn ??? ? ?a ? ?? ? 1? 1? ? ? ? 1? ?? ? ?a ? ?? ? ? ? ? 1? , B ? ? 1? ,? ? ? 1? ?? ? ?? ? ??? ?? ? ?a ? ?? ? ? 1? ? 1? ? 1?对于每个已知点，都可列出以上方程，在 ?? 2 ? min条件下，解出各ai，再 按（9-105）式内插求出待定点的高程异常，从而求出正常高。 (2)曲面样条拟合法 曲面样条拟合法是基于无限大平板小挠度方程的数学模型，设 点的与点的坐标x、y存在如下样条关系:? ( x, y ) ? a0 ? a1 x ? a2 y ? ? Fi ri 2 ln ri 2 ?i ?1 m?Fi ?1mi? ? xi Fi ? ? yi Fi ? 0i ?1 i ?1mm? ? ? ? ? ?(9 ? 89)其中:? ? yi2 )]? i ?1 ? m ? a1 ? ?2? Bixi ? i ?1 ? m ? a0 ? ?2? Biyi ? i ?1 ? ? Fi ? Pi /(16? D ) ? r 2 ? ( x ? xi ) 2 ? ( y ? yi ) 2 ? i ? ? ? a0 ?? [ Ai ? Bi ( xm2 i（9－90） xi、yi为已知点的坐标，x、y 为待求点的坐标，Ai、Bi为待定系数，Pi为点的负载，D为刚度。对于每一个公共点都可以列出一个 ? ( x, y ) 方程，对于n个公共点 列出n+3个方程，求解出n+3个未知系数a0、a1、a2、F1、F2、…、 Fn。求解方程组（9-90）时，至少应有三个公共点。 4.移动曲面法 分三步来说明该法。 (1)为了讨论的方便引入解析坐标，设内插点的坐标 (x j , y j ) (j=1,2,…， m)，相应坐标系的数据点为 ( x , y ) ,(i=1,2,…，n)，对于内插点 (x 'j , y 'j ) 作：' 'i ix & ? xi ? x? ? ' ? y & ? yi ? y j ? ?' j(9 ? 91)的变换，形成新的坐标( x&, y& )为移动坐标。如图（9-4），移动坐标是为了简化计算而引入的，下面可以看 到使用移动坐标的好处。 (2)任一数据点? w(d ) ? w ? ?( xi , y假设距离d的递减函数: i)?xi2 2 ? ? x 'j ? ? ? yi ? yi' ? ? ?(9 ? 92)将 w(d )作为权函数，权的引入是为了在移动时根据内插点到 数据点的距离给出各数据点的不同的影响程度，两点越近影响越 大。它并不像测量平差中的权是由误差定义的。目前在DTM中广 泛使用的权函数有:w(d ) ? exp(?d 2 a2 ) (9 ? 93)w(d ) ? 1 ?1 ? d 2 a 2 ?(9 ? 94)式中a为常数，可由试验给定，一般应取数据点平均间距的两倍为宜。 DTM大量内插计算表明，权的组织上入内插精度会有明显改善。 当然权函数的选取，可根据具体的发问选取不同形式，或同时利用 其他的一些信息，以利于内插。因此它使得移动法比多项式内插， 样条函数内插更加灵活。 (3)一般以某一内插点 ( x 'j , y 'j ) ，若数据点 ( xi , yi ) 中满足:?x ? x ? ?? y ? y ?i ' j 2 i' 2 i? R2(9 ? 95)可用这些数据点参加内插，则称以 移动窗口曲面内插。( x 'j , y 'j )为圆心，半径为R的圆形移动曲面法的基本原理与多项式拟合是类似的，为了进一步理解移 动法原理，下面我们以移动多项式进行公式推导。 设移动到第j个内插点 ( x 'j , y 'j ) 时，欲利用落入该点移动窗口内的 m个数据点 ( xi , yi )上的测值 ?i(i ? 1,2,?, m) ，以下列多项式:? i ? a0 ? a1x ? a 2 y ? a3 xy(9 ? 96)计算第j个内插点函数值。 在m个数据点上建立如下“误差方程”:Vi ? a0 ? a1 xi ? a2 yi ? a3 xi yi ? ? i , P i (9 ? 97) 式中:Pi ? w(d i ) ? w? ? ??xi2 2 ? xi' ? ? ? y i ? y i' ? ? ? ?，（i=1,2,?,m）令:Aj? x1 ? 1 ? ? ?? ? ? 1y1 x2 ? xmx1 y2 ? ym? x2 y 2 ? ? ? ? ? xm y m ? y1Lj ? (?1,? 2 ,?,? m )TP j ? diag( P , P ,?, P ) 1 2 m(9 ? 98)X j ? (a0 , a1 , a2 , a3 )T应用最小二乘原理:V T PV ? ? P Vi 2 ? min i (9 ? 99)可得法方程:( AT Pj Aj ) X j ? AT Pj Lj j j (9 ?100) ( AT Pj Aj )非奇异，因此: jX j ? ( AT Pj Aj )?1 AT Pj Lj j j (9 ?101)由（9-101）式求出Xj后代入（9-96）式可得:? j ? a0 ? a1x j ? a2 y j ? a3 x j y j(9 ?102)当使用移动坐标时:x& ? 0 ? j ? ? & y j ? 0? ? (9 ? 103)代入（9-102）式可得:? j ? a0(9 ?104)由此看出使用移动坐标时，移动多项式的常数项即为内插点的内插 值，这样就给计算带来了很大的方便。解算时可把a0排在最后，然后形 成下角阵，最后的元素便等于a0。使用移动坐标的另一个好处是由于计 & & 算中心化，移动坐标 ( x j , y j ) 相对较小，相对较小，形成的法方程矩阵各 元素大小相差不太悬殊，对改善法方程的数值稳定性和提高解算精度有 一定作用。形成的法方程矩阵各元素大小相差不太悬殊，对改善法方程 的数值稳定性和提高解算精度有一定作用。 移动曲面法另一个优点是它可以给出可靠的精度估计信息。 将每个数据点当作内插点，用周围的数据点按移动法计算该点 的内插值，这样对每个数据点来说，由于本点未参加内插得出 的误差具有类似于真误差的性质，所以，最后精度评定比较客 观，可信程度高。令 ? i 为i点观测值，?i 为点内插，观测值减计算值得: ?? ?i ? ? i ? ? (9 ?105)则可用下式进行精度估计:m? ? ??? ?ii/n(9 ? 106)式中n为数据点总个数。 因此，在一定情况下移动法也可用来探测粗差，如果观 测值与预测值差异较大则表明观测值可能为粗差。移动曲面 法在计算时，通常采用契比雪夫多项式为移动多项式。 设点的 ? 与平面坐标x、y的函数关系? ( x, y ) 可表示成 如下契比雪夫多项式函数: ? ? x, y ? ? ?? AmnTm ( x)Tm ( y )m ?0 n ?0mn(9 ? 107)式中Amn为拟合系数，Tm(x)、Tn(y)为变量，分别为x和y直到m 和n次的契比雪夫多项式。m 2 n ? m ? n ? 1? ! Tm ( x) ? ? ? ?1? n!!(m ? 2n)! (2 x)m?2n 2 n ?0m(9 ? 108且有以下递推公式:Tm?1 ( x) ? 2xTm ( x) ? Tm?1 ( x) (9 ?109)n时: 当观测值个数k&m?? i ? ?? AmnTm ( xi )Tn ( yi ), Pim?0 n ?0 m n(9 ? 110)式中i=1,2,…,k；p1,p2, …,Pk为权函数。 5.地形改正方法设点的高程异常 ? 可表示成:? ? ? 0 ?? r(9 ?111)? 其中 ? 为高程异常长波部分， r 为短波部分。 ? ? 0 可按本章前几节方法求出， r 为地形改正。0按莫洛金斯原理有:?r ? T r(9 ?112)其中T为地形起伏对地面点挠动位的影响，r为地面正常重力值。而:T ? G * ? ?? ?? h ? hr ? r 0 ? d? ? ? ??G*? ? 3 3 ?? ?? h ? hr ? / r0 ?d? ? 6 ?1/ 299 ? 113)其中:2 2 r0 ? ?? x ? x p ? ? ? y ? y p ? ? ? ? ? ?(9 ? 114) G――引力常数， ρ――地球质量密度， hr――参考面的高程（平均高程面）， x,y――高程格网点的坐标， xp,yp――待求点的坐标。 实际计算时，利用测区地形图，用1km×1km格网化，得 测区DTM，再按上式计算。 当测区无法得到DTM时，可采用测区GPS点观测的大地 高差来格网化。这样也能有效地提高山区GPS水准的精度。具体 计算方法是： 第一步：对测区进行1km×1km格网化，求出各格网点的坐标 。 第二步：内插出无GPS点格网近似大地高。 ?r 第三步：按上面公式进行求解T和 。 6.多项式曲面拟合精度评定为了能客观地评定GPS水准计算的精度，在布设几何水准联测 点时，适当多联测几个GPS点，其点位也应均匀地分布全网， 以作外部检核用。 (1)内符合精度根据参与拟合计算已知点的值 ? i 与拟合值 ? i' ，用求拟合残差 V i ， 按下式计算GPS水准拟合计算的内符合精度 ? ：? ? ? ?VV ? ? n ? 1?(9 ? 115)式中n为V的个数。(2)外符合精度根据核检点 ? 与拟合值 ? i' 之差，按下式计算GPS水准的外符合精度M:M ? ? ?VV ? ? n ? 1? (9 ? 116)其中n为检核点数。(3)GPS水准精度评定①根据检核点至已知点的距离L（单位：公里），按下表计算检 核点拟合残差的限值，以此来评定GPS水准所能达到的精度。 ②用GPS水准求出的GPS间的正常高程差，在已知点间组成附合 或闭合高程导线，按计算的闭合差W与表9-1中允许残差比较，来 衡量GPS水准达到的精度。表9-1 GPS水准限差等 级 允许残差（mm）三等几何水准测量? 12 L四等几何水准测量? 20 L普通几何水准测量? 30 L(4)外围点的精度估算各种拟合模型都不宜外推，但在实际工作中，测区的GPS点不 可能全部都包含在已知点连成的几何图形内。对这些外围点， GPS水准计算时只能外推，外推点的残差V按下式来估算：V ? a ? cD (9 ?117) 其中:2 2 ?c ? ? ?? ? ? ? DV ? ? D?V N ? ?? D ? ? ? D ? n ? ? ? ?? ? ?a ? ?V n ? C ? D n ? ? ?(9 ? 118)D是待求点至最近已知点的距离（单位为公里），系数a,c可根 据测区部分外围检核点按式（9-139）计算出。按（9-138）式计算 出残差V，根据表9-1估算精度。当希望外围点达到某一精度，确定 V值，按（9-139）式反求出D，可为布设联测几何水准点方案时参 考。9.5.3 GPS三角高程GPS三角高程是在GPS点上加测各GPS点间的高度角（或天顶 距），利用GPS求出的边长，按三角高程测量公式计算GPS点间的 高差，从而求出GPS点的正常高（或正高）的一种方法。 除以上三种方法外，还有求转换参数法和整体平差法。 求转换参数法的原理是：当一测区内，有一定数量点平面坐标和 高程已知，按坐标转换原理，求出参考椭球面与似大地水准面（或 大地水准面）之间的平移和旋转参数，把这些参数加入GPS网的平 差，以已知点高程约束下，通过平差，在求出各GPS点平面坐标同 时，求出点的正常高（或正高）。有文献报导，在平原地区，这种 方法求出的正常高（或正高）精度可达 5ppm? D 整体平差法的原理是：把测区内重力、水准、三角高程观测的天顶距等一并进行联合平差，求定点的三维坐标。这种方法的精度取决于已知高程点的分布及其精度。9.5.4 应用实例 在局部地区，如某一城市或地区的GPS网中，用几何水准联测部分GPS点 的正常高，用数值拟合的方法求出测区的似大地水准面，计算出未联测 几何水准GPS点的高程异常，从而求出这些GPS点的正常高。这是目前GPS 作业中，应用最广泛的GPS水准方法。 下面给出某区利用多项式曲面拟合法计算GPS点正常高的实例。 为了探讨在一个GPS网中实测几何水准点的最佳点数和点位的最佳布 设，我们对图9-5所示的GPS试验网，选用8种（图9-6）不同点数及点 位布设进行比较。 图9-6中8种方案，其中1、2、3为不同点数的试验，3、4、5、6、7、 8为不同点位布设的试验，计算结果（计算时采用多项式曲面拟合模 型）列于表9-2。表9-2 水准点布设方案计算结果比较表方 案 结点数 插值点数 μ （已知）(mm)1 10 18 ±1. 7 14.3 25 32 7 21 4.0 15.1 27 43 6 22 0 12.4 25 34 6 22 9.6 18.2 32 55 6 22 14.1 19.6 26 46 6 22 0 20.4 46 77 6 22 0 153.0 399 128 6 22 0 22.6 27 6μ （待求）(mm) Vmax(mm) V&20mm个数 从表9-2可看出，方案3为最佳。对这样的试验网（共29个点）， 实测水准点6个即足够。6个水准点的分布，在测区一侧精度为最低（方案7），布设在中部或周围有改善，周围和蹭相结合布设的效果较好，其中经构成几何图形强度较强的方案3布设为最好。 结合有关文献的讨论，可以得出以下一般布设原则：测区中联测几何水准点的点数，视测区的大小，测区似大地水准变化情况而定。一般地区以每20~30km2联测一个几何水准点为宜（或联 测GPS总点数的1/5），平原地区可少一些，山区应多一些。一个局部GPS网中最小联测几何水准的点数，不能少于选用计算模型中未知参数的个数。 联测几何水准点的点位，应均匀地布设于测区。测区周围应有几何水准联测点，由这些已知点连成的多边形，应包围整个测区。因拟合计算不宜外推，否则会发生振荡。 若测区有明显地几种趋势地形，对地形突变部位的GPS点，应联 测几何水准。 9.5.5提高GPS水准精度的措施从理论研究和实践经验可知，提高GPS水准精度应注意以下几个方面：1.提高大地高（差）测定的精度大地高（差）测定的精度是影响GPS水准精度的主要因素 之一。因此，要提高GPS水准的精度，必须有效地提高大地高（ 差）测定的精度，其措施主要有： (1)提高局部GPS网基线解算的起算点坐标的精度 研究表明，当起算点坐标有10m误差时，对其他GPS点的高 程会产生10mm的误差。因此，应尽量采用国家A、B级GPS网点为 局部GPS网的起算点。 (2)改善GPS星历的精度 有关文献分析表明，用精密星历比用广播星历可提高精 度34%。美国实施SA政策后，我们应建立自己的测轨系统。 (3)选用双频GPS接收机。 (4)观测时应选择最佳的卫星分布(5)减弱多路径误差和对流层延迟误差(6)大于10km的GPS网点应实测气象参数 实践表明，当边长大于10km，两端点气压为7mbar，气温差为2℃相对湿度为4%，此时用实测气象参数与取平均气象参数对基线处理的边长仅产生1mm误差，对大地高差产生0.1m误差。2.提高联测几何水准的精度据分析，采用四等几何水准联测的误差，约占GPS水准总误差的30%。因此，尽量采用三等几何水准来联测GPS点。对有特殊应用 的GPS网，用二等精密水准来联测，以利有效地提高GPS水准的精 度。 3.提高转换参数的精度 提高转换参数精度的方法是利用我国已有的VLBI和SLR站的 地心坐标转换参数，或利用国家A、B级GPS网点来推算转换参数。 但这一项误差在GPS水准中是次要的。 4.提高拟合计算的精度 提高拟合计算的精度办法有： (1) 根据测区似大地水准面变化情况，合理地布设已知点，并选 定足够的已知点。 (2) 根据不同测区，选用合适的拟合模型。对高差大于100m的测 区，一般要加地形改正。 (3) 对含有不同趋势地区的大测区，可采取分区计算的办法。 (4) 计算时，坐标取到米或10m，但高程异常应取到毫米。计算 结果应由计算机绘出测区高程异常等值线图，以便分析测区高程异 常变化情况，提高拟合计算精度。 从以上分析和国内外GPS水准实践情况看，在局部GPS网中，采用拟 合法进行计算，GPS水准高程的内符合精度一般可达2ppm.D左右。 对于测区面积不大的平坦地我，特别是测区内高程异常的变化有规 律的地区，公共点分布均匀的情况下，多项目式曲面拟合法能够达 到比较理想的精度。只要用三等几何水准联测已知点，点位分布合 理，点数足够，GPS水准可代替四等几何水准；在山区，只要施加 地形改正，也可达到四等几何水准的精度。 思考题1. 为什么要进行数据预处理？处理过程是怎样的？ 2. 数据的标准化的目的是什么？有哪些项目需要标准化？ 3. 双差观测值误差方程式是如何组成的？在两个测站同时观测Nj颗卫星，Nt个历元其相应的误差方程式组有何特点？4. 双差基线解的方差与协方差有何特点？权如何确定？ 5. 基线解算结果如何分析？6. GPS定位结果转换为国家坐标系的过程是怎样的？7. GPS高程是什么含义？求各未知点的GPS高程有几种计 算方法？ 谢谢！
中国矿业大学环境与测绘学院gps测量原理与利用 第九章GPS测量数据处理--中国矿业大学环境与测绘学院gps测量原理与利用 第九章GPS测量数据处理完整文章阅读,提供各年级学习资料,英语学习,日语学习,交流学习方法,学习心得和学习计划,以及各行专业文献和工作范文的学习网,坚持"用户需要的,就是我们关注的!"关注高价值的实用信息，争取提供更多实用下载资源。"
更多搜索：
All rights reserved Powered by
copyright ©right 。文档资料库内容来自网络，如有侵犯请联系客服。
