Uniform spherical grids via equal area projection from the cube to the sphere
|
|
- Silas George
- 6 years ago
- Views:
Transcription
1 Uniform spherical grids via equal area projection from the cube to the sphere Daniela Roşca Gerlind Plonka April 4, 0 Abstract We construct an area preserving map from the cube to the unit sphere S, both centered at the origin. More precisely, each face F i of the cube is first projected to a curved square S i of the same area, and then each S i is projected onto the sphere by inverse Lambert azimuthal equal area projection, ith respect to the points situated at the intersection of the coordinate axes ith S. This map is then used to construct uniform and refinable grids on a sphere, starting from any grid on a square. Key ords: equal area projection, uniform spherical grid, refinable grid, hierarchical grid. Introduction The generation of suitable finite subsets of the unit sphere is closely related to all kinds of data approximation of the sphere. Generic approximation problems on the unit sphere include the question of ho to find the best element from a defined function set that interpolates or approximates given function values on points of the sphere. Here the question that arises, is ho the given function values should be distributed on the sphere in order to achieve optimal approximation results. In many applications, especially in geosciences and astronomy, but also in medical imaging and computer vision, one requires simple, uniform and refinable grids on the sphere. One simple method to construct such grids is to transfer existing planar grids. A complete description of all knon spherical projections from a sphere or parts of a sphere to the plane, used in cartography, is realized in [, 8]. Hoever, most grid constructions given so far do not provide an equal area partition. But this property is crucial for a series of applications including statistical computations and avelet constructions, since non-equal area partitions can generate severe distortions at large distances. The existence of partitions of S into regions of equal area and small diameter has been already used by Alexander [], ho derives loer bounds for the maximum sum of distances beteen points on the sphere. Based on the construction by Zhou [], Leopardi derives a recursive zonal equal area sphere partitioning algorithm for the unit sphere S d embedded in R n+, see [4]. The constructed partition in [4] for S consists of polar cups and rectilinear regions that are arranged in zonal collars. Besides the problem that e have to deal ith different kinds of areas, the obtained partition is not suitable for various applications here one needs to avoid that vertices of spherical rectangles lie on edges of neighboring rectangles. Technical University of Cluj-Napoca, Department of Mathematics, str. Memorandumului 8, RO-4004 Cluj- Napoca, Romania. Daniela.Rosca@math.utcluj.ro University of Göttingen, Institute for Numerical and Applied Mathematics, Lotzestr. 6-8, Göttingen, Germany. plonka@math.uni-goettingen.de
2 Other constructions for equal area partitions of S used in astronomy include the HEALPix grid [3], providing a hierarchical equal area iso-latitude pixelation, and the icosahedron-based method by Tegmark [0], see also [9]. In [7], an equal area global partition method based upon circle edges is presented. Starting ith a spherical triangulation, obtained e.g. by projecting the faces of an icosahedron to the sphere, a subdivision method is proposed to partition each spherical triangle into four equal area subtriangles. In [5], one of the authors suggested a ne area preserving projection method based on a mapping of the square onto a disc in a first step, folloed 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 [6]. In this paper e construct an area preserving map from a cube to the unit sphere S. Thus, any grid on the cube can be transported to the sphere. The ne construction admits to transform ell-knon techniques from the square to (parts of) the sphere. For example, the ell-knon tensor-product avelet transform for data analysis and denoising can be simply transferred to the sphere, here distortions are negligible because of the equal area property. Further, since arbitrary grids on the cube resp. on the square can no simply be transported to the sphere, e believe that this construction may achieve an essential impact for different applications in geosciences. Since e give explicit formulas both for the map from the cube to the sphere, and from the sphere to the cube, the method is easy to implement. Our ne construction scheme consists of to steps. In the first step, e construct in Section 3 a bijection T from each face F i of the cube, onto a curved square S i. In Section 4 e construct the inverse T. In the second step, e combine T ith the inverse Lambert azimuthal projection, in order to map each face F i of the cube onto a subset F i of the sphere, such that 6 i= F i = S. Finally, e present some examples of the obtained spherical grids. Preliminaries Consider the unit sphere S centered at the origin O and the cube K centered at O, ith the same area. Thus, the edge of the cube has the length a = π/3. Denote = a/ = π/6. We cut the sphere ith the six diagonal planes z = ±x, z = ±y, y = ±x and obtain the curves in Figure. Let us focus on one of these curves C, given by the equations x + y + z =, z = [ x > 0, y 3, 3 ]. Using the parameter t [ 3, 3 ], the curve C is described by (x, y, z) = ( ( t ), t, ( t )). The Lambert azimuthal equal-area projection of the sphere x + y + z = onto the plane z = is given by ( ) (x L, y L ) = + z x, + z y. () Hence, the Lambert projection of the curve C onto the plane z = has the equations t x L =, () + t y L = z L =, t + t, (3)
3 Figure : The curves of intersection of the cube ith the diagonal plans. ith t [ / 3, / 3]. The calculations sho that x L + yl = t, (4) y L t = x. (5) L t We denote by S a the curved square in the tangent plane z =, formed by the curved edge of equations ()-(3) and three other curved edges obtained in the same manner by intersecting the sphere ith the half-planes z = x > 0, z = y > 0 and z = y > 0, respectively (Figure ). Because of the area-preserving property of the Lambert projection, the area A(S a ) of the obtained curved square S a is exactly π 3 = a. With the help of this projection, e have simplified the problem of finding an area preserving map from the cube to the sphere S to a problem of finding a to-dimensional map from a square, i.e., a face of the cube, to the curved square S a. 3 Mapping a square onto a curved square In this section e derive an area preserving bijection T : R R, hich maps the square S a = {(x, y) R, x a, y a} onto the curved square S a that has been constructed in Section. Here, e say that T is area-preserving, if it has the property A(D) = A(T (D)), for every domain D R. (6) Since e consider only the to-dimensional problem in this section, e ill ork in the (x, y)- plane for simplicity, i.e., the curved boundary C a of S a is given by () and (3) ith the notation x = x L and y = y L, and analogously for the three other curves. We focus for the moment on the first octant of the (x, y)-plane I = {(x, y) R, 0 < y x} ith origin O = (0, 0) and take a point M = (x M, y M ) = (x M, mx M ) I, here m is a parameter ith 0 m, see Figure. The map T ill be defined in such a ay that each half-line d m of equation y = m x (0 m ) is mapped onto the half-line d φ(m) of equation y = φ(m)x, here φ : [0, ] R is a differentiable function that satisfies φ(0) = 0, φ() = and 0 φ(m). (7) 3
4 M N P Q O Figure : The square S a, the curved square S a and the action of the transform T. Points N and P are the images of M and Q. This condition ensures that the diagonal of the square S a and the axis x are invariant under φ. We denote by (x N, y N ) the coordinates of the point N = T (M). Let Q be the intersection of OM ith the square S a (see Figure ). The point Q has the coordinates (x Q, y Q ) = (, m), here = a = π 6, and the line ON has the equation y = φ(m)x. Further, let the point P = (x P, y P ) = P (x P, φ(m)x P ) be the intersection of ON ith the curved square S a. Thus, the coordinates of P satisfy the equations (4) and (5) ith some t P [0, 3 ], and from (5) e have φ(m) = tp φ(m), hence t t P = P + φ (m). Replacing t P in () and (3), e obtain the coordinates of P in the form x P = y P = + φ (m), + +φ (m) φ(m) + φ (m). + +φ (m) Some simple calculations yield the distances OM = x M + m, OQ = + m, ON = x N + φ (m), ( / OP =. + φ (m)) We ant to determine the map T such that ON OP = OM OQ. 4
5 From the above calculations e then obtain x N = x + φ (m) y N = φ(m)x N = x + φ (m), φ(m) + φ (m) + φ (m). With this assertion, the map T is no completely described by means of the function φ, and e obtain that T maps the point (x, y) I onto the point (X, Y ) given by X = x Y = x + φ ( y x ) + φ ( y (8) x ), φ( y x ) + φ ( y x ) + φ ( y (9) x ). Next e impose the area preserving property by a suitable determination of φ. For this purpose, e define the function φ such that the Jacobian of T is. After simplification, the Jacobian rites as J(T ) = det ( X x Y x X y Y y ) = ( φ y x) + φ ( y x ) + + φ ( ). y For solving the equation J(T ) = e substitute v := y x, and thus, in the considered case 0 < y x e have v (0, ]. Hence, ith the simplified notation φ = φ(v), e get Integration gives φ + φ + + φ =. arctan φ arctan φ + φ = v + C. The condition φ(0) = 0 yields C = 0. Next, in order to determine φ e use the formula x arctan a arctan b = arctan a b + ab a, b R, ab > and e further obtain φ ( + φ ) + φ + φ = tan v. (0) To simplify this term, e introduce the notation γ = tan v = tan π v, and equality (0) yields + φ (φ γ) = φ(γφ + ). () From the requirements (7) e can deduce that, for (x, y) I ith y > 0, both sides of equality () are positive. Thus, () is equivalent ith ( + φ )(φ γ) = φ (γφ + ) (φ + )(( γ )φ 4φγ + γ ) = 0, hich gives φ, (v) = γ ± γ + γ γ. 5
6 We denote α = v = y x. Thus γ = tan α and some calculations sho that the function φ becomes ( y ) sin α φ, (v) = φ, =, ith α = πy x cos α ± x. If e take into account the condition φ() = and the equality sin π = cos π, the only convenient solution is ( y sin α φ =, x) cos α ith α = πy x. Finally, e need to replace φ in the equations (8)-(9). Some straightforard calculations sho that φ + = cos α cos α, φ + = 3 cos α cos α, φ φ + = sin α 3 cos α. Hence, for (x, y) I, the formulas for the desired area preserving map T (x, y) = (X, Y ) are given by X = /4 x cos α, cos α Y = /4 x sin α, ith α = yπ π cos α x, = 6. Similar arguments for the other seven octants sho that the area preserving map T : R R hich maps squares into curved squares is defined as follos: For y x, For x y, (x, y) (X, Y ) = (x, y) (X, Y ) = /4 x /4 y Figure 3 shos to grids on the curved square. 4 The inverse map cos yπ x, /4 x cos yπ x sin xπ y, /4 y cos xπ y sin yπ x ; () cos yπ x cos xπ y. (3) cos xπ y To make the area preserving map T applicable in practice, e need also to derive a closed simple form for the inverse mapping T. Let us consider first the case hen x y > 0, hen the map T is given by formula (). With the notation η = cos yπ x, (4) 6
7 Figure 3: To grids on the curved square, ith seven and eight equidistant horizontal and vertical lines. e have In particular, X = x ( η ) x, Y ( η ) = η. η X + Y = X + Y = X Y = x 3 η, η x ( η), (5) x η(3η ). η We define Hence, e obtain the equation in η ith the solutions B := X + Y X + Y = ( η) 3 η. η ( B)η + 3B = 0 η, = ( B) ± B(B ). In our case, η (cos π, ), therefore the only convenient solution is η = ( B) + B(B ). Then, from it follos that B = Y (X + Y ) and B = X X + Y η = Y + X X + Y. (6) (X + Y ) For simplicity e introduce the notation := Y/X. Thus, e have [0, ] and ( ) π y + = arccos η = arccos x ( + ) +. ( + ) 7
8 No e use the identity for a = arccos(ab + a b ) = arccos a arccos b, 0 a b, and b = (+ ) + π y x = arccos and obtain ( + ) arccos Other possible expressions are π y + x = arctan arctan = arctan arctan + +. = arctan Y X arctan Y X + Y. (7) For the calculation of x from X and Y, e use the second equality in (5) and η = + + (+ ) in order to find x = X + Y = X η = X Finally, from (4) and (7) e find y = x π x arccos η = π ( + ) η = X + ( + + ) = X + Y ( X + Y + X). ( arctan arctan ) + = X 4 + π + + (arctan arctan 4 = X + Y X + ( X + Y arctan Y X arctan ) + ) Y. X + Y Using similar arguments for the other seven octants, e obtain that the inverse T : R R has the expressions For Y X, ( (X, Y ) (x, y) = sign(x) 4 X + Y X + X + Y, For X Y, 4 X + Y X + X + Y (sign(x) arctan Y X arctan Y X +Y ) ) ; (X, Y ) (x, y)=( 4 X + Y Y + X + Y (sign(y ) arctan X Y arctan X X ), +Y sign(y ) 4 X + Y Y + ) X + Y. 8
9 5 Mapping the curved squares onto the sphere The complete mapping from the cube to the sphere S is no described in to steps. In the first step, each face F i of the cube K ill be mapped onto a domain F i, bounded by a curved square, using the transformation T. In the second step, each F i ill be mapped onto F i S by the inverse Lambert azimuthal projection, ith respect to the center of F i. Obviously 6 F i = i= and 6 F i = S. i= We denote by L (0,0,) the Lambert azimuthal area preserving projection ith respect to the North Pole N = (0, 0, ). Remember that L (0,0,) maps a point (x, y, z) S onto the point (X L, Y L, ) situated in the tangent plane at N, as given in (). Applying the inverse Lambert projection L (0,0,), the point (X, Y, ), situated in the tangent plane to S at the pole N, maps onto (x L, y L, z L ) S given by x L = X + Y X, (8) 4 y L = X + Y Y, (9) 4 z L = X + Y. (0) Thus, the application T L (0,0,) maps the upper face of the cube projects onto F S. The formulas for T L (0,0,)(x, y) are given by (8), (9) and (0), here X, Y are calculated ith formulas (), (3). One can easily obtain similar formulas for T L (0,0, ), T L (0,,0), T L (0,,0), T L (,0,0) and T L (,0,0). Figure 4 shos some grids of the sphere, here a regular partition of the faces of the cube into 49, 64, and 44 equal area squares has been applied. Acknoledgement The ork has been funded by the grant PL 70/4- of the German Research Foundation (DFG). This is gratefully acknoledged. References [] R. Alexander, On the sume of distances beteen N points on the sphere, Acta Mathematica 3 (97), [] E.W. Grafared, F.W. Krumm, Map Projections, Cartographic Information Systems, Springer-Verlag, Berlin, 006. [3] K.M. Górski, B.D. Wandelt, E. Hivon, A.J. Banday, B.D. Wandelt, F.K. Hansen, M. Reinecke, M. Bartelmann, HEALPix: A frameork for high-resolution discretization and fast analysis of data distributed on the sphere, The Astrophysical Journal, 6, (005), p [4] P. Leopardi, A partition of the unit sphere into regions of equal area and small diameter, Electron. Trans. on Numer. Anal. 5 (006), [5] D. Roşca, Ne uniform grids on the sphere, Astronomy & Astrophysics, 50, A63 (00). [6] D. Roşca, Uniform and refinable grids on elliptic domains and on some surfaces of revolution, Appl. Math. Comput, accepted doi:0.06/j.amc
10 Figure 4: Grids on the sphere. 0
11 [7] L. Song, A.J. Kimerling, K. Sahr, Developing an equal area global grid by small circle subdivision, in Discrete Global Grids, M. Goodchild and A.J. Kimerling, eds., National Center for Geographic Information & Analysis, Santa Barbara, CA, USA, 00. [8] J.P. Snyder, Flattening the Earth, University of Chicago Press, 990. [9] N.A. Teanby, An icosahedron-based method for even binning of globally distributed remote sensing data, Computers & Geosciences 3(9) (006), [0] M. Tegmark, An icosahedron-based method for pixelizing the celestial sphere, The Astrophysical Journal 470 (996), L8-L84. [] Y.M. Zhou, Arrangements of points on the sphere, PhD thesis, Mathematics, Tampa, FL, 995.
Chapter 6. Curves and Surfaces. 6.1 Graphs as Surfaces
Chapter 6 Curves and Surfaces In Chapter 2 a plane is defined as the zero set of a linear function in R 3. It is expected a surface is the zero set of a differentiable function in R n. To motivate, graphs
More informationFall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics
Fall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics Lecture 5 August 31 2016 Topics: Polar coordinate system Conversion of polar coordinates to 2-D
More informationCalculus III. Math 233 Spring In-term exam April 11th. Suggested solutions
Calculus III Math Spring 7 In-term exam April th. Suggested solutions This exam contains sixteen problems numbered through 6. Problems 5 are multiple choice problems, which each count 5% of your total
More informationA dynamic programming algorithm for perceptually consistent stereo
A dynamic programming algorithm for perceptually consistent stereo The Harvard community has made this article openly available. Please share ho this access benefits you. Your story matters. Citation Accessed
More informationVectors and the Geometry of Space
Vectors and the Geometry of Space In Figure 11.43, consider the line L through the point P(x 1, y 1, z 1 ) and parallel to the vector. The vector v is a direction vector for the line L, and a, b, and c
More informationThe rhealpix Discrete Global Grid System
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS The rhealpix Discrete Global Grid System To cite this article: R G Gibb 2016 IOP Conf. Ser.: Earth Environ. Sci. 34 012012 View
More informationMeasuring Lengths The First Fundamental Form
Differential Geometry Lia Vas Measuring Lengths The First Fundamental Form Patching up the Coordinate Patches. Recall that a proper coordinate patch of a surface is given by parametric equations x = (x(u,
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Homework 1 - Solutions 3. 2 Homework 2 - Solutions 13
MATH 32B-2 (8) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Homework - Solutions 3 2 Homework 2 - Solutions 3 3 Homework 3 - Solutions 9 MATH 32B-2 (8) (L) G. Liu / (TA) A. Zhou Calculus
More informationCoordinate Transformations in Advanced Calculus
Coordinate Transformations in Advanced Calculus by Sacha Nandlall T.A. for MATH 264, McGill University Email: sacha.nandlall@mail.mcgill.ca Website: http://www.resanova.com/teaching/calculus/ Fall 2006,
More informationJanuary 30, 2019 LECTURE 2: FUNCTIONS OF SEVERAL VARIABLES.
January 30, 2019 LECTURE 2: FUNCTIONS OF SEVERAL VARIABLES 110211 HONORS MULTIVARIABLE CALCULUS PROFESSOR RICHARD BROWN Synopsis Today we begin the course in earnest in Chapter 2, although, again like
More informationDifferential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder]
Differential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder] Preliminaries Recall: Given a smooth function f:r R, the function
More informationMoore Catholic High School Math Department
Moore Catholic High School Math Department Geometry Vocabulary The following is a list of terms and properties which are necessary for success in a Geometry class. You will be tested on these terms during
More informationCS-9645 Introduction to Computer Vision Techniques Winter 2019
Table of Contents Projective Geometry... 1 Definitions...1 Axioms of Projective Geometry... Ideal Points...3 Geometric Interpretation... 3 Fundamental Transformations of Projective Geometry... 4 The D
More information16.6. Parametric Surfaces. Parametric Surfaces. Parametric Surfaces. Vector Calculus. Parametric Surfaces and Their Areas
16 Vector Calculus 16.6 and Their Areas Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. and Their Areas Here we use vector functions to describe more general
More informationFRACTAL 3D MODELING OF ASTEROIDS USING WAVELETS ON ARBITRARY MESHES 1 INTRODUCTION: FRACTAL SURFACES IN NATURE
1 FRACTAL D MODELING OF ASTEROIDS USING WAVELETS ON ARBITRARY MESHES André Jalobeanu Automated Learning Group, USRA/RIACS NASA Ames Research Center MS 269-4, Moffett Field CA 9405-1000, USA ajalobea@email.arc.nasa.gov
More informationAffine Transformations Computer Graphics Scott D. Anderson
Affine Transformations Computer Graphics Scott D. Anderson 1 Linear Combinations To understand the poer of an affine transformation, it s helpful to understand the idea of a linear combination. If e have
More informationB.Stat / B.Math. Entrance Examination 2017
B.Stat / B.Math. Entrance Examination 017 BOOKLET NO. TEST CODE : UGA Forenoon Questions : 0 Time : hours Write your Name, Registration Number, Test Centre, Test Code and the Number of this Booklet in
More informationJim Lambers MAT 169 Fall Semester Lecture 33 Notes
Jim Lambers MAT 169 Fall Semester 2009-10 Lecture 33 Notes These notes correspond to Section 9.3 in the text. Polar Coordinates Throughout this course, we have denoted a point in the plane by an ordered
More informationA TESSELLATION FOR ALGEBRAIC SURFACES IN CP 3
A TESSELLATION FOR ALGEBRAIC SURFACES IN CP 3 ANDREW J. HANSON AND JI-PING SHA In this paper we present a systematic and explicit algorithm for tessellating the algebraic surfaces (real 4-manifolds) F
More informationMATERIAL FOR A MASTERCLASS ON HYPERBOLIC GEOMETRY. Timeline. 10 minutes Exercise session: Introducing curved spaces
MATERIAL FOR A MASTERCLASS ON HYPERBOLIC GEOMETRY Timeline 10 minutes Introduction and History 10 minutes Exercise session: Introducing curved spaces 5 minutes Talk: spherical lines and polygons 15 minutes
More informationEuler s Theorem. Brett Chenoweth. February 26, 2013
Euler s Theorem Brett Chenoweth February 26, 2013 1 Introduction This summer I have spent six weeks of my holidays working on a research project funded by the AMSI. The title of my project was Euler s
More informationAnswer Key: Three-Dimensional Cross Sections
Geometry A Unit Answer Key: Three-Dimensional Cross Sections Name Date Objectives In this lesson, you will: visualize three-dimensional objects from different perspectives be able to create a projection
More informationWe have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance.
Solid geometry We have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance. First, note that everything we have proven for the
More information= f (a, b) + (hf x + kf y ) (a,b) +
Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals
More informationConics, Parametric Equations, and Polar Coordinates. Copyright Cengage Learning. All rights reserved.
10 Conics, Parametric Equations, and Polar Coordinates Copyright Cengage Learning. All rights reserved. 10.5 Area and Arc Length in Polar Coordinates Copyright Cengage Learning. All rights reserved. Objectives
More informationStructure in Quaternions Corresponding to the 4-Dimensional Tetrahedron
Structure in Quaternions Corresponding to the 4-Dimensional Tetrahedron AJ Friend School of Mathematics Georgia Institute of Technology Atlanta, GA Advised by: Adrian Ocneanu Department of Mathematics
More informationLifting Transform, Voronoi, Delaunay, Convex Hulls
Lifting Transform, Voronoi, Delaunay, Convex Hulls Subhash Suri Department of Computer Science University of California Santa Barbara, CA 93106 1 Lifting Transform (A combination of Pless notes and my
More informationMA 243 Calculus III Fall Assignment 1. Reading assignments are found in James Stewart s Calculus (Early Transcendentals)
MA 43 Calculus III Fall 8 Dr. E. Jacobs Assignments Reading assignments are found in James Stewart s Calculus (Early Transcendentals) Assignment. Spheres and Other Surfaces Read. -. and.6 Section./Problems
More informationChapter 11. Parametric Equations And Polar Coordinates
Instructor: Prof. Dr. Ayman H. Sakka Chapter 11 Parametric Equations And Polar Coordinates In this chapter we study new ways to define curves in the plane, give geometric definitions of parabolas, ellipses,
More informationEfficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves
Efficient Storage and Processing of Adaptive Triangular Grids using Sierpinski Curves Csaba Attila Vigh, Dr. Michael Bader Department of Informatics, TU München JASS 2006, course 2: Numerical Simulation:
More informationConics, Parametric Equations, and Polar Coordinates. Copyright Cengage Learning. All rights reserved.
10 Conics, Parametric Equations, and Polar Coordinates Copyright Cengage Learning. All rights reserved. 10.5 Area and Arc Length in Polar Coordinates Copyright Cengage Learning. All rights reserved. Objectives
More informationSection Parametrized Surfaces and Surface Integrals. (I) Parametrizing Surfaces (II) Surface Area (III) Scalar Surface Integrals
Section 16.4 Parametrized Surfaces and Surface Integrals (I) Parametrizing Surfaces (II) Surface Area (III) Scalar Surface Integrals MATH 127 (Section 16.4) Parametrized Surfaces and Surface Integrals
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationLecture 1: Turtle Graphics. the turtle and the crane and the swallow observe the time of their coming; Jeremiah 8:7
Lecture 1: Turtle Graphics the turtle and the crane and the sallo observe the time of their coming; Jeremiah 8:7 1. Turtle Graphics Motion generates geometry. The turtle is a handy paradigm for investigating
More informationQuadratic and cubic b-splines by generalizing higher-order voronoi diagrams
Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Yuanxin Liu and Jack Snoeyink Joshua Levine April 18, 2007 Computer Science and Engineering, The Ohio State University 1 / 24
More informationAssignment 1. Prolog to Problem 1. Two cylinders. ü Visualization. Problems by Branko Curgus
Assignment In[]:= Problems by Branko Curgus SetOptions $FrontEndSession, Magnification Prolog to Problem. Two cylinders In[]:= This is a tribute to a problem that I was assigned as an undergraduate student
More information2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into
2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into the viewport of the current application window. A pixel
More informationTrigonometric Functions. Copyright Cengage Learning. All rights reserved.
4 Trigonometric Functions Copyright Cengage Learning. All rights reserved. 4.7 Inverse Trigonometric Functions Copyright Cengage Learning. All rights reserved. What You Should Learn Evaluate and graph
More informationHomework 8. Due: Tuesday, March 31st, 2009
MATH 55 Applied Honors Calculus III Winter 9 Homework 8 Due: Tuesday, March 3st, 9 Section 6.5, pg. 54: 7, 3. Section 6.6, pg. 58:, 3. Section 6.7, pg. 66: 3, 5, 47. Section 6.8, pg. 73: 33, 38. Section
More information1. A model for spherical/elliptic geometry
Math 3329-Uniform Geometries Lecture 13 Let s take an overview of what we know. 1. A model for spherical/elliptic geometry Hilbert s axiom system goes with Euclidean geometry. This is the same geometry
More informationGlossary of dictionary terms in the AP geometry units
Glossary of dictionary terms in the AP geometry units affine linear equation: an equation in which both sides are sums of terms that are either a number times y or a number times x or just a number [SlL2-D5]
More informationGergonne and Nagel Points for Simplices in the n-dimensional Space
Journal for Geometry and Graphics Volume 4 (2000), No. 2, 119 127. Gergonne and Nagel Points for Simplices in the n-dimensional Space Edwin Koźniewsi 1, Renata A. Górsa 2 1 Institute of Civil Engineering,
More information7. The Gauss-Bonnet theorem
7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed
More informationUse of n-vector for Radar Applications
Use of n-vector for Radar Applications Nina Ødegaard, Kenneth Gade Norwegian Defence Research Establishment Kjeller, NORWAY email: Nina.Odegaard@ffi.no Kenneth.Gade@ffi.no Abstract: This paper aims to
More informationChapter 15: Functions of Several Variables
Chapter 15: Functions of Several Variables Section 15.1 Elementary Examples a. Notation: Two Variables b. Example c. Notation: Three Variables d. Functions of Several Variables e. Examples from the Sciences
More informationSpatial Indexing of Global Geographical Data with HTM
Spatial Indexing of Global Geographical Data ith HTM Zhenhua Lv Yingjie Hu Haidong Zhong Bailang Yu Jianping Wu* Key Laboratory of Geographic Information Science Ministry of Education East China Normal
More informationA1:Orthogonal Coordinate Systems
A1:Orthogonal Coordinate Systems A1.1 General Change of Variables Suppose that we express x and y as a function of two other variables u and by the equations We say that these equations are defining a
More informationDouble Integrals, Iterated Integrals, Cross-sections
Chapter 14 Multiple Integrals 1 ouble Integrals, Iterated Integrals, Cross-sections 2 ouble Integrals over more general regions, efinition, Evaluation of ouble Integrals, Properties of ouble Integrals
More informationParallel Monte Carlo Sampling Scheme for Sphere and Hemisphere
Parallel Monte Carlo Sampling Scheme for Sphere and Hemisphere I.T. Dimov 1,A.A.Penzov 2, and S.S. Stoilova 3 1 Institute for Parallel Processing, Bulgarian Academy of Sciences Acad. G. Bonchev Str., bl.
More informationMathematics in Orbit
Mathematics in Orbit Dan Kalman American University Slides and refs at www.dankalman.net Outline Basics: 3D geospacial models Keyhole Problem: Related Rates! GPS: space-time triangulation Sensor Diagnosis:
More information2 Geometry Solutions
2 Geometry Solutions jacques@ucsd.edu Here is give problems and solutions in increasing order of difficulty. 2.1 Easier problems Problem 1. What is the minimum number of hyperplanar slices to make a d-dimensional
More information5/27/12. Objectives. Plane Curves and Parametric Equations. Sketch the graph of a curve given by a set of parametric equations.
Objectives Sketch the graph of a curve given by a set of parametric equations. Eliminate the parameter in a set of parametric equations. Find a set of parametric equations to represent a curve. Understand
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationV = 2πx(1 x) dx. x 2 dx. 3 x3 0
Wednesday, September 3, 215 Page 462 Problem 1 Problem. Use the shell method to set up and evaluate the integral that gives the volume of the solid generated by revolving the region (y = x, y =, x = 2)
More information751 Problem Set I JWR. Due Sep 28, 2004
751 Problem Set I JWR Due Sep 28, 2004 Exercise 1. For any space X define an equivalence relation by x y iff here is a path γ : I X with γ(0) = x and γ(1) = y. The equivalence classes are called the path
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationWHAT YOU SHOULD LEARN
GRAPHS OF EQUATIONS WHAT YOU SHOULD LEARN Sketch graphs of equations. Find x- and y-intercepts of graphs of equations. Use symmetry to sketch graphs of equations. Find equations of and sketch graphs of
More informationGrad operator, triple and line integrals. Notice: this material must not be used as a substitute for attending the lectures
Grad operator, triple and line integrals Notice: this material must not be used as a substitute for attending the lectures 1 .1 The grad operator Let f(x 1, x,..., x n ) be a function of the n variables
More informationWorksheet 3.5: Triple Integrals in Spherical Coordinates. Warm-Up: Spherical Coordinates (ρ, φ, θ)
Boise State Math 275 (Ultman) Worksheet 3.5: Triple Integrals in Spherical Coordinates From the Toolbox (what you need from previous classes) Know what the volume element dv represents. Be able to find
More informationTopic: Orientation, Surfaces, and Euler characteristic
Topic: Orientation, Surfaces, and Euler characteristic The material in these notes is motivated by Chapter 2 of Cromwell. A source I used for smooth manifolds is do Carmo s Riemannian Geometry. Ideas of
More informationMaths Year 11 Mock Revision list
Maths Year 11 Mock Revision list F = Foundation Tier = Foundation and igher Tier = igher Tier Number Tier Topic know and use the word integer and the equality and inequality symbols use fractions, decimals
More informationform are graphed in Cartesian coordinates, and are graphed in Cartesian coordinates.
Plot 3D Introduction Plot 3D graphs objects in three dimensions. It has five basic modes: 1. Cartesian mode, where surfaces defined by equations of the form are graphed in Cartesian coordinates, 2. cylindrical
More informationMATHEMATICS 105 Plane Trigonometry
Chapter I THE TRIGONOMETRIC FUNCTIONS MATHEMATICS 105 Plane Trigonometry INTRODUCTION The word trigonometry literally means triangle measurement. It is concerned with the measurement of the parts, sides,
More informationChapter 18. Geometric Operations
Chapter 18 Geometric Operations To this point, the image processing operations have computed the gray value (digital count) of the output image pixel based on the gray values of one or more input pixels;
More informationChapter 1. Turtle Graphics. 1.1 Turtle Graphics. The turtle and the crane and the swallow observe the time of their coming Jeremiah 8:7
Goldman/An Integrated Introduction to Computer Graphics and Geometric Modeling K10188_C001 Revise Proof page 3 26.3.2009 7:54am Compositor Name: VAmoudavally Chapter 1 Turtle Graphics The turtle and the
More informationIn this chapter, we will investigate what have become the standard applications of the integral:
Chapter 8 Overview: Applications of Integrals Calculus, like most mathematical fields, began with trying to solve everyday problems. The theory and operations were formalized later. As early as 70 BC,
More informationNormals of subdivision surfaces and their control polyhedra
Computer Aided Geometric Design 24 (27 112 116 www.elsevier.com/locate/cagd Normals of subdivision surfaces and their control polyhedra I. Ginkel a,j.peters b,,g.umlauf a a University of Kaiserslautern,
More informationMATHEMATICS FOR ENGINEERING TRIGONOMETRY
MATHEMATICS FOR ENGINEERING TRIGONOMETRY TUTORIAL SOME MORE RULES OF TRIGONOMETRY This is the one of a series of basic tutorials in mathematics aimed at beginners or anyone wanting to refresh themselves
More information9-1 GCSE Maths. GCSE Mathematics has a Foundation tier (Grades 1 5) and a Higher tier (Grades 4 9).
9-1 GCSE Maths GCSE Mathematics has a Foundation tier (Grades 1 5) and a Higher tier (Grades 4 9). In each tier, there are three exams taken at the end of Year 11. Any topic may be assessed on each of
More informationComputational Geometry
Computational Geometry 600.658 Convexity A set S is convex if for any two points p, q S the line segment pq S. S p S q Not convex Convex? Convexity A set S is convex if it is the intersection of (possibly
More informationData Representation in Visualisation
Data Representation in Visualisation Visualisation Lecture 4 Taku Komura Institute for Perception, Action & Behaviour School of Informatics Taku Komura Data Representation 1 Data Representation We have
More informationTHE BASIC TOPOLOGY ISSUES ON SPHERICAL SURFACE QUTERNARY TRANGULAR MESH
THE BASIC TOPOLOGY ISSUES ON SPHERICAL SURFACE QUTERNARY TRANGULAR MESH HOU Miaole a, ZHAO Xuesheng b, c, CHEN Jun c a Department of Surveying & Spatial Science, BICEA (Beijing), 1 Zhanlanguan Road, Beijing,
More informationGY403 Structural Geology Laboratory. Example Orthographic Solutions to Apparent Dip, True Dip, 3-Point, 3 Intersecting Planes Problems
GY403 Structural Geology Laboratory Example Orthographic Solutions to Apparent Dip, True Dip, 3-Point, 3 and Intersecting Planes Problems Orthographic Example of Apparent Dip Problem Orthographic solutions
More information16.6 Parametric Surfaces and Their Areas
SECTION 6.6 PARAMETRIC SURFACES AND THEIR AREAS i j k (b) From (a), v = w r = =( ) i +( ) j +( ) k = i + j i j k (c) curl v = v = = () () i + ( ) () j + () ( ) k =[ ( )] k = k =w 9. For any continuous
More informationMath 170, Section 002 Spring 2012 Practice Exam 2 with Solutions
Math 170, Section 002 Spring 2012 Practice Exam 2 with Solutions Contents 1 Problems 2 2 Solution key 10 3 Solutions 11 1 1 Problems Question 1: A right triangle has hypothenuse of length 25 in and an
More informationarxiv: v1 [math.co] 17 Jan 2014
Regular matchstick graphs Sascha Kurz Fakultät für Mathematik, Physik und Informatik, Universität Bayreuth, Germany Rom Pinchasi Mathematics Dept., Technion Israel Institute of Technology, Haifa 2000,
More informationDTM Based on an Ellipsoidal Squares
DM Based on an Ellipsoidal Squares KRZYSZOF NAUS Institute of Navigation and Hydrography Polish Naval Academy Śmidoicza 69, 8-3 Gdynia POLAND k.naus@am.gdynia.pl Abstract: - he paper presents the description
More informationMATH 234. Excercises on Integration in Several Variables. I. Double Integrals
MATH 234 Excercises on Integration in everal Variables I. Double Integrals Problem 1. D = {(x, y) : y x 1, 0 y 1}. Compute D ex3 da. Problem 2. Find the volume of the solid bounded above by the plane 3x
More informationMath Boot Camp: Coordinate Systems
Math Boot Camp: Coordinate Systems You can skip this boot camp if you can answer the following question: Staying on a sphere of radius R, what is the shortest distance between the point (0, 0, R) on the
More informationChapter 6 Some Applications of the Integral
Chapter 6 Some Applications of the Integral More on Area More on Area Integrating the vertical separation gives Riemann Sums of the form More on Area Example Find the area A of the set shaded in Figure
More informationA Flavor of Topology. Shireen Elhabian and Aly A. Farag University of Louisville January 2010
A Flavor of Topology Shireen Elhabian and Aly A. Farag University of Louisville January 2010 In 1670 s I believe that we need another analysis properly geometric or linear, which treats place directly
More information9. Three Dimensional Object Representations
9. Three Dimensional Object Representations Methods: Polygon and Quadric surfaces: For simple Euclidean objects Spline surfaces and construction: For curved surfaces Procedural methods: Eg. Fractals, Particle
More informationFinal Exam 1:15-3:15 pm Thursday, December 13, 2018
Final Exam 1:15-3:15 pm Thursday, December 13, 2018 Instructions: Answer all questions in the space provided (or attach additional pages as needed). You are permitted to use pencils/pens, one cheat sheet
More informationVisualisation Pipeline : The Virtual Camera
Visualisation Pipeline : The Virtual Camera The Graphics Pipeline 3D Pipeline The Virtual Camera The Camera is defined by using a parallelepiped as a view volume with two of the walls used as the near
More informationf xx (x, y) = 6 + 6x f xy (x, y) = 0 f yy (x, y) = y In general, the quantity that we re interested in is
1. Let f(x, y) = 5 + 3x 2 + 3y 2 + 2y 3 + x 3. (a) Final all critical points of f. (b) Use the second derivatives test to classify the critical points you found in (a) as a local maximum, local minimum,
More informationMATH 19520/51 Class 15
MATH 19520/51 Class 15 Minh-Tam Trinh University of Chicago 2017-11-01 1 Change of variables in two dimensions. 2 Double integrals via change of variables. Change of Variables Slogan: An n-variable substitution
More informationApplications of Triple Integrals
Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 27 / 45
: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ) Chapter 10: Parametric Equations
More informationPS Computational Geometry Homework Assignment Sheet I (Due 16-March-2018)
Homework Assignment Sheet I (Due 16-March-2018) Assignment 1 Let f, g : N R with f(n) := 8n + 4 and g(n) := 1 5 n log 2 n. Prove explicitly that f O(g) and f o(g). Assignment 2 How can you generalize the
More information( ) Derivation of Polar Reduction Formula for a Calculator Robert Bernecky April, 2018 ( )
Derivation of Polar Reduction Formula for a Calculator Robert Bernecky April, 2018 1 Problem Statement The polar reduction formula takes an observer's assumed position (lat, lon), and a body's celestial
More informationMATH 2023 Multivariable Calculus
MATH 2023 Multivariable Calculus Problem Sets Note: Problems with asterisks represent supplementary informations. You may want to read their solutions if you like, but you don t need to work on them. Set
More informationSPECIAL TECHNIQUES-II
SPECIAL TECHNIQUES-II Lecture 19: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Method of Images for a spherical conductor Example :A dipole near aconducting sphere The
More informationAQA GCSE Maths - Higher Self-Assessment Checklist
AQA GCSE Maths - Higher Self-Assessment Checklist Number 1 Use place value when calculating with decimals. 1 Order positive and negative integers and decimals using the symbols =,, , and. 1 Round to
More informationHow to Project Spherical Conics into the Plane
How to Project pherical Conics into the Plane Yoichi Maeda maeda@keyakiccu-tokaiacjp Department of Mathematics Tokai University Japan Abstract: In this paper, we will introduce a method how to draw the
More informationG 6i try. On the Number of Minimal 1-Steiner Trees* Discrete Comput Geom 12:29-34 (1994)
Discrete Comput Geom 12:29-34 (1994) G 6i try 9 1994 Springer-Verlag New York Inc. On the Number of Minimal 1-Steiner Trees* B. Aronov, 1 M. Bern, 2 and D. Eppstein 3 Computer Science Department, Polytechnic
More informationRELATIONSHIP BETWEEN ASTRONOMIC COORDINATES φ, λ, AND GEODETIC COORDINATES φ G, λ G,
RELTIOSHIP BETWEE STROOMIC COORDITES φ, λ, D EODETIC COORDITES φ, λ, h H In geodesy it is important to know the relationships between observed quantities such as horizontal directions (or azimuths) and
More informationStructured Light Based Depth Edge Detection for Object Shape Recovery
Structured Light Based Depth Edge Detection for Object Shape Recovery Cheolhon Kim Jiyoung Park Juneho Yi Matthe Turk School of Information and Communication Engineering Sungkyunkan University, Korea Biometrics
More information4.2 Description of surfaces by spherical harmonic functions
Chapter 4. Parametrization of closed curves and surfaces Im[z] Im[z] Translation Im[z] Im[z] Rotation Scale Starting point Re[z] Re[z] Re[z] Re[z] a b c d Figure 4.: Normalization steps of Fourier coefficients;
More informationTHE BASIC TOPOLOGY MODEL OF SPHERICAL SURFACE DIGITAL SPACE
HE BASIC OPOLOGY MODEL OF SPHERICAL SURFACE DIGIAL SPACE HOU Miao-le 1 1 Department of Surveying & Land Science, China University of Mining and echnology (Beijing), D11 Xueyuan Road, Beijing, China, 100083
More informationMoore Catholic High School Math Department
Moore Catholic High School Math Department Geometry Vocabulary The following is a list of terms and properties which are necessary for success in a Geometry class. You will be tested on these terms during
More information