Uniform and refinable grids on elliptic domains and on some surfaces of revolution
Introduction
Uniform and refinable grids (UR) are useful in many applications, like construction of multiresolution analysis and wavelets, or for solving numerically partial differential equations. While on a rectangle or on other polygonal domains the construction of UR grids is trivial, it is not immediate on an elliptic domain or on a disc.
In this paper we construct a bijection from to , which maps rectangles onto ellipses and preserves area. This allows us to transport a rectangular grid to an elliptic grid, preserving the area of the cells. In particular, any uniform1 rectangular grid is mapped into a uniform elliptic grid. A refinement process is needed when a grid is not fine enough to solve a problem accurately. A uniform refinement consists in dividing a cell into a given number of smaller cells with the same area. With the procedure described here, any uniform refinement of a rectangular grid leads to a uniform refinement of the corresponding elliptic grid.
In the particular case of the disc, such a bijection was constructed in a previous paper [2] and combined with Lambert azimuthal projection, helped us to construct uniform grids on the sphere. Also, the inverse projection, mapping discs onto squares, can be used to map a hemisphere of the Earth onto a square domain, more suitable for computer manipulations. The main advantage over the existing map projections is the area preserving property, desirable in statistical mapping.
We mention that the problem of finding an analytic function that maps a disc onto the interior of any convex polygon was first solved independently by E.B. Christoffel (1867) and H.A. Schwarz (1869). Their map is a complex integral which is conformal (i.e. preserves angles), but it does not preserves areas. In fact, it is known that no conformal map can preserve areas. Another drawback is that it has no explicit formula, requiring solving a system of nonlinear equations.
Some well-known projections mapping a hemi-sphere onto a square are Pierce quincuncial (1879 - conformal except four points), Guyou (1886 - conformal), Adams (1925 - conformal), Collignon (1865 - area preserving). An area preserving projection from a sphere to a rectangle is Lambert Cylindrical projection. A complete description of all known map projections used in cartography can be found in [1], [4].
Section snippets
Construction of an area preserving bijection in
Consider the ellipse of semi-axes a and b, a, b > 0, of equationand the rectangle with edges 2L1 and 2L2, defined asThe domains enclosed by and will be denoted by and , respectively.
We will construct a bijection which maps each rectangle onto the ellipse and has the area preserving propertyHere denotes the area of D. Thus,
Uniform grids on surfaces of revolution
We have presented in [3] a method for constructing a projection from a surface of revolution onto a plane perpendicular to the rotation axis, which preserves area. The method works for a surface of revolution generated by a plane curve of equation z = φ(x), with φ a nonnegative and increasing piecewise smooth function, defined on the interval I = [0, β) or I = [0, ∞). This projection allows us to construct a continuous wavelet transform and a multiresolution analysis, in the case when the curve φ has
Acknowledgements
I thank to the anonymous referees for their constructive suggestions for improving the paper. The work has been co-funded by the Sectorial Operational Programme Human Resources Development 2007-2013 of the Romanian Ministry of Labor, Family and Social Protection through the Financial AgreementPOSDRU/89/1.5/S/62557.
References (4)
Wavelet analysis on some surfaces of revolution via area preserving projection
Appl. Comput. Harmon. Anal.
(2011)- et al.
Map Projections Cartographic Information Systems
(2006)
Cited by (6)
Uniform spherical grids via equal area projection from the cube to the sphere
2011, Journal of Computational and Applied MathematicsCitation Excerpt :In [10], one of the authors suggested a new area preserving projection method based on a mapping of the square onto a disc in a first step, followed by a lifting to the sphere by the inverse Lambert projection. This idea can also be generalized to construct uniform and refinable grids on elliptic domains and on some surfaces of revolution, see [11]. Fig. 3 shows two grids on the curved square.
Generalized orthonormal piecewise constant wavelet bases with 1-4 splitting
2012, International Journal of Wavelets, Multiresolution and Information ProcessingA Square Equal-Area Map Projection with Low Angular Distortion, Minimal Cusps, and Closed-Form Solutions
2021, ACM Transactions on Spatial Algorithms and SystemsDistortion analysis of approximate equal-area grids based on octahedron
2016, Wuhan Daxue Xuebao (Xinxi Kexue Ban)/Geomatics and Information Science of Wuhan UniversityAn Area Preserving Projection from the Regular Octahedron to the Sphere
2012, Results in Mathematics