MEASURING SURFACE CURRENTS USING IR CAMERAS. Background. Optical Current Meter 06/10/2010. J.Paul Rinehimer ESS522
|
|
- Jordan Sutton
- 5 years ago
- Views:
Transcription
1 J.Paul Rinehimer ESS5 Optical Current Meter 6/1/1 MEASURING SURFACE CURRENTS USING IR CAMERAS Background Most in-situ current measurement techniques are based on sending acoustic pulses and measuring the Doppler Shift of pulses that are reflected back to the transducer. In shallow water (below cm) the signal to noise ratio becomes small and the measurements unreliable. Fig. 1 shows velocity data from an Aquadopp ADCP (Acoustic Doppler Current Profilier), an instrument able to profile the velocity in the entire water column. As the tide falls and the instrument is in less than 1 cm of water, no reliable velocity measurements are obtained. Additionally, velocity measurements near the water surface are often contaminated by high noise levels as the acoustic pulses reflect from the water surface, limiting an ADCP s ability to obtain accurate surface measurements. Figure 1: Example Aquadopp ADCP data from Willapa Channel. Top panel shows along-channel current velocities with depth in m/s. Bottom panel shows salinity and temperature at the aquadopp (AQD) and within the seabed nearby. The Optical Current Meter (OCM) [1] is a data processing technique to obtain surface current measurements from video data. This technique uses a D Fourier Transform to convert a sequence of pixel intensities into an energy spectrum in frequency-wave number space and then into velocity-wave number space. A model of the velocity spectrum is fit to the data and surface velocity estimates and statistics can be calculated. This technique is applicable to any video imagery with a signal that corresponds to surface velocities and can be used to calculate 1 / 9
2 surface currents in shallow water and over larger spatial domains than a traditional one-point measurement. An imaging tower was deployed in Willapa Bay, WA with the camera field of view directed towards a tidal flat channel. Both an infrared (IR) and a visual-band camera (Fig. ) were used on the tower, along with a number of in situ instruments located in the camera field of view. The IR instrumentation was being used to measure the surface water temperature as it leaves the tidal channel in order to determine its source (pore water or drainage of surface water off the flats). An additional advantage of this technique is that a single measurement (IR imagery) can provide multiple parameters of interest in the field (velocity, water temperature, sediment temperature/type). Figure : Sample images from the Willapa Bay imaging tower. Visible light image (left) and IR image (right). The f-stop is open too far for the visible and thus the image is over-exposed and washed out. Darker is colder in the IR imagery. Note the cold buoy used as a marker on the tidal flat and the bubbles in the visible imagery. Dark, cold spots on lower right of IR are footprints from the equipment deployment on the previous day. The main channel is to the bottom of the image with the water in the monitored channel flowing towards the main channel. Before the OCM technique can be applied, the raw IR images (Fig. ) have to be rectified to a real-world coordinate system. The details would probably fill an entire course, but the rectification basically applies a matrix transformation to the raw image data transforming the image intensities at pixels into intensities with physical locations (Fig. 3). For this exercise, a subset of the data was chosen, rectified, and interpolated to a meshed grid. The particular data set chosen was a sequence of IR images along a channel where air-filled, and hence cool, bubbles were moving down the channel. Using this data set, we will obtain surface velocities via the OCM method. / 9 ESS5
3 1 6 y (m, north) y (m, north) x (m, east) x (m, east) 8 1 Figure 3: Rectified IR image. The left panel shows the full field of view while the right panel shows a closeup of the channel. Camera is at origin (,). The x and y axis directions have been reversed so that south is towards the top of the figure. Note how resolution decreases and distortion increase away from the camera. The white pattern in the far field is due to a Matlab aliasing the image when it is written. The red line is the location of the timestack. Method The optical current meter works by transforming a timestack of pixel intensities (Fig. 4) via the Fourier Transform. A timestack consists of the intensities of a single line of pixels as it varies with time. The line is chosen such that its primary axis is in the direction of the velocity you wish to measure. Looking at the example timestack (Fig. 4) the bubble streaks are apparent as regular features with a well-defined slope corresponding to their velocity. The OCM method first calculates the Fourier Transform Z (f,k) of the pixel intensities z(t, x) after multiplying it by a Bartlett taper B(t, x) (the Bartlett taper goes from at the edges to 1 in the center). Z (f,k) = B(t, x)z(t, x)e i πf t e i πkx dtdx (1) Then energy spectral density S(f,k) is then computed: S(f,k) = Z (f,k)z (f,k) () where denotes complex conjugation, producing an energy spectrum in frequency-wavenumber (f k) space. Examining this spectrum (Fig. 5), one can see that there is a large peak of energy along a narrow band with a well defined slope. The slope of this peak corresponds to the primary velocity apparent in the signal, but there are a few more steps need to accurately quantify the measurement. ESS5 3 / 9
4 .5 Distance (m) Time (s) Figure 4: Example timestack. Colormap has been inverted so that white is cool and black is warm. Note the long white streaks corresponding to dark bubbles crossing our section. 9 8 Wave Number (m 1 ) v =.66 m/s Frequency (Hz) Figure 5: Log-energy spectral density log(s( f, k)) of the timestack. The peak corresponds to the line with slope v = f /k =.66m/s. 4 / 9 ESS5
5 Next, the spectra is converted to velocity-wavenumber (v k) space following the transformation: S(k, v) = S(f,k) k (3) Multiplication by k is required in order to preserve the variance of the signal. The narrow, sloped band of energy has now been converted into a vertical band in S(v,k) (Fig. 6). A large amount of noise near low wave numbers is also apparent along with the large blank regions in the upperleft and upper-right corresponding to areas beyond the Nyquist frequency (f nyq and wave number (k nyq ). 9 8 Wavenumber (m 1 ) Velocity (m/s Figure 6: Log-energy spectral density log(s(v, k)) of the timestack. In v, k space, the velocity peak is a sharp vertical band. The S(k, v) spectrum can now be integrated to get a one-dimensional velocity spectrum S(v). If we choose the k limits carefully, we can eliminate a lot of the low-frequency noise while still maintaining much of our spectrum. For this exercise, I have chosen k min =1. Then we can compute S(v) by knyq S(v) = S(v,k)dk (4) k min and thus obtain a single spectrum of velocity (Fig. 7). In this example, the background noise is relatively small, so the maximum velocity peak is not far from the mean of the spectrum, however, the long tail towards negative velocities influences the value somewhat. Also, there is a small energy peak near velocity. This peak is due to background noise and spectral leakage and can often be much larger, and in some spectra may dwarf the actual velocity peak. ESS5 5 / 9
6 . Velocity Spectrum, S(v) Mean v =.51 Max v = S (v)/σs (v) Velocity (m/s) Figure 7: Spectrum of S(v). Note how velocity estimates using the mean and the max of the spectra are different due to the noise throughout all v values. In order to account for the background noise field, a model of the velocity spectrum can be developed: S m (v) = S b (v) + S noi se (v) (5) where S m (v) is the model spectrum, S b (v) is the spectrum due to the moving bubbles, and S noi se (v) is a white-noise spectrum (equal at all frequencies). If the bubble spectrum is modeled as a gaussian with mean v, standard deviation σ b, and a signal amplitude of A b then S b (v) = A b exp [ ( ) ] v v while a white noise spectrum with amplitude A noi se will follow f nyq A noi se v f nyq v S noi se (v) = k nyq knyq A noi se v > f (7) nyq k nyq this gives us a set of 4 unknowns v,σ b, A b, and A noi se that can be fit to the model and allow us to derive a useful statistics such as confidence intervals (Fig. 8). Fitting these equations requires non-linear curve fitting and optimization techniques and will not be covered here. Finally, by subwindowing a long timestack, a velocity timeseries can be produced (Fig. 9). σ b (6) 6 / 9 ESS5
7 1.8 S(v) S model (v).6 S(v) Velocity (m/s) Figure 8: Comparsion of S(v) and modelled S model (v) spectra. Black dashed lines indicate ±σ b Velocity (m/s) Time (s) Figure 9: Timeseries of estimated velocities with ±95% confidence intervals. ESS5 7 / 9
8 Conclusions The OCM technique provides robust, remote way to measure the velocity of surface features. The method works well for shallow water environments where in situ observations are sparse and difficult to obtain, as well as with data sets ill-suited for PIV due to the lack of image texture. Improvements to the technique include extending the analysis to a 3-D timestack, and calculating velocities along an orthogonal component, in order to rectify both the velocity magnitude and direction. Ease of setup, fast data processing, and low-cost make the OCM technique a useful tool to gain current information where traditional instrumentation fail. References [1] Chickadel, CC, RA Holman, and MH Freilich. 3 An optical technique for the measurement of longshore currents. Journal of Geophysical Research 18: C11, doi: 1.19/3JC / 9 ESS5
Workhorse ADCP Multi- Directional Wave Gauge Primer
Acoustic Doppler Solutions Workhorse ADCP Multi- Directional Wave Gauge Primer Brandon Strong October, 2000 Principles of ADCP Wave Measurement The basic principle behind wave the measurement, is that
More informationReconstruction of Images Distorted by Water Waves
Reconstruction of Images Distorted by Water Waves Arturo Donate and Eraldo Ribeiro Computer Vision Group Outline of the talk Introduction Analysis Background Method Experiments Conclusions Future Work
More informationMinimizing Noise and Bias in 3D DIC. Correlated Solutions, Inc.
Minimizing Noise and Bias in 3D DIC Correlated Solutions, Inc. Overview Overview of Noise and Bias Digital Image Correlation Background/Tracking Function Minimizing Noise Focus Contrast/Lighting Glare
More informationTh ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data
Th ELI1 12 Joint Crossline Reconstruction and 3D Deghosting of Shallow Seismic Events from Multimeasurement Streamer Data Y.I. Kamil* (Schlumberger), M. Vassallo (Schlumberger), W. Brouwer (Schlumberger),
More informationComputer Vision 2. SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung. Computer Vision 2 Dr. Benjamin Guthier
Computer Vision 2 SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung Computer Vision 2 Dr. Benjamin Guthier 1. IMAGE PROCESSING Computer Vision 2 Dr. Benjamin Guthier Content of this Chapter Non-linear
More informationOptimised corrections for finite-difference modelling in two dimensions
Optimized corrections for 2D FD modelling Optimised corrections for finite-difference modelling in two dimensions Peter M. Manning and Gary F. Margrave ABSTRACT Finite-difference two-dimensional correction
More informationAquadopp Profiler 1 MHz
Small and compact, with up to 25 m current profiling range; option for PUV wave measurements The Aquadopp Profiler is a highly versatile Acoustic Doppler Current Profiler (ADCP) available in four profiling
More informationMetroPro Surface Texture Parameters
MetroPro Surface Texture Parameters Contents ROUGHNESS PARAMETERS...1 R a, R q, R y, R t, R p, R v, R tm, R z, H, R ku, R 3z, SR z, SR z X, SR z Y, ISO Flatness WAVINESS PARAMETERS...4 W a, W q, W y HYBRID
More informationAn imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults
Earth Planets Space, 52, 3 11, 2000 An imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults Takumi Murakoshi 1, Hiroshi Takenaka
More informationRange Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation
Obviously, this is a very slow process and not suitable for dynamic scenes. To speed things up, we can use a laser that projects a vertical line of light onto the scene. This laser rotates around its vertical
More informationSensor Modalities. Sensor modality: Different modalities:
Sensor Modalities Sensor modality: Sensors which measure same form of energy and process it in similar ways Modality refers to the raw input used by the sensors Different modalities: Sound Pressure Temperature
More informationNAME :... Signature :... Desk no. :... Question Answer
Written test Tuesday 19th of December 2000. Aids allowed : All usual aids Weighting : All questions are equally weighted. NAME :................................................... Signature :...................................................
More information9 length of contour = no. of horizontal and vertical components + ( 2 no. of diagonal components) diameter of boundary B
8. Boundary Descriptor 8.. Some Simple Descriptors length of contour : simplest descriptor - chain-coded curve 9 length of contour no. of horiontal and vertical components ( no. of diagonal components
More informationDigital Image Processing. Image Enhancement in the Frequency Domain
Digital Image Processing Image Enhancement in the Frequency Domain Topics Frequency Domain Enhancements Fourier Transform Convolution High Pass Filtering in Frequency Domain Low Pass Filtering in Frequency
More informationDepth Range Accuracy for Plenoptic Cameras
Depth Range Accuracy for Plenoptic Cameras Nuno Barroso Monteiro Institute for Systems and Robotics, University of Lisbon, Portugal Institute for Systems and Robotics, University of Coimbra, Portugal Simão
More informationAn Intuitive Explanation of Fourier Theory
An Intuitive Explanation of Fourier Theory Steven Lehar slehar@cns.bu.edu Fourier theory is pretty complicated mathematically. But there are some beautifully simple holistic concepts behind Fourier theory
More information1.1.1 Orientation Coordinate Systems
1.1.1 Orientation 1.1.1.1 Coordinate Systems The velocity measurement is a vector in the direction of the transducer beam, which we refer to as beam coordinates. Beam coordinates can be converted to a
More informationRelationship between Fourier Space and Image Space. Academic Resource Center
Relationship between Fourier Space and Image Space Academic Resource Center Presentation Outline What is an image? Noise Why do we transform images? What is the Fourier Transform? Examples of images in
More informationPS User Guide Series Dispersion Image Generation
PS User Guide Series 2015 Dispersion Image Generation Prepared By Choon B. Park, Ph.D. January 2015 Table of Contents Page 1. Overview 2 2. Importing Input Data 3 3. Main Dialog 4 4. Frequency/Velocity
More informationEECS 556 Image Processing W 09. Image enhancement. Smoothing and noise removal Sharpening filters
EECS 556 Image Processing W 09 Image enhancement Smoothing and noise removal Sharpening filters What is image processing? Image processing is the application of 2D signal processing methods to images Image
More informationChapter 37. Wave Optics
Chapter 37 Wave Optics Wave Optics Wave optics is a study concerned with phenomena that cannot be adequately explained by geometric (ray) optics. Sometimes called physical optics These phenomena include:
More informationspecular diffuse reflection.
Lesson 8 Light and Optics The Nature of Light Properties of Light: Reflection Refraction Interference Diffraction Polarization Dispersion and Prisms Total Internal Reflection Huygens s Principle The Nature
More informationLecture 8 Object Descriptors
Lecture 8 Object Descriptors Azadeh Fakhrzadeh Centre for Image Analysis Swedish University of Agricultural Sciences Uppsala University 2 Reading instructions Chapter 11.1 11.4 in G-W Azadeh Fakhrzadeh
More informationPlane Wave Imaging Using Phased Array Arno Volker 1
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16409 Plane Wave Imaging Using Phased Array
More informationTIME-FREQUENCY SPECTRA OF MUSIC
TIME-FREQUENCY SPECTRA OF MUSIC James S. Walker Department of Mathematics University of Wisconsin Eau Claire walkerjs@uwec.edu Amanda J. Potts Department of Mathematics University of Nebraska Lincoln Abstract
More informationFFT and Spectrum Analyzer
Objective: FFT and Spectrum Analyzer Understand the frequency domain and some of it's uses Understand the Discrete Fourier Transform Preparation: (i) PC and data acquisition system (ii)fft software Background:
More informationIntroduction to Raman spectroscopy measurement data processing using Igor Pro
Introduction to Raman spectroscopy measurement data processing using Igor Pro This introduction is intended to minimally guide beginners to processing Raman spectroscopy measurement data, which includes
More informationARCHAEOLOGICAL 3D MAPPING USING A LOW-COST POSITIONING SYSTEM BASED ON ACOUSTIC WAVES
ARCHAEOLOGICAL 3D MAPPING USING A LOW-COST POSITIONING SYSTEM BASED ON ACOUSTIC WAVES P. Guidorzi a, E. Vecchietti b, M. Garai a a Department of Industrial Engineering (DIN), Viale Risorgimento 2, 40136
More informationScene Matching on Imagery
Scene Matching on Imagery There are a plethora of algorithms in existence for automatic scene matching, each with particular strengths and weaknesses SAR scenic matching for interferometry applications
More informationTh LHR5 08 Multi-modal Surface Wave Inversion and Application to North Sea OBN Data
Th LHR5 08 Multi-modal Surface Wave Inversion and pplication to North Sea ON Data S. Hou (CGG), D. Zheng (CGG), X.G. Miao* (CGG) & R.R. Haacke (CGG) SUMMRY Surface-wave inversion (SWI) for S-wave velocity
More informationChapter 37. Interference of Light Waves
Chapter 37 Interference of Light Waves Wave Optics Wave optics is a study concerned with phenomena that cannot be adequately explained by geometric (ray) optics These phenomena include: Interference Diffraction
More informationPart 1: Calculating amplitude spectra and seismic wavelets
Geology 554 Environmental and Exploration Geophysics II Generating synthetic seismograms The simple in-class exercise undertaken last time illustrates the relationship between wavelet properties, interval
More informationAccurately measuring 2D position using a composed moiré grid pattern and DTFT
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 Accurately measuring 2D position using a composed moiré grid pattern and DTFT S. Van
More informationOptics: Laser Light Show Student Advanced Version
Optics: Laser Light Show Student Advanced Version In this lab, you will explore the behavior of light. You will observe reflection and refraction of a laser beam in jello, and use a diffraction pattern
More informationMotion Analysis. Motion analysis. Now we will talk about. Differential Motion Analysis. Motion analysis. Difference Pictures
Now we will talk about Motion Analysis Motion analysis Motion analysis is dealing with three main groups of motionrelated problems: Motion detection Moving object detection and location. Derivation of
More informationTime in ms. Chan. Summary
Marine source signature estimation with dual near-field hydrophones Rob Telling*, Sergio Grion Stuart Denny & R. Gareth Williams, Shearwater GeoServices Summary We derive marine seismic signatures using
More informationNew Results on the Omega-K Algorithm for Processing Synthetic Aperture Radar Data
New Results on the Omega-K Algorithm for Processing Synthetic Aperture Radar Data Matthew A. Tolman and David G. Long Electrical and Computer Engineering Dept. Brigham Young University, 459 CB, Provo,
More informationSUMMARY: DISTINCTIVE IMAGE FEATURES FROM SCALE- INVARIANT KEYPOINTS
SUMMARY: DISTINCTIVE IMAGE FEATURES FROM SCALE- INVARIANT KEYPOINTS Cognitive Robotics Original: David G. Lowe, 004 Summary: Coen van Leeuwen, s1460919 Abstract: This article presents a method to extract
More informationCoherent partial stacking by offset continuation of 2-D prestack data
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 124 Coherent partial stacking by offset continuation of 2-D prestack data Nizar Chemingui and Biondo Biondi 1 ABSTRACT Previously, we introduced
More informationAnalysis of different interpolation methods for uphole data using Surfer software
10 th Biennial International Conference & Exposition P 203 Analysis of different interpolation methods for uphole data using Surfer software P. Vohat *, V. Gupta,, T. K. Bordoloi, H. Naswa, G. Singh and
More informationHigh Resolution Imaging by Wave Equation Reflectivity Inversion
High Resolution Imaging by Wave Equation Reflectivity Inversion A. Valenciano* (Petroleum Geo-Services), S. Lu (Petroleum Geo-Services), N. Chemingui (Petroleum Geo-Services) & J. Yang (University of Houston)
More informationWhen is anti-aliasing needed in Kirchhoff migration? a
When is anti-aliasing needed in Kirchhoff migration? a a Published in SEP report, 80, 467-476 (1994) Dimitri Bevc and David E. Lumley ABSTRACT We present criteria to determine when numerical integration
More informationBack Projection of the Okhotsk earthquake on the vertical plane. To test the candidate geometry of a vertical fault, we perform the BP analysis on the
Supplementary information of The 213 Okhotsk Deep-Focus Earthquake: Rupture Beyond the Metastable Olivine Wedge and Thermally-Controlled Rise Time near the Edge of a Slab Back Projection of the Okhotsk
More informationFourier analysis of low-resolution satellite images of cloud
New Zealand Journal of Geology and Geophysics, 1991, Vol. 34: 549-553 0028-8306/91/3404-0549 $2.50/0 Crown copyright 1991 549 Note Fourier analysis of low-resolution satellite images of cloud S. G. BRADLEY
More informationDownloaded 09/01/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Random Noise Suppression Using Normalized Convolution Filter Fangyu Li*, Bo Zhang, Kurt J. Marfurt, The University of Oklahoma; Isaac Hall, Star Geophysics Inc.; Summary Random noise in seismic data hampers
More informationSingle-point current meter designed for very long-term deployments
Single-point current meter designed for very long-term deployments With all the features and capabilities of the standard Aquadopp, the deepwater current meter has been used and proven by oceanographers
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Review of Motion Modelling and Estimation Introduction to Motion Modelling & Estimation Forward Motion Backward Motion Block Motion Estimation Motion
More informationEstimating the wavelength composition of scene illumination from image data is an
Chapter 3 The Principle and Improvement for AWB in DSC 3.1 Introduction Estimating the wavelength composition of scene illumination from image data is an important topics in color engineering. Solutions
More informationHigh-performance mean currents and turbulence, wave height and direction
High-performance mean currents and turbulence, wave height and direction The ADCP is the optimal tool for turbulence measurements. With a maximum sampling frequency of 16 Hz, it gives the scientific community
More informationAssignment 2. Due Feb 3, 2012
EE225E/BIOE265 Spring 2012 Principles of MRI Miki Lustig Assignment 2 Due Feb 3, 2012 1. Read Nishimura Ch. 3 2. Non-Uniform Sampling. A student has an assignment to monitor the level of Hetch-Hetchi reservoir
More informationPop Quiz 1 [10 mins]
Pop Quiz 1 [10 mins] 1. An audio signal makes 250 cycles in its span (or has a frequency of 250Hz). How many samples do you need, at a minimum, to sample it correctly? [1] 2. If the number of bits is reduced,
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 informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More informationMotion Estimation. There are three main types (or applications) of motion estimation:
Members: D91922016 朱威達 R93922010 林聖凱 R93922044 謝俊瑋 Motion Estimation There are three main types (or applications) of motion estimation: Parametric motion (image alignment) The main idea of parametric motion
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationInternational Symposium Non-Destructive Testing in Civil Engineering (NDT-CE) September 15-17, 2015, Berlin, Germany
More Info at Open Access Database www.ndt.net/?id=18355 Effect of Surface Unevenness on In Situ Measurements and Theoretical Simulation in Non-Contact Surface Wave Measurements Using a Rolling Microphone
More informationRemote sensing In situ T, S and tracers Velocity Observing systems. Course url: Talley SIO 210 (2013)
SIO 210: I. Observational methods and II. Data analysis (combined single lecture) Fall 2013 I. Observations Reading: DPO 6.1, S16.1, S16.4, S16.5, S16.9 Remote sensing In situ T, S and tracers Velocity
More informationREMI and SPAC: A Comparison
REMI and SPAC: A Comparison Don Zhao Geogiga Technology Corp. Outline Overview Comparison by Synthetic Data Linear Array for Omni-Directional Sources Field Setup of Triangular Array for SPAC Conclusion
More informationIntroduction to Digital Image Processing
Fall 2005 Image Enhancement in the Spatial Domain: Histograms, Arithmetic/Logic Operators, Basics of Spatial Filtering, Smoothing Spatial Filters Tuesday, February 7 2006, Overview (1): Before We Begin
More informationWinADCP User s Guide. P/N (November 2001) RD Instruments Acoustic Doppler Solutions
WinADCP User s Guide P/N 957-6158-00 (November 2001) RD Instruments Acoustic Doppler Solutions Table of Contents 1 Introduction... 1 1.1 System Requirements...1 1.2 Software Installation...2 2 Using WinADCP...
More informationLimitations of Rayleigh Dispersion Curve Obtained from Simulated Ambient Vibrations
Limitations of Rayleigh Dispersion Curve Obtained from Simulated Ambient Vibrations KEN TOKESHI and CARLOS CUADRA Department of Architecture and Environment Systems Akita Prefectural University 84-4 Ebinokuchi,
More informationGeometric Rectification of Remote Sensing Images
Geometric Rectification of Remote Sensing Images Airborne TerrestriaL Applications Sensor (ATLAS) Nine flight paths were recorded over the city of Providence. 1 True color ATLAS image (bands 4, 2, 1 in
More informationOriginal grey level r Fig.1
Point Processing: In point processing, we work with single pixels i.e. T is 1 x 1 operator. It means that the new value f(x, y) depends on the operator T and the present f(x, y). Some of the common examples
More informationImporting and processing a DGGE gel image
BioNumerics Tutorial: Importing and processing a DGGE gel image 1 Aim Comprehensive tools for the processing of electrophoresis fingerprints, both from slab gels and capillary sequencers are incorporated
More informationExample Applications of A Stochastic Ground Motion Simulation Methodology in Structural Engineering
Example Applications of A Stochastic Ground Motion Simulation Methodology in Structural Engineering S. Rezaeian & N. Luco U.S. Geological Survey, Golden, CO, USA ABSTRACT: Example engineering applications
More informationSUMMARY. denoise the original data at each iteration. This can be
A comparison of D reconstruction methods Aaron Stanton*, Nadia Kreimer, David Bonar, Mostafa Naghizadeh, and Mauricio Sacchi, Department of Physics, University of Alberta SUMMARY A comparison is made between
More informationRefraction Full-waveform Inversion in a Shallow Water Environment
Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui
More informationDetermining satellite rotation rates for unresolved targets using temporal variations in spectral signatures
Determining satellite rotation rates for unresolved targets using temporal variations in spectral signatures Joseph Coughlin Stinger Ghaffarian Technologies Colorado Springs, CO joe.coughlin@sgt-inc.com
More informationChapter 2: Wave Optics
Chapter : Wave Optics P-1. We can write a plane wave with the z axis taken in the direction of the wave vector k as u(,) r t Acos tkzarg( A) As c /, T 1/ and k / we can rewrite the plane wave as t z u(,)
More informationAnswer Notes Marks 1 (a) MP1. pitch is frequency; allow it for pitch
(a) MP. pitch is frequency; allow it for pitch MP. any one of: whether sound/note sounds high or low; high sound has high frequency ORA; ignore references to amplitude, wavelength allow vibrates more often
More informationA tutorial on optimizing time domain finite-difference schemes: "Beyond Holberg"
Stanford Exploration Project, Report 129, May 6, 2007, pages 33 43 A tutorial on optimizing time domain finite-difference schemes: "Beyond Holberg" John T. Etgen ABSTRACT Time-domain finite-difference
More informationFINAL REPORT: Sonar Detection of Buried Targets at Subcritical Grazing Angles: APL-UW Component. Approved ft r Public Release Distribution Unlimited
FINAL REPORT: Sonar Detection of Buried Targets at Subcritical Grazing Angles: APL-UW Component DISTRIBUTION STATEMENT A Approved ft r Public Release Distribution Unlimited Kevin L. Williams, Eric I. Thorsos,
More informationVERY LARGE TELESCOPE 3D Visualization Tool Cookbook
European Organisation for Astronomical Research in the Southern Hemisphere VERY LARGE TELESCOPE 3D Visualization Tool Cookbook VLT-SPE-ESO-19500-5652 Issue 1.0 10 July 2012 Prepared: Mark Westmoquette
More informationImage Restoration and Reconstruction
Image Restoration and Reconstruction Image restoration Objective process to improve an image, as opposed to the subjective process of image enhancement Enhancement uses heuristics to improve the image
More informationGeogiga Seismic Pro 8.0 Release Notes
Geogiga Seismic Pro 8.0 Release Notes Copyright 2015, All rights reserved. Table of Contents Introduction...1 Part 1 General Enhancements...3 Trace Display...4 Color Section Display...6 Curve Plotting...8
More informationCRISM (Compact Reconnaissance Imaging Spectrometer for Mars) on MRO. Calibration Upgrade, version 2 to 3
CRISM (Compact Reconnaissance Imaging Spectrometer for Mars) on MRO Calibration Upgrade, version 2 to 3 Dave Humm Applied Physics Laboratory, Laurel, MD 20723 18 March 2012 1 Calibration Overview 2 Simplified
More informationFast Camera Image Analysis. B. A. Grierson
Fast Camera Image Analysis B. A. Grierson 9.28.2008 Introduction A high speed camera is used to image the visible light in. The camera is begins recording at the specified trigger time, and records until
More informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
More informationFive Dimensional Interpolation:exploring different Fourier operators
Five Dimensional Interpolation:exploring different Fourier operators Daniel Trad CREWES-University of Calgary Summary Five-Dimensional interpolation has become a very popular method to pre-condition data
More informationImage Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi
Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi 1. Introduction The choice of a particular transform in a given application depends on the amount of
More informationfraction of Nyquist
differentiator 4 2.1.2.3.4.5.6.7.8.9 1 1 1/integrator 5.1.2.3.4.5.6.7.8.9 1 1 gain.5.1.2.3.4.5.6.7.8.9 1 fraction of Nyquist Figure 1. (top) Transfer functions of differential operators (dotted ideal derivative,
More informationSampling and Reconstruction
Sampling and Reconstruction Sampling and Reconstruction Sampling and Spatial Resolution Spatial Aliasing Problem: Spatial aliasing is insufficient sampling of data along the space axis, which occurs because
More information1.12 Optimal Filters (Wiener Filters)
Random Data 75 1.12 Optimal Filters (Wiener Filters) In calculating turbulent spectra we frequently encounter a noise tail just as the spectrum roles down the inertial subrange (-5/3 slope region) toward
More informationFull waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC)
of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) Summary (FWI) is a promising technology for next generation model building. However, it faces with many
More informationThe Argonaut-XR System and the Multi-Cell Feature
The Argonaut-XR System and the Multi-Cell Feature As a current profiling system, the Argonaut-XR offers exceptional value for near shore deployments. Its small size, rugged build quality, and flexible
More informationExperimental and Numerical Analysis of Near Wall Flow at the Intake Valve and its Influence on Large-Scale Fluctuations
Experimental and Numerical Analysis of Near Wall Flow at the Intake Valve and its Influence on Large-Scale Fluctuations Frank Hartmann, Stefan Buhl, Florian Gleiß, Christian Hasse Philipp Barth, Martin
More informationENHANCEMENT OF DIFFUSERS BRDF ACCURACY
ENHANCEMENT OF DIFFUSERS BRDF ACCURACY Grégory Bazalgette Courrèges-Lacoste (1), Hedser van Brug (1) and Gerard Otter (1) (1) TNO Science and Industry, Opto-Mechanical Instrumentation Space, P.O.Box 155,
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Aliasing Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 16, 2014 Bee Bednar (Panorama Technologies) Seismic
More informationScanner Parameter Estimation Using Bilevel Scans of Star Charts
ICDAR, Seattle WA September Scanner Parameter Estimation Using Bilevel Scans of Star Charts Elisa H. Barney Smith Electrical and Computer Engineering Department Boise State University, Boise, Idaho 8375
More informationNew demands on the modeling of the continental shelf require the need for 3D graphics connected to seismic interpretation of seabed sediments.
A New Approach to Mapping of the Continental Shelf Olav Egeland Manager of GeoCap AS, Norway. Website: www.geocap.no Email: olav.egeland@geocap.no Abstract New demands on the modeling of the continental
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 informationPrincipal Component Image Interpretation A Logical and Statistical Approach
Principal Component Image Interpretation A Logical and Statistical Approach Md Shahid Latif M.Tech Student, Department of Remote Sensing, Birla Institute of Technology, Mesra Ranchi, Jharkhand-835215 Abstract
More information4. Refraction. glass, air, Perspex and water.
Mr. C. Grima 11 1. Rays and Beams A ray of light is a narrow beam of parallel light, which can be represented by a line with an arrow on it, in diagrams. A group of rays makes up a beam of light. In laboratory
More informationGeometric Correction
CEE 6150: Digital Image Processing Geometric Correction 1 Sources of Distortion Sensor Characteristics optical distortion aspect ratio non-linear mirror velocity detector geometry & scanning sequence Viewing
More informationUniversity of Technology Building & Construction Department / Remote Sensing & GIS lecture
5. Corrections 5.1 Introduction 5.2 Radiometric Correction 5.3 Geometric corrections 5.3.1 Systematic distortions 5.3.2 Nonsystematic distortions 5.4 Image Rectification 5.5 Ground Control Points (GCPs)
More informationAdvances in radial trace domain coherent noise attenuation
Advances in radial trace domain coherent noise attenuation ABSTRACT David C. Henley* CREWES, Department of Geology and Geophysics University of Calgary, Calgary, AB henley@crewes.org The radial trace transform,
More informationImplementing non-stationary filtering in time and in Fourier domain
Stanford Exploration Project, Report 111, June 9, 2002, pages 347 361 Implementing non-stationary filtering in time and in Fourier domain Gabriel Alvarez 1 ABSTRACT Non-stationary filtering of seismic
More informationComments on wavefield propagation using Reverse-time and Downward continuation
Comments on wavefield propagation using Reverse-time and Downward continuation John C. Bancroft ABSTRACT Each iteration a of Full-waveform inversion requires the migration of the difference between the
More informationThe 2014 IEEE GRSS Data Fusion Contest
The 2014 IEEE GRSS Data Fusion Contest Description of the datasets Presented to Image Analysis and Data Fusion Technical Committee IEEE Geoscience and Remote Sensing Society (GRSS) February 17 th, 2014
More informationAP* Optics Free Response Questions
AP* Optics Free Response Questions 1978 Q5 MIRRORS An object 6 centimeters high is placed 30 centimeters from a concave mirror of focal length 10 centimeters as shown above. (a) On the diagram above, locate
More information