GPS测量成果在约束平差或联合平差前属于大地坐标系的WGS-84坐标系，然而我国建设工程测量多采用参心坐标系，故需将GPS点的地心坐标转换成参心坐标才能被使用^[1-2]。一般采用的方法是利用已知控制点作为约束条件进行约束平差，把GPS网强制符合到控制点所在坐标系中^[3]。由于这些已知控制点的坐标可能存在误差或粗差，如果不进行有效性检核控制处理，随意将这些控制点用作平差计算或坐标转换的约束条件，将得不到高精度的最终成果^[4]。因此想获得高质量的定位成果，必须对已知控制点进行有效性挑选，剔除较大误差或粗差的控制点。有效性控制点指能用于网平差且其成果符合GPS规范要求的已知控制点。

解算未知控制点的坐标时，需要至少3个已知点，但并不是控制点越多越好。控制点的精度、个数及其分布对平差结果有一定的影响，控制点选择不当，将会引起GPS网的扭曲和变形，变形超过一定程度，就会明显降低GPS网的精度，是不能接受的。因此，合理地选择GPS网平差的控制点，对保证GPS控制网的精度起着至关重要的作用，也是GPS数据处理中解决的关键问题。本文以建立江苏连云港市灌南全县范围内的GPS控制网为例进行分析。

一、 GPS网平差的有效控制点选择原则

在GPS网布设中，已知控制点一般有若干个，为了得到高精度高质量的GPS网平差结果，必须选择出GPS网平差有效的起算控制点，其选择原则如下：

1) 精度：在进行GPS网数据处理后得到参心坐标或所需要的坐标时，将已知控制点坐标代入进行约束平差，从而得到最终的平差结果。由此可知，已知控制点的精度是外符合精度，决定了最终的成果质量。由于受到地面变形、人为因素的影响或建立这些控制点时的测量成果质量都会影响这些控制点的坐标，故需要在测量前对这些控制进行检核。

2) 控制点数量：根据公共控制点的选取与分布要求所建立的GPS网的成果与现有成果相似度最好时，则起算的公共控制点数量越多越好。反之，通常情况下一般选3~5个起算控制点为宜，用于实现基准的转换及必要的检核，这种处理方式既保证了新老坐标成果的一致性，也保持了原有GPS网的精度^[5]。

3) 控制点分布：挑选有效控制点位置时，要纵观全局，切勿盲目选取。起算控制点不能位于网的一侧或某个角落，要平均分布于整个测区。两起算控制点间距离增大时，方位角较差和相对较差均减小，精度就会提高，当距离继续增大出现控制点各自分布在控制网边界时，会突然增大相对较差。选择位于网形的中央区域的起算公共控制点，而不能使它们分布在角落^[4-6]。如果GPS网呈带状，则不应使起算控制点位于网的一端。

二、 GPS网平差控制点有效性选择原理 1. 相似变换模型

相似变换模型(平面四参数模型)用于点在不同平面直角坐标系间的转换。设原坐标系下的平面直角坐标为(x₁，y₁)，目标坐标系下的平面直角坐标为(x₂，y₂)。平面坐标转换涉及的4个参数分别为：2个平移参数，1个旋转参数 ，1个尺度参数 ^[4]。其坐标转换模型如下

式中，x₀、y₀为平移参数，坐标单位均为m；α为旋转参数，单位为弧度；k 为尺度参数。

将上式改写为

令m=kcos α，n=ksin α，可得，α=arctan(n/m)，式(2)可转化为

将式(3)写成矩阵形式为

式中，；X=[x₀ y₀ m n]^T。有

则误差方程为

若两坐标系中有n个控制点时，则可列出误差方程为

根据最小二乘原理，由观测方程列出误差方程进而组成法方程，根据法方程，得出计算四参数的公式为

式中，P为单位矩阵。求解转换参数，最后将转换参数回代入式(1)即可完成坐标转换^[7]。

相似变换坐标转换方法是一种降维的坐标转换方法，即由三维空间的坐标转换转化为二维平面的坐标转换，避免了由于已知点高程系统不一致而引起的误差。在区域较小的情况下，应用二维四参数坐标转换将不同坐标系下的直角坐标进行转换^[8]。

2. 最优控制点提取步骤

1) 将控制点坐标换算成高斯平面坐标或空间坐标。

2) 将全部控制点利用相似变换模型求取转换参数，再将求取的转换参数回代，求出全部控制点的转换坐标，进而计算各控制点的坐标残差。

3) 分析残差值较大(大于2倍残差中误差)的点的有效性，检查点位和成果是否正确，对残差较大特别是大于2倍残差中误差的点进行剔除后重新确定坐标转换控制点。

4) 重复上述步骤2)和3)，回代残差到一定范围内和起算公共控制点在3~5个为止^[3]。

