Robust 3D gravity gradient inversion by planting anomalous densities
|
|
- Oscar Long
- 5 years ago
- Views:
Transcription
1 Robust 3D gravity gradient inversion by planting anomalous densities Leonardo Uieda Valéria C. F. Barbosa Observatório Nacional September, 2011
2 Outline
3 Outline Forward Problem
4 Outline Forward Problem Inverse Problem
5 Outline Forward Problem Inverse Problem Planting Algorithm Inspired by René (1986)
6 Outline Forward Problem Inverse Problem Planting Algorithm Inspired by René (1986) Synthetic Data
7 Outline Forward Problem Inverse Problem Planting Algorithm Inspired by René (1986) Synthetic Data Real Data Quadrilátero Ferrífero, Brazil
8 Forward problem
9 Observed g αβ g αβ
10 Observed g αβ g αβ
11 Observed g αβ g αβ Anomalous density
12 Observed g αβ g Want to model this αβ Anomalous density
13 Interpretative model
14 Interpretative model Right rectangular prism Δ ρ= p j
15 Prisms with p j =0 not shown
16 Prisms with p j =0 not shown [] p1 p= p 2 pm
17 Predicted g αβ d Prisms with p j =0 not shown αβ [] p1 p= p 2 pm
18 Predicted g αβ d αβ M d = p j a αβ j=1 Prisms with p j =0 not shown [] p1 p= p 2 pm αβ j
19 Predicted g αβ d αβ M d = p j a αβ j=1 αβ j Contribution of jth prism Prisms with p j =0 not shown [] p1 p= p 2 pm
20 More components: xx d xy d xz d yy d yz d zz d
21 More components: xx d xy d xz d yy d yz d zz d d
22 More components: xx d xy d xz d yy d yz d zz d M d = p j a j j=1
23 More components: xx d xy d xz d yy d yz d zz d M d = p j a j = A p j=1
24 More components: xx d xy d xz d yy d yz d zz d M d = p j a j = A p j=1 Jacobian (sensitivity) matrix
25 More components: xx d xy d xz d yy d yz d zz d M d = p j a j = A p j=1 Column vector of A
26 Forward problem: p M d= p j a j j=1 d
27 Inverse problem: p? g
28 Inverse problem
29 Minimize difference between Residual vector g r= g d and d
30 Minimize difference between Residual vector r= g d Data misfit function: ( ϕ( p)= r 2 = N (gi d i ) i=1 g ) and d
31 Minimize difference between Residual vector r= g d Data misfit function: ( ϕ( p)= r 2 = ℓ2 norm of r N (gi d i ) i=1 g ) and d
32 Minimize difference between Residual vector r= g d Data misfit function: ( ϕ( p)= r 2 = N (gi d i ) i=1 ℓ2 norm of r Least squares fit g ) and d
33 Minimize difference between Residual vector g and d r= g d Data misfit function: ( ϕ( p)= r 2 = N (gi d i ) i=1 ℓ2 norm of r Least squares fit ) N ϕ( p)= r 1= gi d i i=1 ℓ1 norm of r
34 Minimize difference between Residual vector g and d r= g d Data misfit function: ( ϕ( p)= r 2 = N (gi d i ) i=1 ℓ2 norm of r Least squares fit ) N ϕ( p)= r 1= gi d i i=1 ℓ1 norm of r Robust fit
35 ill posed problem non existent non unique non stable
36 constraints ill posed problem non existent non unique non stable
37 constraints ill posed problem well posed problem non existent exist non unique unique non stable stable
38 Constraints: 1. Compact
39 Constraints: 1. Compact no holes inside
40 Constraints: 1. Compact no holes inside 2. Concentrated around seeds
41 Constraints: 1. Compact no holes inside 2. Concentrated around seeds User specified prisms Given density contrasts ρs Any # of density contrasts
42 Constraints: 1. Compact no holes inside 2. Concentrated around seeds User specified prisms Given density contrasts ρs Any # of density contrasts 3. Only p j =0 or p j =ρs
43 Constraints: 1. Compact no holes inside 2. Concentrated around seeds User specified prisms Given density contrasts ρs Any # of density contrasts 3. Only p j =0 or p j =ρs 4. p j =ρs of closest seed
44 Well posed problem: Minimize goal function Γ( p)=ϕ( p)+μ θ( p)
45 Well posed problem: Minimize goal function Γ( p)=ϕ( p)+μ θ( p) Data misfit function
46 Well posed problem: Minimize goal function Γ( p)=ϕ( p)+μ θ( p) Regularizing parameter (Tradeoff between fit and regularization)
47 Well posed problem: Minimize goal function Γ( p)=ϕ( p)+μ θ( p) Regularizing function M θ( p)= j=1 pj p j +ϵ l β j
48 Well posed problem: Minimize goal function Γ( p)=ϕ( p)+μ θ( p) Regularizing function M Similar to Silva Dias et al. (2009) θ( p)= j=1 pj p j +ϵ l β j
49 Well posed problem: Minimize goal function Γ( p)=ϕ( p)+μ θ( p) Regularizing function M Similar to Silva Dias et al. (2009) θ( p)= j=1 pj p j +ϵ l β j Distance between jth prism and seed
50 Well posed problem: Minimize goal function Γ( p)=ϕ( p)+μ θ( p) Regularizing function M Similar to Silva Dias et al. (2009) θ( p)= j=1 pj p j +ϵ l β j Distance between jth prism and seed Imposes: Compactness Concentration around seeds
51 Constraints: 1. Compact 2. Concentrated around seeds 3. Only p j =0 or p j =ρs 4. p j =ρs of closest seed Regularization
52 Constraints: 1. Compact 2. Concentrated around seeds 3. Only p j =0 or p j =ρs 4. p j =ρs of closest seed Regularization Algorithm Based on René (1986)
53 Planting Algorithm
54 Setup: g = observed data
55 Setup: g = observed data Define interpretative model Interpretative model
56 Setup: g = observed data Define interpretative model All parameters zero Interpretative model
57 Setup: g = observed data Define interpretative model All parameters zero N S seeds Interpretative model
58 Setup: g = observed data Define interpretative model All parameters zero N S seeds Include seeds Prisms with p j =0 not shown
59 Setup: g = observed data Define interpretative model All parameters zero N S seeds Include seeds d = predicted data Compute initial residuals (0) r = g d (0) Prisms with p j =0 not shown
60 Setup: g = observed data Define interpretative model All parameters zero N S seeds Include seeds d = predicted data Compute initial residuals (0) r = g d (0) Predicted by seeds Prisms with p j =0 not shown
61 Setup: g = observed data Define interpretative model All parameters zero N S seeds d = predicted data Include seeds Compute initial residuals NS ( ) r (0)= g s=1 ρs a j S Prisms with p j =0 not shown
62 Setup: g = observed data Define interpretative model All parameters zero N S seeds d = predicted data Include seeds Compute initial residuals NS ( ) r (0)= g s=1 ρs a j S Find neighbors of seeds Neighbors Prisms with p j =0 not shown
63 Growth: Try accretion to sth seed: Prisms with p j =0 not shown
64 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit 2. Smallest goal function Prisms with p j =0 not shown
65 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit 2. Smallest goal function j = chosen p j =ρs (New elements) j Prisms with p j =0 not shown
66 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit 2. Smallest goal function p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj j Prisms with p j =0 not shown
67 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit 2. Smallest goal function p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj j Contribution of j Prisms with p j =0 not shown
68 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit 2. Smallest goal function p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj None found = no accretion j Variable sizes Prisms with p j =0 not shown
69 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit NS 2. Smallest goal function p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj None found = no accretion Prisms with p j =0 not shown
70 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit NS 2. Smallest goal function p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj None found = no accretion j Prisms with p j =0 not shown
71 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit NS 2. Smallest goal function p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj None found = no accretion At least one seed grow? j Prisms with p j =0 not shown
72 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit 2. Smallest goal function NS p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj None found = no accretion At least one seed grow? j Yes Prisms with p j =0 not shown
73 Growth: Try accretion to sth seed: Choose neighbor: 1. Reduce data misfit 2. Smallest goal function NS p j =ρs (New elements) j = chosen Update residuals r ( new) =r (old ) pj aj None found = no accretion At least one seed grow? Yes j No Done! Prisms with p j =0 not shown
74 Advantages: Compact & non smooth Any number of sources Any number of different density contrasts No large equation system Search limited to neighbors
75 Remember equations: Initial residual (0) ( r = g Update residual vector NS ρs a j s=1 S ) r (new) =r (old ) pj aj
76 Remember equations: Initial residual (0) ( r = g Update residual vector NS ρs a j s=1 S ) r (new) =r No matrix multiplication (only vector +) (old ) pj aj
77 Remember equations: Initial residual (0) ( r = g Update residual vector NS ρs a j s=1 S ) r (new) =r (old ) pj aj No matrix multiplication (only vector +) Only need some columns of A
78 Remember equations: Initial residual (0) ( r = g Update residual vector NS ρs a j s=1 S ) r (new) =r (old ) pj aj No matrix multiplication (only vector +) Only need some columns of A Calculate only when needed
79 Remember equations: Initial residual (0) ( r = g Update residual vector NS ρs a j s=1 S ) r (new) =r (old ) pj aj No matrix multiplication (only vector +) Only need some columns of A Calculate only when needed & delete after update
80 Remember equations: Initial residual (0) ( r = g Update residual vector NS ρs a j s=1 S ) r (new) =r (old ) pj aj No matrix multiplication (only vector +) Only need some columns of A Calculate only when needed & delete after update Lazy evaluation
81 Advantages: Compact & non smooth Any number of sources Any number of different density contrasts No large equation system Search limited to neighbors
82 Advantages: Compact & non smooth Any number of sources Any number of different density contrasts No large equation system Search limited to neighbors No matrix multiplication (only vector +) Lazy evaluation of Jacobian
83 Advantages: Compact & non smooth Any number of sources Any number of different density contrasts No large equation system Search limited to neighbors No matrix multiplication (only vector +) Lazy evaluation of Jacobian Fast inversion + low memory usage
84 Synthetic Data
85 Data set: 3 components 51 x 51 points 2601 points/component 7803 measurements 5 Eötvös noise
86 Model:
87 Model: 11 prisms
88 Model: 11 prisms 4 outcropping
89 Model: 11 prisms 4 outcropping
90 Model: 11 prisms 4 outcropping
91 Strongly interfering effects
92 Strongly interfering effects What if only interested in these?
93 Common scenario
94 Common scenario May not have prior information Density contrast Approximate depth
95 Common scenario May not have prior information Density contrast Approximate depth No way to provide seeds
96 Common scenario May not have prior information Density contrast Approximate depth No way to provide seeds Difficult to isolate effect of targets
97 Robust procedure:
98 Robust procedure: Seeds only for targets
99 Robust procedure: Seeds only for targets
100 Robust procedure: Seeds only for targets ℓ1 norm to ignore non targeted
101 Robust procedure: Seeds only for targets ℓ1 norm to ignore non targeted
102 Inversion: 13 seeds 7,803 data
103 Inversion: 13 seeds 7,803 data
104 Inversion: 13 seeds 7,803 data
105 Inversion: 13 seeds 7,803 data
106 Inversion: 13 seeds 7,803 data
107 Inversion: 13 seeds 7,803 data 37,500 prisms
108 Inversion: 13 seeds 7,803 data 37,500 prisms Only prisms with zero density contrast not shown
109 Inversion: 13 seeds 7,803 data 37,500 prisms Only prisms with zero density contrast not shown
110 Inversion: 13 seeds 7,803 data 37,500 prisms Only prisms with zero density contrast not shown
111 Inversion: 13 seeds 7,803 data 37,500 prisms Only prisms with zero density contrast not shown
112 Inversion: 13 seeds 7,803 data 37,500 prisms Only prisms with zero density contrast not shown
113 Inversion: 13 seeds 7,803 data 37,500 prisms Recover shape of targets Only prisms with zero density contrast not shown
114 Inversion: 13 seeds 7,803 data Recover shape of targets Total time = 2.2 minutes (on laptop) 37,500 prisms Only prisms with zero density contrast not shown
115 Inversion: 13 seeds 7,803 data 37,500 prisms Predicted data in contours
116 Inversion: 13 seeds 7,803 data 37,500 prisms Predicted data in contours Effect of true targeted sources
117 Real Data
118 Data: 3 components FTG survey Quadrilátero Ferrífero, Brazil
119 Data: 3 components Targets: FTG survey Iron ore bodies Quadrilátero Ferrífero, Brazil BIFs of Cauê Formation
120 Data: 3 components Targets: FTG survey Iron ore bodies Quadrilátero Ferrífero, Brazil BIFs of Cauê Formation
121 Data: Seeds for iron ore: Density contrast 1.0 g/cm3 Depth 200 m
122 Predicted Observed Inversion: 46 seeds 13,746 data
123 Inversion: 46 seeds 13,746 data 164,892 prisms
124 Inversion: 46 seeds 13,746 data 164,892 prisms
125 Inversion: 46 seeds 13,746 data 164,892 prisms
126 Inversion: 46 seeds 13,746 data 164,892 prisms
127 Inversion: 46 seeds 13,746 data 164,892 prisms
128 Inversion: 46 seeds 13,746 data 164,892 prisms
129 Inversion: 46 seeds 13,746 data 164,892 prisms Only prisms with zero density contrast not shown
130 Inversion: 46 seeds 13,746 data 164,892 prisms Only prisms with zero density contrast not shown
131 Inversion: 46 seeds 13,746 data 164,892 prisms Only prisms with zero density contrast not shown
132 Inversion: 46 seeds 13,746 data 164,892 prisms
133 Inversion: 46 seeds 13,746 data 164,892 prisms
134 Inversion: 46 seeds 13,746 data 164,892 prisms Agree with previous interpretations (Martinez et al., 2010)
135 Inversion: 46 seeds 13,746 data 164,892 prisms Agree with previous interpretations (Martinez et al., 2010) Total time = 14 minutes (on laptop)
136 Conclusions
137 Conclusions New 3D gravity gradient inversion Multiple sources Interfering gravitational effects Non targeted sources No matrix multiplications No linear systems Lazy evaluation of Jacobian matrix
138 Conclusions Estimates geometry Given density contrasts Ideal for: Sharp contacts Well constrained physical properties Ore bodies Intrusive rocks Salt domes
139 Thank you
Computation of the gravity gradient tensor due to topographic masses using tesseroids
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
More informationIntroduction to Geophysical Inversion
Introduction to Geophysical Inversion Goals Understand the non-uniqueness in geophysical interpretations Understand the concepts of inversion. Basic workflow for solving inversion problems. Some important
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 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 informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
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 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 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 informationIntroduction to Geophysical Modelling and Inversion
Introduction to Geophysical Modelling and Inversion James Reid GEOPHYSICAL INVERSION FOR MINERAL EXPLORERS ASEG-WA, SEPTEMBER 2014 @ 2014 Mira Geoscience Ltd. Forward modelling vs. inversion Forward Modelling:
More informationApplication of wavelet theory to the analysis of gravity data. P. Hornby, F. Boschetti* and F. Horowitz, Division of Exploration and Mining, CSIRO,
Application of wavelet theory to the analysis of gravity data. P. Hornby, F. Boschetti* and F. Horowitz, Division of Exploration and Mining, CSIRO, Australia. Summary. The fundamental equations of potential
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 informationPerception Viewed as an Inverse Problem
Perception Viewed as an Inverse Problem Fechnerian Causal Chain of Events - an evaluation Fechner s study of outer psychophysics assumes that the percept is a result of a causal chain of events: Distal
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
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 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 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 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 informationLocally Weighted Learning for Control. Alexander Skoglund Machine Learning Course AASS, June 2005
Locally Weighted Learning for Control Alexander Skoglund Machine Learning Course AASS, June 2005 Outline Locally Weighted Learning, Christopher G. Atkeson et. al. in Artificial Intelligence Review, 11:11-73,1997
More informationMultiview Stereo COSC450. Lecture 8
Multiview Stereo COSC450 Lecture 8 Stereo Vision So Far Stereo and epipolar geometry Fundamental matrix captures geometry 8-point algorithm Essential matrix with calibrated cameras 5-point algorithm Intersect
More informationB15 Enhancement of Linear Features from Gravity Anomalies by Using Curvature Gradient Tensor Matrix
B5 Enhancement of Linear Features from Gravity Anomalies by Usin Curvature Gradient Tensor Matrix B. Oruç* (Kocaeli University) SUMMARY In this study, a new ede enhancement technique based on the eienvalues
More informationFiltering Images. Contents
Image Processing and Data Visualization with MATLAB Filtering Images Hansrudi Noser June 8-9, 010 UZH, Multimedia and Robotics Summer School Noise Smoothing Filters Sigmoid Filters Gradient Filters Contents
More informationNormalized cuts and image segmentation
Normalized cuts and image segmentation Department of EE University of Washington Yeping Su Xiaodan Song Normalized Cuts and Image Segmentation, IEEE Trans. PAMI, August 2000 5/20/2003 1 Outline 1. Image
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 information26257 Nonlinear Inverse Modeling of Magnetic Anomalies due to Thin Sheets and Cylinders Using Occam s Method
26257 Nonlinear Inverse Modeling of Anomalies due to Thin Sheets and Cylinders Using Occam s Method R. Ghanati* (University of Tehran, Insitute of Geophysics), H.A. Ghari (University of Tehran, Insitute
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
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 informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
More informationCalibrating an Overhead Video Camera
Calibrating an Overhead Video Camera Raul Rojas Freie Universität Berlin, Takustraße 9, 495 Berlin, Germany http://www.fu-fighters.de Abstract. In this section we discuss how to calibrate an overhead video
More informationCS4733 Class Notes, Computer Vision
CS4733 Class Notes, Computer Vision Sources for online computer vision tutorials and demos - http://www.dai.ed.ac.uk/hipr and Computer Vision resources online - http://www.dai.ed.ac.uk/cvonline Vision
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 information18.02 Final Exam. y = 0
No books, notes or calculators. 5 problems, 50 points. 8.0 Final Exam Useful formula: cos (θ) = ( + cos(θ)) Problem. (0 points) a) (5 pts.) Find the equation in the form Ax + By + z = D of the plane P
More informationMonocular Human Motion Capture with a Mixture of Regressors. Ankur Agarwal and Bill Triggs GRAVIR-INRIA-CNRS, Grenoble, France
Monocular Human Motion Capture with a Mixture of Regressors Ankur Agarwal and Bill Triggs GRAVIR-INRIA-CNRS, Grenoble, France IEEE Workshop on Vision for Human-Computer Interaction, 21 June 2005 Visual
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 informationPS wave AVO aspects on processing, inversion, and interpretation
PS wave AVO aspects on processing, inversion, and interpretation Yong Xu, Paradigm Geophysical Canada Summary Using PS wave AVO inversion, density contrasts can be obtained from converted wave data. The
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 10 130221 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Canny Edge Detector Hough Transform Feature-Based
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 information3D modeling of the Quest Projects Geophysical Datasets. Nigel Phillips
3D modeling of the Quest Projects Geophysical Datasets Nigel Phillips Advanced Geophysical Interpretation Centre Undercover Exploration workshop KEG-25 April 2012 Mineral Physical Properties: density sus.
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 informationDense Image-based Motion Estimation Algorithms & Optical Flow
Dense mage-based Motion Estimation Algorithms & Optical Flow Video A video is a sequence of frames captured at different times The video data is a function of v time (t) v space (x,y) ntroduction to motion
More informationFEMAP Tutorial 2. Figure 1: Bar with defined dimensions
FEMAP Tutorial 2 Consider a cantilevered beam with a vertical force applied to the right end. We will utilize the same geometry as the previous tutorial that considered an axial loading. Thus, this tutorial
More informationFitting: The Hough transform
Fitting: The Hough transform Voting schemes Let each feature vote for all the models that are compatible with it Hopefully the noise features will not vote consistently for any single model Missing data
More informationModel-based segmentation and recognition from range data
Model-based segmentation and recognition from range data Jan Boehm Institute for Photogrammetry Universität Stuttgart Germany Keywords: range image, segmentation, object recognition, CAD ABSTRACT This
More informationConcept of Curve Fitting Difference with Interpolation
Curve Fitting Content Concept of Curve Fitting Difference with Interpolation Estimation of Linear Parameters by Least Squares Curve Fitting by Polynomial Least Squares Estimation of Non-linear Parameters
More informationPerception. Autonomous Mobile Robots. Sensors Vision Uncertainties, Line extraction from laser scans. Autonomous Systems Lab. Zürich.
Autonomous Mobile Robots Localization "Position" Global Map Cognition Environment Model Local Map Path Perception Real World Environment Motion Control Perception Sensors Vision Uncertainties, Line extraction
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 informationTensor Based Approaches for LVA Field Inference
Tensor Based Approaches for LVA Field Inference Maksuda Lillah and Jeff Boisvert The importance of locally varying anisotropy (LVA) in model construction can be significant; however, it is often ignored
More informationMath Exam III Review
Math 213 - Exam III Review Peter A. Perry University of Kentucky April 10, 2019 Homework Exam III is tonight at 5 PM Exam III will cover 15.1 15.3, 15.6 15.9, 16.1 16.2, and identifying conservative vector
More informationMulticomponent f-x seismic random noise attenuation via vector autoregressive operators
Multicomponent f-x seismic random noise attenuation via vector autoregressive operators Mostafa Naghizadeh and Mauricio Sacchi ABSTRACT We propose an extension of the traditional frequency-space (f-x)
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 informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationLecture 15: Segmentation (Edge Based, Hough Transform)
Lecture 15: Segmentation (Edge Based, Hough Transform) c Bryan S. Morse, Brigham Young University, 1998 000 Last modified on February 3, 000 at :00 PM Contents 15.1 Introduction..............................................
More informationPatch-based Object Recognition. Basic Idea
Patch-based Object Recognition 1! Basic Idea Determine interest points in image Determine local image properties around interest points Use local image properties for object classification Example: Interest
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 informationNoise in FTG Data and its Comparison With Conventional Gravity Surveys
Noise in FTG Data and its Comparison With Conventional Gravity Surveys G. Barnes, J. umley, P. Houghton, R. Gleave ARKeX td. Cambridge, U.K. Summary Employing various assumptions and approximations, we
More informationIntegral equation method for anisotropic inversion of towed streamer EM data: theory and application for the TWOP survey
Integral equation method for anisotropic inversion of towed streamer EM data: theory and application for the TWOP survey Michael S. Zhdanov 1,2, Masashi Endo 1, Daeung Yoon 1,2, Johan Mattsson 3, and Jonathan
More informationAPPROXIMATING PDE s IN L 1
APPROXIMATING PDE s IN L 1 Veselin Dobrev Jean-Luc Guermond Bojan Popov Department of Mathematics Texas A&M University NONLINEAR APPROXIMATION TECHNIQUES USING L 1 Texas A&M May 16-18, 2008 Outline 1 Outline
More informationEstimation of Depth and Shape Factor of Buried Structure From Residual Gravity Anomaly Data
Australian Journal of Basic and Applied Sciences, 5(): 2-25, 2 ISSN 99-878 Estimation of Depth and Shape Factor of Buried Structure From Residual Gravity Anomaly Data MojtabaBabaee, Ahmad Alvandi, 2 Hosein
More informationFeature Matching and Robust Fitting
Feature Matching and Robust Fitting Computer Vision CS 143, Brown Read Szeliski 4.1 James Hays Acknowledgment: Many slides from Derek Hoiem and Grauman&Leibe 2008 AAAI Tutorial Project 2 questions? This
More informationMAT203 OVERVIEW OF CONTENTS AND SAMPLE PROBLEMS
MAT203 OVERVIEW OF CONTENTS AND SAMPLE PROBLEMS MAT203 covers essentially the same material as MAT201, but is more in depth and theoretical. Exam problems are often more sophisticated in scope and difficulty
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 informationGENERAL AUTOMATED FLAW DETECTION SCHEME FOR NDE X-RAY IMAGES
GENERAL AUTOMATED FLAW DETECTION SCHEME FOR NDE X-RAY IMAGES Karl W. Ulmer and John P. Basart Center for Nondestructive Evaluation Department of Electrical and Computer Engineering Iowa State University
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture XV (04.02.08) Contents: Function Minimization (see E. Lohrmann & V. Blobel) Optimization Problem Set of n independent variables Sometimes in addition some constraints
More informationP257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies
P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies. azemi* (University of Alberta) & H.R. Siahkoohi (University of Tehran) SUMMARY Petrophysical reservoir properties,
More informationFunction Approximation. Pieter Abbeel UC Berkeley EECS
Function Approximation Pieter Abbeel UC Berkeley EECS Outline Value iteration with function approximation Linear programming with function approximation Value Iteration Algorithm: Start with for all s.
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 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 informationMotion and Tracking. Andrea Torsello DAIS Università Ca Foscari via Torino 155, Mestre (VE)
Motion and Tracking Andrea Torsello DAIS Università Ca Foscari via Torino 155, 30172 Mestre (VE) Motion Segmentation Segment the video into multiple coherently moving objects Motion and Perceptual Organization
More informationA spectral boundary element method
Boundary Elements XXVII 165 A spectral boundary element method A. Calaon, R. Adey & J. Baynham Wessex Institute of Technology, Southampton, UK Abstract The Boundary Element Method (BEM) is not local and
More informationIntroduction to seismic body waves Tomography
Introduction to seismic body waves Tomography Presented by Han Zhang Prepared by group work of Justin Wilgus and Han Zhang Outlines 1. Velocity derived from wave equation 2. Two commonly used reference
More informationIMAGE PROCESSING >FILTERS AND EDGE DETECTION FOR COLOR IMAGES UTRECHT UNIVERSITY RONALD POPPE
IMAGE PROCESSING >FILTERS AND EDGE DETECTION FOR COLOR IMAGES UTRECHT UNIVERSITY RONALD POPPE OUTLINE Filters for color images Edge detection for color images Canny edge detection FILTERS FOR COLOR IMAGES
More informationLearning Two-View Stereo Matching
Learning Two-View Stereo Matching Jianxiong Xiao Jingni Chen Dit-Yan Yeung Long Quan Department of Computer Science and Engineering The Hong Kong University of Science and Technology The 10th European
More informationELEG Compressive Sensing and Sparse Signal Representations
ELEG 867 - Compressive Sensing and Sparse Signal Representations Gonzalo R. Arce Depart. of Electrical and Computer Engineering University of Delaware Fall 211 Compressive Sensing G. Arce Fall, 211 1 /
More informationThe points (2, 2, 1) and (0, 1, 2) are graphed below in 3-space:
Three-Dimensional Coordinate Systems The plane is a two-dimensional coordinate system in the sense that any point in the plane can be uniquely described using two coordinates (usually x and y, but we have
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationLarge Mesh Deformation Using the Volumetric Graph Laplacian
Large Mesh Deformation Using the Volumetric Graph Laplacian Kun Zhou1 Jin Huang2 John Snyder3 Xinguo Liu1 Hujun Bao2 Baining Guo1 Heung-Yeung Shum1 1 Microsoft Research Asia 2 Zhejiang University 3 Microsoft
More information11/10/2011 small set, B, to probe the image under study for each SE, define origo & pixels in SE
Mathematical Morphology Sonka 13.1-13.6 Ida-Maria Sintorn ida@cb.uu.se Today s lecture SE, morphological transformations inary MM Gray-level MM Applications Geodesic transformations Morphology-form and
More informationStabilization of the Euler deconvolution algorithm by means of a two steps regularization approach
Stabilization of the Euler deconvolution algorithm by means of a two steps regularization approach R. Pašteka ( 1,2 ), D. Kušnirák ( 1 ), H.-J. Götze ( 2 ) ( 1 )Department of Applied Geophysics, Comenius
More informationDiffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012
Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring 12 1. Introduction. This paper discusses research conducted in order to investigate
More informationStability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow
Stability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow F. Haider 1, B. Courbet 1, J.P. Croisille 2 1 Département de Simulation Numérique des Ecoulements
More informationRobotics. Lecture 5: Monte Carlo Localisation. See course website for up to date information.
Robotics Lecture 5: Monte Carlo Localisation See course website http://www.doc.ic.ac.uk/~ajd/robotics/ for up to date information. Andrew Davison Department of Computing Imperial College London Review:
More informationLecture 8: Jointly distributed random variables
Lecture : Jointly distributed random variables Random Vectors and Joint Probability Distributions Definition: Random Vector. An n-dimensional random vector, denoted as Z = (Z, Z,, Z n ), is a function
More information8(x 2) + 21(y 1) + 6(z 3) = 0 8x + 21y + 6z = 55.
MATH 24 -Review for Final Exam. Let f(x, y, z) x 2 yz + y 3 z x 2 + z, and a (2,, 3). Note: f (2xyz 2x, x 2 z + 3y 2 z, x 2 y + y 3 + ) f(a) (8, 2, 6) (a) Find all stationary points (if any) of f. et f.
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 informationOutcomes List for Math Multivariable Calculus (9 th edition of text) Spring
Outcomes List for Math 200-200935 Multivariable Calculus (9 th edition of text) Spring 2009-2010 The purpose of the Outcomes List is to give you a concrete summary of the material you should know, and
More informationBioimage Informatics
Bioimage Informatics Lecture 14, Spring 2012 Bioimage Data Analysis (IV) Image Segmentation (part 3) Lecture 14 March 07, 2012 1 Outline Review: intensity thresholding based image segmentation Morphological
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 information302 CHAPTER 3. FUNCTIONS OF SEVERAL VARIABLES. 4. Function of several variables, their domain. 6. Limit of a function of several variables
302 CHAPTER 3. FUNCTIONS OF SEVERAL VARIABLES 3.8 Chapter Review 3.8.1 Concepts to Know You should have an understanding of, and be able to explain the concepts listed below. 1. Boundary and interior points
More informationLecture 7: Most Common Edge Detectors
#1 Lecture 7: Most Common Edge Detectors Saad Bedros sbedros@umn.edu Edge Detection Goal: Identify sudden changes (discontinuities) in an image Intuitively, most semantic and shape information from the
More informationEdges and Binary Images
CS 699: Intro to Computer Vision Edges and Binary Images Prof. Adriana Kovashka University of Pittsburgh September 5, 205 Plan for today Edge detection Binary image analysis Homework Due on 9/22, :59pm
More informationLinear Analysis of an Arch Dam
Linear Analysis of an Arch Dam Outline 1 Description 2 Finite Element Model 2.1 Coincidence Tolerance 2.2 Units 2.3 Import CAD File 2.4 Properties 2.5 Boundary Conditions 2.6 Loads 2.7 Meshing 3 Linear
More informationPractical applications of uniqueness theorems in gravimetry: Part I Constructing sound interpretation methods
GEOPHYSICS, VOL. 67, NO. 3 (MAY-JUNE 2002); P. 788 794, 6 FIGS., 1 TABLE. 10.1190/1.1484522 Practical applications of uniqueness theorems in gravimetry: Part I Constructing sound interpretation methods
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 informationf for Directional Derivatives and Gradient The gradient vector is calculated using partial derivatives of the function f(x,y).
Directional Derivatives and Gradient The gradient vector is calculated using partial derivatives of the function f(x,y). For a function f(x,y), the gradient vector, denoted as f (pronounced grad f ) is
More informationEfficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers
Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou
More informationDS Machine Learning and Data Mining I. Alina Oprea Associate Professor, CCIS Northeastern University
DS 4400 Machine Learning and Data Mining I Alina Oprea Associate Professor, CCIS Northeastern University September 20 2018 Review Solution for multiple linear regression can be computed in closed form
More informationA Method of Automated Landmark Generation for Automated 3D PDM Construction
A Method of Automated Landmark Generation for Automated 3D PDM Construction A. D. Brett and C. J. Taylor Department of Medical Biophysics University of Manchester Manchester M13 9PT, Uk adb@sv1.smb.man.ac.uk
More informationEdge Detection (with a sidelight introduction to linear, associative operators). Images
Images (we will, eventually, come back to imaging geometry. But, now that we know how images come from the world, we will examine operations on images). Edge Detection (with a sidelight introduction to
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 informationCS 231A Computer Vision (Winter 2014) Problem Set 3
CS 231A Computer Vision (Winter 2014) Problem Set 3 Due: Feb. 18 th, 2015 (11:59pm) 1 Single Object Recognition Via SIFT (45 points) In his 2004 SIFT paper, David Lowe demonstrates impressive object recognition
More information3D Computer Vision. Dense 3D Reconstruction II. Prof. Didier Stricker. Christiano Gava
3D Computer Vision Dense 3D Reconstruction II Prof. Didier Stricker Christiano Gava Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de
More information