THE UNIVERSITY OF CALGARY
|
|
- Hollie Norton
- 6 years ago
- Views:
Transcription
1 Important Notice This copy may be used only for the purposes of research and private study, and any use of the copy for a purpose other than research or private study may require the authorization of the copyright owner of the work in question. Responsibility regarding questions of copyright that may arise in the use of this copy is assumed by the recipient.
2 Formulation
3 THE UNIVERSITY OF CALGARY FACULTY
4 Abstract One approach to local least-squares polynomial approximation of data in two independent variables consists of centring a square subarea on each point in a data set and computing the coefficients that define the best-fit polynomial over the subset defined
5 Approximating
6 Acknowledgements I am grateful to my supervisor, Dr. Brown, for his guidance through the course of this work. His urging to explore questions that were of foremost interest to me made this project
7 Table of Contents Approval Page... ii Abstract...
8 Introduction Regional and residual anomaly maps Long-wavelength anomaly maps...
9 Table 2.1. Terms included List of Tables
10 List
11 Figure 3.7. Transfer function of the ^-component gradient operator for a 5x5 window, and a second-order polynomial Figure 3.8. Transfer function of the ^-component gradient operator for a 7x7 window, and a first- (a), second- (b) and third-order (c) polynomial... 73,74,75 Figure 3.9. Transfer function of the ^-component gradient operator for a 9x9 (a) and an 11x11
12 approximating polynomial 2
13 Chapter 1 Introduction 1 There are two main objectives of this thesis. The first is to use least-squares polynomial approximation methods
14 3 the local least-squares polynomial approach are first derived and studied; implications for
15 Chapter
16 addressed 5
17 r=0 5=0
18 0» -My)I < e (2.2) is satisfied over the interval defined by x and y. In addition, if f(x,y) is 2-n periodic, continuous everywhere,
19 inhibited by this restriction. That is, polynomials of the form of equations 2.1 can approximate data with infinitely many wavelengths regardless of window size, the inherent difference being that polynomial approximation requires no assumptions regarding the periodicity of the signal. Thus Fourier series expansion does not lend itself
20 coefficients can be considered to be the regional. This regional can either be subtracted from 9
21 10 that removing residual features via local least-squares approximation is equivalent to linear digital filtering (Wood
22 been given by Davis (1986). Hence only a cursory discussion is presented here. It is important 11
23 12 as = O, (2.5 for
24 13 features from data in two independent variables has been given by Hayes (1970). However in Hayes's discussion, the computational procedure is not stated in terms of a convolution. Hence, after summarizing the relevant elements of Hayes's review, this procedure,
25 14 where t+w E X tpr i+w E PW t-i-w (2.9) t+w B f-i-w
26 15 t+w j+\v E EEE-^2^ r=i-w u=j-w r s "" ** w (2.13) -. /' Evaluating equation 2.13 at X=X 1 and y=jj gives an expression for the best-fit polynomial
27 16 This Z ns( x t*yj) ~ z (*i»)0 ~ J^ 53 ^tn z (*f v a)- ( 2-16 ) t=-~ «= -
28 17 even functions, and orthogonal polynomials of odd orders are odd functions. That is, from equation 2.9, p r (x)
29 the window size 18
30 Wavenumber-domain response of local least-squares filters 19
31 20 window size =* 0.3 =
32 k y (cycles/grid interval)
33 exhibit reliable filtering characteristics. 22
34 domain coordinate axes. 23 The spectrum in figure 2.1b shows that the amplitude spectrum is symmetric about
35 , O window size = 0.3 TJ 0-2 fr
36 25 figures 2.3 is the suite of all possible filters that may be generated for a 9x9 window. Increasing the polynomial order increases the cutoff wavenumber of the filter. As well, these spectra demonstrate that
37 o.o window size = 9x9 polynomial order = 0/ k x (cycles/grid interval) (a) Figure 2.3a-d. Amplitude spectra
38 ky (cycles/grid interval) to
39 ^l 1 8-.«e. ^ ^g- H u 5? ^* ttl p b \ o p N> p O4 to OO k y (cycles/grid interval) P b p p pppppp bo ki b> In > Giho^ P P k> ^- O a S O /""^* f.' * < * S P b
40 ky (cycles/grid Interval) K)
41 , window size = 7x7 polynomial order = 2/ TJ , k x (cycles/grid interval) Figure 2.4. Amplitude spectrum of a filter computed using a 7x7 window and a secondorder polynomial.
42 31
43 32 ee f order = wavenumber (cycles/grid interval) 0.5 (a) Figure 2.5. Bandwidth
44 I 33
45 window size. This is an important criterion for selecting window size and polynomial order 34
46 35 19-i order = 3 order wavelength (grid intervals) (a) Figure 2.6a. Cutoff wavelengths (a) along the frequency axes, and (b) along k x =k r A first-order best-fit line
47 order = 3 order = 1 O wavelength (grid intervals) (b)
48 37 r=0 (2.26) where
49 38 Study of such filters is complicated by the fact that they introduce phase distortion. Hence both amplitude and phase spectra must be studied when evaluating the performance
50 window size =
51 k y (cycles/grid interval)
52 window size = 41
53 42 o.o k x (cycles/grid interval) (contour units (b)
54 43 (2.28) B=-/ 5=0 l
55 44 difficulties encountered when the infinite operator is truncated to give a useable operator.
56 Chapter 45
57 46 ~^% (3-2a) and direction 0 = tan' 1 -^. (3.2b) ^ The present study focuses on computing and applying the magnitude of the horizontal gradient. Equations 3.1 and 3.2a indicate the theoretical differences between acceleration and gradient magnitudes.
58 A review 47
59 Gradient-magnitude data have been computed by this technique for several 48 published studies. For example, Cordell (1979) inferred that the maxima of such data located graben-bounding faults
60 operators 49
61 50 M M AB Al,.-_ r=-_ 4rCv+y) (3.6b) where S^
62 It will 51
63 52 filtering data were previously discussed. These advantages hold true for gradient computations.
64 The method of nonorthogonal polynomials 53 As before, the data are assumed to be on a rectangular grid, with dimensions max by p&y where Ax and Ay are the grid intervals in the x and v directions. The moving window is specified by {*,., y }.J] through {x i+w, y j+w )}, where i = 1,2,...,w, y = 1,2,...,«, w = 1,2,...,2&+1, k < (p-l)/2 and k < (m-l)/2. Again, let Z(x,y) represent the best-fit polynomial in the window; then, Z(XyV) n n rw/ \ V^ "T^ T S C\~l\
65 r=0 5=0 '- 1 (3-8b) Imposing the L 2 norm is one method of determining the best-fit polynomial for the data, 54 where de = O (3.9 da rj for all r = 0,1,2,... «, and 5 = 0,1,2,...«Solving i+w J+\v E = E E {ztw,)-^,)) 2. f=i-w u-j-w (3JO)
66 and 55 Ul \
67 *> %% ax i i (3.12a) 56
68 57 (3.14a) and
69 58 function over a symmetrical interval, and thus is equal to zero. The product under summation in the numerator of equation 3.14b can be seen to be odd by recalling the property
70 59 There are four properties of the impulse response of gradient-component operators that are not shared by the regional filters. These are here discussed.
71 60 function reduces to a one-dimensional operator. Using the notation of equations 3.12, non-zero coefficients occur for m equal to zero for ^-component gradient operators, and
72 band 61
73 ^I * I! C> t1 i;? idea t : C-C C i c 1 ^ S i»1npontenit 5" d Ler ' It 1 fil tor » I i < > t << h I C 5 t i ^; > I ' 0.1 ^ 5 \ S < *. «1 «< t «I«3 1 ^ S 9 «5 0.0 OO I OJ ki(cycles/grid interval) (a) 0.4 OJ I T I I 2» Ui IS ideal y-compooeac gradient filter s O4 0.4 I l I I
74 t I 63 11S'W 56'N C.I mgal/km O 50 I km I Figure 3.2. Horizontal-gradient magnitude of the Bouguer data shown in figure 4.1, computed
75 functional frequency range. The following three incarnations of this ideal case attempt 64
76 I I I 0.4 O m n m n 0 O n?n n " 1 In 0.3 TJ
77 frequency half of the passband, and thus should give reliable results for studies 66 focusing
78 k x (cycles/grid interval) 0.4
79 o.o k x (cycles/grid interval) Figure 3.5. Amplitude spectrum
80 band-pass filter designated by 0.04 and 0.08 cycles per grid interval, and a low-pass filter designated 69
81 k x (cycles/grid interval) Figure 3.6. Amplitude spectrum of y-component gradient operator from Kis (1983), calculated
82 T~ window size =
83 pass virtually 72
84 ,
85 , , i,, -, window size =* 7x7 polynomial order = k x (cycles/grid interval) (b)
86 p Oi p b k y (cycles/grid interval) i - a K) <P 3 P A " Ir p o
87 window size = 9x9 polynomial order = Tf 0.3 S3I
88 window size =* 11x11 polynomial order =
89 These data were obtained by empirically estimating the limits of the passband for a range of window sizes and polynomial orders. Because these filters are band-pass with respect to one frequency axis, and low-pass with respect to the other, a complete specification
90 33-i first-order polynomial CO high-pass wavenumbcr ^ low-pass wavenumber I I I I I I I I I I I I l I I I I T I I I I I
91 80 33-i third-order polynomial OT high-pass wavenumber "%k low-pass wavenumber k y (cycles/grid interval) (b)
92 81 passband widths I I third-order polynomial first-order polynomial k y (cycles/grid interval) (c)
93 82 component operators. First, these plots provide an estimate of the frequency content of data filtered by a weight function of a given polynomial order and window size. This serves as a guide for selecting parameters according to the objectives of the study, as well as indicating the frequencies represented on gradient data. In addition they indicate
94 Chapter 4 - Wavelength Filtering of Bouguer Gravity Data 83 Introduction
95 84 116'W t W5 /) VV Y(( u 56'N C.I. = 4 mgal km 50 Figure 4.1. Bouguer gravity
96 85 residuals has been computed. Filtering was performed in the frequency domain with a sixth-order Butterworth filter (Kanasewich, 1981, p. 265) using software developed by Eaton (1988).
97 X.
98 116* W Residual X, = 62.5 km Regional 116'W I 56'N C.l. = 2 mgal C.I.
99 116'W Residual 56'N Regional 116'W X, = 125 km I 56'N C.I. = 2 mgal C.I.
100
101 C.I.
102 Long-wavelength anomaly maps 91 Kane and Godson (1985) studied attenuation caused by filtering, and concluded that anomaly maps containing wavelengths predominantly shorter than
103 area is probably not in isostatic equilibrium. These authors surmised that this anomaly is 92
104 Relative gravity highs 93
105 94 On the map shown in figure 4.2h the most prominent feature is an approximately 170 km x 60 km north-south (in the northern portion) and southwestnortheast trending
106 basement rocks corresponding 95
107 adjacent pluton 96
108 transformation with horizontal-gradient component filters. Frequency filtering 97
109 Chapter 98
110 Regional/residual separation using polynomial approximation 99 In chapter
111 implemented 100
112 (a) (b) 56'N i data-set size=181xl81 order =11 11 rw 116"W I 56'N C.I.= 2 mgal ton BO C.I.- 2 mgal Figure 5.1. Data from figure
113 116'W 116* W (a) (b) km SO window size = 57x57 order = 2/3 I 56'N C.l.«= 2 C.I.- 2 rngaj Figure 5.2. Data from figure
114 data-set size
115 window size = 55x55 order = 0/ 'W C.I.=
116 Regional maps computed using polynomial approximation 105 reproduces No one regional field computed using polynomial approximation wholly
117 approximation, it is quite likely an artefact. The fact that it does not appear when the data 106
118 Finally, 107
119 content 108
120 109 various combinations of parameters. However, irrespective of the behaviour of the filter, transforming an anomaly to its gradient-magnitude counterpart inherently alters the frequency content. This occurs because the magnitude of the gradient of an input waveform has wavelengths half as long, provided these wavelengths fall within the passband of the gradient-component operators. This
121 110 N- Figure 5.5. A monochromatic waveform, and its gradient magnitude. A single cycle of the input waveform spans two cycles of the gradient magnitude.
122 instance, gradient-component operators based on a 31x31 window and a first-order Ill polynomial enhance wavelengths approximately between
123 'W 56'N C.I mgal/km, km 50 Figure 5.6. Horizontal-gradient magnitude data computed using a 5x5 window, and a third-order polynomial.
124 'W t I 56*N C.I mgal/km O 50 I km Figure 5.7. Horizontal-gradient magnitude data computed using a 17x17 window, and a second-order polynomial.
125 'W t I 56'N C.I mgal/km. km 50 Figure 5.8. Horizontal-gradient magnitude data computed using
126 'W I 56* N C.I.B 0.05 mgal/km O ^ 50 km I Figure 5.9. Horizontal-gradient magnitude data computed using a 51x51 window, and a first-order polynomial.
127 in 116
128 overestimating 117
129 5.7, 118
130 portion 119
131 120 transmission of considerably higher frequencies along the directions of the coordinate axes. This degrades features that are oblique to the axes. Polynomials defined by equation 2.1b
132 121 the gradient-component operators. This makes the horizontal-gradient magnitude technique unsuitable for studying long-wavelength crustal features that are masked by higher-amplitude shorter-wavelength features.
133 Chapter 122
134 proportional 123
135 in the frequency domain. The former are restricted because the passbands cannot be varied, whereas 124
136 References 125 Abdelrahman, E.M., Bayoumi, A. L, Abdelhady Y.E., Gobashy, M.M., and El-Araby, H.M., 1989, Gravity interpretation using correlation factors between successive least-squares residual anomalies: Geophysics,
137 Chandler, V.W., Koski, IS., Hinze, WJ., and Braile, L.W., 1981, Analysis of multisource gravity 126
138 Grauch, V.J.S., 127
139 128 Kis, K., 1983, Determination
140 Stewart, R.R., 1985, Median filtering: Review 129
141 Zurflech, E.G., 1967, Application of two-dimensional linear wavelength filtering: Geophysics, 130
Geologic applications of magnetic data and using enhancements for contact mapping
Geologic applications of magnetic data and using enhancements for contact mapping M. Pilkington ( 1 ), P. Keating ( 1 ) ( 1 ) Geological Survey of Canada, Ottawa, Canada Summary Magnetic field data are
More informationLecture 27 Frequency Response 2
Lecture 27 Frequency Response 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/12 1 Application of Ideal Filters Suppose we can generate a square wave with a fundamental period
More informationPotential Field investigation of Williston Basin Basement
Potential Field investigation of Williston Basin Basement Jiakang Li* and Igor Morozov, University of Saskatchewan, Canada Summary Major faults and domains are identified in the Precambrian basement of
More informationOUTLINE. Many of us secretly dream of six months without gravity COURSE DESCRIPTION
GEOL 481.3 OUTLINE POTENTIAL FIELD METHODS GEOL 481.3 email: jim.merriam@usask.ca POTENTIAL FIELD METHODS Many of us secretly dream of six months without gravity Allan Fotheringham COURSE DESCRIPTION This
More informationDepth and shape solutions from residual gravity anomalies due to simple geometric structures using a statistical approach
Contributions to Geophysics and Geodesy Vol. 47/2, 2017 (113 132) Depth and shape solutions from residual gravity anomalies due to simple geometric structures using a statistical approach El-Sayed ABDELRAHMAN,
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More information1. FIR Filter Design
ELEN E4810: Digital Signal Processing Topic 9: Filter Design: FIR 1. Windowed Impulse Response 2. Window Shapes 3. Design by Iterative Optimization 1 1. FIR Filter Design! FIR filters! no poles (just zeros)!
More informationIterative 3D gravity inversion with integration of seismologic data
BOLLETTINO DI GEOFISICA TEORICA ED APPLICATA VOL. 40, N. 3-4, pp. 469-475; SEP.-DEC. 1999 Iterative 3D gravity inversion with integration of seismologic data C. BRAITENBERG and M. ZADRO Department of Earth
More informationA SPECTRAL ATTENUATION MODEL FOR JAPAN USING DIGITAL STRONG MOTION RECORDS OF JMA87 TYPE
A SPECTRAL ATTENUATION MODEL FOR JAPAN USING DIGITAL STRONG MOTION RECORDS OF JMA87 TYPE Shuji KOBAYASHI 1, Tetsuichi TAKAHASHI 2, Shin'ichi MATSUZAKI 3, Masafumi MORI 4, Yoshimitsu FUKUSHIMA 5, John X
More informationChapter 5: Application of Filters to Potential Field Gradient Tensor Data
Chapter 5: Filters 98 Chapter 5: Application of Filters to Potential Field Gradient Tensor Data 5.1 Introduction One of the objectives of this research is to investigate the use of spatial filters on potential
More informationGravity Support for Hydrocarbon Exploration at the Prospect Level
Journal of Emerging Trends in Engineering and Applied Sciences (JETEAS) 2 (1): 1-6 Scholarlink Research Institute Journals, 2011 (ISSN: 2141-7016) jeteas.scholarlinkresearch.org Journal of Emerging Trends
More informationChapter 7: Filter Design 7.1 Practical Filter Terminology
hapter 7: Filter Design 7. Practical Filter Terminology Analog and digital filters and their designs constitute one of the major emphasis areas in signal processing and communication systems. This is due
More informationAvailable online Journal of Scientific and Engineering Research, 2016, 3(2):1-7. Research Article
Available online www.jsaer.com, 2016, 3(2):1-7 Research Article ISSN: 2394-2630 CODEN(USA): JSERBR Assessment of the Reliability of Magnetic Method to Delineate Geologic Features in a Basement Complex:
More informationC5 Magnetic exploration methods data analysis techniques
C5 Magnetic exploration methods data analysis techniques C5.1 Data processing and corrections After magnetic field data have been collected a number of corrections are applied to simplify the interpretation.
More informationWork with Gravity Data in GM-SYS Profile
Work with Gravity Data in GM-SYS Profile 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 informationChapter 8. Conclusions and Further Work Overview of Findings
8.1. Overview of Findings The aim of this thesis was to analyse, theoretically and experimentally, the signal and noise aspects of digital magnetic recording on longitudinal thin-film disk media. This
More informationFourier Transform for Continuous Functions
Fourier Transform for Continuous Functions Central goal: representing a signal by a set of orthogonal bases that are corresponding to frequencies or spectrum. Fourier series allows to find the spectrum
More informationA. Relationship of DSP to other Fields.
1 I. Introduction 8/27/2015 A. Relationship of DSP to other Fields. Common topics to all these fields: transfer function and impulse response, Fourierrelated transforms, convolution theorem. 2 y(t) = h(
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More informationLecture 16: Bessel s Inequality, Parseval s Theorem, Energy convergence
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. ot to be copied, used, or revised without explicit written permission from the copyright owner. ecture 6: Bessel s Inequality,
More informationEdge Detection in Gravity Field of the Gheshm Sedimentary Basin
Int J Min & Geo-Eng (IJMGE), Vol. 47, No. 1, Jun. 013, pp. 41-50 Edge Detection in Gravity Field of the Gheshm Sedimentary Basin Seyed Ali Akbar Hosseini 1*, Faramarz Doulati Ardejani, Seyed Hashem Tabatabaie
More informationEmpirical Mean and Variance!
Global Image Properties! Global image properties refer to an image as a whole rather than components. Computation of global image properties is often required for image enhancement, preceding image analysis.!
More informationReduction to the Pole of Magnetic Anomalies Using Analytic Signal. A.H. Ansari and K. Alamdar
World Applied Sciences Journal 7 (4): 405-409, 2009 SSN 1818-4952 DOS Publications, 2009 Reduction to the Pole of Magnetic Anomalies Using Analytic Signal A.H. Ansari K. Alamdar Department of Mining Metallurgical
More informationContinuous Fourier transform of a Gaussian Function
Continuous Fourier transform of a Gaussian Function Gaussian function: e t2 /(2σ 2 ) The CFT of a Gaussian function is also a Gaussian function (i.e., time domain is Gaussian, then the frequency domain
More informationRotational Raman Spectroscopy
Rotational Raman Spectroscopy If EM radiation falls upon an atom or molecule, it may be absorbed if the energy of the radiation corresponds to the separation of two energy levels of the atoms or molecules.
More informationSound & Vibration Magazine March, Fundamentals of the Discrete Fourier Transform
Fundamentals of the Discrete Fourier Transform Mark H. Richardson Hewlett Packard Corporation Santa Clara, California The Fourier transform is a mathematical procedure that was discovered by a French mathematician
More informationELEG 5173L Digital Signal Processing Ch. 5 Digital Filters
Department of Electrical Engineering University of Aransas ELEG 573L Digital Signal Processing Ch. 5 Digital Filters Dr. Jingxian Wu wuj@uar.edu OUTLINE 2 FIR and IIR Filters Filter Structures Analog Filters
More informationThe Approximation Problem
EE 508 Lecture 3 The Approximation Problem Classical Approximating Functions - Thompson and Bessel Approximations Review from Last Time Elliptic Filters Can be thought of as an extension of the CC approach
More information176 5 t h Fl oo r. 337 P o ly me r Ma te ri al s
A g la di ou s F. L. 462 E l ec tr on ic D ev el op me nt A i ng er A.W.S. 371 C. A. M. A l ex an de r 236 A d mi ni st ra ti on R. H. (M rs ) A n dr ew s P. V. 326 O p ti ca l Tr an sm is si on A p ps
More informationLABORATORY 3 FINITE IMPULSE RESPONSE FILTERS
LABORATORY 3 FINITE IMPULSE RESPONSE FILTERS 3.. Introduction A digital filter is a discrete system, used with the purpose of changing the amplitude and/or phase spectrum of a signal. The systems (filters)
More informationA Note on the Use of the Second Vertical Derivative (SVD) of Gravity Data with Reference to Indonesian Cases
J. Eng. Technol. Sci., Vol. 50, No. 1, 018, 17-139 17 A Note on the Use of the Second Vertical Derivative (SVD) of Gravity Data with Reference to Indonesian Cases Prihadi Sumintadireja 1, Darharta Dahrin
More informationGeophysical Investigations of the Precambrian Basement of the Williston Basin in south-eastern Saskatchewan and south-western Manitoba
Final Project Report: Geophysical Investigations of the Precambrian Basement of the Williston Basin in south-eastern Saskatchewan and south-western Manitoba Jiakang Li and Igor Morozov University of Saskatchewan
More information9.4 Enhancing the SNR of Digitized Signals
9.4 Enhancing the SNR of Digitized Signals stepping and averaging compared to ensemble averaging creating and using Fourier transform digital filters removal of Johnson noise and signal distortion using
More informationA Simple Iterative Approach to Shape and Depth Estimation from SP Residual Anomalies
JKAU: Earth Sci., vol. 14, A pp.1-19 Simple Iterative... (2002-2003 A.D. / 1423-1424 A.H.) 1 A Simple Iterative Approach to Shape and Depth Estimation from SP Residual Anomalies E.M. ABDELRAHMAN, T.M.
More informationDigital Baseband Systems. Reference: Digital Communications John G. Proakis
Digital Baseband Systems Reference: Digital Communications John G. Proais Baseband Pulse Transmission Baseband digital signals - signals whose spectrum extend down to or near zero frequency. Model of the
More informationELEG 3124 SYSTEMS AND SIGNALS Ch. 5 Fourier Transform
Department of Electrical Engineering University of Arkansas ELEG 3124 SYSTEMS AND SIGNALS Ch. 5 Fourier Transform Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Introduction Fourier Transform Properties of Fourier
More informationCONSTRAINED PERCOLATION ON Z 2
CONSTRAINED PERCOLATION ON Z 2 ZHONGYANG LI Abstract. We study a constrained percolation process on Z 2, and prove the almost sure nonexistence of infinite clusters and contours for a large class of probability
More informationAeromagnetic Data Interpretation for Geostructural Analysis of Ibadan, Southwestern Nigeria
Aeromagnetic Data Interpretation for Geostructural Analysis of Ibadan, Southwestern Nigeria 1 Sunmonu, L.A., 2* Olasunkanmi, N. K. and 3 Alagbe, O. A. 1 Pure and Applied Physics Department, Ladoke Akintola
More informationEE123 Digital Signal Processing. M. Lustig, EECS UC Berkeley
EE123 Digital Signal Processing Today Last time: DTFT - Ch 2 Today: Continue DTFT Z-Transform Ch. 3 Properties of the DTFT cont. Time-Freq Shifting/modulation: M. Lustig, EE123 UCB M. Lustig, EE123 UCB
More informationAcoustic holography. LMS Test.Lab. Rev 12A
Acoustic holography LMS Test.Lab Rev 12A Copyright LMS International 2012 Table of Contents Chapter 1 Introduction... 5 Chapter 2... 7 Section 2.1 Temporal and spatial frequency... 7 Section 2.2 Time
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing 6. DSP - Recursive (IIR) Digital Filters Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher: John Wiley
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationFourier Series Representation of
Fourier Series Representation of Periodic Signals Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline The response of LIT system
More informationConversion from Linear to Circular Polarization and Stokes Parameters in FPGA. Koyel Das, Alan Roy, Gino Tuccari, Reinhard Keller
Conversion from Linear to Circular Polarization and Stokes Parameters in FPGA Koyel Das, Alan Roy, Gino Tuccari, Reinhard Keller Purpose 1. Conventionally, for the formation of circular polarization, analogue
More informationr 2. (5) This agrees of course with the traditional multipole way of writing it: for an axial point dipole at the coordinate origin
1. A Digression on Dipoles From the perspective of potential theory, the dipole represents the point source obtained by allowing two point charges of opposite sign to approach each other, in such a way
More informationREDUNDANCY ANALYSIS AN ALTERNATIVE FOR CANONICAL CORRELATION ANALYSIS ARNOLD L. VAN DEN WOLLENBERG UNIVERSITY OF NIJMEGEN
PSYCHOMETRIKA-VOL. 42, NO, 2 JUNE, 1977 REDUNDANCY ANALYSIS AN ALTERNATIVE FOR CANONICAL CORRELATION ANALYSIS ARNOLD L. VAN DEN WOLLENBERG UNIVERSITY OF NIJMEGEN A component method is presented maximizing
More information2.2 Gravity surveys. Gravity survey
2.2 Gravity surveys Gravity survey The effect of latitude The effect of elevation The Bouguer effect Topographic effect The effect of tides Summary of corrections Gravity in boreholes Gravity survey In
More informationPower Spectral Density of Digital Modulation Schemes
Digital Communication, Continuation Course Power Spectral Density of Digital Modulation Schemes Mikael Olofsson Emil Björnson Department of Electrical Engineering ISY) Linköping University, SE-581 83 Linköping,
More information3D Geometry of the Xade Complex inferred from Gravity and Magnetic Data
Geophysical Case Histories 3D Geometry of the Xade Complex inferred from Gravity and Magnetic Data 1. British Geological Survey, Edinburgh, United Kingdom Paper 92 Pouliquen, G. [1], Key, R. [1] ABSTRACT
More informationLINEAR-PHASE FIR FILTERS DESIGN
LINEAR-PHASE FIR FILTERS DESIGN Prof. Siripong Potisuk inimum-phase Filters A digital filter is a minimum-phase filter if and only if all of its zeros lie inside or on the unit circle; otherwise, it is
More information(2) D. Mastellone, M. Fedi, V. Paoletti DISTAR, Department of Earth Sciences, Environment and Resources, University of Naples Federico II, Italy
Accurate evaluation of edges and dip of faults and contacts through the Volume Upward Continuation (VUC) of gravity data D. Mastellone, M. Fedi, V. Paoletti DISTAR, Department of Earth Sciences, Environment
More informationSome notes on processing: causal vs. acausal low-cut filters version 1.0. David M. Boore. Introduction
File: c:\filter\notes on processing.tex Some notes on processing: causal vs. acausal low-cut filters version. David M. Boore Introduction These are some informal notes showing results of some procedures
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationThe failure of the sounding assumption in electroseismic investigations
The failure of the sounding assumption in electroseismic investigations F.D. Fourie, J.F. Botha Institute for Groundwater Studies, University of the Free State, PO Box 339, Bloemfontein, 9300, South Africa
More informationDiscrete-time Signals and Systems in
Discrete-time Signals and Systems in the Frequency Domain Chapter 3, Sections 3.1-39 3.9 Chapter 4, Sections 4.8-4.9 Dr. Iyad Jafar Outline Introduction The Continuous-Time FourierTransform (CTFT) The
More informationTime-Gate Filtering in Geophysical Response Estimation
Terrigena Ltd. Time-Gate Filtering in Geophysical Response Estimation Author: John E. E. Kingman Terrigena Technical White Paper Number 2005FEB01.0104(b) 2005 - FEB - 01 Copyright Terrigena, Ltd., 2005
More informationTheory and Problems of Signals and Systems
SCHAUM'S OUTLINES OF Theory and Problems of Signals and Systems HWEI P. HSU is Professor of Electrical Engineering at Fairleigh Dickinson University. He received his B.S. from National Taiwan University
More informationFourier Series Example
Fourier Series Example Let us compute the Fourier series for the function on the interval [ π,π]. f(x) = x f is an odd function, so the a n are zero, and thus the Fourier series will be of the form f(x)
More informationINTRODUCTION TO APPLIED GEOPHYSICS
INTRODUCTION TO APPLIED GEOPHYSICS EXPLORING THE SHALL0W SUBSURFACE H. Robert Burger Anne F. Sheehan Craig H.Jones VERSITY OF COLORADO VERSITY OF COLORADO W. W. NORTON & COMPANY NEW YORK LONDON Contents
More informationLecture 3 - Design of Digital Filters
Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can
More informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture : Design of Digital IIR Filters (Part I) Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner (Univ.
More informationDesign of IIR filters
Design of IIR filters Standard methods of design of digital infinite impulse response (IIR) filters usually consist of three steps, namely: 1 design of a continuous-time (CT) prototype low-pass filter;
More informationADVANCED PROBLEMS AND SOLUTIONS Edited by RAYMOND E.WHITNEY Lock Haven State College, Lock Haven, Pennsylvania
ADVANCED PROBLEMS AND SOLUTIONS Edited by RAYMOND E.WHITNEY Lock Haven State College, Lock Haven, Pennsylvania Send all communications concerning Advanced Problems and Solutions to Raymond E. Whitney,
More informationAeromagnetic map of the Death Valley ground-water model area, Nevada and California
Aeromagnetic map of the Death Valley ground-water model area, Nevada and California By D.A. Ponce and R.J. Blakely Prepared in cooperation with the Nevada Operations Office National Nuclear Security Administration
More informationDepth determination of 2-D SP anomaly source using Energy Spectrum method and its advantages
P-187 Depth determination of 2-D SP anomaly source using Energy Spectrum method and its advantages Summary Baisakhi Das*, NGRI, Prof. B.N.P. Agarwal, Indian School of Mines Self potential (SP) method is
More informationFirst scalar magnetic anomaly map from CHAMP satellite data indicates weak lithospheric field
GEOPHYSICAL RESEARCH LETTERS, VOL. 29, NO. 14, 10.1029/2001GL013685, 2002 First scalar magnetic anomaly map from CHAMP satellite data indicates weak lithospheric field S. Maus, M. Rother, R. Holme, H.
More informationMITOCW watch?v=jtj3v Rx7E
MITOCW watch?v=jtj3v Rx7E The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To
More informationSignals and Systems Laboratory with MATLAB
Signals and Systems Laboratory with MATLAB Alex Palamides Anastasia Veloni @ CRC Press Taylor &. Francis Group Boca Raton London NewYork CRC Press is an imprint of the Taylor & Francis Group, an informa
More informationGravity data reduction
Gravity data reduction REDUCTION: raw data à gravity anomaly data Temporal corrections tides and instrument drift Spatial corrections latitude and elevation GRS67 = gravity variation with latitude at sea
More informationWALL D*TA PRINT-OUT. 3 nmth ZONE
T A B L E A 4. 3 G L A S S D A T A P R I N T - O U T H T C L».>qth» H e ig h t n u «b»r C L A S S D A T A P R I N T O U T it************************************ 1*q o v»rh # n g recm oi*ion*l orient n
More informationComputer Vision & Digital Image Processing. Periodicity of the Fourier transform
Computer Vision & Digital Image Processing Fourier Transform Properties, the Laplacian, Convolution and Correlation Dr. D. J. Jackson Lecture 9- Periodicity of the Fourier transform The discrete Fourier
More information5. LIGHT MICROSCOPY Abbe s theory of imaging
5. LIGHT MICROSCOPY. We use Fourier optics to describe coherent image formation, imaging obtained by illuminating the specimen with spatially coherent light. We define resolution, contrast, and phase-sensitive
More informationContinuous Wavelet Transform Reconstruction Factors for Selected Wavelets
Continuous Wavelet Transform Reconstruction Factors for Selected Wavelets General Background This report expands on certain aspects of the analytical strategy for the Continuous Wavelet Transform (CWT)
More informationCHAPTER 9 ELECTROMAGNETIC WAVES
CHAPTER 9 ELECTROMAGNETIC WAVES Outlines 1. Waves in one dimension 2. Electromagnetic Waves in Vacuum 3. Electromagnetic waves in Matter 4. Absorption and Dispersion 5. Guided Waves 2 Skip 9.1.1 and 9.1.2
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
IT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete all 2008 or information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. assachusetts
More informationLong-period Ground Motion Characteristics of the Osaka Sedimentary Basin during the 2011 Great Tohoku Earthquake
Long-period Ground Motion Characteristics of the Osaka Sedimentary Basin during the 2011 Great Tohoku Earthquake K. Sato, K. Asano & T. Iwata Disaster Prevention Research Institute, Kyoto University, Japan
More informationDistortion Effects of Faults on Gravity Worm Strings Robin O Leary
Distortion Effects of Faults on Gravity Worm Strings Robin O Leary Problem Can upward continued gravity anomaly worm strings be used as a tool to determine the dip direction of an offsetting structure
More informationOptimum Ordering and Pole-Zero Pairing of the Cascade Form IIR. Digital Filter
Optimum Ordering and Pole-Zero Pairing of the Cascade Form IIR Digital Filter There are many possible cascade realiations of a higher order IIR transfer function obtained by different pole-ero pairings
More informationModule 1: Signals & System
Module 1: Signals & System Lecture 6: Basic Signals in Detail Basic Signals in detail We now introduce formally some of the basic signals namely 1) The Unit Impulse function. 2) The Unit Step function
More informationDigital Signal Processing:
Digital Signal Processing: Mathematical and algorithmic manipulation of discretized and quantized or naturally digital signals in order to extract the most relevant and pertinent information that is carried
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationSPEECH ANALYSIS AND SYNTHESIS
16 Chapter 2 SPEECH ANALYSIS AND SYNTHESIS 2.1 INTRODUCTION: Speech signal analysis is used to characterize the spectral information of an input speech signal. Speech signal analysis [52-53] techniques
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 34 Outline 1 Infinite Order Finite Differences
More informationCHAPTER 29 COMPUTER MODELLING OF DIFFRACTION OF WIND WAVES ABSTRACT
CHAPTER 29 COMPUTER MODELLING OF DIFFRACTION OF WIND WAVES by Shou-shan Fan- 1 and L E Borgman 2 ABSTRACT A digital computer model for diffraction of wind waves behind a breakwater is developed The model
More informationFast wavefield extrapolation by phase-shift in the nonuniform Gabor domain
Gabor wavefield extrapolation Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureux ABSTRACT Wavefield extrapolation for
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationLab 8: Gravity and Isostasy (35 points)
it's not the most important thing in your life right now. But what is important is gravity. Arnold Schwarzenegger as Colonel John Matrix, Commando (Check out this classic of American cinema!) Lab 8: Gravity
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationPhysics Assessment Unit A2 3
New Specification Centre Number 71 Candidate Number ADVANCED General Certificate of Education 2010 Physics Assessment Unit A2 3 Practical Techniques (Internal Assessment) Session 1 [AY231] AY231 WEDNESDAY
More informationFast wavefield extrapolation by phase-shift in the nonuniform Gabor domain
Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain Jeff P. Grossman* and Gary F. Margrave Geology & Geophysics, University of Calgary 2431 22A Street NW, Calgary, AB, T2M 3X8 grossman@geo.ucalgary.ca
More informationEfficient Algorithms for Pulse Parameter Estimation, Pulse Peak Localization And Pileup Reduction in Gamma Ray Spectroscopy M.W.Raad 1, L.
Efficient Algorithms for Pulse Parameter Estimation, Pulse Peak Localization And Pileup Reduction in Gamma Ray Spectroscopy M.W.Raad 1, L. Cheded 2 1 Computer Engineering Department, 2 Systems Engineering
More informationUpdating the Chiou and YoungsNGAModel: Regionalization of Anelastic Attenuation
Updating the Chiou and YoungsNGAModel: Regionalization of Anelastic Attenuation B. Chiou California Department of Transportation R.R. Youngs AMEC Environment & Infrastructure SUMMARY: (10 pt) Ground motion
More information5.3 Rotational Raman Spectroscopy General Introduction
5.3 Rotational Raman Spectroscopy 5.3.1 General Introduction When EM radiation falls on atoms or molecules, it may be absorbed or scattered. If λis unchanged, the process is referred as Rayleigh scattering.
More informationmm. COPY HO. ARDC HARD COPY MICROFICHE COPY.._^ OF _X- TECHNICAL LIBRARY DOC. NO. i> COPY NO IA.4 $. / ^ US A F. AIR FORCE BALLISTIC MISSILE DIVISION
DOC. NO. i> COPY NO < COPY.._^ OF _X- HARD COPY MICROFICHE 3; $. / ^ IA.4 TECHNICAL LIBRARY AIR FORCE BALLISTIC MISSILE DIVISION HO. ARDC US A F. if i \ mm. COPY / - / >' ; ""?>* woo romm on i J n i t
More informationUNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter?
UNIT - III PART A. Mention the important features of the IIR filters? i) The physically realizable IIR filters does not have linear phase. ii) The IIR filter specification includes the desired characteristics
More informationMathematics for Engineers II. lectures. Power series, Fourier series
Power series, Fourier series Power series, Taylor series It is a well-known fact, that 1 + x + x 2 + x 3 + = n=0 x n = 1 1 x if 1 < x < 1. On the left hand side of the equation there is sum containing
More informationSummary so far. Geological structures Earthquakes and their mechanisms Continuous versus block-like behavior Link with dynamics?
Summary so far Geodetic measurements velocities velocity gradient tensor (spatial derivatives of velocity) Velocity gradient tensor = strain rate (sym.) + rotation rate (antisym.) Strain rate tensor can
More informationIndex. p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96
p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96 B 1,94-96 M,94-96 B oro!' 94-96 BIro!' 94-96 I/r, 79 2D linear system, 56 2D FFT, 119 2D Fourier transform, 1, 12, 18,91 2D sinc, 107, 112
More informationHRAM FAULT INTERPRETATION USING MAGPROBE DEPTH ESTIMATES AND NON- TRADITIONAL FILTERING
HRAM FAULT INTERPRETATION USING MAGPROBE DEPTH ESTIMATES AND NON- TRADITIONAL FILTERING Serguei A. Goussev 1, Robert A. Charters 1, Hassan H. Hassan 1, John W. Peirce 1 and James A. Genereux 2 Canadian
More informationGravity Methods (VII) more wrap up
Environmental and Exploration Geophysics II Gravity Methods (VII) more wrap up tom.h.wilson tom.wilson@mail.wvu.edu Department of Geology and Geography West Virginia University Morgantown, WV 0.4 0.35
More information