椭球面三角形外心的地图代数解法  

Determination of Circumcenter of Triangle on Ellipsoidal Surface Based on Map Algebra

作　　者：蒋会平[1,2,3] 谭树东[4] 胡海[1] 

机构地区：[1]武汉大学资源与环境科学学院,湖北武汉430079 [2]中国科学院地理科学与资源研究所资源与环境信息系统国家重点实验室,北京100101 [3]中国科学院大学,北京100049 [4]国家海洋信息中心,天津300171

出　　处：《测绘学报》2016年第2期241-249,共9页Acta Geodaetica et Cartographica Sinica

基　　金：国家自然科学基金(41271443;41471328);国家863计划(2009AA12Z224)~~

摘　　要：椭球面三角形外心到3个相邻顶点的大地线距离都相等。面向椭球面空间的外心大地坐标的求解对于椭球面Voronoi图的生成和椭球面Delaunay三角网的构造具有重要作用。利用基于地图代数理论的矢栅结合方法,首先基于地图代数测地变换建立高精度椭球面空间距离场,再通过边界跟踪配对确定外心所在的栅格范围,最后通过数值计算内插生成初始等距点并不断逼近外心的精确大地坐标。试验结果表明,采用本文方法求解的椭球面三角形外心大地坐标,在103~104 km跨度内其定位误差小于0.001m,且算法非常适用于海量空间数据的高精度快速计算。The geodesic distances from the circumcenter to 3vertexes of the triangle on ellipsoidal surface are equal.The ellipsoid-oriented determination of circumcenter of triangle on ellipsoidal surface is applicable when it comes to generation of the Voronoi diagram and construction of the Delaunay triangulation net on the ellipsoidal earth,which can be considered as a solution of significance in computation of geometries and spatial analysis on the ellipsoid.Based on the idea of combining the raster and vector methods and the theory of map algebra,the working process can be described as below：firstly,initiate the geographical distance transformation and create the distance field with a high degree of accuracy;secondly,conduct boundary tracking and matching and then determinate the range of grids where the circumcenter of triangle locates;thirdly,interpolate the initial equidistant point;finally,approximate the circumcenter of triangle on earth ellipsoidal surface by means of numeric calculation.The positioning error of this algorithm is controlled less than 0.001 m within several thousand kilometers range of span.As regards the method proposed in the present paper,its computational efficiency is O（m）where mis the number of pixels in the image,i.e.,grid resolution.In conclusion,this algorithm can be considered as both ellipsoid-oriented and not content-related,which is especially appropriate for complex geocomputation globally.

关 键 词：地图代数 测地变换 椭球面 三角形外心 矢栅结合 

分 类 号：P208[天文地球—地图制图学与地理信息工程]