Computation of the gravity gradient tensor due to topographic masses using tesseroids
|
|
- Cora Hutchinson
- 5 years ago
- Views:
Transcription
1 Computation of the gravity gradient tensor due to topographic masses using tesseroids Leonardo Uieda 1 Naomi Ussami 2 Carla F Braitenberg 3 1. Observatorio Nacional, Rio de Janeiro, Brazil 2. Universidade de São Paulo, São Paulo, Brazil 3. University of Trieste, Trieste, Italy. August 9, 2010
2 Outline The Gravity Gradient Tensor (GGT) What is a tesseroid Why use tesseroids Numerical issues Modeling topography with tesseroids Topographic effect in the Paraná Basin region Further applications Concluding remarks
3 Gravity Gradient Tensor
4 Gravity Gradient Tensor Hessian matrix of gravitational potential
5 Gravity Gradient Tensor Hessian matrix of gravitational potential 2 V x 2 g xx g xy g xz GGT = g yx g yy g yz = 2 V g zx g zy g zz y x 2 V z x 2 V x y 2 V y 2 2 V z y 2 V x z 2 V y z 2 V z 2
6 Gravity Gradient Tensor Hessian matrix of gravitational potential 2 V x 2 g xx g xy g xz GGT = g yx g yy g yz = 2 V g zx g zy g zz y x 2 V z x Volume integrals 2 V x y 2 V y 2 2 V z y 2 V x z 2 V y z 2 V z 2
7 Gravity Gradient Tensor Hessian matrix of gravitational potential 2 V x 2 g xx g xy g xz GGT = g yx g yy g yz = 2 V g zx g zy g zz y x 2 V z x Volume integrals 2 V x y 2 V y 2 2 V z y ˆ g ij (x, y, z) = Kernel(x, y, z, x, y, z ) dω Ω 2 V x z 2 V y z 2 V z 2
8 Gravity Gradient Tensor Can discretize volume Ω using:
9 Gravity Gradient Tensor Can discretize volume Ω using: Rectangular prisms
10 Gravity Gradient Tensor Can discretize volume Ω using: Rectangular prisms Tesseroids (spherical prisms)
11 What is a tesseroid?
12 What is a tesseroid? Z Tesseroid r X λ φ Y
13 What is a tesseroid? Delimited by: Z 2 meridians r X λ 1 λ φ Y
14 What is a tesseroid? Delimited by: Z 2 meridians r X λ 2 λ φ Y
15 What is a tesseroid? Delimited by: Z 2 meridians 2 parallels φ 1 r X λ φ Y
16 What is a tesseroid? Delimited by: Z 2 meridians 2 parallels φ 2 r X λ φ Y
17 What is a tesseroid? Delimited by: Z 2 meridians 2 parallels 2 concentric spheres r 1 r X λ φ Y
18 What is a tesseroid? Delimited by: Z 2 meridians 2 parallels 2 concentric spheres r 2 r X λ φ Y
19 Why use tesseroids?
20 Why use tesseroids? Earth Core Matle Crust
21 Why use tesseroids? Earth Core Matle Crust
22 Why use tesseroids? Want to model the geologic body Observation Point Geologic body
23 Why use tesseroids? Flat Earth Observation Point
24 Why use tesseroids? Flat Earth + Rectangular Prisms Observation Point
25 Why use tesseroids? Good for small regions (Rule of thumb: < 2500 km) Flat Earth + Rectangular Prisms Observation Point
26 Why use tesseroids? Good for small regions (Rule of thumb: < 2500 km) Flat Earth + Rectangular Prisms Observation Point and close observation point
27 Why use tesseroids? Good for small regions (Rule of thumb: < 2500 km) Flat Earth + Rectangular Prisms Observation Point and close observation point Not very accurate for larger regions
28 Why use tesseroids? Spherical Earth Observation Point
29 Why use tesseroids? Spherical Earth + Rectangular Prisms Observation Point
30 Why use tesseroids? Spherical Earth + Rectangular Prisms Observation Point
31 Why use tesseroids? Usually accurate enough (if mass of prisms = mass of tesseroids) Spherical Earth + Rectangular Prisms Observation Point
32 Why use tesseroids? Usually accurate enough (if mass of prisms = mass of tesseroids) Involves many coordinate changes Spherical Earth + Rectangular Prisms Observation Point
33 Why use tesseroids? Usually accurate enough (if mass of prisms = mass of tesseroids) Involves many coordinate changes Spherical Earth + Rectangular Prisms Observation Point Computationally slow
34 Why use tesseroids? Spherical Earth Observation Point
35 Why use tesseroids? Spherical Earth + Tesseroids Observation Point
36 Why use tesseroids? Spherical Earth + Tesseroids Observation Point
37 Why use tesseroids? As accurate as Spherical Earth + rectangular prisms Spherical Earth + Tesseroids Observation Point
38 Why use tesseroids? As accurate as Spherical Earth + rectangular prisms But faster Spherical Earth + Tesseroids Observation Point
39 Why use tesseroids? As accurate as Spherical Earth + rectangular prisms But faster Spherical Earth + Tesseroids Observation Point As shown in Wild-Pfeiffer (2008)
40 Why use tesseroids? As accurate as Spherical Earth + rectangular prisms But faster Spherical Earth + Tesseroids Observation Point As shown in Wild-Pfeiffer (2008) Some numerical problems
41 Numerical issues
42 Numerical issues Gravity Gradient Tensor (GGT) volume integrals solved:
43 Numerical issues Gravity Gradient Tensor (GGT) volume integrals solved: Analytically in the radial direction
44 Numerical issues Gravity Gradient Tensor (GGT) volume integrals solved: Analytically in the radial direction Numerically over the surface of the sphere
45 Numerical issues Gravity Gradient Tensor (GGT) volume integrals solved: Analytically in the radial direction Numerically over the surface of the sphere Using the Gauss-Legendre Quadrature (GLQ)
46 Numerical issues At 250 km height with Gauss-Legendre Quadrature (GLQ) order 2
47 Numerical issues At 50 km height with Gauss-Legendre Quadrature (GLQ) order 2
48 Numerical issues At 50 km height with Gauss-Legendre Quadrature (GLQ) order 10
49 Numerical issues General rule:
50 Numerical issues General rule: Distance to computation point > Distance between nodes
51 Numerical issues General rule: Distance to computation point > Distance between nodes Increase number of nodes
52 Numerical issues General rule: Distance to computation point > Distance between nodes Increase number of nodes Divide the tesseroid in smaller parts
53 Modeling topography with tesseroids
54 Modeling topography with tesseroids Computer program: Tesseroids
55 Modeling topography with tesseroids Computer program: Tesseroids Python programming language
56 Modeling topography with tesseroids Computer program: Tesseroids Python programming language Open Source (GNU GPL License)
57 Modeling topography with tesseroids Computer program: Tesseroids Python programming language Open Source (GNU GPL License) Project hosted on Google Code
58 Modeling topography with tesseroids Computer program: Tesseroids Python programming language Open Source (GNU GPL License) Project hosted on Google Code
59 Modeling topography with tesseroids Computer program: Tesseroids Python programming language Open Source (GNU GPL License) Project hosted on Google Code Under development:
60 Modeling topography with tesseroids Computer program: Tesseroids Python programming language Open Source (GNU GPL License) Project hosted on Google Code Under development: Optimizations using C coded modules
61 Modeling topography with tesseroids To model topography:
62 Modeling topography with tesseroids To model topography: Digital Elevation Model (DEM) Tesseroid model
63 Modeling topography with tesseroids To model topography: Digital Elevation Model (DEM) Tesseroid model 1 Grid Point = 1 Tesseroid
64 Modeling topography with tesseroids To model topography: Digital Elevation Model (DEM) Tesseroid model 1 Grid Point = 1 Tesseroid Top centered on grid point
65 Modeling topography with tesseroids To model topography: Digital Elevation Model (DEM) Tesseroid model 1 Grid Point = 1 Tesseroid Top centered on grid point Bottom at reference surface
66 Topographic effect in the Paraná Basin region
67 Topographic effect in the Paraná Basin region Digital Elevation Model (DEM) Grid:
68 Topographic effect in the Paraná Basin region Digital Elevation Model (DEM) Grid: ETOPO1
69 Topographic effect in the Paraná Basin region Digital Elevation Model (DEM) Grid: ETOPO1 10 x 10 Grid
70 Topographic effect in the Paraná Basin region Digital Elevation Model (DEM) Grid: ETOPO1 10 x 10 Grid ~ 23,000 Tesseroids
71 Topographic effect in the Paraná Basin region Digital Elevation Model (DEM) Grid: ETOPO1 10 x 10 Grid ~ 23,000 Tesseroids Density = 2.67 g cm 3
72 Topographic effect in the Paraná Basin region Digital Elevation Model (DEM) Grid: ETOPO1 10 x 10 Grid ~ 23,000 Tesseroids Density = 2.67 g cm 3 Computation height = 250 km
73 Topographic effect in the Paraná Basin region 5 S 10 S 15 S 20 S 25 S 30 S 35 S 40 S 75 W 70 W 65 W 60 W 55 W 50 W 45 W 40 W 35 W Height [m]
74 Topographic effect in the Paraná Basin region Height of 250 km
75 Topographic effect in the Paraná Basin region Topographic effect in the region has the same order of magnitude as a km tesseroid (10 0 Eötvös)
76 Topographic effect in the Paraná Basin region Topographic effect in the region has the same order of magnitude as a km tesseroid (10 0 Eötvös) Need to take topography into account when modeling (even at 250 km altitudes)
77 Further applications
78 Further applications Satellite gravity data = global coverage
79 Further applications Satellite gravity data = global coverage + Tesseroid modeling:
80 Further applications Satellite gravity data = global coverage + Tesseroid modeling: Regional/global inversion for density (Mantle)
81 Further applications Satellite gravity data = global coverage + Tesseroid modeling: Regional/global inversion for density (Mantle) Regional/global inversion for relief of an interface (Moho)
82 Further applications Satellite gravity data = global coverage + Tesseroid modeling: Regional/global inversion for density (Mantle) Regional/global inversion for relief of an interface (Moho) Joint inversion with seismic tomography
83 Concluding remarks
84 Concluding remarks Developed a computational tool for large-scale gravity modeling with tesseroids
85 Concluding remarks Developed a computational tool for large-scale gravity modeling with tesseroids Better use tesseroids than rectangular prisms for large regions
86 Concluding remarks Developed a computational tool for large-scale gravity modeling with tesseroids Better use tesseroids than rectangular prisms for large regions Take topographic effect into consideration when modeling density anomalies within the Earth
87 Concluding remarks Developed a computational tool for large-scale gravity modeling with tesseroids Better use tesseroids than rectangular prisms for large regions Take topographic effect into consideration when modeling density anomalies within the Earth Possible application: tesseroids in regional/global gravity inversion
88 Thank you
89 References WILD-PFEIFFER, F. A comparison of different mass elements for use in gravity gradiometry. Journal of Geodesy, v. 82 (10), p , 2008.
Robust 3D gravity gradient inversion by planting anomalous densities
Robust 3D gravity gradient inversion by planting anomalous densities Leonardo Uieda Valéria C. F. Barbosa Observatório Nacional September, 2011 Outline Outline Forward Problem Outline Forward Problem Inverse
More informationprismatic discretization of the digital elevation model, and their associated volume integration problems. Summary
A new method of terrain correcting airborne gravity gradiometry data using 3D Cauchy-type integrals Michael S. Zhdanov*, University of Utah and TechnoImaging, Glenn A. Wilson, TechnoImaging, and Xiaojun
More informationDownloaded 18 Jul 2011 to Redistribution subject to SEG license or copyright; see Terms of Use at
Rapid gravity and gravity gradiometry terrain correction via adaptive quadtree mesh discretization Kristofer Davis, M. Andy Kass, and Yaoguo Li, Center for Gravity, Electrical and Magnetic Studies, Colorado
More informationSection 4.2 selected answers Math 131 Multivariate Calculus D Joyce, Spring 2014
4. Determine the nature of the critical points of Section 4. selected answers Math 11 Multivariate Calculus D Joyce, Spring 014 Exercises from section 4.: 6, 1 16.. Determine the nature of the critical
More informationDifferent topographic reduction techniques for GOCE gravity data
University of Kiel Different topographic reduction techniques for GOCE gravity data N. Koether W. Szwillius H.-J. Goetze Motivation Introduction to Bouguer Anomalies Overview about algorithms Tritop Polyhedral
More informationEMIGMA V9.x Premium Series April 8, 2015
EMIGMA V9.x Premium Series April 8, 2015 EMIGMA for Gravity EMIGMA for Gravity license is a comprehensive package that offers a wide array of processing, visualization and interpretation tools. The package
More informationEfficient 3D Gravity and Magnetic Modeling
Efficient 3D Gravity and Magnetic Modeling X. Li Fugro Gravity & Magnetic Services Inc., Houston, Texas, USA Summary There are many different spatial-domain algorithms for 3D gravity and magnetic forward
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 8: Joint Geophysical Inversions Wed, April 20, 2011
GEOPHYS 4: Near Surface Geophysical Imaging Class 8: Joint Geophysical Inversions Wed, April, 11 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical,
More informationCUDA Parallel Algorithms for Forward and Inverse Structural Gravity Problems
CUDA Parallel Algorithms for Forward and Inverse Structural Gravity Problems Alexander Tsidaev Bulashevich Institute of Geophysics, Yekaterinburg, Russia pdc@tsidaev.ru Abstract. This paper describes usage
More informationObjectives for Terrain Week
The Terrain Correction Chuck Connor, Laura Connor Potential Fields Geophysics: Terrain Week for Terrain Week Learn about the terrain correction The inner terrain correction Learn about The outer terrain
More information(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1
Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same
More informationEvaluation of precise terrain effects using high-resolution digital elevation models
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb005639, 2009 Evaluation of precise terrain effects using high-resolution digital elevation models Dimitrios Tsoulis, 1 Pavel Novák, 2,3 and
More informationIntroduction to Computer Graphics. Modeling (3) April 27, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (3) April 27, 2017 Kenshi Takayama Solid modeling 2 Solid models Thin shapes represented by single polygons Unorientable Clear definition of inside & outside
More informationDOWN PLUNGE CROSS SECTIONS
GG303 Lab 7 10/6/10 1 DOWN PLUNGE CROSS SECTIONS I Main Topics A Cylindrical folds B Downplunge cross-section views C Apparent dip II Cylindrical folds A Surface of a cylindrical fold is parallel to a
More informationFundamentals of Structural Geology Exercise: concepts from chapter 2
0B Reading: Fundamentals of Structural Geology, Ch 2 1) Develop a MATLAB script that plots the spherical datum (Fig. 2.1a) with unit radius as a wire-frame diagram using lines of constant latitude and
More information) in the k-th subbox. The mass of the k-th subbox is M k δ(x k, y k, z k ) V k. Thus,
1 Triple Integrals Mass problem. Find the mass M of a solid whose density (the mass per unit volume) is a continuous nonnegative function δ(x, y, z). 1. Divide the box enclosing into subboxes, and exclude
More informationSpectral-Element and Adjoint Methods in Seismology
Spectral-Element and Adjoint Methods in Seismology Jeroen Tromp Min Chen, Vala Hjorleifsdottir, Dimitri Komatitsch, Swami Krishnan, Qinya Liu, Alessia Maggi, Anne Sieminski, Carl Tape & Ying Zhou Introduction
More information21-256: Lagrange multipliers
21-256: Lagrange multipliers Clive Newstead, Thursday 12th June 2014 Lagrange multipliers give us a means of optimizing multivariate functions subject to a number of constraints on their variables. Problems
More informationSummer 2017 MATH Suggested Solution to Exercise Find the tangent hyperplane passing the given point P on each of the graphs: (a)
Smmer 2017 MATH2010 1 Sggested Soltion to Exercise 6 1 Find the tangent hyperplane passing the given point P on each of the graphs: (a) z = x 2 y 2 ; y = z log x z P (2, 3, 5), P (1, 1, 1), (c) w = sin(x
More informationHow to Use GOCE Level 2 Products
How to Use GOCE Level 2 Products Thomas Gruber 1), Reiner Rummel 1), Radboud Koop 2) 1) Institute of Astronomical and Physical Geodesy, Technical University Munich 2) Netherlands Institute for Space Research
More information(1) Given the following system of linear equations, which depends on a parameter a R, 3x y + 5z = 2 4x + y + (a 2 14)z = a + 2
(1 Given the following system of linear equations, which depends on a parameter a R, x + 2y 3z = 4 3x y + 5z = 2 4x + y + (a 2 14z = a + 2 (a Classify the system of equations depending on the values of
More informationGeophysics 224 B2. Gravity anomalies of some simple shapes. B2.1 Buried sphere
Geophysics 4 B. Gravity anomalies of some simple shapes B.1 Buried sphere Gravity measurements are made on a surface profile across a buried sphere. The sphere has an excess mass M S and the centre is
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 informationVector Calculus: Understanding the Cross Product
University of Babylon College of Engineering Mechanical Engineering Dept. Subject : Mathematics III Class : 2 nd year - first semester Date: / 10 / 2016 2016 \ 2017 Vector Calculus: Understanding the Cross
More informationHierarchical Hybrid Grids
Hierarchical Hybrid Grids IDK Summer School 2012 Björn Gmeiner, Ulrich Rüde July, 2012 Contents Mantle convection Hierarchical Hybrid Grids Smoothers Geometric approximation Performance modeling 2 Mantle
More informationUnstructured grid modelling to create 3-D Earth models that unify geological and geophysical information
Unstructured grid modelling to create 3-D Earth models that unify geological and geophysical information Peter Lelièvre, Angela Carter-McAuslan, Cassandra Tycholiz, Colin Farquharson and Charles Hurich
More informationChristian Hirt 1,2 *, Elisabeth Reußner 1, Moritz Rexer 1,2 and Michael Kuhn 3. Technische Universität München, Germany
Citation: Hirt, C., E. Reußner, M. Rexer and M. Kuhn (2016), Topographic gravity modelling for global Bouguer maps to degree 2,160: Validation of spectral and spatial domain forward modelling techniques
More informationMath 113 Calculus III Final Exam Practice Problems Spring 2003
Math 113 Calculus III Final Exam Practice Problems Spring 23 1. Let g(x, y, z) = 2x 2 + y 2 + 4z 2. (a) Describe the shapes of the level surfaces of g. (b) In three different graphs, sketch the three cross
More informationRemoving Noise due to Aircraft Motion from Gravity Gradient Data. Wolobah B. Sali
Removing Noise due to Aircraft Motion from Gravity Gradient Data Wolobah B. Sali Declaration I declare that this thesis is my own, unaided work. It is being submitted for the degree of Master of science
More informationMath 233. Lagrange Multipliers Basics
Math 233. Lagrange Multipliers Basics Optimization problems of the form to optimize a function f(x, y, z) over a constraint g(x, y, z) = k can often be conveniently solved using the method of Lagrange
More informationGravity Data Correction in QCTool
Gravity Data Correction in QCTool January 20, 2017 info@qc-tool.com 2 1. Importing Data You must import to QCTool measured gravity data from ASCII files. Then you need to import the data from CSV files
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 informationNavigation coordinate systems
Lecture 3 Navigation coordinate systems Topic items: 1. Basic Coordinate Systems. 2. Plane Cartesian Coordinate Systems. 3. Polar Coordinate Systems. 4. Earth-Based Locational Reference Systems. 5. Reference
More informationAdding Magnetic Data to a GM-SYS Profile Model
Adding Magnetic Data to a GM-SYS Profile Model In GM-SYS Profile, a "Station" is a location at which an anomaly component is calculated and, optionally, was measured. In order for GM-SYS Profile to calculate
More informationOasis montaj How-To Guide. VOXI Earth Modelling - Running an AGG Unconstrained Inversion
Oasis montaj How-To Guide VOXI Earth Modelling - Running an AGG Unconstrained Inversion The software described in this manual is furnished under license and may only be used or copied in accordance with
More informationA few multilinear algebraic definitions I Inner product A, B = i,j,k a ijk b ijk Frobenius norm Contracted tensor products A F = A, A C = A, B 2,3 = A
Krylov-type methods and perturbation analysis Berkant Savas Department of Mathematics Linköping University Workshop on Tensor Approximation in High Dimension Hausdorff Institute for Mathematics, Universität
More informationSpectral analysis of the Earth s topographic potential via 2D DFT a new data based degree variance model to degree 90,000
Citation: Rexer, M. and C. Hirt (2015), Spectral analysis of the Earth s topographic potential via 2D DFT a new data based degree variance model to degree 90,000. Journal of Geodesy 89(9), 887 909, doi:
More informationChapter 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 informationQuaternion properties: addition. Introduction to quaternions. Quaternion properties: multiplication. Derivation of multiplication
Introduction to quaternions Definition: A quaternion q consists of a scalar part s, s, and a vector part v ( xyz,,, v 3 : q where, [ s, v q [ s, ( xyz,, q s+ ix + jy + kz i 2 j 2 k 2 1 ij ji k k Quaternion
More information3.3 Optimizing Functions of Several Variables 3.4 Lagrange Multipliers
3.3 Optimizing Functions of Several Variables 3.4 Lagrange Multipliers Prof. Tesler Math 20C Fall 2018 Prof. Tesler 3.3 3.4 Optimization Math 20C / Fall 2018 1 / 56 Optimizing y = f (x) In Math 20A, we
More informationHw 4 Due Feb 22. D(fg) x y z (
Hw 4 Due Feb 22 2.2 Exercise 7,8,10,12,15,18,28,35,36,46 2.3 Exercise 3,11,39,40,47(b) 2.4 Exercise 6,7 Use both the direct method and product rule to calculate where f(x, y, z) = 3x, g(x, y, z) = ( 1
More information2 Second Derivatives. As we have seen, a function f (x, y) of two variables has four different partial derivatives: f xx. f yx. f x y.
2 Second Derivatives As we have seen, a function f (x, y) of two variables has four different partial derivatives: (x, y), (x, y), f yx (x, y), (x, y) It is convenient to gather all four of these into
More informationMath 2130 Practice Problems Sec Name. Change the Cartesian integral to an equivalent polar integral, and then evaluate.
Math 10 Practice Problems Sec 1.-1. Name Change the Cartesian integral to an equivalent polar integral, and then evaluate. 1) 5 5 - x dy dx -5 0 A) 5 B) C) 15 D) 5 ) 0 0-8 - 6 - x (8 + ln 9) A) 1 1 + x
More informationEffects of Hypothetical Complex Mass-Density Distributions on Geoidal Height
Effects of Hypothetical Complex Mass-Density Distributions on Geoidal Height 51 Robert Kingdon, Petr Vaníček, and Marcelo Santos Abstract Geoid computation according to the Stokes-Helmert scheme requires
More informationOPENRAVE TUTORIAL. Robot Autonomy Spring 2014
OPENRAVE TUTORIAL Robot Autonomy Spring 2014 OPENRAVE Stands for: Open Robotics Automation Virtual Environment All-in-one Robotics Package Contains: Kinematics (forward, inverse, velocity, etc.) Collision
More informationExercise Lab: Where is the Himalaya eroding? Using GIS/DEM analysis to reconstruct surfaces, incision, and erosion
Exercise Lab: Where is the Himalaya eroding? Using GIS/DEM analysis to reconstruct surfaces, incision, and erosion 1) Start ArcMap and ensure that the 3D Analyst and the Spatial Analyst are loaded and
More informationCOORDINATE TRANSFORMATION. Lecture 6
COORDINATE TRANSFORMATION Lecture 6 SGU 1053 SURVEY COMPUTATION 1 Introduction Geomatic professional are mostly confronted in their work with transformations from one two/three-dimensional coordinate system
More informationAnimation Curves and Splines 2
Animation Curves and Splines 2 Animation Homework Set up Thursday a simple avatar E.g. cube/sphere (or square/circle if 2D) Specify some key frames (positions/orientations) Associate Animation a time with
More information7.3 FIRST INVERSE ANTECEDENT BAF(BAF) OPERATION
7.3 FIRST INVERSE ANTECEDENT BAFBAF OPERATION The inverse of the problem of the Direct BAFBAF Operation of the precedent Sections, could be also considered as a problem of Boolean Equations. Thus, if we
More informationTLS Parameters, Workflows and Field Methods
TLS Parameters, Workflows and Field Methods Marianne Okal, UNAVCO June 20 th, 2014 How a Lidar instrument works (Recap) Transmits laser signals and measures the reflected light to create 3D point clouds.
More informationIn-plane principal stress output in DIANA
analys: linear static. class: large. constr: suppor. elemen: hx24l solid tp18l. load: edge elemen force node. materi: elasti isotro. option: direct. result: cauchy displa princi stress total. In-plane
More informationAn introduction to the GOCE error variance covariance products: A user s perspective. Rory Bingham Newcastle University
An introduction to the GOCE error variance covariance products: A user s perspective Rory Bingham Newcastle University rory.bingham@ncl.ac.uk The ESA GOCE Virtual Archive The variance-covariance matrices
More informationNUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS
- 465 - NUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS LIU, Z. Q. SUN, J. G. * SUN, H. * LIU, M. C. * GAO, Z. H. College for Geoexploration Science
More informationLithological and surface geometry joint inversions using multi-objective global optimization methods
Lithological and surface geometry joint inversions using multi-objective global optimization methods Peter G. Lelièvre 1, Rodrigo Bijani and Colin G. Farquharson 1 plelievre@mun.ca http://www.esd.mun.ca/~peter/home.html
More informationA Wavelet-Based Method for Simulation of Seismic Wave Propagation
AGU Fall Meeting, A Wavelet-Based Method for Simulation of Seismic Wave Propagation Tae-Kyung Hong & B.L.N. Kennett Research School of Earth Sciences The Australian National University Abstract Seismic
More informationA simple example. Assume we want to find the change in the rotation angles to get the end effector to G. Effect of changing s
CENG 732 Computer Animation This week Inverse Kinematics (continued) Rigid Body Simulation Bodies in free fall Bodies in contact Spring 2006-2007 Week 5 Inverse Kinematics Physically Based Rigid Body Simulation
More informationSensor fusion for motion processing and visualization
Sensor fusion for motion processing and visualization Ali Baharev, PhD TÁMOP 4.2.2 Szenzorhálózat alapú adatgyűjtés és információfeldolgozás workshop April 1, 2011 Budapest, Hungary What we have - Shimmer
More informationThere are 10 problems, with a total of 150 points possible. (a) Find the tangent plane to the surface S at the point ( 2, 1, 2).
Instructions Answer each of the questions on your own paper, and be sure to show your work so that partial credit can be adequately assessed. Put your name on each page of your paper. You may use a scientific
More informationBernstein-Bezier Splines on the Unit Sphere. Victoria Baramidze. Department of Mathematics. Western Illinois University
Bernstein-Bezier Splines on the Unit Sphere Victoria Baramidze Department of Mathematics Western Illinois University ABSTRACT I will introduce scattered data fitting problems on the sphere and discuss
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2017) 209, 406 424 Advance Access publication 2017 January 11 GJI Gravity, geodesy and tides doi: 10.1093/gji/ggx010 A high-order 3-D spectral-element
More informationTLS Parameters, Workflows and Field Methods
TLS Parameters, Workflows and Field Methods Marianne Okal, UNAVCO GSA, October 20 th, 2017 How a Lidar instrument works (Recap) Transmits laser signals and measures the reflected light to create 3D point
More informationTesting Joint Inversion Code with Geologically Realistic Models
Testing Joint Inversion Code with Geologically Realistic Models Angela Carter-McAuslan, Peter Lelièvre, Gary Blades, Colin Farquharson GAC-MAC Joint Annual Meeting May 27 th, 2012 Background Joint Inversion
More informationPolSARpro v4.03 Forest Applications
PolSARpro v4.03 Forest Applications Laurent Ferro-Famil Lecture on polarimetric SAR Theory and applications to agriculture & vegetation Thursday 19 April, morning Pol-InSAR Tutorial Forest Application
More informationTLS Parameters, Workflows and Field Methods
TLS Parameters, Workflows and Field Methods Marianne Okal, UNAVCO GSA, September 23 rd, 2016 How a Lidar instrument works (Recap) Transmits laser signals and measures the reflected light to create 3D point
More informationWorksheet 3.4: Triple Integrals in Cylindrical Coordinates. Warm-Up: Cylindrical Volume Element d V
Boise State Math 275 (Ultman) Worksheet 3.4: Triple Integrals in Cylindrical Coordinates From the Toolbox (what you need from previous classes) Know what the volume element dv represents. Be able to find
More informationAnimating orientation. CS 448D: Character Animation Prof. Vladlen Koltun Stanford University
Animating orientation CS 448D: Character Animation Prof. Vladlen Koltun Stanford University Orientation in the plane θ (cos θ, sin θ) ) R θ ( x y = sin θ ( cos θ sin θ )( x y ) cos θ Refresher: Homogenous
More informationExperimental evaluation of static stiffness of a spatial translational parallel manipulator.
Experimental evaluation of static stiffness of a spatial translational parallel manipulator. CORRAL, J. (1), PINTO, Ch. (), ALTUZARRA, O. (), PETUA, V. (), DEL POZO, D. (1), LÓPEZ, J.M. (1) (1) Robotier
More informationSUMMARY. method to synthetic datasets is discussed in the present paper.
Geophysical modeling through simultaneous Joint Inversion of Seismic, Gravity and Magnetotelluric data Michele De Stefano (1), Daniele Colombo (1) WesternGeco EM - Geosystem, via Clericetti 42/A, 20133
More informationProgram Families in Scientific Computing
Slide 1 of 20 Mesh Program Families in Scientific Computing Spencer Smith, John McCutchan and Fang Cao Department of Computing and Software McMaster University The 7 th OOPSLA Workshop on Domain-Specific
More informationLiDAR data overview. Dr. Keiko Saito Global Facility for Disaster Reduction and Recovery (GFDRR)
LiDAR data overview Dr. Keiko Saito Global Facility for Disaster Reduction and Recovery (GFDRR) LiDAR (Light Detecting And Ranging) 3D height profile Laser emitted from sensor onboard aircraft to measure
More informationERTH 4121 Gravity and Magnetic Exploration Session 5. Gravity forward modelling
ERTH 4121 Gravity and Magnetic Exploration Session 5 Gravity forward modelling Lecture schedule (subject to change) Minimum 10 x 3 hour lecture sessions: 1:30pm Tuesdays Aug 2 : 1. Introduction to gravity
More information22. LECTURE 22. I can define critical points. I know the difference between local and absolute minimums/maximums.
. LECTURE Objectives I can define critical points. I know the difference between local and absolute minimums/maximums. In many physical problems, we re interested in finding the values (x, y) that maximize
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 informationMath 233. Lagrange Multipliers Basics
Math 33. Lagrange Multipliers Basics Optimization problems of the form to optimize a function f(x, y, z) over a constraint g(x, y, z) = k can often be conveniently solved using the method of Lagrange multipliers:
More informationGravity Methods (VII) wrap up
Environmental and Exploration Geophysics II Gravity Methods (VII) wrap up tom.h.wilson tom.wilson@mail.wvu.edu Department of Geology and Geography West Virginia University Morgantown, WV Items on the list
More informationDIGITAL TERRAIN MODELS
DIGITAL TERRAIN MODELS 1 Digital Terrain Models Dr. Mohsen Mostafa Hassan Badawy Remote Sensing Center GENERAL: A Digital Terrain Models (DTM) is defined as the digital representation of the spatial distribution
More informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationINTERNATIONAL JOURNAL OF GEOMATICS AND GEOSCIENCES Volume 5, No 2, 2014
INTERNATIONAL JOURNAL OF GEOMATICS AND GEOSCIENCES Volume 5, No 2, 2014 Copyright by the authors - Licensee IPA- Under Creative Commons license 3.0 Research article ISSN 0976 4380 Edge detection process
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 informationChapter 15 Notes, Stewart 7e
Contents 15.2 Iterated Integrals..................................... 2 15.3 Double Integrals over General Regions......................... 5 15.4 Double Integrals in Polar Coordinates..........................
More informationBasic landscape analysis
Chapter 4 Basic landscape analysis Tangible Landscape allows us to explore the spatial patterns of topographic parameters and their relation to basic surface geometry. To analyze the topography of a landscape
More informationReal Geodetic Map (Map without Projection) Abstract Keywords: 1. Introduction
Real ( without Projection) Ahmad Shaker 1 Abdullah Saad 1 Abdurrahman Arafa 2* 1.Surveying Dep., Shoubra Faculty of Engineering, Benha University, Egypt 2.Manager of Surveying Dep. in Horse Company. Egypt
More informationf xx + f yy = F (x, y)
Application of the 2D finite element method to Laplace (Poisson) equation; f xx + f yy = F (x, y) M. R. Hadizadeh Computer Club, Department of Physics and Astronomy, Ohio University 4 Nov. 2013 Domain
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
Calculus III-Final review Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Find the corresponding position vector. 1) Define the points P = (-,
More informationApplying Geodetic Coordinate Reference Systems in Building Information Modeling (BIM)
Applying Geodetic Coordinate Reference Systems in Building Information Modeling (BIM) Presented at the FIG Working Week 2017, May 29 - June 2, 2017 in Helsinki, Finland Robert Kaden und Christian Clemen
More informationMultiple Integrals. max x i 0
Multiple Integrals 1 Double Integrals Definite integrals appear when one solves Area problem. Find the area A of the region bounded above by the curve y = f(x), below by the x-axis, and on the sides by
More informationDownloaded 09/16/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Time-domain incomplete Gauss-Newton full-waveform inversion of Gulf of Mexico data Abdullah AlTheyab*, Xin Wang, Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 9// to
More informationGEOGRAPHIC INFORMATION SYSTEMS Lecture 24: Spatial Analyst Continued
GEOGRAPHIC INFORMATION SYSTEMS Lecture 24: Spatial Analyst Continued Spatial Analyst - Spatial Analyst is an ArcGIS extension designed to work with raster data - in lecture I went through a series of demonstrations
More information16. LECTURE 16. I understand how to find the rate of change in any direction. I understand in what direction the maximum rate of change happens.
6. LETURE 6 Objectives I understand how to find the rate of change in any direction. I understand in what direction the maximum rate of change happens. So far, we ve learned the definition of the gradient
More informationReality Check: Processing LiDAR Data. A story of data, more data and some more data
Reality Check: Processing LiDAR Data A story of data, more data and some more data Red River of the North Red River of the North Red River of the North Red River of the North Introduction and Background
More informationChoose the best answer. For 1 3, use the figure. For 4 5, use the figure. Which face is parallel to face DEF? BCEF
Choose the best answer. For, use the figure. Page For, use the figure. How many pairs of parallel faces does the figure have? Which face is parallel to face DEF? BCEF ACED ABFD ABC What is the value of
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 informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
More informationN.J.P.L.S. An Introduction to LiDAR Concepts and Applications
N.J.P.L.S. An Introduction to LiDAR Concepts and Applications Presentation Outline LIDAR Data Capture Advantages of Lidar Technology Basics Intensity and Multiple Returns Lidar Accuracy Airborne Laser
More informationI. An Intro to ArcMap Version 9.3 and 10. 1) Arc Map is basically a build your own Google map
I. An Intro to ArcMap Version 9.3 and 10 What is Arc Map? 1) Arc Map is basically a build your own Google map a. Display and manage geo-spatial data (maps, images, points that have a geographic location)
More informationHigh-Precision Positioning Unit 2.2 Student Exercise: Calculating Topographic Change
High-Precision Positioning Unit 2.2 Student Exercise: Calculating Topographic Change Ian Lauer and Ben Crosby (Idaho State University) Change is an inevitable part of our natural world and varies as a
More informationCalibration of NFR models with interpreted well-test k.h data. Michel Garcia
Calibration of NFR models with interpreted well-test k.h data Michel Garcia Calibration with interpreted well-test k.h data Intermediate step between Reservoir characterization Static model conditioned
More informationEOSC 454 Lab #3. 3D Magnetics. Date: Feb 3, Due: Feb 10, 2009.
EOSC 454 Lab #3 3D Magnetics Date: Feb 3, 2009. Due: Feb 10, 2009. 1 Introduction In this exercise you will perform both forward models and inversions of magnetic data. You will compare the response of
More informationf(x, y, z)dv = As remarked above, triple integrals can be evaluated as iterated integrals.
7.5 Triple Integrals These are just like double integrals, but with another integration to perform. Although this is conceptually a simple extension of the idea, in practice it can get very complicated.
More information