区域大地水准面精化理论与应用
|
摘要
随着GPS技术的发展，由GPS定位获取的大地高直接解求厘米级精度的正高或正常高已成为一种趋势，其实现可以减轻外业工作量，大大提高工作效率，从而带来直接的经济效益。而实现这一目标的前提条件是确定一个达到足够精度的(似)大地水准面。确定厘米级的区域(似)大地水准面是当今大地测量学研究的热点，它将为基础测绘、数字中国地理空间基础框架、区域沉降监测、地学研究、交通、水利、电力等多学科研究．与应用提供必要的测绘服务。
本文首先介绍了大地水准面精化的意义和国内外大地水准面精化发展动态，系统地研究了精化区域大地水准面的原理和各种方法，并分别对其中的重力法和几何法进行了详细的讨论和研究。在重力法中，研究了重力法精化大地水准面的理论和主要技术、(似)大地水准面的计算模型、区域厘米级大地水准面精化的实施方法。在几何法中，介绍了天文水准、卫星测高和GPS／水准三种方法，并着重研究了GPS／水准拟合法：常数拟合、曲线拟合(多项式曲线拟合、三次样条曲线拟合和Akima法曲线拟合)和曲面拟合(多项式曲面拟合、移动曲面拟合、多面函数拟合和样条函数拟合)。此外，还研究了地球重力场模型和拟合相结合的方法。结合测区GPS／水准资料，采用多项式曲面拟合模型、移动曲面拟合模型和多面函数模型，用VB自编程序做了拟合试验和精度分析，得出相关结论，取得良好效果，拟合的区域似大地水准面可应用于GPS大地高与正常高之间的转换。之后，介绍了地球重力场模型理论，研究了如何利用地球重力场模型和GPS／水准数据，采用移去～恢复方法计算区域似大地水准面，进而求解正常高的方法，并结合实例将其与常规拟合的结果进行比较，对其在生产应用中的可行性进行了分析。
最后，分析了大地水准面精化的误差来源，并总结了不同情况下区域大地水准面精化的方法和改进办法。
With the recent advance in GPS technique, direct transformation from GPS-derived ellipsoidal heights to orthometric heights or normal heights in centimeter level becomes inevitable. Its implement can produce direct economic benefits by tremendously reducing laborious field tasks so as to increase efficiency. However, the precondition of realizing this goal is to determinate an accurate geoid or quasi-geoid. Currently, regional geoid or quasi-geoid determination with centimeter level is the focus of geodesy research, it will provide mufti-disciplinary research and application such as base map up date, digital China geo-spatial framework, regional monitoring on settlement, earth science research, traffic, water conservancy and electric power with essential mapping service.
To begin with, the article analyses the importance of refining regional geoid and its developments in China and abroad. Then it studies systematically the basic theories and principles of refining regional geoid and digs into the methods of gravity and geometer respectively. In the gravity method, it describes the theories and techniques in determining a precise regional geoid, the computing model of geoid, means of determining regional centimeter rank accurate geoid. In the geometrical method, it introduces three methods, astronomy leveling, satellite altimetry and GPS/leveling. It especially describes GPS/leveling fitting methods: figure fitting, curve fitting (multinomial curve fitting, cubic sp line curve fitting and Akima curve fitting) and curved surface fitting (multinomial curved surface fitting, moving curved surface fitting, Hardy function interpolation fitting and SP function fitting). Besides, it also describes the method of integrating fitting and global geopotential model. Combining measuring region GPS/leveling data thesis adopts the multinomial surface fitting model and move curve surface fitting model and Hardy function interpolation model, author has developed a VB program to do the test and analyze its accuracy and draws several satisfied conclusions. The fitting regional quasi-gieod can be applied into the shift between GPS height and normal height. Then, it introduces the theory of global geopotential model and studies how to determine regional quasi-geoid with remove～restore method, GPS/leveling data and global geopotential model, gaining the normal height at last. Based on the comparison, its feasibility is also analyzed.
At last, it analyzes the error resources of regional geoid refining and summarizes how to choose a proper method in different conditions.
本文首先介绍了大地水准面精化的意义和国内外大地水准面精化发展动态，系统地研究了精化区域大地水准面的原理和各种方法，并分别对其中的重力法和几何法进行了详细的讨论和研究。在重力法中，研究了重力法精化大地水准面的理论和主要技术、(似)大地水准面的计算模型、区域厘米级大地水准面精化的实施方法。在几何法中，介绍了天文水准、卫星测高和GPS／水准三种方法，并着重研究了GPS／水准拟合法：常数拟合、曲线拟合(多项式曲线拟合、三次样条曲线拟合和Akima法曲线拟合)和曲面拟合(多项式曲面拟合、移动曲面拟合、多面函数拟合和样条函数拟合)。此外，还研究了地球重力场模型和拟合相结合的方法。结合测区GPS／水准资料，采用多项式曲面拟合模型、移动曲面拟合模型和多面函数模型，用VB自编程序做了拟合试验和精度分析，得出相关结论，取得良好效果，拟合的区域似大地水准面可应用于GPS大地高与正常高之间的转换。之后，介绍了地球重力场模型理论，研究了如何利用地球重力场模型和GPS／水准数据，采用移去～恢复方法计算区域似大地水准面，进而求解正常高的方法，并结合实例将其与常规拟合的结果进行比较，对其在生产应用中的可行性进行了分析。
最后，分析了大地水准面精化的误差来源，并总结了不同情况下区域大地水准面精化的方法和改进办法。
With the recent advance in GPS technique, direct transformation from GPS-derived ellipsoidal heights to orthometric heights or normal heights in centimeter level becomes inevitable. Its implement can produce direct economic benefits by tremendously reducing laborious field tasks so as to increase efficiency. However, the precondition of realizing this goal is to determinate an accurate geoid or quasi-geoid. Currently, regional geoid or quasi-geoid determination with centimeter level is the focus of geodesy research, it will provide mufti-disciplinary research and application such as base map up date, digital China geo-spatial framework, regional monitoring on settlement, earth science research, traffic, water conservancy and electric power with essential mapping service.
To begin with, the article analyses the importance of refining regional geoid and its developments in China and abroad. Then it studies systematically the basic theories and principles of refining regional geoid and digs into the methods of gravity and geometer respectively. In the gravity method, it describes the theories and techniques in determining a precise regional geoid, the computing model of geoid, means of determining regional centimeter rank accurate geoid. In the geometrical method, it introduces three methods, astronomy leveling, satellite altimetry and GPS/leveling. It especially describes GPS/leveling fitting methods: figure fitting, curve fitting (multinomial curve fitting, cubic sp line curve fitting and Akima curve fitting) and curved surface fitting (multinomial curved surface fitting, moving curved surface fitting, Hardy function interpolation fitting and SP function fitting). Besides, it also describes the method of integrating fitting and global geopotential model. Combining measuring region GPS/leveling data thesis adopts the multinomial surface fitting model and move curve surface fitting model and Hardy function interpolation model, author has developed a VB program to do the test and analyze its accuracy and draws several satisfied conclusions. The fitting regional quasi-gieod can be applied into the shift between GPS height and normal height. Then, it introduces the theory of global geopotential model and studies how to determine regional quasi-geoid with remove～restore method, GPS/leveling data and global geopotential model, gaining the normal height at last. Based on the comparison, its feasibility is also analyzed.
At last, it analyzes the error resources of regional geoid refining and summarizes how to choose a proper method in different conditions.
引文
[1] 师芸，区域大地水准面的精化方法研究，[硕士论文]，西安科技大学，2004，1～2
[2] Denker H, Behrend D, Torge W, The European Gravimetfic Quasi-geoid EGG95[J], 1996, IAG Bulletin of D'information 77, IGES Bulletin No. 4(Special Issue), 3～11
[3] Zur Verwendung, Europaischen Gravimetrischen Quasigeoids EGG97 in Deutschland[J], Manuscripat Geodaetica, 1999, 18(5), 154～166
[4] 李建成，陈俊勇，宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，6～9
[5] 宁津生，罗志才，李建成，我国省市级太地水准面精化的现状及技术模式[J]，大地测量与地球动力学，2004，Vol．24，No．1，5～7
[6] 邬建耀，关于福建地区GPS网的建立及大地水准面的精化[J]，福建建设科技，2001(3)，34～36
[7] 陈俊勇，李建成，全国及部分省市地区高精度高分辨率似大地水准面的研究和实施[J]，测绘通报，2005(5)，1～4
[8] 李骏元，程传录，郭春喜，大地水准面确定的几种方法和比较[J]，测绘技术装备，2005，Vol．7，No．2，15～16
[9] 雷晓霞，基于重力与GPS水准组合法的大地水准面精化研究，[硕士论文]，长安大学，2005，7～8
[10] Weikko A Heiskanen & Helmut Moritz，Physical Geodey，物理大地测量学(中译本)，北京，测绘出版社，1967，137～138
[11] 管泽霖，地球形状及外部重力场，北京，测绘出版社，1980，3：55～39
[12] Ardalan A A. Grafarend E, Ellipsoidal Geoidal Undulations (Ellipsoidal Bruns Formula) [J], Journal of Geodesy, 2001, 75(9/10), 544～552
[13] 李建成，陈俊勇，宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，17～18
[14] 管泽霖，宁津生，地球形状及外部重力场，北京，测绘出版社，1981，57～62
[15] Bjerhammar A, Discrete physical geodesy, Dept. of Geodetic Science and Surveying, 1987, The Ohio State University, Columbus, 125～129
[16] 许厚泽，陆仲元，中国重力场与大地水准面，解放军出版社，1997，63～67
[17] Brovar, On the Solutions of Monodensky's Boundary Value problom[J], Bull. Geoid, 1964(72), 112～115
[18] 马新莹，局部大地水准面精化中几个问题的探讨，[硕士论文]，武汉大学，2005，13～15
[19] 陈俊勇，李建成，我国新一代高精度高分辨率大地水准面的研究和实施[J]，武汉大学学报，2001，30(2)，95～100
[20] 李春华，成都市似大地水准面精化分析研究，[硕士论文]，西南交通大学，2004，27～28
[21] 郭春喜，伍寿兵，王惠民等，区域厘米级大地水准面的确定[J]，测绘通报，2000，19(2)，26～27
[22] 李建成、陈俊勇、宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，115～130
[23] 陈俊勇，李建成，宁津生等，我国大陆高精度高分辨率大地水准面的研究和实施[J]，测绘学报，2001，30(2)，3～6
[24] Smith D A, Milbert D G, The GEOID96 High Resolution Geoid Height Model[J], Journal of Geodesy, 1999(7), 219～236
[25] 武汉测绘科技大学测量平差教研室，测量平差基础，测绘出版社，1993，65～66
[26] 李晓桓，GPS水准拟合模型的优选[J]’测绘通报，2003(7)，11～14
[27] 李春华，成都市似大地水准面精化分析研究，[硕士论文]，西南交通大学，2001，30～31
[28] 徐绍铨，李振洪，吴云孙，GPS高程拟合系统研制[J]，武汉测绘科技大学学报，1999(12)，Vol．24，No．4，336～338
[29] Wemzel HG, Determination of Active Units with Different Kinematic Property and Their Activity Pattern Based on the Data from GPS Remeasurements[J], Acta Seismologica Sinica, 2001(1), 2～28
[30] 施宝湘，施一民，逐点分别拟合及网格技术在GPS水准中的应用研究[J]，工程勘 察，2006(2)，50～51
[31] 金时华，多面函数拟合法转换GPS高程[J]，测绘与空间地理信息，2005(12)，44～47
[32] 解海中，GPS测定正常高方法研究[J]，测绘通报，1999(7)，20～21
[33] 乔仰文，GPS高程转换的若干问题的研究[J]，测绘通报，1999(11)，17～19
[34] 黄莜蓉，周光文，用多面函数法内插GPS点的高程异常[J]，工程勘察，1995(4)，36～39
[35] Lemoine F G, The development of the Joint NASA GSFC and the National Imagery and Mapping, Agency(NIMA) geopotential model EDM96, 1998, 30～34
[36] L I J. A, Formula for Computing the Gravity Disturbance from the Second Radial Derivative of the Disturbing Po2 tential[J], Journal of Geodesy, 2002(76), 226～231
[37] 李建成，陈俊勇，宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，73～74
[38] 李建成，全国及部分省市地区高精度高分辨率似大地水准面的研究和实施[J]，测绘通报，2005(5)，1～4
[39] 管泽霖，李建成，晁定波等，WZD94中国大地水准面研究[J]，武汉测绘科技大学学报，1994，19(4)，292～297
[40] Luo Z C, Chen Y Q, Evaluation of Geopotential Models EGM96, WDM94 and GPM98CR in Hong Kong and Shen Zhen[J], Journal of Geospatial Engineering, 2002, 4(1), 21～30
[41] 魏子卿，王刚，用地球位模型和GPS／水准数据确定我国大陆似大地水准面[J]，测绘学报，2003，32(1)，1～5
[2] Denker H, Behrend D, Torge W, The European Gravimetfic Quasi-geoid EGG95[J], 1996, IAG Bulletin of D'information 77, IGES Bulletin No. 4(Special Issue), 3～11
[3] Zur Verwendung, Europaischen Gravimetrischen Quasigeoids EGG97 in Deutschland[J], Manuscripat Geodaetica, 1999, 18(5), 154～166
[4] 李建成，陈俊勇，宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，6～9
[5] 宁津生，罗志才，李建成，我国省市级太地水准面精化的现状及技术模式[J]，大地测量与地球动力学，2004，Vol．24，No．1，5～7
[6] 邬建耀，关于福建地区GPS网的建立及大地水准面的精化[J]，福建建设科技，2001(3)，34～36
[7] 陈俊勇，李建成，全国及部分省市地区高精度高分辨率似大地水准面的研究和实施[J]，测绘通报，2005(5)，1～4
[8] 李骏元，程传录，郭春喜，大地水准面确定的几种方法和比较[J]，测绘技术装备，2005，Vol．7，No．2，15～16
[9] 雷晓霞，基于重力与GPS水准组合法的大地水准面精化研究，[硕士论文]，长安大学，2005，7～8
[10] Weikko A Heiskanen & Helmut Moritz，Physical Geodey，物理大地测量学(中译本)，北京，测绘出版社，1967，137～138
[11] 管泽霖，地球形状及外部重力场，北京，测绘出版社，1980，3：55～39
[12] Ardalan A A. Grafarend E, Ellipsoidal Geoidal Undulations (Ellipsoidal Bruns Formula) [J], Journal of Geodesy, 2001, 75(9/10), 544～552
[13] 李建成，陈俊勇，宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，17～18
[14] 管泽霖，宁津生，地球形状及外部重力场，北京，测绘出版社，1981，57～62
[15] Bjerhammar A, Discrete physical geodesy, Dept. of Geodetic Science and Surveying, 1987, The Ohio State University, Columbus, 125～129
[16] 许厚泽，陆仲元，中国重力场与大地水准面，解放军出版社，1997，63～67
[17] Brovar, On the Solutions of Monodensky's Boundary Value problom[J], Bull. Geoid, 1964(72), 112～115
[18] 马新莹，局部大地水准面精化中几个问题的探讨，[硕士论文]，武汉大学，2005，13～15
[19] 陈俊勇，李建成，我国新一代高精度高分辨率大地水准面的研究和实施[J]，武汉大学学报，2001，30(2)，95～100
[20] 李春华，成都市似大地水准面精化分析研究，[硕士论文]，西南交通大学，2004，27～28
[21] 郭春喜，伍寿兵，王惠民等，区域厘米级大地水准面的确定[J]，测绘通报，2000，19(2)，26～27
[22] 李建成、陈俊勇、宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，115～130
[23] 陈俊勇，李建成，宁津生等，我国大陆高精度高分辨率大地水准面的研究和实施[J]，测绘学报，2001，30(2)，3～6
[24] Smith D A, Milbert D G, The GEOID96 High Resolution Geoid Height Model[J], Journal of Geodesy, 1999(7), 219～236
[25] 武汉测绘科技大学测量平差教研室，测量平差基础，测绘出版社，1993，65～66
[26] 李晓桓，GPS水准拟合模型的优选[J]’测绘通报，2003(7)，11～14
[27] 李春华，成都市似大地水准面精化分析研究，[硕士论文]，西南交通大学，2001，30～31
[28] 徐绍铨，李振洪，吴云孙，GPS高程拟合系统研制[J]，武汉测绘科技大学学报，1999(12)，Vol．24，No．4，336～338
[29] Wemzel HG, Determination of Active Units with Different Kinematic Property and Their Activity Pattern Based on the Data from GPS Remeasurements[J], Acta Seismologica Sinica, 2001(1), 2～28
[30] 施宝湘，施一民，逐点分别拟合及网格技术在GPS水准中的应用研究[J]，工程勘 察，2006(2)，50～51
[31] 金时华，多面函数拟合法转换GPS高程[J]，测绘与空间地理信息，2005(12)，44～47
[32] 解海中，GPS测定正常高方法研究[J]，测绘通报，1999(7)，20～21
[33] 乔仰文，GPS高程转换的若干问题的研究[J]，测绘通报，1999(11)，17～19
[34] 黄莜蓉，周光文，用多面函数法内插GPS点的高程异常[J]，工程勘察，1995(4)，36～39
[35] Lemoine F G, The development of the Joint NASA GSFC and the National Imagery and Mapping, Agency(NIMA) geopotential model EDM96, 1998, 30～34
[36] L I J. A, Formula for Computing the Gravity Disturbance from the Second Radial Derivative of the Disturbing Po2 tential[J], Journal of Geodesy, 2002(76), 226～231
[37] 李建成，陈俊勇，宁津生等，地球重力场逼近理论与中国2000似大地水准面的确定，武汉大学出版社，2003，73～74
[38] 李建成，全国及部分省市地区高精度高分辨率似大地水准面的研究和实施[J]，测绘通报，2005(5)，1～4
[39] 管泽霖，李建成，晁定波等，WZD94中国大地水准面研究[J]，武汉测绘科技大学学报，1994，19(4)，292～297
[40] Luo Z C, Chen Y Q, Evaluation of Geopotential Models EGM96, WDM94 and GPM98CR in Hong Kong and Shen Zhen[J], Journal of Geospatial Engineering, 2002, 4(1), 21～30
[41] 魏子卿，王刚，用地球位模型和GPS／水准数据确定我国大陆似大地水准面[J]，测绘学报，2003，32(1)，1～5