Array signal processing for radio-astronomy imaging and future radio telescopes
|
|
- Giles Burke
- 5 years ago
- Views:
Transcription
1 Array signal processing for radio-astronomy imaging and future radio telescopes School of Engineering Bar-Ilan University Amir Leshem Faculty of EEMCS Delft University of Technology Joint wor with Alle-Jan van der Veen, Chen Ben-David, Ronny Levanda, A. Boonstra Some data and simulations were provided by S. van der Toll, S. Wijnholds, Huib Intema, Huub Rottgering and T. Clare Page 1
2 What is parametric imaging? Page 2
3 References C. Ben David and A. Leshem. Parametric high resolution techniques for radio astronomical image formation. IEEE JSTSP October A. Leshem and A-J. van der Veen. Radioastronomical imaging in the presence of strong radio interference. IEEE trans. on Information Theory. Special issue on information theoretic imaging, no. 5, pages August A. Leshem, A-J. van der Veen and A. J. Boonstra. Multichannel interference mitigation techniques in radio-astronomy. The Astrophysical Journal supplements. Volume 131, , November Page 3
4 IEEE Journal of selected topics in Signal Processing October 2008 Page 4
5 Presentation overview The matrix measurement equation Planar array Non co-planar array Polarization Calibration* Spatially varying calibration* Some ideas related to parametric imaging LS-MVI Acceleration and hardware implementation Joint imaging and calibration through SDP L1 optimization* Examples Conclusions Page 5
6 High dynamic range 1 dimensional images for simplicity Dirty images only! Theoretical dynamic range within the image 1,000,000:1 Weaest source is only 10x noise RMSE! Sources are extended off the grid One source is moved through FOV Page 6
7 High dynamic range - Array beam Page 7
8 High dynamic range - sources Page 8
9 High dynamic range - Dirty images Three 1-D dirty images: Classical, AAR and AAR with 3% RMSE calibration error on each array element in each direction Page 9
10 Imaging The correlation process T = 1 p Signals are staced into vectors x( t) x ( t),..., x ( t) Location of antenna i at time t is r ( t) Baseline between antennas i, j is r( t) r ( t) i i j Page 10
11 Imaging - the correlation process 1 R x x x x 10 s/ T s * *, E ( t) ( t ) ω = ω ω = ω( t + nts) ω( t + nts) N n= 0 ( ω ) R corresponds to the visibility V ri( t), rj( t) at frequency, i, j ω 2π Page 11
12 Coordinate system Correlations are measured only for baselines r( t) r ( t) i j Page 12
13 Imaging equation Planar array ( ) ( ) 2 2 jul vm Vf u, v A ( l, m) I f ( l, m) e π + dldm = We limit ourselves to single frequency. Alm (, ) is the "nown" amplitude response I ( l, m) is the intensity at location ( l, m). Vf of the antennas. f ( u v), is the visibility function. Page 13
14 Imaging equation: Non-coplanar array V f u, v, w ( ) 1 l = 1 m 2 2 (,, ) (, ) ( ) 2 2 jul vm wn A l m n I f l m e π + + dldm Page 14
15 The discrete measurement equation T j ( ) V ( ), ( ) ( ) l ij = V i t j t = I l e s ri rj r r s d l= 1 Page 15
16 Matrix formulation of the imaging equation Recall T j ( ) T jsr i R V = I s e = e I s e ( ) ( ) i, j d d s r r T jsr ( ) = 1 = 1 In "our" notation this translates to a where ij R = A BA H i j j I( s1 ) 0 A = a ( s1),..., a ( sd), B= 0 I( s ) d ( s) T jsr1 ( t e ) = T jsr p ( t ) e Page 16
17 The non co-planar array case Page 17
18 Matrix formulation of the imaging equation H H R = Γ( t ) ABAΓ( t ) + R t ( ) nn Page 18 Γ( t ) is a self-calibration slowly time varying pxp matrix Rnn( t ) is the noise covariance matrix. Typically AWGN but when interference exists it dominates
19 Polarimetric imaging When imaging polarized sources each element receives the two polarizations. The steering vector at epoch is now replaced by two orthogonal vectors describing circularly left and right polarization components a (), s a () s, L, R For all s, s we have a ( s ) a ( s ) 1 2, L 1, R 2 For a quasi mono-chromatic source the polarization can be described by the parameters b11 b12 b21 b 22 which are the coherence between right and left circular polarization components. Page 19
20 Polarimetric imaging The sources cross polarization coefficients b, b, b, b are connected to the Stoes parameters I,Q,U,V through the matrices b I b j 0 Q = b j 0 U b V This implies that it will be sufficient to deconvolve the parameters b, b, b, b and then recover the Stoes parameter Page 20
21 Polarimetric imaging The matrix form of the (calibrated) measurement equation now becomes (assuming N point sources) A = a, L( s1 ), a, R( s1 ),..., a, L( sn), a, R( sn) 1 1 b11 b b21 b22 B = N N b11 b12 N N b21 b 22 R = A BA H + σ 2 I, =1,...,K Page 21
22 Self calibration with polarized sources 2 R = Γ ABA Γ + σ I H H, =1,...,K ( ) Γ = diag Γ,..., Γ Γ,1, p i, ( γ ) ( ) 11 γ12 i, i, = ( γ ) ( ) 21 γ i, 22 i, Page 22
23 Deconvolution Page 23
24 Deconvolution in the (u,v) plain Clar proposed to estimate several CLEAN component in a single dirty image and subtract them Cotton-Schwab proposed to perform the subtraction in the (u,v) domain. This allows for non grid points We will always consider the Cotton-Schwab approach Perform non-grid estimation whenever possible We will discuss extension of the Clar approach Page 24
25 Deconvolution Classical Fourier imaging H I( s) = a ( s) R a ( s) D K = 1 Page 25 Super-resolution MVDR based imaging Pseudo-spectrum MVDR criterion ' wˆ ( s) = arg min I ( s) such tha t Solution w D ( ) ' H I( s) = w ( s) R w ( s) D H w ( s) a ( s) = 1 K = 1 1 wˆ ( s) = β R a s where β = K ' 1 I D( s) = a s R a s 1 ( ) ( ) H = 1 1 a s R a s 1 ( ) ( ) H
26 The adapted angular response We would lie to maintain the Gaussian receiver noise isotropic To that end we add to the MVDR equations the requirement w = 1 Still we want a solution of the form w ( ) α 1 s = R a ( s) We obtain that R a () s w = a s R a s 1 H 2 () () and the angular spectrum becomes a () s R a () s a () s R a () s K H 1 '' D s = H 2 = 1 I ( ) Page 26
27 Enhancing the image using independence Since the covariance matrices at different epochs represent indepndent observations it is reasonable to assume that the overall covariance of the data is given by MVDR R = R 1 = 1 ( ) R K T T T = 1 K Choosing a vector w for R yields w w,..., w where w = R a s 1 In this case the MVDR dirty image becomes: I ( s) = K 1 H 1 a ( s) R a ( s) Page 27
28 The AAR dirty image I ( s) AAR = K = 1 K = 1 a s R a s ( ) 1 ( ) H a s R a s ( ) 2 ( ) H Page 28
29 Dirty images (a) Classical (b) MVDR (c) AAR Page 29
30 Bandwidth synthesis Extension to bandwidth synthesis is different than the classical imaging due to non-linear structure Spectral index of sources can be included as extra parameter (not shown below) I ( s) AAR = ω K ω = 1 K = 1 a s R a s ( ) 1 ( ) H, ω, ω, ω a s R a s ( ) 2 ( ) H, ω, ω, ω Page 30
31 Classical beamforming based CLEAN H max w ( s) R w ( s) = 1 ( ) H ( ) σ 2 R = b a s a s + I n n n n= 1 H I( s) = a ( s) R a ( s) D s K N K = 1 sˆ = arg max I ( s) ˆ α = max I ( sˆ) H R = R γαˆ a ( s) a ( s) D s D Page 31
32 MVI CLEAN imaging sˆ = arg max I ( s) ˆ α = max I ( sˆ) AAR H R = R γαˆ a ( s) a ( s) s AAR Page 32
33 Least Squares power estimator to remove bias λ K 2 ( n) ( n) ( ) H ( ) = arg min n n R αa( s ) a ( s ) α =1 F Independently used in the context of multi-resolution CLEAN by Bhatnagar & Cornwell, A&A, 2004 Page 33
34 The LS-MVI algorithm Measure R,..., R Set R (0) 1 Until convergence criterion Source location estimator Power estimator Update K = R γ =0.1 R λ ( n) s = arg max I ( s) ( n) (n) AAR s K 2 ( n) ( n) ( n) H ( n) = arg min R αas ( ) a ( s ) α =1 F R R a( s ) a ( s ) ( n+ 1) ( n) ( n) ( n) H ( n) = γλ Page 34
35 Enforcing the non-negative definite constraint All source powers are non-negative and residual covariance is also non-negative 2 ( n) H ( n) R σ I αa( s ) a ( s ) 0 α 0 Solution is simple: Solve the unconstrained problem and use bisection on α finding sequentially over α that is good for all Page 35
36 Other enhancements Clar type algorithm using SDP Cotton-Schwab with off-grid estimation Page 36
37 Complexity analysis Main bottlenec: Recomputing the dirty image at each stage Complexity of naive implemntation is M 2 where M is the number of pixels in the image, p number of antennas p 2 2 K per iteration We can reduce this to M 2 pk per iteration 2 Still high compared to gridded CLEAN with 2M log 2 M operations for recomputing the dirty image, but managable. These techniques are easy to implement using parallel processors! Dedicated accelerators might be a viable option! Can use coarse grid and use local optimization to obtain accurate location estimate! We do not need intermediate images just accurate parameter estimates! Page 37
38 Robust beamforming and self calibration Assume that a C ( s) is the nominal value. ( s) is not completely nown. is a positive definite matrix describing the uncertainty at epoch a H ( a s a s ) C ( a s a s ) () () () () 1 The robust beamforming problem is given by [ ρ aˆ s aˆ s ], ( ),..., ( ) = arg max 1 subject to K 2 H H ( a a s ) C ( a a s ) 1,..., a () () 1 ρ, a R σ I ρa a 0 ρ 0 K ρ Page 38
39 Equivalent SDP formulation [ ρ aˆ s aˆ s ], ( ),..., ( ) = arg min 1 subject to ρ =1/ τ K τ, a 1,..., a K ( ()) C a a s H ( a a () s ) 1 τ 2 R σ I 0 a H 0 a τ 0 Page 39
40 Robust beamforming based self calibration K N ( ) ( ) arg min n n ˆ ( ) ( ) Γ = 1 n= 1 Γˆ = a s Γa s 2 Page 40
41 Extensions Simplified Maximum lielihood (Leshem and van der Veen 2000) Maximum lielihood using EM (Lanterman) Robust beamforming techniques, e.g., diagonal loading Generalized LS-MVI using parametric source templates (templets: shapelets, curvelets, ridgelets, wavelets ) Global optimization techniques: L1 and others Performance analysis Sensitivity to modeling assumptions Resolution limits Page 41
42 Simulated experiments -East west array with 10 elements and longest baseline 1000λ. 720 covariance matrices - Fully calibrated array Page 42
43 Resolution example Page 43
44 Extended sources example Page 44
45 Extended sources example (a) CLEAN (b) LS_MVI Page 45
46 Conclusions Exciting science can be done with future instruments Science strongly depends on technology Great challenges: Calibration is extremely challenging! Image processing can tae us beyond the equipment!! Real time adaptive interference mitigation is crucial Architectures and algorithms for extremely large arrays Comm 2-4 sensors Radar - thousands of sensors Radio astronomy 10, ,000,000 elements Page 46
arxiv: v1 [astro-ph.im] 16 Apr 2009
Closed form solution of the maximum entropy equations with application to fast radio astronomical image formation arxiv:0904.2545v1 [astro-ph.im] 16 Apr 2009 Amir Leshem 1 School of Engineering, Bar-Ilan
More informationParametric High Resolution Techniques for Radio Astronomical Imaging Chen Ben-David and Amir Leshem, Senior Member, IEEE
670 IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING, VOL. 2, NO. 5, OCTOBER 2008 Parametric High Resolution Techniques for Radio Astronomical Imaging Chen Ben-David and Amir Leshem, Senior Member,
More informationParametric high resolution techniques for radio astronomical imaging
Parametric high resolution techniques for radio astronomical imaging 1 arxiv:0807.4852v1 [astro-ph] 30 Jul 2008 Chen Ben-David and Amir Leshem (Senior Member) School of Engineering Bar-Ilan University
More informationDeconvolving Primary Beam Patterns from SKA Images
SKA memo 103, 14 aug 2008 Deconvolving Primary Beam Patterns from SKA Images Melvyn Wright & Stuartt Corder University of California, Berkeley, & Caltech, Pasadena, CA. ABSTRACT In this memo we present
More informationImaging using GPU. V-K Veligatla, Kapteyn Institute P. Labropoulos, ASTRON and Kapteyn Institute L. Koopmans, Kapteyn Institute
Imaging using GPU V-K Veligatla, Kapteyn Institute P. Labropoulos, ASTRON and Kapteyn Institute L. Koopmans, Kapteyn Institute Introduction What is a GPU? Why another Imager? Large amount of data to be
More informationPrinciples of Interferometry. Hans-Rainer Klöckner IMPRS Black Board Lectures 2014
Principles of Interferometry Hans-Rainer Klöckner IMPRS Black Board Lectures 2014 acknowledgement Mike Garrett lectures James Di Francesco crash course lectures NAASC Lecture 5 calibration image reconstruction
More informationWide-Field Imaging: I
Wide-Field Imaging: I S. Bhatnagar NRAO, Socorro Twelfth Synthesis Imaging Workshop 010 June 8-15 Wide-field imaging What do we mean by wide-field imaging W-Term: D Fourier transform approximation breaks
More informationCompressed Sensing: Extending CLEAN and NNLS
Compressed Sensing: Extending CLEAN and NNLS Ludwig Schwardt SKA South Africa (KAT Project) Calibration & Imaging Workshop Socorro, NM, USA 31 March 2009 Outline 1 Compressed Sensing (CS) Introduction
More informationLow-frequency radio astronomy and wide-field imaging
Low-frequency radio astronomy and wide-field imaging James Miller-Jones (NRAO Charlottesville/Curtin University) ITN 215212: Black Hole Universe Many slides taken from NRAO Synthesis Imaging Workshop (Tracy
More informationA ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT
Progress In Electromagnetics Research Letters, Vol. 16, 53 60, 2010 A ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT Y. P. Liu and Q. Wan School of Electronic Engineering University of Electronic
More informationNew calibration sources for very long baseline interferometry in the 1.4-GHz band
New calibration sources for very long baseline interferometry in the 1.4-GHz band M K Hailemariam 1,2, M F Bietenholz 2, A de Witt 2, R S Booth 1 1 Department of Physics, University of Pretoria, South
More informationImaging, Deconvolution & Image Analysis I. Theory. (IRAM/Obs. de Paris) 7 th IRAM Millimeter Interferometry School Oct. 4 - Oct.
Imaging, Deconvolution & Image Analysis I. Theory Jérôme PETY (IRAM/Obs. de Paris) 7 th IRAM Millimeter Interferometry School Oct. 4 - Oct. 8 2010, Grenoble Scientific Analysis of a mm Interferometer Output
More informationPerformance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna
Performance Analysis for Strong Interference Remove of Fast Moving Target in Linear Array Antenna Kwan Hyeong Lee Dept. Electriacal Electronic & Communicaton, Daejin University, 1007 Ho Guk ro, Pochen,Gyeonggi,
More informationOverview of Technical Approaches
Overview of Technical Approaches Robust receivers Edit / excise / blank - frequency / time domain - might lose/corrupt astronomy signal Cancel / subtract / null - identify / characterize / subtract - frequency
More informationADAPTIVE ANTENNAS. SPATIAL BF
ADAPTIVE ANTENNAS SPATIAL BF 1 1-Spatial reference BF -Spatial reference beamforming may not use of embedded training sequences. Instead, the directions of arrival (DoA) of the impinging waves are used
More information2-D SENSOR POSITION PERTURBATION ANALYSIS: EQUIVALENCE TO AWGN ON ARRAY OUTPUTS. Volkan Cevher, James H. McClellan
2-D SENSOR POSITION PERTURBATION ANALYSIS: EQUIVALENCE TO AWGN ON ARRAY OUTPUTS Volkan Cevher, James H McClellan Georgia Institute of Technology Atlanta, GA 30332-0250 cevher@ieeeorg, jimmcclellan@ecegatechedu
More informationApplications of Robust Optimization in Signal Processing: Beamforming and Power Control Fall 2012
Applications of Robust Optimization in Signal Processing: Beamforg and Power Control Fall 2012 Instructor: Farid Alizadeh Scribe: Shunqiao Sun 12/09/2012 1 Overview In this presentation, we study the applications
More informationDetection and Localization of Tones and Pulses using an Uncalibrated Array
Detection and Localization of Tones and Pulses using an Uncalibrated Array Steven W. Ellingson January 24, 2002 Contents 1 Introduction 2 2 Traditional Method (BF) 2 3 Proposed Method Version 1 (FXE) 3
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationA Crash Course in Radio Astronomy and Interferometry: 2. Aperture Synthesis
A Crash Course in Radio Astronomy and Interferometry: 2. Aperture Synthesis James Di Francesco National Research Council of Canada North American ALMA Regional Center Victoria (thanks to S. Dougherty,
More informationRadio Astronomical Image Formation using Constrained Least Squares and Krylov Subspaces
Radio Astronomical Image Formation using Constrained Least Squares and Krylov Subspaces A. Mouri Sardarabadi, Amir Leshem 2, Alle-Jan van der Veen arxiv:55.348v [astro-ph.im] 6 May 25 Abstract Image formation
More informationQuantifying and mitigating correlated noise between formed beams on the ASKAP Phased Array Feeds
Quantifying and mitigating correlated noise between formed beams on the ASKAP Phased Array Feeds Ian Heywood,PaoloSerra,AidanHotan,DavidMcConnell CSIRO Astronomy & Space Science Abstract Pairs of beams
More informationSingle-User MIMO systems: Introduction, capacity results, and MIMO beamforming
Single-User MIMO systems: Introduction, capacity results, and MIMO beamforming Master Universitario en Ingeniería de Telecomunicación I. Santamaría Universidad de Cantabria Contents Introduction Multiplexing,
More informationImage restoration. An example in astronomy
Image restoration Convex approaches: penalties and constraints An example in astronomy Jean-François Giovannelli Groupe Signal Image Laboratoire de l Intégration du Matériau au Système Univ. Bordeaux CNRS
More informationLecture 9: Indirect Imaging 2. Two-Element Interferometer. Van Cittert-Zernike Theorem. Aperture Synthesis Imaging. Outline
Lecture 9: Indirect Imaging 2 Outline 1 Two-Element Interferometer 2 Van Cittert-Zernike Theorem 3 Aperture Synthesis Imaging Cygnus A at 6 cm Image courtesy of NRAO/AUI Very Large Array (VLA), New Mexico,
More informationRobust Adaptive Beamforming Based on Interference Covariance Matrix Reconstruction and Steering Vector Estimation
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 60, NO. 7, JULY 0 388 Robust Adaptive Beamforming Based on Interference Covariance Matrix Reconstruction and Steering Vector Estimation Yujie Gu and Amir Leshem
More informationRobust Capon Beamforming
Robust Capon Beamforming Yi Jiang Petre Stoica Zhisong Wang Jian Li University of Florida Uppsala University University of Florida University of Florida March 11, 2003 ASAP Workshop 2003 1 Outline Standard
More informationAdaptive beamforming for uniform linear arrays with unknown mutual coupling. IEEE Antennas and Wireless Propagation Letters.
Title Adaptive beamforming for uniform linear arrays with unknown mutual coupling Author(s) Liao, B; Chan, SC Citation IEEE Antennas And Wireless Propagation Letters, 2012, v. 11, p. 464-467 Issued Date
More informationHigh-performance computing and the Square Kilometre Array (SKA) Chris Broekema (ASTRON) Compute platform lead SKA Science Data Processor
High-performance computing and the Square Kilometre Array (SKA) Chris Broekema (ASTRON) Compute platform lead SKA Science Data Processor Radio Astronomy and Computer Science Both very young sciences (1950s)
More informationJ. Liang School of Automation & Information Engineering Xi an University of Technology, China
Progress In Electromagnetics Research C, Vol. 18, 245 255, 211 A NOVEL DIAGONAL LOADING METHOD FOR ROBUST ADAPTIVE BEAMFORMING W. Wang and R. Wu Tianjin Key Lab for Advanced Signal Processing Civil Aviation
More informationTOWARDS A NEW GENERATION OF RADIO ASTRONOMICAL INSTRUMENTS: signal processing for large distributed arrays
Proposal for a special session at ICASSP 2005: TOWARDS A NEW GENERATION OF RADIO ASTRONOMICAL INSTRUMENTS: signal processing for large distributed arrays This session was organized at the request of the
More informationOPTIMAL ADAPTIVE TRANSMIT BEAMFORMING FOR COGNITIVE MIMO SONAR IN A SHALLOW WATER WAVEGUIDE
OPTIMAL ADAPTIVE TRANSMIT BEAMFORMING FOR COGNITIVE MIMO SONAR IN A SHALLOW WATER WAVEGUIDE Nathan Sharaga School of EE Tel-Aviv University Tel-Aviv, Israel natyshr@gmail.com Joseph Tabrikian Dept. of
More informationRadio Interferometry and Aperture Synthesis
Radio Interferometry and Aperture Synthesis Phil gave a detailed picture of homodyne interferometry Have to combine the light beams physically for interference Imposes many stringent conditions on the
More information1 EM algorithm: updating the mixing proportions {π k } ik are the posterior probabilities at the qth iteration of EM.
Université du Sud Toulon - Var Master Informatique Probabilistic Learning and Data Analysis TD: Model-based clustering by Faicel CHAMROUKHI Solution The aim of this practical wor is to show how the Classification
More informationThe data model. Panos Labropoulos. Kapteyn Instituut
The data model Panos Labropoulos " Challenges The actual 21-cm signal is several times weaker than the foregrounds and sits below the telescope noise Foregrounds: galactic and extragalactic (polarization?)
More informationSIGNAL PROCESSING FOR RADIO ASTRONOMICAL ARRAYS
SIGAL PROCESSIG FOR RADIO ASTROOMICAL ARRAYS Alle-Jan van der Veen, Amir Leshem 2 and Albert-Jan Boonstra 3 Delft Univ. of Technology, Dept. Electrical Eng./DIMES, Mekelweg 4, 2628 CD Delft, The etherlands
More informationTECHNICAL REPORT NO. 86 fewer points to average out the noise. The Keck interferometry uses a single snapshot" mode of operation. This presents a furt
CHARA Technical Report No. 86 1 August 2000 Imaging and Fourier Coverage: Mapping with Depleted Arrays P.G. Tuthill and J.D. Monnier 1. INTRODUCTION In the consideration of the design of a sparse-pupil
More informationSpatial Filtering of RF Interference in Radio Astronomy Using a Reference Antenna Array
432 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 64, NO. 2, JANUARY 15, 2016 Spatial Filtering of RF Interference in Radio Astronomy Using a Reference Antenna Array Ahmad Mouri Sardarabadi, Alle-Jan van
More informationMassive MIMO: Signal Structure, Efficient Processing, and Open Problems II
Massive MIMO: Signal Structure, Efficient Processing, and Open Problems II Mahdi Barzegar Communications and Information Theory Group (CommIT) Technische Universität Berlin Heisenberg Communications and
More informationFARADAY ROTATION MEASURE SYNTHESIS OF UGC 10288, NGC 4845, NGC 3044
FARADAY ROTATION MEASURE SYNTHESIS OF UGC 10288, NGC 4845, NGC 3044 PATRICK KAMIENESKI DYLAN PARÉ KENDALL SULLIVAN PI: DANIEL WANG July 21, 2016 Madison, WI OUTLINE Overview of RM Synthesis: Benefits Basic
More informationThe Square Kilometre Array & High speed data recording
The Square Kilometre Array & High speed Astrophysics Group, Cavendish Laboratory, University of Cambridge http://www.mrao.cam.ac.uk/ bn204/ CRISP Annual Meeting FORTH, Crete, Greece March 2012 Outline
More informationThe Australia Telescope. The Australia Telescope National Facility. Why is it a National Facility? Who uses the AT? Ray Norris CSIRO ATNF
The Australia Telescope National Facility The Australia Telescope Ray Norris CSIRO ATNF Why is it a National Facility? Funded by the federal government (through CSIRO) Provides radio-astronomical facilities
More informationArray signal processing for radio astronomy
Array signal processing for radio astronomy Alle-Jan van der Veen 1, Amir Leshem 2 and Albert-Jan Boonstra 3 1 Delft Univ. of Technology, Dept. Electrical Eng./DIMES, Mekelweg 4, 2628 CD Delft, The Netherlands;
More informationStatistical inversion of the LOFAR Epoch of Reionization experiment data model
Statistical inversion of the LOFAR Epoch of Reionization experiment data model ASTRON, Oude Hoogeveensedijk 4, 7991 PD, Dwingeloo, the Netherlands Kapteyn Astronomical Institute, Landleven 12, 9747 AD,
More informationI. Introduction. Index Terms Multiuser MIMO, feedback, precoding, beamforming, codebook, quantization, OFDM, OFDMA.
Zero-Forcing Beamforming Codebook Design for MU- MIMO OFDM Systems Erdem Bala, Member, IEEE, yle Jung-Lin Pan, Member, IEEE, Robert Olesen, Member, IEEE, Donald Grieco, Senior Member, IEEE InterDigital
More informationCLEAN Based on Spatial Source Coherence
CLEN Based on Spatial Source Coherence Pieter Sijtsma Nationaal Lucht- en Ruimtevaartlaboratorium National erospace Laboratory NLR coustic array measurements in closed wind tunnel test section irbus 340
More information4674 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 58, NO. 9, SEPTEMBER 2010
4674 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 58, NO. 9, SEPTEMBER 2010 A Low Complexity Blind Estimator of Narrowband Polynomial Phase Signals Alon Amar, Member, IEEE, Amir Leshem, Senior Member,
More informationASKAP. and phased array feeds in astronomy. David McConnell CASS: ASKAP Commissioning and Early Science 16 November 2017
ASKAP and phased array feeds in astronomy David McConnell CASS: ASKAP Commissioning and Early Science 16 November 2017 Image credit: Alex Cherney / terrastro.com 1 Credits ASKAP Commissioning & Early Science
More informationInterferometry of Solar System Objects
Interferometry of Solar System Objects Bryan Butler National Radio Astronomy Observatory Atacama Large Millimeter/submillimeter Array Expanded Very Large Array Robert C. Byrd Green Bank Telescope Very
More informationNon-Closing Offsets on the VLA. R. C. Walker National Radio Astronomy Observatory Charlottesville VA.
VLA SCIENTIFIC MEMORANDUM NO. 152 Non-Closing Offsets on the VLA R. C. Walker National Radio Astronomy Observatory Charlottesville VA. March 1984 Recent efforts to obtain very high dynamic range in VLA
More informationELEC E7210: Communication Theory. Lecture 10: MIMO systems
ELEC E7210: Communication Theory Lecture 10: MIMO systems Matrix Definitions, Operations, and Properties (1) NxM matrix a rectangular array of elements a A. an 11 1....... a a 1M. NM B D C E ermitian transpose
More information12.4 Known Channel (Water-Filling Solution)
ECEn 665: Antennas and Propagation for Wireless Communications 54 2.4 Known Channel (Water-Filling Solution) The channel scenarios we have looed at above represent special cases for which the capacity
More informationCorrelation functions in optics and quantum optics, 4
Correlation functions in optics and quantum optics, 4 State Key Laboratory of Precision Spectroscopy East China Normal University, Shanghai, China Luis A. Orozco www.jqi.umd.edu The slides of the course
More informationLeast Squares. Ken Kreutz-Delgado (Nuno Vasconcelos) ECE 175A Winter UCSD
Least Squares Ken Kreutz-Delgado (Nuno Vasconcelos) ECE 75A Winter 0 - UCSD (Unweighted) Least Squares Assume linearity in the unnown, deterministic model parameters Scalar, additive noise model: y f (
More informationIntroduction to Interferometry
Introduction to Interferometry Ciro Pappalardo RadioNet has received funding from the European Union s Horizon 2020 research and innovation programme under grant agreement No 730562 Radioastronomy H.Hertz
More informationBeamforming. A brief introduction. Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University
Beamforming A brief introduction Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University March 2008 References Barry D. Van Veen and Kevin Buckley, Beamforming:
More informationRobust Adaptive Beamforming Based on Low-Complexity Shrinkage-Based Mismatch Estimation
1 Robust Adaptive Beamforming Based on Low-Complexity Shrinkage-Based Mismatch Estimation Hang Ruan and Rodrigo C. de Lamare arxiv:1311.2331v1 [cs.it] 11 Nov 213 Abstract In this work, we propose a low-complexity
More informationMerging Single-Dish Data
Lister Staveley-Smith Narrabri Synthesis Workshop 15 May 2003 1 Merging Single-Dish Data!Bibliography!The short-spacing problem!examples Images Power spectrum!fourier (or UV) plane.v. sky (or image) plane
More informationWavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis
Wavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis Jeffrey S. Morris University of Texas, MD Anderson Cancer Center Joint wor with Marina Vannucci, Philip J. Brown,
More informationAn introduction to closure phases
An introduction to closure phases Michelson Summer Workshop Frontiers of Interferometry: Stars, disks, terrestrial planets Pasadena, USA, July 24 th -28 th 2006 C.A.Haniff Astrophysics Group, Department
More informationAn Uplink-Downlink Duality for Cloud Radio Access Network
An Uplin-Downlin Duality for Cloud Radio Access Networ Liang Liu, Prati Patil, and Wei Yu Department of Electrical and Computer Engineering University of Toronto, Toronto, ON, 5S 3G4, Canada Emails: lianguotliu@utorontoca,
More informationDesign Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
CENTER FOR COMPUTER RESEARCH IN MUSIC AND ACOUSTICS DEPARTMENT OF MUSIC, STANFORD UNIVERSITY REPORT NO. STAN-M-4 Design Criteria for the Quadratically Interpolated FFT Method (I): Bias due to Interpolation
More informationSelf-Calibration and Biconvex Compressive Sensing
Self-Calibration and Biconvex Compressive Sensing Shuyang Ling Department of Mathematics, UC Davis July 12, 2017 Shuyang Ling (UC Davis) SIAM Annual Meeting, 2017, Pittsburgh July 12, 2017 1 / 22 Acknowledgements
More informationDevelopment of Radio Astronomy at the Bosscha Observatory
Proceedings of the Conference of the Indonesia Astronomy and Astrophysics, 29-31 October 2009 Premadi et al., Eds. c HAI 2010 Development of Radio Astronomy at the Bosscha Observatory T. Hidayat 1, M.
More informationA New High-Resolution and Stable MV-SVD Algorithm for Coherent Signals Detection
Progress In Electromagnetics Research M, Vol. 35, 163 171, 2014 A New High-Resolution and Stable MV-SVD Algorithm for Coherent Signals Detection Basma Eldosouky, Amr H. Hussein *, and Salah Khamis Abstract
More informationTop Le# side TEST Right side bo.om
Top bo.om e# side TEST Right side Correlation functions in optics and quantum optics, 4 University of Science and Technology of China Hefei, China Luis A. Orozco www.jqi.umd.edu The slides of the course
More informationVirtual Array Processing for Active Radar and Sonar Sensing
SCHARF AND PEZESHKI: VIRTUAL ARRAY PROCESSING FOR ACTIVE SENSING Virtual Array Processing for Active Radar and Sonar Sensing Louis L. Scharf and Ali Pezeshki Abstract In this paper, we describe how an
More informationDATA ASSIMILATION FOR FLOOD FORECASTING
DATA ASSIMILATION FOR FLOOD FORECASTING Arnold Heemin Delft University of Technology 09/16/14 1 Data assimilation is the incorporation of measurement into a numerical model to improve the model results
More informationarxiv:astro-ph/ v1 17 May 2000
MULTICHANNEL INTERFERENCE MITIGATION TECHNIQUES IN RADIO ASTRONOMY Amir Leshem, Alle-Jan van der Veen arxiv:astro-ph/0005359v1 17 May 2000 Information Technology and Systems Delft University of Technology
More informationRobust Space-Time Adaptive Processing Using Projection Statistics
Robust Space-Time Adaptive Processing Using Projection Statistics André P. des Rosiers 1, Gregory N. Schoenig 2, Lamine Mili 3 1: Adaptive Processing Section, Radar Division United States Naval Research
More informationIvan Valtchanov Herschel Science Centre European Space Astronomy Centre (ESAC) ESA. ESAC,20-21 Sep 2007 Ivan Valtchanov, Herschel Science Centre
SPIRE Observing Strategies Ivan Valtchanov Herschel Science Centre European Space Astronomy Centre (ESAC) ESA Outline SPIRE quick overview Observing with SPIRE Astronomical Observation Templates (AOT)
More informationOn the Optimality of Multiuser Zero-Forcing Precoding in MIMO Broadcast Channels
On the Optimality of Multiuser Zero-Forcing Precoding in MIMO Broadcast Channels Saeed Kaviani and Witold A. Krzymień University of Alberta / TRLabs, Edmonton, Alberta, Canada T6G 2V4 E-mail: {saeed,wa}@ece.ualberta.ca
More informationNon-Imaging Data Analysis
Outline 2 Non-Imaging Data Analysis Greg Taylor Based on the original lecture by T.J. Pearson Introduction Inspecting visibility data Model fitting Some applications Superluminal motion Gamma-ray bursts
More informationLecture 7 MIMO Communica2ons
Wireless Communications Lecture 7 MIMO Communica2ons Prof. Chun-Hung Liu Dept. of Electrical and Computer Engineering National Chiao Tung University Fall 2014 1 Outline MIMO Communications (Chapter 10
More informationNoise covariance model for time-series InSAR analysis. Piyush Agram, Mark Simons
Noise covariance model for time-series InSAR analysis Piyush Agram, Mark Simons Motivation Noise covariance in individual interferograms has been well studied. (Hanssen, 2001) Extend the noise model to
More informationOverview of Beamforming
Overview of Beamforming Arye Nehorai Preston M. Green Department of Electrical and Systems Engineering Washington University in St. Louis March 14, 2012 CSSIP Lab 1 Outline Introduction Spatial and temporal
More informationCombining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation
UIUC CSL Mar. 24 Combining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation Yuejie Chi Department of ECE and BMI Ohio State University Joint work with Yuxin Chen (Stanford).
More informationDOA Estimation of Uncorrelated and Coherent Signals in Multipath Environment Using ULA Antennas
DOA Estimation of Uncorrelated and Coherent Signals in Multipath Environment Using ULA Antennas U.Somalatha 1 T.V.S.Gowtham Prasad 2 T. Ravi Kumar Naidu PG Student, Dept. of ECE, SVEC, Tirupati, Andhra
More informationSKA Science Data Processing or MAKING SENSE OF SENSOR BIG DATA
SKA Science Data Processing or MAKING SENSE OF SENSOR BIG DATA Bojan Nikolic & SDP Consortium Astrophysics Group, Cavendish Laboratory Project Engineer/Architecture Lead SKA Science Data Processor What
More informationHIGH RESOLUTION DOA ESTIMATION IN FULLY COHERENT ENVIRONMENTS. S. N. Shahi, M. Emadi, and K. Sadeghi Sharif University of Technology Iran
Progress In Electromagnetics Research C, Vol. 5, 35 48, 28 HIGH RESOLUTION DOA ESTIMATION IN FULLY COHERENT ENVIRONMENTS S. N. Shahi, M. Emadi, and K. Sadeghi Sharif University of Technology Iran Abstract
More informationTHE estimation of covariance matrices is a crucial component
1 A Subspace Method for Array Covariance Matrix Estimation Mostafa Rahmani and George K. Atia, Member, IEEE, arxiv:1411.0622v1 [cs.na] 20 Oct 2014 Abstract This paper introduces a subspace method for the
More informationNear Optimal Adaptive Robust Beamforming
Near Optimal Adaptive Robust Beamforming The performance degradation in traditional adaptive beamformers can be attributed to the imprecise knowledge of the array steering vector and inaccurate estimation
More informationCalibration of Radio Interferometers Using a Sparse DoA Estimation Framework
Calibration of Radio Interferometers Using a Sparse DoA Estimation Framework M. BROSSARD, M. N. EL KORSO, M. ESAVENTO, R. BOYER and. LARZABAL SATIE, UMR 809, École Normale Supérieure de Cachan, Cachan,
More informationLinear and Nonlinear Iterative Multiuser Detection
1 Linear and Nonlinear Iterative Multiuser Detection Alex Grant and Lars Rasmussen University of South Australia October 2011 Outline 1 Introduction 2 System Model 3 Multiuser Detection 4 Interference
More informationNew Coherence and RIP Analysis for Weak. Orthogonal Matching Pursuit
New Coherence and RIP Analysis for Wea 1 Orthogonal Matching Pursuit Mingrui Yang, Member, IEEE, and Fran de Hoog arxiv:1405.3354v1 [cs.it] 14 May 2014 Abstract In this paper we define a new coherence
More informationFuture radio galaxy surveys
Future radio galaxy surveys Phil Bull JPL/Caltech Quick overview Radio telescopes are now becoming sensitive enough to perform surveys of 107 109 galaxies out to high z 2 main types of survey from the
More informationFAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE
Progress In Electromagnetics Research C, Vol. 6, 13 20, 2009 FAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE Y. Wu School of Computer Science and Engineering Wuhan Institute of Technology
More informationPrecision calibration of radio interferometers using redundant baselines
Precision calibration of radio interferometers using redundant baselines The MIT Faculty has made this article openly available Please share how this access benefits you Your story matters Citation As
More informationASKAP Commissioning Update #7 April-July Latest results from the BETA test array
ASKAP Commissioning Update #7 April-July 2014 Welcome to the seventh edition of the ASKAP Commissioning Update. This includes all the latest results from commissioning of the six-antenna Boolardy Engineering
More informationSolar System Objects. Bryan Butler National Radio Astronomy Observatory
Solar System Objects Bryan Butler National Radio Astronomy Observatory Atacama Large Millimeter/submillimeter Array Expanded Very Large Array Robert C. Byrd Green Bank Telescope Very Long Baseline Array
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationSome Synthesis Telescope imaging algorithms to remove nonisoplanatic and other nasty artifacts
ASTRONOMY & ASTROPHYSICS MAY I 1999, PAGE 603 SUPPLEMENT SERIES Astron. Astrophys. Suppl. Ser. 136, 603 614 (1999) Some Synthesis Telescope imaging algorithms to remove nonisoplanatic and other nasty artifacts
More informationJoint Azimuth, Elevation and Time of Arrival Estimation of 3-D Point Sources
ISCCSP 8, Malta, -4 March 8 93 Joint Azimuth, Elevation and Time of Arrival Estimation of 3-D Point Sources Insaf Jaafar Route de Raoued Km 35, 83 El Ghazela, Ariana, Tunisia Email: insafjaafar@infcomrnutn
More informationAdaptive Filters. un [ ] yn [ ] w. yn n wun k. - Adaptive filter (FIR): yn n n w nun k. (1) Identification. Unknown System + (2) Inverse modeling
Adaptive Filters - Statistical digital signal processing: in many problems of interest, the signals exhibit some inherent variability plus additive noise we use probabilistic laws to model the statistical
More informationAdaptive Noise Cancellation
Adaptive Noise Cancellation P. Comon and V. Zarzoso January 5, 2010 1 Introduction In numerous application areas, including biomedical engineering, radar, sonar and digital communications, the goal is
More informationPolynomial Matrix Formulation-Based Capon Beamformer
Alzin, Ahmed and Coutts, Fraser K. and Corr, Jamie and Weiss, Stephan and Proudler, Ian K. and Chambers, Jonathon A. (26) Polynomial matrix formulation-based Capon beamformer. In: th IMA International
More informationImage Processing in Astrophysics
AIM-CEA Saclay, France Image Processing in Astrophysics Sandrine Pires sandrine.pires@cea.fr NDPI 2011 Image Processing : Goals Image processing is used once the image acquisition is done by the telescope
More informationRecent progress with the BYU/NRAO phased array feed
Recent progress with the BYU/NRAO phased array feed Brian D. Jeffs, Karl F. Warnick, Jonathan Landon, Michael Elmer Department of Electrical and Computer Engineering Brigham Young University, Provo, UT
More informationOverview of technical approaches to radio frequency interference mitigation
RADIO SCIENCE, VOL. 40,, doi:10.1029/2004rs003160, 2005 Overview of technical approaches to radio frequency interference mitigation F. H. Briggs Research School of Astronomy and Astrophysics, Australian
More informationDesign of Adaptive Filtering Algorithm for Relative Navigation
Design of Adaptive Filtering Algorithm for Relative Navigation Je Young Lee, Hee Sung Kim, Kwang Ho Choi, Joonhoo Lim, Sung Jin Kang, Sebum Chun, and Hyung Keun Lee Abstract Recently, relative navigation
More information