Array Shape Tracking Using Active Sonar Reverberation
|
|
- Rolf Freeman
- 6 years ago
- Views:
Transcription
1 Lincoln Laboratory ASAP-2003 Worshop Array Shape Tracing Using Active Sonar Reverberation Vijay Varadarajan and Jeffrey Kroli Due University Department of Electrical and Computer Engineering Durham, NC Supported by the ONR 321US 6.1 Program 1
2 Array Shape Tracing from Active Sonar Clutter OBJECTIVE: To estimate and trac the shape of a towed distorted-linear array using reverberation from active sonar pings as a distributed source of opportunity. BACKGROUND: Array shape estimation using heading and depth sensors has been previously developed, e.g. Gray et al.(1993), by applying Kalman filtering with a state equation derived from spatio-temporal discretization of the simplified Paidoussis equation (water-pulley model) for array motion. Acoustic sources of opportunity have been employed for array shape estimation by fitting transverse element displacements to measuring inter-element phase differences e.g. Owsley (1980). For mid-frequency active sonar arrays the above methods may be precluded since it may be infeasible to instrument the array with a sufficient number of heading sensors and strong point sources of opportunity may not always be available in the presence of strong reverberation. Clutter has previously been used for element gain and phase calibration of a uniform linear airborne radar array (Robey, Fuhrmann, and Krisch, 1994) but not for array shape estimation. We propose constrained maximum lielihood array shape estimation from clutter (ASEC) which uses an array shape-dependent model of spatially-distributed, Doppler-spread reverberation. Array shape parameters are traced within and across pings by using ML ASEC heading estimates as input to a Kalman filter whose state equation incorporates a dynamical model for array motion. 2
3 Motivation for Array Shape Estimation Detection Detection of of slowly-moving slowly-moving targets targets with with Doppler-sensitive Doppler-sensitive waveforms waveforms is is limited limited by by beamformer beamformer sidelobes. sidelobes. Uncompensated Uncompensated array array distortion distortion causes causes increased increased clutter clutter ran ran and andtarget masing. masing. Simulation Simulation example example angle-doppler angle-doppler spectrum spectrum from from perfectly perfectly compensated compensated array array (left) (left) vs. vs. uncompensated uncompensated distorted distortedarray (right). (right). Note Note target target mased mased by by reverberation reverberation due due to to beamformer beamformer sidelobes sidelobes of of uncompensated uncompensated array array Target 3
4 Array Shape Estimation from Clutter Concept Constrained Constrained ML ML ASEC ASEC assimilates assimilates clutter clutter data data with with a a few few heading heading sensor sensor outputs outputs to to obtain obtain array array shape shape parameters parameters from from each each sub-pulse sub-pulse length length space-time space-time snapshot. snapshot. Kalman Kalman tracer tracer filters filters sub-cpi sub-cpi shape shape coefficients coefficients using using a a model model for for array array dynamics. dynamics. Example SWAC-3 Reverberation Return Sub-CPI space-time clutter snapshots r Heading sensors g r r 1 Constrained ML ASEC Tow-cable Excitation u µ Kalman Tracer with Array Dynamics estimated heading vector Filtered array shape 4
5 Reverberation Received at a Distorted Array th The clutter data from the n sensor located at ( x, y ) at time t mt for a distorted array moving with velocity v rnm, l l (, ) l e a n n m r along the x direction is modeled as: 2 2va j (sin cos x cos cos y ) j2 sin cos mt l n l n l r where, is the complex Gaussian scatter amplitude from clutter at azimuth,, and multipath elevation,. l Sensor coordinates can be expressed in terms of heading and inter-element spacing n as x d cos, y d sin, 0 n N 1, x,y 0,0. n i n i i 1 i 1 n 0 0 d 5
6 Space-Time Reverberation Model The space-time data snapshot at time t j max consisting of clutter from all azimuths,, 1 i N and elevation angles,1 j N, can be written as i 1 N 1 T where µ, Vµ v ( 1, 1, µ ) v ( N, N, µ ) is the clutter MN 1 2 steering matrix, represents unnown scattering, and noise e has covariance I. r V µ e th i N j V µ v i j µ b ij a i j µ i j The ( 1) column of is (,, ),, which represents the return from a single clutter patch at location,. e MN Array shape parameters µ are assumed constant over the sub-cpi. Array shape parameters µ estimated by fitting the low ran ( MN) clutter subspace of V µ to the observed received space-time snapshot, r. 6
7 Clutter Eigenvalues for a Distorted Linear Array Eigenvalues Eigenvaluesof of the the asymptotic asymptotic space-time space-time covariance covariance matrix matrix for for a a modestly modestly distorted distorted array array versus versus ideal ideal uniform uniform linear linear array array for for N=32 N=32 and and M=5. M=5. Observe Observe that that ran ran inflation inflation with with 00 db dbbaclobes baclobesagrees agrees with with 2-D 2-D Brennan s Brennan s rule rule [Varadarajan [Varadarajanand andkroli, 2002]. 2002] Eigen Values (db) D Brennan's rule 2-D Brennan's rule -300 Linear Array Distorted Array (No-Multipath) Distorted Array (Multipath) Eigen Value Index 7
8 Maximum Lielihood Array Shape Estimation The N -1 array shape headings can be parameterized by a low dimensional subspace using N 1 ( L N ) arbitrary but nown shape basis as µ =. The lielihood function for for the model reduces to minimizing the projection: ˆ argmin f ( ) P ( ) r 2 where P ( ) I V ( ) V ( ) H 1 orthogonal complement of the clutter subspace. H V ( ) V ( ) is the projection matrix onto the Evaluation of log( f( )) (right) for L 2 demonstrates left-right ambiguity in shape estimate corresponding to mirrored solutions about the array axis
9 Constrained ML ASEC Algorithm The left-right shape ambiguity can be resolved with nowledge of the position or heading of a single off-axis sensor which can be expressed as a linear constraint on the shape basis coefficients as T c g. The ML ASEC estimate can be obtained iteratively using a gradient projection approach (e.g. Frost (1972)): j 1 j j ˆ ˆ c P c f' ˆ j T 1 th where c P c ccc ( ) g is the projection of the current solution (j iteration) onto the constraint j subspace, Pc I-Pc, f( ) -1 f' H H H = =-2 r P ( ) V i( ) V ( ) V ( ) V ( ) r, 0< 1. i i j ˆ V ( ) The matrix V ( ) can be computed analytically by assuming the temporal component of the i i steering matrix is independent of array shape distortion over a sub-cpi, and an analytic form for V ( ) obtained by judicious sampling of the clutter wavenumber spectrum. The ran of V ( ) is assumed constant and can be chosen so that it does not change over the set of possible array distortions. 9
10 Incorporating Array Dynamics into ASEC A Kalman filter can be used to trac the array shape coefficients across multiple sub-cpi intervals during a ping. The state vector µ can be related to µ using the water-pulley model: 1 µ Fµ u v 1 1 H 1 N 2 2 v1 and represents white state noise with v I where F (1 ) I L is the state transition matrix, u 0 is the tow-cable driving term, 1. It can be shown that the displacement velocity along the array is determined by 2. An observation equation can be defined using the MLE z µ ˆ from each sub-cpi: z µ v 2 2 where v2 is the measurement noise with covariance v I. 2 The predicted MMSE array shape is used to initialize the next ASEC MLE search. 10
11 Array Dynamical Model Validation Using TB-29 Data Real Real heading heading sensor sensor data data from from a a TB-29 TB-29 array array (left) (left) vs. vs. time time (sec), (sec), courtesy courtesy of of Bruce Bruce Newhall Newhall at at APL/JHU, APL/JHU, indicates indicates water water pulley pulley model model valid valid for for mild mild maneuvers. maneuvers. Transverse Transverse distortion distortion over over CPI CPI (right) (right) validates validates assumption assumption of of rigid rigid short-time short-time motion. motion. 11
12 TB-29 Heading Sensor Eigenbasis Empirical Empirical orthogonal orthogonal basis basis vectors vectors (left) (left) for for headings headings derived derived from from 66 heading heading sensor sensor TB-29 TB-29 data data (interpolated (interpolated to to sensors) sensors) demonstrate demonstrate characteristic characteristic behavior. behavior. TB-29 TB-29 heading heading sensor sensor eigen-spectrum eigen-spectrum (right) (right) indicates indicates most most of of the the variation variation captured captured by by less less than than 44 modes. modes. 12
13 ASEC versus Heading Only Shape Estimation Illustrative Illustrative simulation simulation comparison comparison of of ASEC ASEC versus versus headings headings only only tracing tracing of of a a mild mild maneuver maneuver as as seen seen in in the the middle middle of of the the array array based based on on TB-29 TB-29 heading heading data data (left). (left). Simulation Simulation assumes assumes ASEC ASEC with with s. s. moving moving window window sub-cpi, sub-cpi, 44 EOF EOF basis, basis, spacetimtime snapshot snapshot dimension dimension and and clutter clutter ran ran of of space- Heading Heading (left) (left) and and error error (right) (right) for for ASEC ASEC and and 22 heading-sensor heading-sensor tracing tracing as as in in MFTA. MFTA. 13
14 RMS Heading Error along the Array verus Time Simulated Simulated RMSE RMSE for for 2-sensor 2-sensor headings-only headings-only tracing tracing error error (left) (left) and and ASEC ASEC (right) (right) over over second second ping ping along along the the array, array, ensemble ensemble averaged averaged over over initial initial condition condition perturbations perturbations and and clutter clutter realizations. realizations. Observe Observe propagation propagation of of error error down down the the array array with with time time using using water-pulley water-pulley model model (left). (left). Error Error largest largest in in the the middle middle of of the the array array since since headings headingsnown at at the the end end and and at at tow-point tow-point offset offset from from first first sensor. sensor. 14
15 Angle-Doppler Spectrum for Distorted Array Impact Impact of of ASEC ASEC tracing tracing on on sonar sonar performance performance illustrated illustrated by by conventional conventional beamforming-doppler beamforming-doppler spectra spectra outputs outputs simulated simulated for for continuous continuous array array motion. motion. Time-evolving Time-evolving simulated simulated array array shapes shapes based based on on scaled scaled TB-29 TB-29 motion motion model model (left). (left). Angle-Doppler Angle-Doppler spectrum spectrum with with db db target target for for perfect perfect array array compensation compensation (right). (right). Target 15
16 ASEC vs. Headings-only Tracing Angle-Doppler Spectra Spectrum Spectrum for for conventional conventional array array tracing tracing for for WP WP model model with with 22 heading heading sensors sensors (left). (left). Spectrum Spectrum for for ASEC ASEC tracing tracing with with 44 basis basis functions, functions, s. s. sub-cpi sub-cpi sliding sliding window, window, and and CNR CNR = db db (right). (right). Target Target visible visible at at zero-doppler zero-doppler and and degree degree bearing. bearing. Target 16
17 ASEC Summary and Future Wor Distortion of nominally linear arrays will often result in a substantial increase in spatial sidelobe levels which can mas slowly-moving targets. ML estimation of array shape is facilitated by fitting a reduced-ran reverberation model to a sliding window of sub-cpi space-time snapshots over the extent of each sonar return. Shape ambiguities can be resolved by efficiently maximizing lielihood subject to a constraint which incorporates measurements from at least a single heading sensor. Constrained ML ASEC heading estimates are used as inputs to a Kalman tracer incorporating an array dynamical model and driven by a tow-cable heading sensor output. Simulations using array motion scaled from real TB-29 heading data suggests that ASEC tracing can facilitate improved array compensation compared to headings-only tracing. Future wor will include ASEC performance evaluation with real 53C/MFTA data. 17
Passive Differential Matched-field Depth Estimation of Moving Acoustic Sources
Lincoln Laboratory ASAP-2001 Workshop Passive Differential Matched-field Depth Estimation of Moving Acoustic Sources Shawn Kraut and Jeffrey Krolik Duke University Department of Electrical and Computer
More informationOptimum Array Processing
Optimum Array Processing Part IV of Detection, Estimation, and Modulation Theory Harry L. Van Trees WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION Preface xix 1 Introduction 1 1.1 Array Processing
More informationRadar Detection Improvement by Integration of Multi- Object Tracking
Radar Detection Improvement by Integration of Multi- Object Tracing Lingmin Meng Research and Technology Center Robert Bosch Corp. Pittsburgh, PA, U.S.A. lingmin.meng@rtc.bosch.com Wolfgang Grimm Research
More informationNew Approaches for EEG Source Localization and Dipole Moment Estimation. Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine
New Approaches for EEG Source Localization and Dipole Moment Estimation Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine Outline Motivation why EEG? Mathematical Model equivalent
More informationm=[a,b,c,d] T together with the a posteriori covariance
zimuthal VO analysis: stabilizing the model parameters Chris Davison*, ndrew Ratcliffe, Sergio Grion (CGGVeritas), Rodney Johnston, Carlos Duque, Musa Maharramov (BP). solved using linear least squares
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 informationRuch (Motion) Rozpoznawanie Obrazów Krzysztof Krawiec Instytut Informatyki, Politechnika Poznańska. Krzysztof Krawiec IDSS
Ruch (Motion) Rozpoznawanie Obrazów Krzysztof Krawiec Instytut Informatyki, Politechnika Poznańska 1 Krzysztof Krawiec IDSS 2 The importance of visual motion Adds entirely new (temporal) dimension to visual
More informationWorkhorse 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 informationResearch Article Cross Beam STAP for Nonstationary Clutter Suppression in Airborne Radar
Antennas and Propagation Volume 213, Article ID 27631, 5 pages http://dx.doi.org/1.1155/213/27631 Research Article Cross Beam STAP for Nonstationary Clutter Suppression in Airborne Radar Yongliang Wang,
More informationAdaptive de-ghosting by kurtosis maximisation. Sergio Grion, Rob Telling, Janet Barnes, Dolphin Geophysical. Summary
Adaptive de-ghosting by kurtosis maximisation Sergio Grion, Rob Telling, Janet Barnes, Dolphin Geophysical Summary We discuss the reasons for adaptive de-ghosting and its advantages. The success of broadband
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 informationL1 REGULARIZED STAP ALGORITHM WITH A GENERALIZED SIDELOBE CANCELER ARCHITECTURE FOR AIRBORNE RADAR
L1 REGULARIZED STAP ALGORITHM WITH A GENERALIZED SIDELOBE CANCELER ARCHITECTURE FOR AIRBORNE RADAR Zhaocheng Yang, Rodrigo C. de Lamare and Xiang Li Communications Research Group Department of Electronics
More informationFeature Tracking and Optical Flow
Feature Tracking and Optical Flow Prof. D. Stricker Doz. G. Bleser Many slides adapted from James Hays, Derek Hoeim, Lana Lazebnik, Silvio Saverse, who 1 in turn adapted slides from Steve Seitz, Rick Szeliski,
More informationAnalysis of Functional MRI Timeseries Data Using Signal Processing Techniques
Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Sea Chen Department of Biomedical Engineering Advisors: Dr. Charles A. Bouman and Dr. Mark J. Lowe S. Chen Final Exam October
More informationMotion Tracking and Event Understanding in Video Sequences
Motion Tracking and Event Understanding in Video Sequences Isaac Cohen Elaine Kang, Jinman Kang Institute for Robotics and Intelligent Systems University of Southern California Los Angeles, CA Objectives!
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 informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: EKF-based SLAM Dr. Kostas Alexis (CSE) These slides have partially relied on the course of C. Stachniss, Robot Mapping - WS 2013/14 Autonomous Robot Challenges Where
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 informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationAn Adaptive Eigenshape Model
An Adaptive Eigenshape Model Adam Baumberg and David Hogg School of Computer Studies University of Leeds, Leeds LS2 9JT, U.K. amb@scs.leeds.ac.uk Abstract There has been a great deal of recent interest
More informationAN acoustic array consists of a number of elements,
APPLICATION NOTE 1 Acoustic camera and beampattern Jørgen Grythe, Norsonic AS, Oslo, Norway Abstract The wavenumber-frequency response of an array describes the response to an arbitrary plane wave both
More informationLessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield?
Lessons in Spatial Sampling Or... Does Anybody Know the Shape of a Wavefield? Norm Cooper - Mustagh Resources Ltd., Calgary, Canada In a homogeneous, isotropic, infinite half space (remember those regimes
More informationGround target tracking with signal adaptive measurement error covariance matrix
Ground target tracing with signal adaptive measurement error covariance matrix Robert Kohlleppel Fraunhofer Institute for High Frequency Physics and Radar Techniques FHR Neuenahrer Str. 0 53343 Wachtberg,
More information2D Video Stabilization for Industrial High-Speed Cameras
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 15, No 7 Special Issue on Information Fusion Sofia 2015 Print ISSN: 1311-9702; Online ISSN: 1314-4081 DOI: 10.1515/cait-2015-0086
More informationMobile Robotics. Mathematics, Models, and Methods. HI Cambridge. Alonzo Kelly. Carnegie Mellon University UNIVERSITY PRESS
Mobile Robotics Mathematics, Models, and Methods Alonzo Kelly Carnegie Mellon University HI Cambridge UNIVERSITY PRESS Contents Preface page xiii 1 Introduction 1 1.1 Applications of Mobile Robots 2 1.2
More informationI. INTRODUCTION. J. Acoust. Soc. Am. 113 (5), May /2003/113(5)/2719/13/$ Acoustical Society of America
Source motion mitigation for adaptive matched field processing a) Lisa M. Zurk, b) Nigel Lee, and James Ward MIT Lincoln Laboratory, 244 Wood Street, Lexington, Massachusetts 02420 Received 30 November
More informationENHANCED RADAR IMAGING VIA SPARSITY REGULARIZED 2D LINEAR PREDICTION
ENHANCED RADAR IMAGING VIA SPARSITY REGULARIZED 2D LINEAR PREDICTION I.Erer 1, K. Sarikaya 1,2, H.Bozkurt 1 1 Department of Electronics and Telecommunications Engineering Electrics and Electronics Faculty,
More informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationCoupling of surface roughness to the performance of computer-generated holograms
Coupling of surface roughness to the performance of computer-generated holograms Ping Zhou* and Jim Burge College of Optical Sciences, University of Arizona, Tucson, Arizona 85721, USA *Corresponding author:
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 informationMotion. 1 Introduction. 2 Optical Flow. Sohaib A Khan. 2.1 Brightness Constancy Equation
Motion Sohaib A Khan 1 Introduction So far, we have dealing with single images of a static scene taken by a fixed camera. Here we will deal with sequence of images taken at different time intervals. Motion
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationA Multi-hypothesis Filter for Passive Tracking of Surface and Sub-Surface Target Localization
A Multi-hypothesis Filter for Passive racing of Surface and Sub-Surface arget Localization Wilbur J., Hyland J., Dobec G. and Sermarini C. Code HS3, NSWCPC West Hwy 98 Panama City Beach, FL 38 joellen.wilbur@navy.mil
More informationA Multiscale Nested Modeling Framework to Simulate the Interaction of Surface Gravity Waves with Nonlinear Internal Gravity Waves
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Multiscale Nested Modeling Framework to Simulate the Interaction of Surface Gravity Waves with Nonlinear Internal Gravity
More informationUnconstrained Beamforming : A Versatile Approach.
Unconstrained Beamforming : A Versatile Approach. Monika Agrawal Centre for Applied Research in Electronics, IIT Delhi October 11, 2005 Abstract Adaptive beamforming is minimizing the output power in constrained
More informationThe Staggered SAR Concept: Imaging a Wide Continuous Swath with High Resolution
The Staggered SAR Concept: Imaging a Wide Continuous Swath with High Resolution Michelangelo Villano *, Gerhard Krieger *, Alberto Moreira * * German Aerospace Center (DLR), Microwaves and Radar Institute
More informationEffects of Weight Approximation Methods on Performance of Digital Beamforming Using Least Mean Squares Algorithm
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331,Volume 6, Issue 3 (May. - Jun. 2013), PP 82-90 Effects of Weight Approximation Methods on Performance
More informationPATENT COUNSEL 1176 HOWELL ST. CODE 00OC, BLDG. 11 NEWPORT, RI 02841
DEPARTMENT OF THE NAVY NAVAL UNDERSEA WARFARE CENTER DIVISION NEWPORT OFFICE OF COUNSEL PHONE: (401) 832-3653 FAX: (401) 832-4432 NEWPORT DSN: 432-3653 Attorney Docket No. 83417 Date: 20 June 2007 The
More informationCertificate in Clinician Performed Ultrasound (CCPU)
Certificate in Clinician Performed Ultrasound (CCPU) Syllabus Physics Tutorial Physics Tutorial Purpose: Training: Assessments: This unit is designed to cover the theoretical and practical curriculum for
More informationV. Zetterberg Amlab Elektronik AB Nettovaegen 11, Jaerfaella, Sweden Phone:
Comparison between whitened generalized cross correlation and adaptive filter for time delay estimation with scattered arrays for passive positioning of moving targets in Baltic Sea shallow waters V. Zetterberg
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationFeature Tracking and Optical Flow
Feature Tracking and Optical Flow Prof. D. Stricker Doz. G. Bleser Many slides adapted from James Hays, Derek Hoeim, Lana Lazebnik, Silvio Saverse, who in turn adapted slides from Steve Seitz, Rick Szeliski,
More informationChristian Doppler Laboratory for Dependable Wireless Connectivity for the Society in Motion Three-Dimensional Beamforming
Christian Doppler Laboratory for Three-Dimensional Beamforming Fjolla Ademaj 15.11.216 Studying 3D channel models Channel models on system-level tools commonly 2-dimensional (2D) 3GPP Spatial Channel Model
More informationLine of Sight Stabilization Primer Table of Contents
Line of Sight Stabilization Primer Table of Contents Preface 1 Chapter 1.0 Introduction 3 Chapter 2.0 LOS Control Architecture and Design 11 2.1 Direct LOS Stabilization 15 2.2 Indirect LOS Stabilization
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationSS4 - Spatial array processing
Doctoral Program on Electrical Engineering and Communications Sérgio M. Jesus (sjesus@ualg.pt) Universidade do Algarve, PT-8005-139 Faro, Portugal www.siplab.fct.ualg.pt 24 February 2010 thanks to Paulo
More informationSIGNAL SEPARATION USING CIRCULAR ARRAYS. I. Jouny
SIGNAL SEPARATION USING CIRCULAR ARRAYS I Jouny Department of Electrical Engineering Lafayette College Easton, PA 18042 ABSTRACT Signal separation algorithms can be utilized to improve channel capacity
More informationMultitarget/Multisensor Tracking using only Range and Doppler Measurements
1 Multitarget/Multisensor Tracing using only Range and Doppler Measurements Ross Deming, Member, IEEE, John Schindler, Life Fellow, IEEE, and Leonid Perlovsy, Senior Member, IEEE Abstract A new approach
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Underwater Acoustics Session 2aUW: Wave Propagation in a Random Medium
More informationDo It Yourself 8. Polarization Coherence Tomography (P.C.T) Training Course
Do It Yourself 8 Polarization Coherence Tomography (P.C.T) Training Course 1 Objectives To provide a self taught introduction to Polarization Coherence Tomography (PCT) processing techniques to enable
More informationAdvanced Processing Techniques and Classification of Full-waveform Airborne Laser...
f j y = f( x) = f ( x) n j= 1 j Advanced Processing Techniques and Classification of Full-waveform Airborne Laser... 89 A summary of the proposed methods is presented below: Stilla et al. propose a method
More informationSEG/New Orleans 2006 Annual Meeting
Accuracy improvement for super-wide angle one-way waves by wavefront reconstruction Ru-Shan Wu* and Xiaofeng Jia, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary To
More informationResidual move-out analysis with 3-D angle-domain common-image gathers
Stanford Exploration Project, Report 115, May 22, 2004, pages 191 199 Residual move-out analysis with 3-D angle-domain common-image gathers Thomas Tisserant and Biondo Biondi 1 ABSTRACT We describe a method
More informationMultiple target detection in video using quadratic multi-frame correlation filtering
Multiple target detection in video using quadratic multi-frame correlation filtering Ryan Kerekes Oak Ridge National Laboratory B. V. K. Vijaya Kumar Carnegie Mellon University March 17, 2008 1 Outline
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 information( ) =cov X Y = W PRINCIPAL COMPONENT ANALYSIS. Eigenvectors of the covariance matrix are the principal components
Review Lecture 14 ! PRINCIPAL COMPONENT ANALYSIS Eigenvectors of the covariance matrix are the principal components 1. =cov X Top K principal components are the eigenvectors with K largest eigenvalues
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationAdaptive radar sensing strategies
Adaptive radar sensing strategies A.Hero Univ. of Michigan Ann Arbor 1st Year Review, AFRL, 09/07 AFOSR MURI Integrated fusion, performance prediction, and sensor management for ATE (PI: R. Moses) Outline
More information10. Cartesian Trajectory Planning for Robot Manipulators
V. Kumar 0. Cartesian rajectory Planning for obot Manipulators 0.. Introduction Given a starting end effector position and orientation and a goal position and orientation we want to generate a smooth trajectory
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 informationFusion of Radar and EO-sensors for Surveillance
of Radar and EO-sensors for Surveillance L.J.H.M. Kester, A. Theil TNO Physics and Electronics Laboratory P.O. Box 96864, 2509 JG The Hague, The Netherlands kester@fel.tno.nl, theil@fel.tno.nl Abstract
More informationMulti-azimuth velocity estimation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 87 Multi-azimuth velocity estimation Robert G. Clapp and Biondo Biondi 1 ABSTRACT It is well known that the inverse problem of estimating interval
More informationExperiments with Edge Detection using One-dimensional Surface Fitting
Experiments with Edge Detection using One-dimensional Surface Fitting Gabor Terei, Jorge Luis Nunes e Silva Brito The Ohio State University, Department of Geodetic Science and Surveying 1958 Neil Avenue,
More informationAnisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS
Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS Summary Anisotropic depth model building using surface seismic data alone is non-unique and
More informationMemorandum. Clint Slatton Prof. Brian Evans Term project idea for Multidimensional Signal Processing (EE381k)
Memorandum From: To: Subject: Date : Clint Slatton Prof. Brian Evans Term project idea for Multidimensional Signal Processing (EE381k) 16-Sep-98 Project title: Minimizing segmentation discontinuities in
More informationStructured light 3D reconstruction
Structured light 3D reconstruction Reconstruction pipeline and industrial applications rodola@dsi.unive.it 11/05/2010 3D Reconstruction 3D reconstruction is the process of capturing the shape and appearance
More informationMotion Estimation for Video Coding Standards
Motion Estimation for Video Coding Standards Prof. Ja-Ling Wu Department of Computer Science and Information Engineering National Taiwan University Introduction of Motion Estimation The goal of video compression
More informationGE Healthcare. Agile Ultrasound. The Next Revolution in Ultrasound Imaging
Agile Ultrasound The Next Revolution in Ultrasound Imaging Abstract Diagnostic use of ultrasound has greatly expanded over the past couple of decades because it offers many advantages as an imaging modality.
More informationComparison of Several Spectral Estimation Methods for Application to ISAR Imaging
Methods for Application to ISAR Imaging A. Karakasiliotis and P. Frangos Division of Information Transmission Systems and Materials Technology School of Electrical and Computer Engineering National Technical
More informationMACHINE VISION AS A METHOD FOR CHARACTERIZING SOLAR TRACKER PERFORMANCE
MACHINE VISION AS A METHOD FOR CHARACTERIZING SOLAR TRACKER PERFORMANCE M. Davis, J. Lawler, J. Coyle, A. Reich, T. Williams GreenMountain Engineering, LLC ABSTRACT This paper describes an approach to
More informationPrincipal Component Analysis (PCA) is a most practicable. statistical technique. Its application plays a major role in many
CHAPTER 3 PRINCIPAL COMPONENT ANALYSIS ON EIGENFACES 2D AND 3D MODEL 3.1 INTRODUCTION Principal Component Analysis (PCA) is a most practicable statistical technique. Its application plays a major role
More informationSelection and Integration of Sensors Alex Spitzer 11/23/14
Selection and Integration of Sensors Alex Spitzer aes368@cornell.edu 11/23/14 Sensors Perception of the outside world Cameras, DVL, Sonar, Pressure Accelerometers, Gyroscopes, Magnetometers Position vs
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationSuccessful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment.
Successful application of joint reflection/refraction tomographic velocity inversion in a shallow water marine environment. Sergey Birdus 1, Dean Criddle 2, Alexey Artyomov 1, Li Li 1, Qingbing Tang 1,
More informationSeismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy
Seismic data Interpolation in the Continuous Wavenumber Domain, Flexibility and Accuracy Ye Zheng Geo-X Exploration Services Inc. Summary Our recently developed algorithm, ASFT (Arbitrarily Sampled Fourier
More informationCommon-angle processing using reflection angle computed by kinematic pre-stack time demigration
Common-angle processing using reflection angle computed by kinematic pre-stack time demigration Didier Lecerf*, Philippe Herrmann, Gilles Lambaré, Jean-Paul Tourré and Sylvian Legleut, CGGVeritas Summary
More informationTh N Deghosting by Echo-deblending SUMMARY. A.J. Berkhout* (Delft University of Technology) & G. Blacquiere (Delft University of Technology)
Th N103 11 Deghosting by Echo-deblending A.J. Berkhout* (Delft University of Technology) & G. Blacquiere (Delft University of Technology) SUMMARY Because of the strong sea surface reflectivity, a marine
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 NOISE REDUCTION BY SPATIAL SLOPE ESTIMATION
: steerable and non-linear filter approach Abdulaziz M. Almuhaidib and M. Nafi Toksöz, Earth Resources Laboratory, MIT SUMMARY We present an approach based on local-slope estimation for the separation
More informationChapter 4. Clustering Core Atoms by Location
Chapter 4. Clustering Core Atoms by Location In this chapter, a process for sampling core atoms in space is developed, so that the analytic techniques in section 3C can be applied to local collections
More informationNote Set 4: Finite Mixture Models and the EM Algorithm
Note Set 4: Finite Mixture Models and the EM Algorithm Padhraic Smyth, Department of Computer Science University of California, Irvine Finite Mixture Models A finite mixture model with K components, for
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 informationScale-invariant visual tracking by particle filtering
Scale-invariant visual tracing by particle filtering Arie Nahmani* a, Allen Tannenbaum a,b a Dept. of Electrical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel b Schools of
More informationTu N Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction
Tu N114 05 Internal Multiple Attenuation on Radial Gathers With Inverse-scattering Series Prediction F. Xavier de Melo* (Schlumberger), C. Kostov (Schlumberger) & J. Wu (Schlumberger) SUMMARY We present
More informationInternal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction
Internal Multiple Attenuation on Radial Gathers With Inverse- Scattering Series Prediction Frederico Xavier de Melo, Clément Kostov, James Wu Schlumberger Summary We present a novel workflow for pre-stack
More informationTerrain correction. Backward geocoding. Terrain correction and ortho-rectification. Why geometric terrain correction? Rüdiger Gens
Terrain correction and ortho-rectification Terrain correction Rüdiger Gens Why geometric terrain correction? Backward geocoding remove effects of side looking geometry of SAR images necessary step to allow
More informationRadar Coherent Backscatter!
Radar Coherent Backscatter! Pixels in a radar image are a complex phasor representation of the coherent backscatter from the resolution element on the ground and the propagation phase delay! Interferometric
More informationAdaptive Doppler centroid estimation algorithm of airborne SAR
Adaptive Doppler centroid estimation algorithm of airborne SAR Jian Yang 1,2a), Chang Liu 1, and Yanfei Wang 1 1 Institute of Electronics, Chinese Academy of Sciences 19 North Sihuan Road, Haidian, Beijing
More informationComputer Experiments: Space Filling Design and Gaussian Process Modeling
Computer Experiments: Space Filling Design and Gaussian Process Modeling Best Practice Authored by: Cory Natoli Sarah Burke, Ph.D. 30 March 2018 The goal of the STAT COE is to assist in developing rigorous,
More informationSPACE-TIME adaptive processing (STAP) techniques have
4182 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 58, NO 8, AUGUST 2010 Reduced-Rank STAP Schemes for Airborne Radar Based on Switched Joint Interpolation, Decimation Filtering Algorithm Rui Fa, Rodrigo
More informationNumerical Robustness. The implementation of adaptive filtering algorithms on a digital computer, which inevitably operates using finite word-lengths,
1. Introduction Adaptive filtering techniques are used in a wide range of applications, including echo cancellation, adaptive equalization, adaptive noise cancellation, and adaptive beamforming. These
More informationA Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation
, pp.162-167 http://dx.doi.org/10.14257/astl.2016.138.33 A Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation Liqiang Hu, Chaofeng He Shijiazhuang Tiedao University,
More informationThe interpolation of a 3-D data set by a pair of 2-D filters
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 275 The interpolation of a 3-D data set by a pair of 2-D filters Matthias Schwab and Jon F. Claerbout 1 ABSTRACT Seismic 3-D field data is
More informationModule 4. K-Space Symmetry. Review. K-Space Review. K-Space Symmetry. Partial or Fractional Echo. Half or Partial Fourier HASTE
MRES 7005 - Fast Imaging Techniques Module 4 K-Space Symmetry Review K-Space Review K-Space Symmetry Partial or Fractional Echo Half or Partial Fourier HASTE Conditions for successful reconstruction Interpolation
More information9.8 Application and issues of SZ phase coding for NEXRAD
9.8 Application and issues of SZ phase coding for NEXRAD J.C. Hubbert, G. Meymaris, S. Ellis and M. Dixon National Center for Atmospheric Research, Boulder CO 1. INTRODUCTION SZ phase coding has proved
More informationRobot Mapping. A Short Introduction to the Bayes Filter and Related Models. Gian Diego Tipaldi, Wolfram Burgard
Robot Mapping A Short Introduction to the Bayes Filter and Related Models Gian Diego Tipaldi, Wolfram Burgard 1 State Estimation Estimate the state of a system given observations and controls Goal: 2 Recursive
More information1/12/2009. Image Elements (Pixels) Image Elements (Pixels) Digital Image. Digital Image =...
PAM3012 Digital Image Processing for Radiographers Image Sampling & Quantization In this lecture Definitions of Spatial l & Gray-level l resolution Perceived Image Quality & Resolution Aliasing & Moire
More informationCFP migration; practical aspects
13 CFP migration; practical aspects Jan Thorbecke 1 13.1 Introduction In the past year the Common Focus Point (CFP) technology has become an important paradigm in the DELPHI research program. The interpretation
More informationA Correlation Test: What were the interferometric observation conditions?
A Correlation Test: What were the interferometric observation conditions? Correlation in Practical Systems For Single-Pass Two-Aperture Interferometer Systems System noise and baseline/volumetric decorrelation
More informationParticle Velocimetry Data from COMSOL Model of Micro-channels
Particle Velocimetry Data from COMSOL Model of Micro-channels P.Mahanti *,1, M.Keebaugh 1, N.Weiss 1, P.Jones 1, M.Hayes 1, T.Taylor 1 Arizona State University, Tempe, Arizona *Corresponding author: GWC
More information