第25第2别 海洋测绘 VoL 26.No.2 2005年1月 HYDROGRAMIC SRVEYNG AND CIARTNG Mar.,2006 论导线测角粗差的产生及其影响 然炳节”石忠科2，夏双喜 (1,生安理工大学水比学究，碳出安101长2齿安市用河施市段管理中心，碳齿齿安？10016) 摘要：针对导线测量中出见的测角租整对粗差的产生原因.有无粗差的坐标关系进行了分析对粗差判据方 法注行了里论推证，实测验还和时论，证实了墙方法具有星论的完整性和实后的可操作性，并对减夕粗差出现的措 提出了速议， 关续词：测量T芝：测角相某，前推计算：后进计京 中图分英号：四11”.4 文献标1识码：A文章第号·167-3044(236)0-00？-04 1前言 位角闭合楚严重超限。对此，有的测量小组是白目 地全线运工：有的小组则速反作业规定仪凭回忆来 用染的存在对测量平差结果会产生硅坏性的影 判断可能出观的帮熟站。 南。关于导线水平角相差站的校出方法，目前有 图解法和公式法”。利用检出方法有效而准确地 2有无测角粗差时点位坐标关系 检出外业测角粗差产生的测站，对指孝外业生产和 图1为一用合导线，A、R(=角、C(=路，D 提高上作效中有重要意文。如，测完一条导线后发 分别为已知点，A、乃、…P、…P分别为新设导线 观方位角闭合差高达几分以上，经校在手薄记录， 点，月，月、月、…，且，、…且。、B分别为前进方间 算均无误，测证明导线有测角用熟存在。若不用科 实就头方向的观测左角，1〔=0、，、兔、… 学而严密的方法去玉别相来站，育目去个线工，会 x,…怎~飞1(=分别为礼剃的水半组离。若按 对工怡效幸和队员积极性产生负面影响，本文针树 前进方向岸算御支导线计算任意点P的坐标为 等线测量水平角相差产生的原因，各点有无粗差时 真伪坐标的数学关系，以及减少粗差出现几率的措 旄进衍讨论. 1 出于测绘外业生产不境恶劣、劳动强度人以及 某些不完善的作业制度等深因，使得测绘人员的工 作周量难以一流稳定在良好状药，工作中较易出现 2/ 用染。一般而言，产生相差的主要原因有：①度盘 式中，T为己知方位 ,=01,2、w中1e 实际配置位置与本该配置位置不相符，而按本该配 置的位置读数的：②观测读数特有祖差的。包括： 印象读数，如将盘右100110紫成18000'10等； 读激口齿不清，如4与0,6与0,7与1的发音含混 等：#位读数，如将69读成9613引读成113等： ③记薄奶所错和记错的：印象记簿，如将盘右180° 0110记成1800010等：④仪墨不经意张动，气泡 图1附合导线的N测图 严重复斜而未发现的：⑤分次测完一条导线时，衔 若按后退方向唐简头方向性算，则·坐闲 按边两嘴点中最少有一点两次位置不一数的；圆定 加上标以区别于前桩坐标为： 向方向的目标点不在正确位置上，或测站上带点等， 无论种情况出观的狂整都不容忽视，会使方 收璃日期：2005-11-01：修▣日期：200602-10 基金度日：陕西省自船科学基发项口003E器). 作者篱介：魏销乾(3，男峡西兴平人对上出数援主妻从事河流动力学发定沙湖服叶究 994-2009 Clm Acadenuc Jownal Elecmrowe Pablsling Houre.Ail nighes rererred hup:/wcnft.ner
第 26卷第 2期 2006年 3月 海 洋 测 绘 HYDROGRAPH IC SURVEYING AND CHARTING Vol126,No12 Mar1, 2006 收稿日期 : 2005211201; 修回日期 : 2006202210 基金项目 : 陕西省自然科学基金项目 (2003E208)。 作者简介 : 魏炳乾 (19632) ,男 ,陕西兴平人 ,博士 ,副教授 ,主要从事河流动力学及泥沙测报研究。 论导线测角粗差的产生及其影响 魏炳乾 1 ,石忠科 2 ,夏双喜 1 (11西安理工大学 水电学院 ,陕西 西安 710048; 21西安市渭 河城市段管理中心 ,陕西 西安 710016) 摘要 : 针对导线测量中出现的测角粗差 ,对粗差的产生原因、有无粗差的坐标关系进行了分析 ,对粗差判据方 法进行了理论推证、实例验证和讨论 ,证实了该方法具有理论的完整性和实际的可操作性 ,并对减少粗差出现的措 施提出了建议。 关键词 : 测量平差 ;测角粗差 ;前推计算 ;后推计算 中图分类号 : P211 + 14 文献标识码 : A 文章编号 : 167123044 (2006) 0220009204 1 前 言 粗差的存在对测量平差结果会产生破坏性的影 响 [ 1 ]。关于导线水平角粗差站的检出方法 ,目前有 图解法和公式法 [ 2～3 ]。利用检出方法有效而准确地 检出外业测角粗差产生的测站 ,对指导外业生产和 提高工作效率有重要意义。如 ,测完一条导线后发 现方位角闭合差高达几分以上 ,经检查手簿记录、计 算均无误 ,则证明导线有测角粗差存在。若不用科 学而严密的方法去甄别粗差站 ,盲目去全线返工 ,会 对工作效率和队员积极性产生负面影响。本文针对 导线测量水平角粗差产生的原因、各点有无粗差时 真伪坐标的数学关系、以及减少粗差出现几率的措 施进行讨论。 由于测绘外业生产环境恶劣、劳动强度大以及 某些不完善的作业制度等原因 ,使得测绘人员的工 作质量难以一直稳定在良好状态 ,工作中较易出现 粗差。一般而言 ,产生粗差的主要原因有 : ①度盘 实际配置位置与本该配置位置不相符 ,而按本该配 置的位置读数的 ; ②观测读数带有粗差的。包括 : 印象读数 ,如将盘右 180°01′10″读成 180°00′10″等 ; 读数口齿不清 ,如 4与 10, 6与 0, 7与 1的发音含混 等 ;错位读数 ,如将 69 读成 96, 131 读成 113 等 ; ③记簿员听错和记错的 ;印象记簿 ,如将盘右 180° 01′10″记成 180°00′10″等 ; ④仪器不经意碰动、气泡 严重倾斜而未发现的 ; ⑤分次测完一条导线时 ,衔 接边两端点中最少有一点两次位置不一致的 ; ⑥定 向方向的目标点不在正确位置上 ,或测站上错点等。 无论哪种情况出现的粗差都不容忽视 ,会使方 位角闭合差严重超限。对此 ,有的测量小组是盲目 地全线返工 ;有的小组则违反作业规定仅凭回忆来 判断可能出现的粗差站。 2 有无测角粗差时点位坐标关系 图 1为一附合导线 , A、B ( = p0 ) 、C ( = pn + 1 )、D 分别为已知点 , P1、P2、…、Pi、…、Pn分别为新设导线 点 , β0、β1、β2、…、βi、…、βn、βn + 1分别为前进方向 (实箭头方向 )的观测左角 , s- 1 ( = 0)、s0、s1、s2、…、 si、…、sn、sn + 1 ( = 0)分别为观测的水平距离。若按 前进方向推算 (即支导线计算 )任意点 Pi的坐标为 : xpi = xB + ∑ i j=0 sj- 1 cos( TAB + ∑ j k =1 (βk - 1 + 180°) ) (1) ypi = yB + ∑ i j =0 sj- 1 sin ( TAB + ∑ j k =1 (βk - 1 + 180°) ) (2) 式中 , TAB为已知方位角 , I = 0、1、2、…n + 1。 图 1 附合导线的观测图 若按后退方向 (虚箭头方向 )推算 , 则 Pi坐标 (加上标 ′以区别于前推坐标 )为 : x′pi = xC + ∑ i j =n +1 sj cos( TDC + ∑ j k =n (180°- βk +1 ) ) (3) ' 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
海洋测始 第6卷 顾及2式，同理可得 是+￡5油+E1”,我》周 式中，为D边的已如方位角. 在不考需器深染条件下。·点的前推与后推 坐标应相蒂，=。=成： 50s7+月：+1809》 112 为+1m+用+180- 后推时，粗差△将对点到B点间所有边的方位角 元+m+10.》 以及所有点的标产生彩有：而对C到向所有边 的方位角及所有点的坐标无彬响，同于前推的方 很设在任意站有测角粗差楼入，即该站观 法即有 测角值为：，+4，那么.的推时X,Y表示行粗差。 当n+1时， 排入后所推坐标).相差凸将对兵到C点间所有边 的方位州以及所有点的标俏产生彩响：而对B点 到间两有边的方位角及所有点的坐标无彰响。即 13) 可写成 当0w制 14) 当m-1即A与Af合到. 7 =车+3mfx+E10.R》 15) 8 当m=,即4与亚介时 16) 当>m0时， 9) 100 当<m≤n+1时，顾及1式.并经配项，和 整化积等量学变换，可得： /17 o) -2 ∑1sfTu- 2 花 18) 者在有测尚拥整△，出7)和17)式、8影和 (18/式可得Aw<点前指与后推的坐标关系为 110 代中的为无相整狗伯： J994-30019 Chuoa Acmdrooe looranl Electrme Puhfisivne ffuwre.AU ngivr rvorwd hn/nn cul w
海 洋 测 绘 第 26卷 y′pi = yC + ∑ i j=n+1 sj sin ( TDC + ∑ j k =n (180°- βk +1 ) ) (4) 式中 , TDC为 DC边的已知方位角。 在不考虑偶然误差条件下 , Pi点的前推与后推 坐标应相等 xpi = xp′i , ypi = yp′i ,或 : xB + ∑ i j =0 sj- 1 cos( TAB + ∑ j k =1 (βk - 1 + 180°) ) = xC + ∑ i j=n +1 sj cos( TDC + ∑ j k = n (180°- βk +1 ) ) (5) yB + ∑ i j =0 sj- 1 sin ( TAB + ∑ j k =1 (βk - 1 + 180°) ) = yC + ∑ i j =n +1 sj sin ( TDC + ∑ j k = n (180°- βk +1 ) ) (6) 假设在任意站 pi有测角粗差 Δ掺入 ,即该站观 测角值为 :βi +Δ,那么 ,前推时 (X、Y表示有粗差 Δ 掺入后所推坐标 ) ,粗差 Δ将对 pi到 C点间所有边 的方位角以及所有点的坐标值产生影响;而对 B 点 到 pi间所有边的方位角及所有点的坐标无影响。即 可写成 : 当 0≤m i时 , X p′m = xC + ∑ m j =n +1 sj cos( TDC + ∑ j k =n (180°- βk +1 ) ) (13) Yp′m = yC + ∑ m j= n +1 sj sin ( TDC + ∑ j k =n (180°- βk +1 ) ) (14) 当 m = i,即 pm与 pi重合时 , Xp′m = xC + ∑ i j= n +1 sj cos( TDC + ∑ j k = n (180°- βk +1 ) ) (15) Yp′m = yC + ∑ i j= n +1 sj sin ( TDC + ∑ j k =n (180°- βk +1 ) ) (16) 当 i >m ≥0时 , X′pm = xC + ∑ m j=n+1 sj cos( TDC + ∑ j k =n (180°- βk+1 ) ) - 2 sin Δ 2 ( sin Δ 2 ∑ m j= i- 1 sj cos( Ti+1→i + ∑ j k = i- 1 (180°- βk+1 ) ) - cos Δ 2 ∑ m j= i- 1 sj sin ( Ti+1→i + ∑ j k = i- 1 (180°- βk+1 ) ) ) (17) Y′pm = yC + ∑ m j=n+1 sjsin ( TDC + ∑ j k =n (180°- βk+1 ) ) - 2 sin Δ 2 ( sin Δ 2 ∑ m j= i- 1 sj sin ( Ti+1→i + ∑ j k = i- 1 (180°- βk+1 ) ) + cos Δ 2 ∑ m j= i- 1 sj cos( Ti+1→i + ∑ j k = i- 1 (180°- βk+1 ) ) ) (18) 若在 pi有测角粗差 Δ,由 ( 7)和 ( 17)式、( 8)和 (18)式可得 pm (m < i)点前推与后推的坐标关系为 (式中的 β为无粗差角值 ) : 10 ' 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
第2用 纯炳花，等论导线测角祖差的产生及其影响 X-.2s号s,go 3应用与讨论 以上推宁是杆不考虑倘然狱差的假定下进行 的.面实际测量的观测误差不可避免，所以判断时， 各点排后推坐标相差最小的测的即为粗兰站。 图2为文款3的算例结果.图中点名不带的位置 19) 就是无粗差时的正确位置。图中实线和虚线分别为 成推和后推的计算路线。 过方 ●7 (20/ 同理.由11和13式、12和14式可得A 侧之点的前带与听推的坐标关系为 ann 种岁1对 s(.- 图2导线计京翰图 21 由图2可知 28 sin(T. 1)成推与后推的坐标结果在D2站处基本 重合，也就是说该站有测角用差存在： +成+180四s 2,前在时Dk2站之后的所有测站在Dx2处 整体向右旋转&而偏离真拉置：后准时D2站以后 的所有站卿DXI和未利在D2处都向左沿 4m (22 后推方向霜旋转△面偏商真位置： 3/无粗差制站X,了坐标的偏离量正是式 以上四式中点，m工-+系，TB+ 19一22/右边的第二实值： 180和三，0亿-，+点不月+1809) 4/在不计码然误差时，D2站的前推与后推 分别为到m站的x的坐标增量，式21、22) 坐标敌.粗差角值可山任一点的前，斤排坐标与粗 与式19、20同样，由于等号右边第2项括号中 整点坐标推求： 的何一量的符号不问，使得拼号内的值不可能同时 5,当粗差结在已知点时.粗望A可以推求， 为零风有在x、的坐标招量都为零时才为零，暴然 不任己知点时.由于角闭合差是粗差d和全学线测 这与1中w矛盾)所以A点的前推点位x.-) 角偶然误差W之和，若按一定的条件仰点位闭合 与后推点位，是不量合的：山当围一时，根 差最小或角闭合差最小求A时，木必能反陕真实 据式5)、6.由9)和13式，10/和16式可 而且这包与规范爱求相饽， 得点的前推与后推坐称才相等，即 为涯免相差日现，测骨时应根据理测员与记簿 0=Xg= 23/ 员的言语况，个体情况等钊定对策。如：对阿拉伯 上式表明，不计偶然保差时，即使在点有测角 数字的读数.为程免1与7,4与10,6与0的祖斯混 粗左凸，n点的的推与后推坐标仍是相等的。这藏 记可规定区分明确的代字：读数时禁止中象读要： 是解析法判断粗差站点的理论依据。 记簿员责任心买强：不能在读爱前记录，记录时若感 199-4-2009 (h Aemueowe Joornal Elecfrenie Pubfistvng ffoscie.AN nglr renerwd.hp /.cul m
第 2期 魏炳乾 ,等 论导线测角粗差的产生及其影响 Xp′m = Xpm - 2 sin Δ 2 ( sin Δ 2 ∑ m j= i- 1 sj cos( Ti+1→i + ∑ j k = i- 1 (180°- βk +1 ) ) - cos Δ 2 ∑ m j= i- 1 sj sin ( Ti+1→i + ∑ j k = i- 1 (180°- βk +1 ) ) ) (19) Yp′m = Ypm - 2 sin Δ 2 ( sin Δ 2 ∑ m j= i- 1 sj sin ( Ti+1→i + ∑ j k = i- 1 (180°- βk +1 ) ) + cos Δ 2 ∑ m j= i- 1 sj cos( Ti+1→i + ∑ j k = i- 1 (180°- βk +1 ) ) ) (20) 同理 ,由 (11)和 (13)式、(12)和 ( 14)式可得 pm (m > i)点的前推与后推的坐标关系为 : Xpm = X p′m - 2 sin Δ 2 ( sin Δ 2 ∑ m j = i+1 sj- 1 cos( Ti- 1→i + ∑ j k = i+1 (βk - 1 + 180°) ) + cos Δ 2 ∑ m j = i+1 sj- 1 sin ( Ti- 1→i + ∑ j k = i+1 (βk - 1 + 180°) ) ) (21) Ypm = Yp′m - 2 sin Δ 2 ( sin Δ 2 ∑ m j = i+1 sj- 1 sin ( Ti- 1→i + ∑ j k = i+1 (βk - 1 + 180°) ) - cos Δ 2 ∑ m j = i+1 sj- 1 cos( Ti- 1→i + ∑ j k = i+1 (βk - 1 + 180°) ) ) (22) 以上四式中 , ∑ m j = iº 1 sj j- 1 sin ( Ti ±1→i + ∑ j k = iº 1 ( º βk ±1 + 180°) )和 ∑ m j = iº 1 sj j- 1 cos( Ti ±1→i + ∑ j k = iº 1 ( ºβk ±1 + 180°) ) 分别为 i到 m 站的 x、y的坐标增量 ,式 ( 21)、( 22) 与式 (19)、(20)同样 ,由于等号右边第 2项括号中 的同一量的符号不同 ,使得括号内的值不可能同时 为零 (只有在 x、y的坐标增量都为零时才为零 ,显然 这与 i≠m 矛盾 ) ,所以 pm点的前推点位 (Xpm , Ypm ) 与后推点位 (Xpm′, Ypm′)是不重合的;但当 m = i时 ,根 据式 (5) 、(6) ,由 ( 9)和 ( 13)式、( 10)和 ( 16)式可 得 pi点的前推与后推坐标才相等 ,即 : X p′i = Xpi , Yp′i = Ypi (23) 上式表明 ,不计偶然误差时 ,即使在 pi点有测角 粗差 Δ, pi点的前推与后推坐标仍是相等的。这就 是解析法判断粗差站点的理论依据。 3 应用与讨论 以上推导是在不考虑偶然误差的假定下进行 的 ,而实际测量的观测误差不可避免 ,所以判断时 , 各点前推、后推坐标相差最小的测站即为粗差站。 图 2为文献 [ 3 ]的算例结果 ,图中点名不带 ′的位置 就是无粗差时的正确位置。图中实线和虚线分别为 前推和后推的计算路线。 图 2 导线计算简图 由图 2可知 : (1) 前推与后推的坐标结果在 DX2站处基本 重合 ,也就是说该站有测角粗差存在; (2) 前推时 DX2站之后的所有测站在 DX2处 整体向右旋转 Δ而偏离真位置;后推时 DX2站以后 的所有测站 (即 DX1′和宋村 ′)在 DX2处都向左 (沿 后推方向看 )旋转 Δ而偏离真位置; (3) 无粗差测站 X、Y 坐标的偏离量正是式 (19) ～(22)右边的第二项值; (4) 在不计偶然误差时 , DX2站的前推与后推 坐标一致 ,粗差角值可由任一点的前、后推坐标与粗 差点坐标推求; (5) 当粗差站在已知点时 ,粗差 Δ可以推求。 不在已知点时 ,由于角闭合差是粗差 Δ和全导线测 角偶然误差 w 之和 ,若按一定的条件 (如点位闭合 差最小或角闭合差最小 )求Δ时 ,未必能反映真实 , 而且这也与规范要求相悖。 为避免粗差出现 ,测量时应根据观测员与记簿 员的言语情况、个体情况等制定对策。如 :对阿拉伯 数字的读数 ,为避免 1与 7、4与 10、6与 0的混听混 记 ,可规定区分明确的代字;读数时禁止印象读数; 记簿员责任心要强;不能在读数前记录;记录时若感 11 ' 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
海洋测绘 第26拾 觉异常，要及利雍认：观测员测目问必须检有气泡： 对各宁戏点的世置，礼测色须做到心中有数，随时 检查点的位置等。 参考文献： 11川就★效.用的定性分析1卫！测染后具与上和， 4站京语 2002.27(1):13 面过理诊推导实除证和讨论可知：当导线方 [2】王成牌周治邮，商充导线测量中角度相整点的两 位角同合差严重短限仅在某测站有测角相差时， 件定位方法叠式注及CD法1」1携矿测拾，201， 可按支导线计算各点的前推和后推坐标，若两套坐 12):22-15 3】等坏德，向帝君.导投署量水平角粗左站的刺断A1 标伯扫势最小者就是相整利站：为了减少相整的出 重庆市测给学会202年学术论义集【C】域花：内南 现几率，作业小组应制定详尽的工作钊度，以杜绝粗 文道大学出饭社，001 差的出现。 Study on the Mistakes of llorintal Angle in A Connecting Traversing and Its Influences WEIBingqian SHI Zhong-ke XM Shuangxi (L nstiute of W ater Kexres and Hydmekctne Engineering, Xi'an Unnersiry of 'Technolgy,XI'an,Shaanxi,710045. 2 Xi'an Uiban Segnent Adm nisirton Cenkr ofWerChan River,Xi'an,Shamnxi,710016) Abstract n thispaper.the nean of horiontal angk mistake uccurring in a coenocting tranvens is discussod firsth Th:relationshp of cooddiates beteen mard and backard direction cakulton is derived ad analyzed Furhemore the utiliry of the method is mestigated usng te reauhs of an exaple Fially e suggestions n moid mistkes are pmposed Key woads adjusment ofmeasunment.anglar guss enor,and eakulaton,lskwand eakulatin 上换第8页) Rescarch on App mxmation Method on Marine Magnctic Anomaly JN Shao-hua',YU Bo,LU Yan-chun',ZHA I Guo-jun',B AN Giang (1.Deparment of Hydpgrphy and Carpgmhpby,Dalan Nmal Acadany,Dalin,L soning 116018, 2 Naval hstitk of Hydmgphk Surveyig and Chartng Tanjin.300061) Abstmet Wim the analysis aed study of genenl numencal app oomaton methods according the chanekristie of he=arine magelic suney,nn电curved s:e当ead,Hayo小hedal fancton metd Shupand mefod and kriging mefod ae aulated and anabad in cae uf differen!intepoblon accunky Simulerously.the inkrpolaton accuncy of he four knds of app uximation methods ax cakubkd as example The resalts of the anulation and example show that the intepolton accuracy of these hods are different wih distrbutng o对wn sarvey point and change ofmageti aoma水For d terent nstaces白is paper ammar oomepondent methods suitble fr marine mngnetic ananal appmomaton ey wond×mnc mayelic amaly,津pxmi道，neny 'y429 Ch Acuee loarf Eleemo企'isge2.是gr6 Amp.ww.c点g
海 洋 测 绘 第 26卷 觉异常 ,要及时确认;观测员测回间必须检查气泡; 对各导线点的位置 ,观测员必须做到心中有数 ,随时 检查点的位置等。 4 结束语 通过理论推导、实例验证和讨论可知 :当导线方 位角闭合差严重超限 ,仅在某一测站有测角粗差时 , 可按支导线计算各点的前推和后推坐标 ,若两套坐 标值相差最小者就是粗差测站;为了减少粗差的出 现几率 ,作业小组应制定详尽的工作制度 ,以杜绝粗 差的出现。 参考文献 : [ 1 ] 姚宜斌. 粗差的定性分析 [ J ]. 测绘信息与工程 , 2002, 27 (1) : 1～3. [ 2 ] 王政梅 ,周治雄 ,高 亮. 导线测量中角度粗差点的两 种定位方法 —公式法及 CAD法 [J ]. 地矿测绘 , 2001, (2) : 22～25. [ 3 ] 魏炳乾 ,向泽君. 导线测量水平角粗差站的判断 [A ]. 重庆市测绘学会 2002年学术论文集 [ C ]. 成都 :西南 交通大学出版社 , 2003. Study on the M istakes of Horizontal Angle in A Connecting Traversing and Its Influences WE IBing2qian 1 , SH I Zhong2ke 2 , X IA Shuang2xi 1 (1. Institute ofWater Resources and Hydro2electric Engineering, Xi’an University of Technology, Xi’an, Shaanxi, 710048; 2. Xi’an U rban Segment Adm inistration Center ofWei2Chan River, Xi’an, Shaanxi, 710016) Ab stra c t: In this paper, the reason of horizontal angle m istake occurring in a connecting traverse is discussed firstly. The relationship of coordinates between forward and backward direction calculation is derived and analyzed. Furthermore the utility of the method is investigated using the results of an examp le. Finally some suggestions to avoid mistakes are p roposed. Key wo rd s: adjustment of measurement; angular gross error; forward calculation; backward calculation (上接第 8页 ) Research on App roximation Method on Marine Magnetic Anomaly J IN Shao2hua 1 , YU Bo 1 ,L IU Yan2chun 1 , ZHA I Guo2jun 2 ,B IAN Gang 1 (1. Department of Hydrography and Cartograhphy , Dalian Naval Academy , Dalian,Liaoning, 116018; 2. Naval Institute of Hydrographic Surveying and Charting , Tianjin , 300061) Ab stra c t: W ith the analysis and study of general numerical app roximation methods, according to the characteristic of the marine magnetic survey, moving curved surface method, Hardy polyhedral function method, Shepard method and kriging method are emulated and analyzed in case of different interpolation accuracy. Simultaneously, the interpolation accuracy of the four kinds of app roximation methods are calculated as examp le. The results of the emulation and examp le show that the interpolation accuracy of these methods are different with distributing of known survey point and change of magnetic anomaly. For different instances, this paper summarizes correspondent methods suitable for marine magnetic anomaly app roximation. Key wo rd s: marine magnetic anomaly; app roximation; interpolation accuracy 12 ' 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net