Multiresolution terrain modeling using level curve information

The core of the presented multiresolution method is an algorithm for removing recursively the level curves according to some error criterion. This allows us to obtain a sequence of approximations of the terrain where the difference between two consecutive approximations is only one curve: the less ''important''. In other words, the input curves are sorted in such a way that the n most relevant curves are contained in the n-th resolution level. For a given curve, the relevance criterion is the error computed using a function interpolating the remaining curves. Hence, to fully formulate the multiresolution algorithm a function interpolating the contour lines of the terrain is necessary. We note that a contour line representation has a higher density in the horizontal direction than in the vertical. To alleviate this problem, a horizontal simplification algorithm for each curve is proposed. Computational efficiency concerns arising from the size of the datasets, such as the computation of the distance from a point to a polygon (with a huge number of vertices), and the point location problem, are addressed. To obtain an efficient implementation of the proposed method, it was necessary to use adequate data structures and computational geometry algorithms, in order to solve several subproblems, for instance: the computation of the distance from a point to a polygon (with a huge amount of vertices), the simplification of the level curves and the point location problem. Finally, we show how to visualize the interpolant corresponding to the n-th resolution level.

[1]  Honglei Zhu,et al.  Triangulated irregular network optimization from contour data using bridge and tunnel edge removal , 2001, Int. J. Geogr. Inf. Sci..

[2]  Hans-Peter Seidel,et al.  A General Framework for Mesh Decimation , 1998, Graphics Interface.

[3]  Alper Üngör O.-Centers: A New Type of Steiner Points for Computing Size-Optimal Quality-Guaranteed Delaunay Triangulations , 2004, LATIN.

[4]  Peter Schröder,et al.  C1-Continuous Terrain Reconstruction from Sparse Contours , 2003, VMV.

[5]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[6]  Eric Keppel,et al.  Approximating Complex Surfaces by Triangulation of Contour Lines , 1975, IBM J. Res. Dev..

[7]  Jean Duchon,et al.  Splines minimizing rotation-invariant semi-norms in Sobolev spaces , 1976, Constructive Theory of Functions of Several Variables.

[8]  A. B. Ekoule,et al.  A triangulation algorithm from arbitrary shaped multiple planar contours , 1991, TOGS.

[9]  Hendrik Speleers,et al.  Local subdivision of Powell-Sabin splines , 2006, Comput. Aided Geom. Des..

[10]  Günther Greiner,et al.  Interpolating and approximating scattered 3D-data with hierarchical tensor product B-splines , 2010 .

[12]  Nira Dyn,et al.  Adaptive Thinning for Terrain Modelling and Image Compression , 2005, Advances in Multiresolution for Geometric Modelling.

[13]  Brian G. Schunck,et al.  A Two-Stage Algorithm for Discontinuity-Preserving Surface Reconstruction , 1992, IEEE Trans. Pattern Anal. Mach. Intell..

[14]  Alper Üngör,et al.  Off-centers: A new type of Steiner points for computing size-optimal quality-guaranteed Delaunay triangulations , 2004, Comput. Geom..

[15]  Luiz Velho,et al.  Multi-resolution triangulations with adaptation to the domain based on physical compression , 2004, Proceedings. 17th Brazilian Symposium on Computer Graphics and Image Processing.

[16]  Andrew H. Gee,et al.  Surface interpolation from sparse cross sections using region correspondence , 2000, IEEE Transactions on Medical Imaging.

[17]  W. Randolph Franklin,et al.  Constructing a dem from grid-based data by computing intermediate contours , 2003, GIS '03.

[18]  M. Floater,et al.  Thinning algorithms for scattered data interpolation , 1998 .

[19]  Jean-Daniel Boissonnat,et al.  Shape reconstruction from planar cross sections , 1988, Comput. Vis. Graph. Image Process..

[20]  R. Schaback,et al.  Numerical Techniques Based on Radial Basis Functions , 2000 .

[21]  Larry L. Schumaker,et al.  Curve and Surface Fitting: Saint-Malo 1999 , 2000 .

[22]  Adhemar Bultheel,et al.  Surface Compression with Hierarchical Powell-Sabin B-splines , 2006, Int. J. Wavelets Multiresolution Inf. Process..

[23]  Bernd Hamann,et al.  Piecewise optimal triangulation for the approximation of scattered data in the plane , 2000, Comput. Aided Geom. Des..