另外，在上述步骤中，不要仅仅追求转换参数的回代精度，要同时兼顾网结构和图形强度。

3. 坐标转换精度估计

根据已知控制点和坐标转换模型计算出原坐标系到目标坐标系的转换参数，然后对坐标转换精度进行估计，精度评判如下：

1) V (残差)=转换坐标-已知坐标 。

2) 平面坐标x、y残差中误差，大地高H残差中误差(见表 1)。坐标转换的精度通过求算转换参数的控制点的中误差体现。

三、 案例分析 1. 测区概况和工作内容

(1) 测区概况

灌南县位于江苏省北部海滨城市连云港市的南翼，地处北纬33°59′—34°27′、东经119°07′—119°48′之间。东濒黄海，四市交界，西接宿迁的沭阳，南邻淮安的涟水，东与盐城的响水相连。县境最大直线距离东西长71 km，南北宽30 km。县域地势南高北低，西高东低，地面高程西南部达5.9 m，东部2.0 m，地面坡降1/18 000。东北东友引河一线高程仅1.6 m。由西南向东北逐渐倾斜，地形西宽东窄。

(2) 工作内容

在灌南县第二次土地调查期间布设的(基础控制测量网点D级GPS点25个)控制网基础上，引入三等以上水准点5个，建立覆盖全县范围内的E级GPS控制网，为全县下一步开展农村宅基地及集体土地建设用地使用权确权登记发证打下坚实的基础。

(3) 测区已有资料及利用

采用委托方提供的全县D级GPS控制点成果表，该套控制点属于1980西安坐标系，中央子午线为120°，1985国家高程基准。经实地调查、踏勘后，有25个控制点符合要求，分布如图 1所示。这25个点现场保存完好、交通便利、视野开阔，周围没有大型障碍物和信号干扰设备，符合GPS观测要求，可以作为E级平面控制测量起算依据。由于该测区海拔较低，属于软土地基，故易受到地面变形影响，为有效建立GPS控制网，需对这25个控制点进行分析与挑选。

2. 数据处理分析

针对江苏连云港市灌南县GPS测量控制网，在选择有效性控制点时，主要根据上述原理进行挑选。该GPS控制网共有25个控制点，将无约束平差得到的三维空间直角坐标转换为高斯平面直角坐标 GPS(x₁，y₁) ，将GPS(x₁，y₁)和目标坐标(x₂，y₂)运用最小二乘原理，根据相似变换模型法计算得4个转换参数，然后将GPS(x₁，y₁)高斯平面直角坐标进行回代计算，得到起算控制点在地方坐标系中的坐标计算值(x₂，y₂) ，然后计算目标坐标(x₂，y₂)与计算值(x₂，y₂)的坐标差V_x、V_y(残差)。通过分析残差及中误差m_p值的大小，剔除精度较差的控制点^[9]。利用相似变换原理，求出4个参数，回代计算控制点坐标，比较计算坐标与目标坐标的差值，然后计算出x₂残差中误差和y₂残差中误差、点位中误差和两倍中误差，以及单位权中误差和点位残差，见表 2。剔除残差大于2倍残差中误差的控制点。如表 2 所示，控制点D250的V_x(3.72 cm)大于2倍m_x(1.8 cm)，D196和D217的V_y大于2倍m_y(2.55 cm)，因此剔除D196、D217和D250 3个控制点，然后对剩下的22个控制点重新计算，最后得到表 3—表 5。

其中，表 3是5个控制点坐标转换残差结果，其m_x为5.6 mm，m_y为10 mm，m_p为11.5 mm；表 4是4个控制点坐标转换残差结果，其m_x为6.6 mm，m_y为9.8 mm，m_p为11.8 mm；表 5是3个控制点坐标转换残差结果，其m_x为4.1 mm，m_y为9.2 mm，m_p为10.9 mm；根据表 3—表 5和控制点分布要求原则，在后续的约束平差中，如果选择5个控制点来进行计算，则5个控制点分别是：D221、D239、D257、D266和D269；如果选择4个控制点来进行计算，则4个控制点分别是：D221、D239、D257和D266；如果选择3个控制点来进行计算，则3个控制点分别是：D221、D257和D266。其相应的有效控制点分布如图 2所示，这些有效控制点平均分布于整个测区，防止了误差的积累。

四、 结束语

本文以连云港市灌南县第二次土地调查GPS测量为例，结合GPS网平差控制点有效性选择原则，利用相似坐标变换模型，经过反复迭代，并根据控制点分布要求原则，在后续的约束平差中，如果选择5个控制点来进行计算，则5个控制点分别是D221、D239、D257、D266和D269；如果选择4个控制点来进行计算，则4个控制点分别是D221、D239、D257和D266；如果选择3个控制点来进行计算，则3个控制点分别是：D221、D257和D266。从结果证明，本文提出的方法对这样的工程案例是有效的，可以作为借鉴与参考。