Triangular grid-based rapid mapping of scattered data
LI Xiao-Dong1,2,3, JIN Sheng1,2, WANG Yang-Ling1,2, ZHANG Jia-Hong4, CHENG Li-Hui1,2
1. School of Geophysics and Information Technology, China University of Geosciences, Beijing 100083, China; 2. Geo-detection and instruments Laboratory, Ministry of Education, China University of Geosciences, Beijing 100083, China; 3. Non-ferrous Metals and Nuclear Industry Geological Exploration Bureau of Guizhou, Geophysical and Geochemical Exploring Group, Dujun 558004, China; 4. China Aero Geophysics Survey & Remote Sensing Center foe Land and Resources, Beijing 100083, China
Abstract:Conventional contour mapping performs interpolation based on a rectangular grid. A linear interpolation method is presented in this paper based on triangular mesh. Triangular mesh can better approximate the boundary of scattered data and the morphology of geophysical field, which makes the contour maps smoother. By searching boundary, triangulated mesh, linear interpolation, search contours, Bezier curves and smooth contours, five steps can be carried out quickly for any scattered data mapping. The actual data mapping results show that the interpolation method is good in that no data extrapolation is needed, the contour map obtained directly reflect the spatial location of scattered data, and the mapping is speeded. The method can therefore greatly improve the efficiency of the actual work.
[1] 刘兆平,杨进,武炜.地球物理数据网格化方法的选取[J].物探与化探,2010,34(1):93-97. [2] 陈欢欢,李星,丁文秀.Surfer8.0等值线绘制中的十二种插值方法[J].工程地球物理学报,2007,4(1):52. [3] Goetze X L A H. Comparison of some gridding methods [J]. The Leading Edge, 1999, (8):898-900. [4] 郭良辉, 孟小红,郭志宏, 等.地球物理不规则分布数据的空间网格化法[J].物探与化探,2005,29(5):438-442. [5] Hardy R L. Multiquadric Equations of Topography and Other Irregular Surfaces [J]. Journal of Geophysics Research,1971,78: 1905-1915. [6] 郭志宏. 一种实用的等值线型数据网格化方法 [J].物探与化探,2001,25(3):203-208. [7] 孟小红, 侯建全, 梁宏英, 等. 离散光滑插值方法在地球物理位场中的快速实现 [J].物探与化探, 2002, 26(4): 302-306. [8] Cormen T H.算法导论[M].潘金贵,顾铁成,译.北京:机械工业出版社,2013. [9] Graham R L.An efficient algorithm for determining the convex hull of a finite planar set[J].Information processing letters, 1972, 1(4): 132-133. [10] Jarvis R A.On the identification of the convex hull of a finite set of points in the plane[J].Information Processing Letters,1973,(2):18-21. [11] 陈涛, 李光耀. 平面离散点集的边界搜索算法 [J].计算机仿真,2004,21(3):21-23. [12] Lee D T,Schacheer B J.Two algorithms for constructing a Delaunay triangulation[J].International Journal of Computer and Information Science,1980,9(3):219-242. [13] Watson D F.Computing the n-dimension Delaunay tessellation with application to Voronoi polygons[J].The computer journal, 1981, 24(2): 167-172. [14] Persson P O, Strang G. A simple mesh generator in MATLAB[J]. SIAM review, 2004, 46(2): 329-345. [15] Shewchuk J R. Triangle: Engineering a 2D quality mesh generator and delaunay triangulator[J]. Applied Computational Geometry Towards Geometric Engineering,1996,1148:124-133. [16] 徐世浙.地球物理中的有限单元法[M].北京:科学出版社,1994:40-42. [17] 柳建新, 蒋鹏飞, 童孝忠, 等. 不完全 LU 分解预处理的 BICGSTAB 算法在大地电磁二维正演模拟中的应用[J].中南大学学报: 自然科学版, 2009, 40(2): 484-491. [18] 苗润忠.光滑的等值线生成算法[J].长春理工大学学报,2004,(1):16-18. [19] 施法中. 计算机辅助几何设计与非均匀有理B样条[M]. 北京: 北京航空航天大学出版社, 1994. [20] ROBERT C B. An Introduction to the Curves and Surfaces of Computer-aided Design [M]. Stanford Linear Accelerator Center: Van Nostrand Reinhold,1991. [21] Mallet J L. Discrete smooth interpolation in geometric modeling [J]. Computer-aided design, 1992, 24(4): 178-191. [22] 王正林,龚纯.精通matlab科学计算[M].北京:电子工业出版社,2007. [23] Farin G. Algorithms for rational Bézier curves[J].Computer-Aided Design,1983,15(2):73-77.