区域似大地水准面精化基本技术规定

区域似大地水准面精化基本技术规定

国家标准 GB/T2370g一2009 区域似大地水准面精化基本技术规定 Basicspeeifeatonsforregionquasrgeoiddetermimation 200g-05-06发布 2009-10-01实施 国家质量监督检验检疫总局 发布 国家标准化管理委员会国家标准
GB/T23709一2009 前 言 本标准的附录A、附录B,附录C,附录D都是资料性附录 本标准由国家测绘局提出 本标准由全国地理信息标准化技术委员会归口 本标准主要起草单位:国家测绘局测绘标准化研究所,国家测绘局大地测量数据处理中心,国家基 础地理信息中心 本标准主要起草人:郭春喜、肖学年、王斌,章磊 业
GB/T2370g一2009 引 言 区域似大地水准面精化的目的是综合利用重力资料、地形资料、重力场模型与GPS/水准成果,采 用物理大地测量理论与方法,应用移去-恢复技术确定区域性精密似大地水准面 通过似大地水准面精 化,利用GPS技术结合高精度高分辨率似大地水准面模型,已成为高程测量的一种方式 本标准依据国家强制性标准GB22021一2008《国家大地测量基本技术规定》中规定的国家将省级、 城市似大地水准面的精度和基本技术指标,结合我国已有的似大地水准面精化工程成果的经验,参考国 内外该领域的研究成果,制定了本标准,用于规范和指导我国的区域似大地水准面精化
GB/T23709一2009 区域似大地水准面精化基本技术规定 范围 本标准规定了我国国家、省级和城市区域似大地水准面精化的精度,布测原则、数据处理、精度检查 等内容 本标准适用于我国国家、省级和城市等区域似大地水准面精化 规范性引用文件 下列文件中的条款通过本标准的引用而成为本标准的条文 凡是注日期的引用文件,其随后所有 的修改单(不包括勘误的内容)或修订版均不适用于本标准,然而,鼓励根据本标准达成协议的各方研究 是否可使用这些文件的最新版本 凡是不注日期的引用文件,其最新版本适用于本标准 GB/T12897 国家一、二等水准测量规范 GB/T12898国家三、四等水准测量规范 GB/T17944 加密重力测量规范 GB/T18314全球定位系统(G;PS)测量规范 CH1002测绘产品检查验收规定 CH1003测绘产品质量评定标准 CH/T1004测绘技术设计规定 术语和定义 下列术语和定义适用于本标准 似大地水准面quasigeoid 从地面点沿正常重力方向至正常高端点所构成的曲面 3.2 似大地水准面分辨率resolutionofquasigeoiad 似大地水准面模型采用的格网尺寸 3.3 高程异常heightanomaly 似大地水准面相对于地球椭球面的高度 3 高程异常控制点controlpointofheightanomaly GPSs/水准点GPSs/levellingpoint 大地高由GPS测定、,正常高由水准测量测定的大地点 5 3 平均重力异常分辨率resolutionofaveraggraxityanomal 平均重力异常数字模型采用的格网 6 3. 重力似大地水准面gravityquasigeoid 先利用地球重力场模型、数字地形模型和地面重力资料获得该地区高分辨率剩余重力异常,然后利
GB/T2370g一2009 用移去-恢复技术确定的区域似大地水准面称为重力似大地水准面,其实质是利用重力观测数据和数字 地形模型改进由地球重力场模型确定的模型似大地水准面 总则 区域似大地水准面的参考基准 大地坐标系;2000国家大地坐标系 高程基准;1985国家高程基准 重力基准:2000国家重力基本网 4.2似大地水准面精度与分辨率 4.2.1似大地水准面的精度由格网平均高程异常相对于本区域内各高程异常控制点的高程异常平均 中误差表示 4.2.2似大地水准面的分辨率由似大地水准面模型采用的等角格网间距表示 4.2.3我国似大地水准面按范围和精度,分为国家似大地水准面、省级似大地水准面和城市似大地水 准面 各级似大地水准面的精度和分辨率应不低于表1的规定 表1各级似大地水准面精度和分辨率 似大地水准面精度 级 似大地水准面分辨率 平地、丘陵地 山地、高山地 士0.3m 士0.6m 国家 15'×15 省级 士0.1m 士0.3m 5'×5 城市 士0.05m 2.5'×2.5 地理区域较小的城市或局部的似大地水准面的精度和分辨率在满足表1 4.2. 中城市似大地水准面 要求的前提下,可根据应用需要设计 4.3格网平均重力异常及分辨率 4.3.1在似大地水准面计算时,应利用已有的或通过实测重力资料确定的格网平均重力异常,作为重 力似大地水准面计算的基础数据 4.3.2格网平均重力异常的分辨率应与似大地水准面分辨率及该区域内重力点的密度相匹配 每个 平均重力异常格网中宜有一个实测重力点,其精度应不低于加密重力点的精度 各级似大地水准面计算采用的格网平均重力异常分辨率应不低于表2的规定 4.3.3 表2各级似大地水准面格网平均重力异常分辨率 平均重力异常分辨率 级 别 平地、丘陵地 山地、高山地 5'×5' 15'×15 国家 2.5'×2.5" 省级 城市 2.5'×2.5' 4.3.4格网平均重力异常的精度以格网平均重力异常的代表误差表示,格网平均重力异常的代表误差 计算如公式(1): =2.7eA 式中: o 格网平均重力异常代表误差,单位为毫伽(10-"ms-'); 平均重力异常格网分辨率,单位为角分('). 平均重力异常代表误差系数,各种地形类别对应平均重力异常代表误差系数按表3规定 执行
GB/T23709一2009 表3地形类别与格网平均重力异常代表误差系数 类 别 地形 重力异常代表误差系数 平原 0.54 0.81 丘陵 1.08 山地 高山地 1.50 4.4数字高程模型(DEM)的分辨率和精度 4.4.1所采用的数字高程模型分辨率应不低于表4的规定 表4各级似大地水准面数字高程模型分辨率 别 级 数字高程模型分辨率 30"X30" 家 国 3"×3" 省 3"×3"”" 城 4.4.2数字高程模型应使用精度不低于国家1:50000比例尺数字高程模型的数据,其格网间距不大 于25m×25m,格网高程中误差不大于表5的要求 表5基础数据格网高程中误差 类 别 格网高程中误差 平原 士4m 丘陵 士7m 山地 士llm 高山地 士19m 4.5高程异常控制点测量精度 用于精化国家似大地水准面的高程异常控制点,其坐标和高程精度应不低于B级GPs网点和 4.5.1 国家二等水准网点的精度 4.5.2用于精化省级似大地水准面的高程异常控制点,其坐标和高程精度应不低于C级GPs网点和 国家三等水准网点的精度 4.5.3用于精化城市似大地水准面的高程异常控制点,其坐标和高程精度应不低于C级GPs网点和 国家三等水准网点的精度 4.6仪器的技术要求 似大地水准面精化所用仪器应经过法定计量检定部门检验合格后,在有效期内使用 其检验方法 和技术要求按相关标准执行 高程异常控制点的布设 5.1技术设计 .1.1技术设计的基本要求 5. 高程异常控制点布测前应进行技术设计,以得到最优的布测方案 技术设计书的格式、内容,要求 与审批程序按照CH/T1004执行 5.1.2技术设计准备 a 根据需要,收集测区范围已有的大地测量成果和资料 b 搜集测区范围内有关的地形图,交通图等资料;必要时,还应搜集有关的地震、地质资料等
GB/T2370g一2009 技术设计前,应对上述资料进行分析研究,必要时进行实地勘察,然后进行图上设计 5.1.3高程异常控制点位布设 a 高程异常控制点应均匀分布于似大地水准面精化区域; b)高程异常控制点应具有代表性,点位分布应顾及平原、丘陵和山地等不同的地形类别区域,点 位在不同地形类别均应占有一定的比例,在可能的情况下,对丘陵和山地等地形变化剧烈地 区应适当加大高程异常控制点分布密度; 各级似大地水准面的高程异常控制点宜利用不低于4.5规定精度的大地控制网点和水准 网点; 相邻高程异常控制点最大间距不宜大于公式(2)计算结果 d=7.19mc-x+ 式中: 相邻高程异常控制网点最大间距,单位为千米(km); 似大地水准面的精度,单位为厘米(em); m1 平均重力异常代表误差系数,按表3规定执行; 平均重力异常格网分辨率,单位为角分('. 5.1.4技术设计完成后应编写并上交技术设计书 5.2高程异常控制点选点与埋石 高程异常控制点的选点与埋石应符合以下要求: 高程异常控制点应满足GB/T18314规定的相应等级大地控制点点位要求,以及GB/T12897 或GB/T12898规定的相应等级的水准点点位要求 新埋设的高程异常控制点,其标石可采用GB/T18314规定的天线墩,其上埋设满足GPS和 水准测量的标志; 当利用已有大地控制点和水准点时,应检查该点的稳定性、可靠性和完好性,符合要求方可 利用; 选点与埋石工作结束后应上交选点工作总结、高程异常控制点选点图、点之记,选点搜集的各 种资料等 外业观测 GPS观测 高程异常控制点GPs测量等级按4.5规定执行,GPS观测的技术要求、观测的作业要求、外业成果 记录整理、外业数据质量检核以及成果验收与上交资料,应符合GB/T 18314的相关规定 6.2水准测量 6.2.1高程异常控制点的正常高程应以国家一等或二等水准点作为起算点,其点位应保存完好、观测 资料与成果齐全且地质条件稳定 6.2.2高程异常控制点水准测量等级按4.5规定执行,水准观测的技术要求和外业成果的记录整理应 符合GB/T12897或GB/T12898的相关规定 6.3重力测量 采用加密重力测量方法实测重力点,其技术要求按GB/T17944执行;采用航空重力测量方法获取 的地面平均重力异常值应满足公式(1)的要求 数据处理 高程异常控制点数据处理 高程异常控制点GPS测量数据处理按照GB/T18314的要求执行
GB/T2370g一2009 7.1.2高程异常控制点水准测量数据处理按照GB/T12897或GB/T12898的要求执行 7.1.3计算高程异常控制点的高程异常,按公式(3)执行 =H一h ws 式中: 高程异常,单位为米(m); gas H大地高,由GPS测量方法获得,单位为米(m); -正常高,由水准测量方法获得,单位为米(m). 7.2似大地水准面计算 7.2.1似大地水准面计算流程 似大地水准面计算流程见图1 参考重力场模型 高程异常控制点数据 格网重力异常 DEM数据 移去-恢复法等 融合技术 重力似大地水准面 拟合后的似大地水准面 图1似大地水准面计算流程 7.2.2收集似大地水准面精化区域的重力资料与数字高程模型资料,并按格网平均重力异常计算要求 对数据进行整理 7.2.3可采用地形均衡重力归算等方法完成重力点重力归算与格网平均重力异常计算 重力归算与 格网平均重力异常计算可参照附录B 7.2.4根据不同情况选择适当的参考重力场模型,采用移去-恢复技术,完成重力似大地水准面计算 重力似大地水准面计算可参照附录C 7.2.5采用融合技术消除或削弱高程异常控制点与对应的重力似大地水准面的不符值,完成与国家高 程系统一致的似大地水准面计算 重力似大地水准面与高程异常控制点的多项式与多面函数拟合法参 见附录D 似大地水准面精度检验 8.1检验点布设原则 a)检验点的点位应分布均匀,在平原、丘陵和山区等不同的地形类别以及有效区域边缘地区均应 布设检验点;应采用未参加似大地水准面计算的实测高程异常点作为检验点 国家似大地水准面检验点的间距不宜超过300km,检验点总数不应少于200个;省级似大地 水准面相邻检验点的间距不宜超过100km,检验点总数不应少于50个;城市似大地水准面相 邻检验点的间距不宜超过30km，检验点总数不应少于20个 检验点与用于区域似大地水准面精化的高程异常控制点间的距离应不小于似大地水准面格网 间距 检验点应满足GPS观测与水准联测条件 dD 在利用旧点作为检验点时,应检查旧点的稳定性、,可靠性和完好性,以及是否满足GPs观测与 水准观测,符合要求方可利用
GB/T2370g一2009 8.2检验点观测 检验点的测量精度应不低于区域似大地水准面精化时高程异常控制点的测量精度 a b)检验点的外业观测与区域似大地水准面精化时高程异常控制点的测量要求一致 8.3检验点数据处理 a GPS数据处理按照GB/T18314的要求执行; b)水准数据处理按照GB/T12897和GB/T12898的要求执行; 按公式(3)计算检验点的实测高程异常; d)利用检验点的大地坐标和拟合后似大地水准面模型计算各检验点的高程异常 似大地水准面精度评定 由似大地水准面模型计算的各检验点高程异常与其实测高程异常不符值计算的中误差,作为似大 地水准面精度 成果检查验收 区域似大地水准面精化成果应采用二级检查一级验收制 9.1 9.2GPS测量、水准测量、重力测量等成果检查验收按CH1002和CH1003执行 成果上交 10 区域似大地水准面精化工作完成后应上交成果,主要包括 技术设计书 a) b) 数据处理方案; GPs观测数据及成果， c) 水准观测数据及成果; d 高程异常控制点成果表; 区域似大地水准面模型成果， g)技术总结 精度检验报告; h 7 检查验收报告
GB/T23709一2009 附 录A 资料性附录 相邻高程异常控制点的最大间距计算 根据似大地水准面精化的精度要求,按公式(2)计算的相邻高程异常控制点的最大间距见表A.1， 供似大地水准面精化设计时参考 表A.1相邻高程异常控制点的最大间距 高程异常控制点的最大间距/lkm 地形类别 m;/m A=15 A=2.5 253 0.3 103 179 69 19 168 0,2 平地 6o 84 0.l 34 17 0.05 30 42 69 0.3 119 68 46 0,2 79 112 丘陵地 0.1 23 40 56 0,05 12 20 28 103 179 253 0,6 0,5 86 149 21m 69 119 168 0.4 山地 0.3 51 89 126 0.2 34 60 84 17 0.1 30 42 0.05 15 21 182 0,6 129 0,5 62 107 152 86 121 0.4 50 37 高山地 0.3 64 91 25 0,2 43 61 12 21 0.1 30 0.05 15 注1，";表示似大地水准面的精度 注2:入表示平均重力异常格网分辨率
GB/T2370g一2009 B 附 录 资料性附录 重力归算与格网平均重力异常计算 B.1重力观测值归算 根据我国地形和重力分辨率情况,在重力归算时一般采用地形均衡异常归算方法 空间重力异常的定义为式(B.1): B.1) g=g一Y十g 式中: 重力观测值,单位为毫伽(10ms-') -2) 椭球面上正常重力值,单位为毫伽(I0ms 空间改正,单位为毫伽(10-"ms一3)y g og=0.3086h一0.72×10'h 式中: -海拔高程,单位为米(m) 布格重力异常的定义为式(B.2) Ag8=g一丫十og十og B.2 =Ag十g 式中: 一2). 层间改正,单位为毫伽(10-"nms 0g 且 og 2rGph=-0.ll19 式中: 地心引力常数 地壳密度 地形均衡重力异常的定义为式(B.3): A以s=从一y十ogi十o6;十a4十os (B.3 AgB十g十gs 式中: -局部地形改正,单位为毫伽(10-ms一3); 0g -均衡改正,单位为毫伽(10-ms') ogs 地形及均衡改正采用四棱柱法计算,其形式为式(B.4): drdyds B.4) g= 3E -句真, .r十y十 积分后变为式(B.5). 3 y 江 og=-Gpxln(y十r)十yln(.r十)十这arctan B.5 工y y之 式(B.5)中,令1=0,;=H一H=Ah就是地形改正g
GB/T23709一2009 式(B.5)中,令1=一 -(H 十T),2=-(H 十T十4.45H),并用补偿密度A代替p,即为均衡改正 ogs,T为均衡深度 B.2格网地形及均衡改正计算 可采用传统积分法与谱方法相结合的组合法完成计算 格网地形改正的计算公式为: B.6 8g =oso+ogmn 式中: 5 ms-); -中央区部分格网地形改正,可采用四棱法计算,单位为毫伽(10 g0 积分区域内中央区以外部分格网地形改正,可采用谱方法计算,单位为毫伽(10-"ms) ogc 格网均衡改正的计算公式为 (B.7 ag 一以丽十Aa 式中: 最小计算区域内(一般取1"X1“区域)公共山根部分对均衡改正的影响(一个计算区域计 g" 算一次),用四棱柱法计算,单位为毫伽(10-ms-'); -公共山根以外区域的影响,采用谱方法计算,单位为毫伽(10-"ms) ogs1 B.3重力异常拟合计算 利用离散点的均衡重力异常值作为已知值,可采用线性移动拟合或Shepard拟合等方法计算格网 均衡异常 B.4格网平均异常的计算实施 在格网平均空间异常计算时,可按以下步骤实施: 利用数字地形模型完成格网点的层间改正、局部地形改正和均衡改正计算; a 完成离散点重力观测值的布格,均衡异常归算; b 利用离散点的均衡重力异常值作为已知值,采用线性移动拟合法计算格网点的均衡异常; 由格网点的均衡异常恢复格网点的空间异常; d 由格网点的空间异常,计算格网平均空间异常
GB/T2370g一2009 g 附 录 资料性附录 重力似大地水准面计算 C.1利用地球重力场模型计算模型似大地水准面 由地球重力场模型计算扰动位的球谐表达式为 学写>c ,cosn十Ssinm)Pcos/0) --- C.1 T(y,0,入=" 式中: -计算点的地心距离、余纬和经度; r,0,入 GM -地心引力常数与地球质量的乘积; 参考椭球长半轴，单位为米(m): e 完全规格化位系数; ,Sm 完全规格化缔合Legendre函数 P(cosO 地球重力场模型展开的最高阶数 Nm 已知位系数,根据Bruns公式来确定模型似大地水准面,即 ->之(c cosmn十Ssin)PcosO)(C.2) gM(r,0,入= 式中: -计算点正常重力 利用地球重力场模型计算模型重力异常 重力异常A与扰动位丁关系为: aT2T C.3) Ag r 将式(c.1)代人式(C.3)可得 S(e了" 一1习(Cmcosna十S,,sinna)P(co-0 g(r,0,A= (C.4 c.3由移去-恢复(remov-restore)技术计算重力似大地水准面 似大地水准面的计算公式 一(Gooh”十K(c.5 = -TGoh Js(e(-Km+T)da-一AKA 元 式中: s(p) 6s" 一4十10s一3(1一2)ln(s十 sin(e/2 h'=0.453一0.018sin(B)十0.087cos(B)cos(L)十0.204cos(B)sin(L N=%一0.3086 TC -地形改正,单位为毫伽(10-ms-=); -计算点高程,单位为米(m); 10
GB/T23709一2009 sms-2 椭球面上的正常重力,单位为毫伽(10 地球的平均曲率半径,单位为米(m); 地球平均正常重力,单位为毫伽(10-" ms; N 剩余空间异常(实际值与模型值的差值),单位为毫伽(10"ms'); Ag da 单位球面a上的面元; 重力异常,单位为毫伽(10-ms-2) Ag G 地心引力常数 地球平均密度 利用重力场模型计算的模型似大地水准面; s(p Stokes函数 -球面角距; B 大地纬度 大地经度 1l
GB/T2370g一2009 附 录 D 资料性附录 重力似大地水准面与高程异常控制网的多项式与多面函数拟合 完成与似大地水准面分辨率一致的规则格网重力似大地水准面计算后,需要利用高程异常控制网 成果作为控制,对重力似大地水准面进行拟合,获得适配于该区域与国家正常高系统定义一致的拟合似 大地水准面 拟合方法如下 D.1高程异常控制点重力似大地水准计算 在完成规则格网重力似大地水准面计算后,可采用Shepard插值法,双线性内插等方法内插各高程 异常控制点的重力似大地水准面 D.2拟合模型 重力似大地水准面与高程异常控制网的拟合时,可根据实际情况选择多项式,多面函数等拟合 方法 三次多项式数学模型见式(D.1)和式(D.2) a )?十 A=a十a9一m -ee(A一入m十a.(9一fm)十a(9一pm m)(入一入m十a.(A一入m )'(入 )(入一入)?十a(入一入)D.1 p- 十a7(4 anp Pm 入m十as9 Pnm D.2 Ag=Swis 6 式中 为多项式拟合系数; aa,ai,as 实测高程异常与重力似大地水准面之差,单位为米(m); A 分别为实测高程异常和重力似大地水准面高,单位为米(m) ns、s 分别为拟合区的中心纬度和经度 Om,人 多面函数方程的一般形式为式(D.3) b D.3 ya,g(.ry,riy, 式中 待定系数; a -r,y的二次函数 gr,y，.Z，,y D.3拟合方案 重力似大地水准面与高程异常控制网的拟合时,一般可采用以下三种拟合方案 全网一次性拟合,求解一套全区域拟合系数 该方案适用于全区高程异常控制点分布均匀,重力似大地水准面与高程异常控制网分辨率和 精度大致一致,且全区范围不是过大(如一般城市). 将全网适当分成若干个子区域,分区求解各区拟合系数 b 该方案适用于全区高程异常控制点的分布有地区性差异,重力似大地水准面与高程异常控制 网的分辨率和精度都有明显的分区特点,且全区范围很大,全区用一个拟合多项式难于描述 这种区域性差异 当采用分区拟合时,对于接边问题的处理,采取对每个分区向邻区扩大一个重叠宽度,所求解 的拟合系数规定只属于本区,这样可以避免或减弱边界两边拟合系数的突变现象,形成平滑 过渡 12
GB/T23709一2009 为避免分区拟合出现的相邻区的接边问题,以一个点为中心,确定一个适当的球冠或方块形拟 合区域,其中分布有必要数量的高程异常控制点,将求解的拟合系数定义为中心点的转换 参数 该方案适用于中等大小区域的拟合 D.4拟合系数求解 由高程异常控制点上的重力似大地水准面高与实测高程异常的差值,组成不符值序列; a b)由不符值序列和相应高程异常控制点的坐标组成多项式拟合“观测方程”,其中未知参数为多 项式系数; 按最小二乘原理求解拟合多项式系数 D.5推算拟合后的似大地水准面格网值 由拟合多项式系数和格网中心点坐标,对重力似大地水准面进行拟合纠正,获得与国家正常高系统 定义一致的格网似大地水准面模型