Image Reconstruction Algorithms for 2D Aperture Synthesis Radiometers
|
|
- Sherilyn Byrd
- 5 years ago
- Views:
Transcription
1 Image Reconstruction Algorithms for 2D Aperture Synthesis Radiometers E. Anterrieu 1 and A. Camps 2 1 Dept. Signal, Image & Instrumentation Laboratoire d Astrophysique de Toulouse-Tarbes Université de Toulouse & CNRS 2 Dept. of Signal Theory and Communications, Universitat Politècnica de Catalunya, Campus Nord, D4-016, Barcelona - IEEC CRAE/UPC, camps@tsc.upc.edu 1
2 Outline of the presentation: 1. Introduction: Aperture Synthesis in Radioastronomy and Earth Observation 2. Fundamentals on Aperture Synthesis Techniques 3. Image Reconstruction Algorithms: Ideal Situation 4. Image Reconstruction Algorithms: Real Situation Basic equation governing an interferometric radiometer Pre-processing of Visibility Samples: Cancellation of Sun and Moon effects Computation of differential visibilities Instrument equation in terms of differential visibilities Solution of Instrument Equation in Terms of Differential Visibilities: direct, iterative 5. Regularization of an ill-posed inverse problem Ill-posed problem Tikhonov / Minimum-norm / Band-limited regularization 6. Projection onto Earth grid and interpolation 7. Apodization and adpative strip-processing 8. Conclusions 2
3 1. Introduction: Aperture Synthesis in Radioastronomy and Earth Observation Aperture synthesis radiometers conceived as a way to increase the angular resolution achievable by real aperture antennas in radioastronomy VLA, New Mexico, Socorro ESTAR MIRAS 3
4 2. Fundamentals on Aperture Synthesis Techniques Channel 1 H 1 (f) H 2 (f) Channel 2 Baseline b 1 (t) b 2 (t) Complex Correlator 1 2 bb * 1 2 T ( ξη, ) = 1 1 * V ( u,v ) = b 1( t ) b2 ( t ) k BB GG 2 B * F ( ξη) ( ξη) ( ξη, ) T n1, F n, T B ΩΩ ξ η 2 ph, rec (u,v): ( uv, ) = ( Δx, Δy) λ 0 Director cosines ( ξ, η) = ( sinθ cosϕ,sinθ sinϕ) = antenna spacing normalized to the wavelength set of (u,v) points depend on the antenna separation and the array geometry set of spatial frequencies where the visibility function V(u,v) is sampled Ideal case: - Identical antenna patterns - Negligible spatial decorrelation - No antenna positioning errors 2D Fourier Transform V ( uv, ) = F T( ξ, η) 4
5 3. Image Reconstruction Algorithms: Ideal Situation (i) Antenna Positions Spatial frequencies (u,v) Periodic extension y [wavelengths] v u x [wavelengths] V ( uv, ) = F T( ξ, η) T ( ξη, ) = 2 F ( ) TB (, ) T n ξη, ξη ph rec Ω 1 ξ η 5
6 3. Image Reconstruction Algorithms: Ideal Situation (ii) In SMOS the alias-free FOV can be enlarged since part of the alias images are the cold sky (including the galaxy!) T B image limited by Earth replicas β η ξ Extension of Alias-Free FOV 6
7 4. Image Reconstruction Algorithms: Real Situation (i) Even in the ideal case: - Antenna spacing > λ/ 3 aliasing - Gibbs phenomenon (ringing) near the sharp transitions (mainly alias borders) In the real case: - Antenna patterns are different - Receivers frequency responses are different ( FWF different) - Antenna positioning errors (u,v,w) real different from (u,v,0) ideal IHFFT cannot be used as image reconstruction method More sophisticated algorithms must be devised But it will be good that the second ones tend to IHFFT in ideal conditions and obviously instrumental errors must be calibrated first! 7
8 4. Image Reconstruction Algorithms: Real Situation (ii) Basic equation governing an interferometric radiometer V 12 ( ξη) δ, Tph,rec * u12ξ + v12η + w12 1 ξ η F 1 ( ξη) ( ξη) np, Fnp 2, r12 1 ξ η f0 1 T = % ΩΩ 1 2 ξ + η 1 fore ( j ( u12 + v12 + w12 1 )) exp 2π ξ η ξ η d ξd η Includes: 1. -T rec (Corbella eqn. term) 2. On-ground measured antenna patterns (amplitude and phase: F np, i ξ, η ) 3. Measured fringe-washing function ( ) from correlations at different time lags r% 12 with correlated noise injected 4. Antenna positioning errors (phase center): (u 12, v 12, w 12 ) 5. Antenna voltage pattern frequency-dependence can be minimized with the following weighted average 6. Extension to the polarimetric case including the co- and cross-polar patterns 7. A second integral must be added with -w 12 to account for antenna backlobes ( ) 8
9 4. Image Reconstruction Algorithms: Real Situation (iii) Pre-processing of Visibility Samples: 1.Computation of auxiliary visibilities to extend the AF-FOV to the periodic repetition of Earth aliases : a) Term corresponding to the physical temperature of the receivers (T ph,rec ): VR ( uv, ) Appears as an offset in the visibilities, except in V(0,0), which is measured by the 3 NIRs. V R ξ η ξ η ( ) ˆq u12 + v12 + w12 n ( ) ( ) Δ 1 -T 1 = Fˆ ξ, η F ξ, η ˆ r exp j2 π( u ξ + v η + w 1 ξ η ) dξdη, % ph,rec p * n ΩΩ ξ + η ξ η f0 This term can also be estimated from the Flat Target Response (FTR) when looking at a homogeneous cold sky at T sky, as: b) V (0,0) = T A : weighted average of the 3 NIR measurements V c) Term corresponding to the sky: V sky (u,v), sky V R -Tph,rec = FTR, T -T sky ph,rec L-band noise map of the sky (cosmic + galactic contributions) Δ ( ξη) 1 T ξ η ξ η Bsky, + + = ˆp ξ η ˆq * ξ η ˆ u12 v12 w12 1 π ξ + η + ξ η ξ η ΩΩ F % n1, Fn 2, 12 r exp j2 ( u12 v12 w12 1 ) d d, ξ + η ξ η 1 2 sky 1 f0 ( ) ( ) ( ) 9
10 4. Image Reconstruction Algorithms: Real Situation (iv) d) Term corresponding to the antenna back lobes, since there are two directions (θ,φ) and (π-θ,φ) that are imaged in the same (ξ,η) point: ( ) ( ) ( ( )) ( ξη) ξ + η ξ η = 1 T, back, T rec ˆ ξ η ˆ* ξ η π ξ + η ξ η ξ η ΩΩ ˆ u % 12 v12 w12 1 V,, ξ + η ξ η xp back Fnp Fnp r e j u v w d d f0 Note: uncertainty in measured antenna patterns from back side is large, comparable to its value. e) Term coming from a constant T B over land (TBD) and a modeled T B over the sea so that differential visibilities are zero-mean: V V land sea Δ { ( ) ( ) ' 0, 0 0 ( ) ( ) ( ) ( ),,,,,, 0 0 A Vsky Vs ea Vback TSun dir VSun dir u,v TSun dirv Sun scatt, TMoon dirv Moon dir, TMoon, dirv Moon, scatt ( 0, 0)} land ( 0,0) T = T - 0, 0,0 00 V, land ( ) ˆq * u12ξ + v12η+ w12 1 ξ η ( )% n ( ) Δ 1 1 = Fˆ ξ, η F ξ, η rˆ exp j2 π( u ξ + v η+ w 1 ξ η ) dξdη, p n ΩΩ ξ + η ξ η f land or iced sea 0 ( ξη) ξ η ξ η ( ) ˆq * + + ( ) ˆ u12 v12 w12 1 % n ( 2 ) Δ 1 T, = Fˆ ξ, η F ξ, η r exp j2 π( u ξ + v η+ w 1 ξ η ) dξdη, Bs, ea p n ΩΩ ξ + η ξ η 1 2 se 1 f a 0 Note: this term will be computed once the Sun (++) and Moon (--) contributions are estimated 10
11 4. Image Reconstruction Algorithms: Real Situation (v) Cancellation of Sun and Moon effects: ˆ B R, -1 - Make a raw T B image using an IHFFT: T = F V (u,v) V ( u v) -The Sun is such a bright source that it completely masks the rest of the image. - Estimate the image of a normalized (T B =1 K) Sun and Moon if visible: ˆ -1 T = F V (u,v) Sun / Moon, dir / ref Sun / Moon, dir / ref - Estimate T B from the Sun as: ˆ T ( ξ,, η, ) T TSun, dir = ˆ T ξ, η B Sun dir Sun dir background ( ) Sun, dir Sun, dir Sun, dir - Compute the contributions to the measured visibilities from the Sun: Number of appearances Most probable value around Sun pixel [Kelvin] V ( u,v ) = T V ( u,v) Sun / Moon, dir / ref Sun / Moon, dir / ref Sun / Moon, dir / ref -T B Moon is much smaller and cannot be easily seen in the raw T B image. T BMoon ~250 K 11
12 Computation of differential visibilities: ( u, v ) back (, v) Vs ky ( u, v) sea (, ) land land (, ) ( ) ( u,v ) T V ( u,v ) T V ( u,v ) ΔV ( u,v) = V (u,v) V V u V u v T V u v R ' T V u,v T V Sun, dir Sun, dir Sun, dir Sun, scatt Moon, dir Moon, dir Moon, dir Moon, scatt, Note: V R (0,0) = 0 and V (0,0) = T A (without the T rec term) Instrument equation in terms of differential visibilities: Δ 1 Δ T B ( ξη, ) ξ + η = ˆp ( ξ η) ˆq * ( ξ η) ˆ u12 v12 V F ( π 12ξ + 12η ) ξ η n Fn r j u v d d ξ η f0 Δ ΩΩ,, % exp 2 ( ), ξ + η Earth Where: T ( ξ, η ) = ˆ ΔT ( ξη, ) + T ( ξ, η ) + T land ( ξ η ) B B sea land, The equation is discretized dec W ( u, v) ΔV ( u, v) = G ΔT ( ξ, η ) Notes: Windowing (W(u,v)) can also be performed in the (ξ,η) domain after image reconstruction Redundant baselines (quasi linear dependent equations) can be weighted-averaged prior to construction of G improves condition number 12
13 If redundant baselines are averaged G-matrix rows correspond to the discretization of the integral equation, sampled at the (ξ,η) points of the (u,v) reciprocal grid (slide #5) Reduces to IHFFT in the ideal case no error amplification (cond. number = 1). Redundant visibilities (same (u,v) point) are averaged to reduce noise the corresponding equations are averaged as well. Solution of Instrument Equation in Terms of Differential Visibilities: Direct: G-matrix pseudo-inverse (used in ESTAR) dec H 1 H ( ξη, ) ( ) (, ) (, ) Δ T = G G G W u v ΔV u v Iterative methods (coded in SEPS can run in standard PC): Extended-CLEAN (UPC) Stabilization + conjugate gradient (UPC): (, ) Δ (, ) = Δ (, ) H H dec G W u v V u v G G T ξ η OK for small arrays Calibration + Sun/Moon cancellation + Differential visibilities + G-matrix formulation + Image reconstruction SMOS L1 Processor 13
14 Simulation examples Original TB FFT retrieved TB G-matrix retrieved TB after visibility decomposition Error: Original TB - FFT retrieved TB Original TB - G-matrix retrieved TB after visibility decomposition [from IGARSS 2006 & TGRS, Jan 2008] 14
15 5. Regularization of an ill-posed inverse problem What is an inverse problem? T G T = V direct problem inverse problem T = G -1 V V What is an ill-posed problem? J. HADAMARD (1902), R. COURANT (1962): A problem satisfying the requirements of existence, uniqueness, and continuity is said to be well-posed. How to regularize an ill-posed inverse problem? J. HADAMARD (1902), R. COURANT (1962): The general principle of regularization is to introduce a priori information in order to compensate for the lack/loss of information in the imaging process. 15
16 Ill-posed problem min V - G T 2 T eigenvalues of G * G G * G T = G * V T r = G + V with G + = (G * G) -1 G * G * G is singular 16
17 Tikhonov regularization min V - G T 2 + μ T 2 T,μ eigenvalues of G * G + μi (G * G + μi) T = G * V T r = G + μ V with G + μ = (G * G + μi) -1 G * μ is the LAGRANGE regularization parameter 17
18 Minimum-norm regularization min T 2 T { V = G T singular values of G Σ i 1 T r = G + 1 T V with Σ T G + = V i U i G = σ i U i V i i 1 σ i (SVD) T r = G + m V with Σ 1 T G + = V i U i i m Σ T m σ G = σ i U i V i i m i m (TSVD) m plays the role of a regularization parameter 18
19 Band-limited regularization min V - G T 2 T { (I PH ) T = 0 singular values of J J * J ^ T H = J * V with J = G U * Z T r = U * Z J + V with J + = (J * J) -1 J * P H plays the role of a regularization parameter 19
20 Examples Observed scenes: 20
21 Examples Retrieved maps: 21
22 6. Projection onto Earth grid and interpolation Retrieved maps are reconstructed at antenna level but they should be projected onto Earth surface for scientific purpose. As a consequence, an interpolation or re-gridding procedure should be performed FOV in antenna frame FOV at ground level 22
23 6. Projection onto Earth grid and interpolation Retrieved maps are reconstructed at antenna level but they should be projected onto Earth surface for scientific purpose. As a consequence, an interpolation or re-gridding procedure should be performed FOV in antenna frame FOV at ground level 23
24 6. Projection onto Earth grid and interpolation Retrieved maps are reconstructed at antenna level but they should be projected onto Earth surface for scientific purpose. As a consequence, an interpolation or re-gridding procedure should be performed FOV in antenna frame FOV at ground level without introducing artifacts! 24
25 6. Projection onto Earth grid and interpolation Resampling of a band-limited image with a bi-linear interpolation is translation dependent, this not the case of the Fourier approach: Fourier interpolation bi-linear interpolation 25
26 6. Projection onto Earth grid and interpolation Resampling of a band-limited image with a bi-linear interpolation is translation dependent, this not the case of the Fourier approach: Fourier interpolation bi-linear interpolation 26
27 7. Apodization and strip-adaptive processing Gibbs oscillations Owing to the finite extent of the experimental aperture and to the associated sharp frequency cut-off, Gibbs oscillations should be expected: d 27
28 7. Apodization and strip-adaptive processing Gibbs oscillations Owing to the finite extent of the experimental aperture and to the associated sharp frequency cut-off, Gibbs oscillations should be expected: λ o 2 3 d 28
29 7. Apodization and strip-adaptive processing Gibbs oscillations Owing to the finite extent of the experimental aperture and to the associated sharp frequency cut-off, Gibbs oscillations should be expected: 29
30 7. Apodization and strip-adaptive processing Gibbs oscillations Owing to the finite extent of the experimental aperture and to the associated sharp frequency cut-off, Gibbs oscillations should be expected: These effects can be attenuated by using an appropriate apodization window. 30
31 7. Apodization and strip-adpative processing Projection onto Earth Geometrical effect: - constant angular resolution at antenna level - variable spatial resolution at ground level 31
32 7. Apodization and strip-adaptive processing Projection onto Earth Strip adaptive processing with multiple apodization windows: - variable angular resolution at antenna level - constant spatial resolution at ground level 32
33 8. Conclusions Image reconstruction algorithm for SMOS presented, including: Foreign sources removal (sky, Sun, Moon) and visibility decomposition (Corbella term, backlobe, land + iced sea & sea) reduces reconstruction errors and bias G-matrix formulation & inversion Projection onto Earth grid Apodization and strip-adaptive processing 33
ORIGIN and REDUCTION of THE RECONSTRUCTION BIAS. Eric ANTERRIEU
6th SMOS WORKSHOP May 15-17, 2006 Technical University of Denmark - Lyngby, DK ORIGIN and REDUCTION of THE RECONSTRUCTION BIAS Eric ANTERRIEU Observatoire Midi-Pyrénées Laboratoire d Astrophysique de Toulouse
More informationSMOS L1 Processor Prototype Phase 3
SMOS L1 Processor Prototype Phase 3 6 th SMOS Workshop Antonio Gutiérrez Peña TUD Denmark 15-17th May 2006 1 SMOS L1 Processor Prototype Phase 3 L1 & ADF Product Definitions TUD Denmark 15-17th May 2006
More informationSSS retrieval from space Comparison study using Aquarius and SMOS data
44 th International Liège Colloquium on Ocean Dynamics 7-11 May 2012 SSS retrieval from space Comparison study using Aquarius and SMOS data Physical Oceanography Department Institute of Marine Sciences
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 informationSea Surface Salinity Retrieval within the ESA Soil Moisture and Ocean Salinity (SMOS) Mission
Sea Surface Salinity Retrieval within the ESA Soil Moisture and Ocean Salinity (SMOS) Mission Roberto Sabia, Adriano Camps, Mercè Vall-llossera, Ramón Villarino, Jorge Miranda, Alessandra Monerris, Miguel
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 informationSMOS L1 Sun BT Validation against on-ground radio-telescope network. Daniele Casella, Raffaele Crapolicchio, Emiliano Capolongo
SMOS L1 Sun BT Validation against on-ground radio-telescope network Daniele Casella, Raffaele Crapolicchio, Emiliano Capolongo Background Focus : SMOS L1b Sun BT Objectives: Evaluate the improvements of
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 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 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 informationAquarius Data Release V2.0 Validation Analysis Gary Lagerloef, Aquarius Principal Investigator H. Kao, ESR And Aquarius Cal/Val Team
Aquarius Data Release V2.0 Validation Analysis Gary Lagerloef, Aquarius Principal Investigator H. Kao, ESR And Aquarius Cal/Val Team Analysis period: Sep 2011-Dec 2012 SMOS-Aquarius Workshop 15-17 April
More informationInterpolation of nonstationary seismic records using a fast non-redundant S-transform
Interpolation of nonstationary seismic records using a fast non-redundant S-transform Mostafa Naghizadeh and Kris A. Innanen CREWES University of Calgary SEG annual meeting San Antonio 22 September 2011
More informationWavefront reconstruction for adaptive optics. Marcos van Dam and Richard Clare W.M. Keck Observatory
Wavefront reconstruction for adaptive optics Marcos van Dam and Richard Clare W.M. Keck Observatory Friendly people We borrowed slides from the following people: Lisa Poyneer Luc Gilles Curt Vogel Corinne
More informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
More informationPHASE UNWRAPPING. Sept. 3, 2007 Lecture D1Lb4 Interferometry: Phase unwrapping Rocca
PHASE UNWRAPPING 1 Phase unwrapping 2 1D Phase unwrapping Problem: given the wrapped phase ψ=w(φ) find the unwrapped one ψ. The wrapping operator: W(φ)=angle(exp(j φ)), gives always solution -π π and is
More informationSMOS L1 Sun BT Valida0on against on-ground radio-telescope network
SMOS L1 Sun BT Valida0on against on-ground radio-telescope network Daniele Casella 1, Raffaele Crapolicchio 1;2, Nicola Compare7 1, Christophe Marqué 3 1 Serco 2 ESA-ESRIN 3 Royal Observatory of Belgium
More informationDeveloping Vicarious Calibration for Microwave Sounding Instruments using Lunar Radiation
CICS Science Meeting, College Park, 2017 Developing Vicarious Calibration for Microwave Sounding Instruments using Lunar Radiation Hu(Tiger) Yang Contributor: Dr. Jun Zhou Nov.08, 2017 huyang@umd.edu Outline
More informationAn Introduction to Radio Astronomy
An Introduction to Radio Astronomy Bernard F. Burke Massachusetts Institute of Technology and Francis Graham-Smith Jodrell Bank, University of Manchester CAMBRIDGE UNIVERSITY PRESS Contents Preface Acknowledgements
More informationOptical interferometry: problems and practice
Outline Optical interferometry: problems and practice Chris Haniff Aims. What is an interferometer? Fundamental differences between optical and radio. Implementation at optical wavelengths. Conclusions.
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 informationRegularization methods for large-scale, ill-posed, linear, discrete, inverse problems
Regularization methods for large-scale, ill-posed, linear, discrete, inverse problems Silvia Gazzola Dipartimento di Matematica - Università di Padova January 10, 2012 Seminario ex-studenti 2 Silvia Gazzola
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
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 information3. Lecture. Fourier Transformation Sampling
3. Lecture Fourier Transformation Sampling Some slides taken from Digital Image Processing: An Algorithmic Introduction using Java, Wilhelm Burger and Mark James Burge Separability ² The 2D DFT can be
More informationAssignment 2 Solutions
EE225E/BIOE265 Spring 2012 Principles of MRI Miki Lustig Assignment 2 Solutions 1. Read Nishimura Ch. 3 2. Non-Uniform Sampling. A student has an assignment to monitor the level of Hetch-Hetchi reservoir
More informationA sparsity-based variational approach for the restoration of SMOS images from L1A data
A sparsity-based variational approach for the restoration of SMOS images from L1A data Javier Preciozzi, Andrès Almansa, Pablo Musé, Sylvain Durand, Ali Khazaal, Bernard Rougé To cite this version: Javier
More informationMoon disc temperature as a function of phase at 1700 MHz
APRAXOS Radio Telescope ETH Zürich, Switzerland Moon disc temperature as a function of phase at 1700 MHz Marina Battaglia Departement Physik, ETHZ, Switzerland, mabattag@student.ethz.ch Created 26.02.2003
More informationArray signal processing for radio-astronomy imaging and future radio telescopes
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
More informationNPP ATMS Instrument On-orbit Performance
NPP ATMS Instrument On-orbit Performance K. Anderson, L. Asai, J. Fuentes, N. George Northrop Grumman Electronic Systems ABSTRACT The first Advanced Technology Microwave Sounder (ATMS) was launched on
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 informationPROOF-OF-CONCEPT DEMONSTRATION OF A MILLIMETRE WAVE IMAGING SOUNDER FOR GEOSTATIONARY EARTH ORBIT
PROOF-OF-CONCEPT DEMONSTRATION OF A MILLIMETRE WAVE IMAGING SOUNDER FOR GEOSTATIONARY EARTH ORBIT Anders Carlström 1, Jacob Christensen 1, Anders Emrich 2, Johan Embretsén 2, Karl-Erik Kempe 2, and Peter
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationRadio sources. P. Charlot Laboratoire d Astrophysique de Bordeaux
Radio sources Laboratoire d Astrophysique de Bordeaux Outline Introduction Continuum and spectral line emission processes The radio sky: galactic and extragalactic History of radioastronomy The first 50
More informationDealing with Noise. Stéphane GUILLOTEAU. Laboratoire d Astrophysique de Bordeaux Observatoire Aquitain des Sciences de l Univers
Dealing with Noise Stéphane GUILLOTEAU Laboratoire d Astrophysique de Bordeaux Observatoire Aquitain des Sciences de l Univers I - Theory & Practice of noise II Low S/N analysis Outline 1. Basic Theory
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 informationORBS and ORCS. Reduction and analysis of SITELLE's data
ORBS and ORCS Reduction and analysis of SITELLE's data omas Martin (Université Laval) Laurent Drissen (Université Laval) Laurie Rousseau-Nepton (Université Laval), Alexandre Alarie (Université Laval),
More informationLecture 9: Speckle Interferometry. Full-Aperture Interferometry. Labeyrie Technique. Knox-Thompson Technique. Bispectrum Technique
Lecture 9: Speckle Interferometry Outline 1 Full-Aperture Interferometry 2 Labeyrie Technique 3 Knox-Thompson Technique 4 Bispectrum Technique 5 Differential Speckle Imaging 6 Phase-Diverse Speckle Imaging
More informationPhase-Referencing and the Atmosphere
Phase-Referencing and the Atmosphere Francoise Delplancke Outline: Basic principle of phase-referencing Atmospheric / astrophysical limitations Phase-referencing requirements: Practical problems: dispersion
More informationAquarius Cal/Val: Salinity Retrieval Algorithm V1.3. Thomas Meissner Frank Wentz
Aquarius Cal/Val: Salinity Retrieval Algorithm V1.3 Thomas Meissner Frank Wentz L2 Algorithm: Overview Aquarius Radiometer Counts Earth + Calibration View Radiometer Calibration Algorithm RFI flagging
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More informationAn Introduction to Radio Astronomy
An Introduction to Radio Astronomy Second edition Bernard F. Burke and Francis Graham-Smith CAMBRIDGE UNIVERSITY PRESS Contents Preface to the second edition page x 1 Introduction 1 1.1 The role of radio
More informationLa vraie faute est celle qu on ne corrige pas Confucius
La vraie faute est celle qu on ne corrige pas Confucius An introduction to optical/ir interferometry Brief summary of main results obtained during the last lecture: V = γ (0,u, v) = 12 S I '(ζ,η)exp{ i2π(
More informationLecture 5. Epipolar Geometry. Professor Silvio Savarese Computational Vision and Geometry Lab. 21-Jan-15. Lecture 5 - Silvio Savarese
Lecture 5 Epipolar Geometry Professor Silvio Savarese Computational Vision and Geometry Lab Silvio Savarese Lecture 5-21-Jan-15 Lecture 5 Epipolar Geometry Why is stereo useful? Epipolar constraints Essential
More informationMTG-IRS processing overview and performances. Dorothee Coppens and Bertrand Theodore
MTG-IRS processing overview and performances Dorothee Coppens and Bertrand Theodore 1 ITSC-XXI conference, Darmstadt, 29 Nov 5 Dec 2017 Outlines 1) Overview of MTG mission and IRS instrument 2) Overview
More informationBINGO simulations and updates on the performance of. the instrument
BINGO simulations and updates on the performance of BINGO telescope the instrument M.-A. Bigot-Sazy BINGO collaboration Paris 21cm Intensity Mapping Workshop June 2014 21cm signal Observed sky Credit:
More informationarxiv: 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 informationAnalyzing Spiral Galaxies Observed in Near-Infrared
Analyzing Spiral Galaxies Observed in Near-Infrared Preben Grosbøl European Southern Observatory Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany Abstract A sample of 54 spiral galaxies was observed
More informationRADIO ASTRONOMY II. QPR No. 83. Academic and Research Staff
II. RADIO ASTRONOMY Academic and Research Staff Prof. A. H. Barrett Prof. L. B. Lenoir Dr. S. H. Zisk Prof. B. F. Burke Prof. D. H. Staelin Patricia P. Crowther Prof. M. Loewenthal E. R. Jensen Graduate
More informationOutline. Microwave Radiometry. Thermal Radiation. Thermal Radiation. Atmospheric Windows. Molecular Radiation Spectra. Dr. Sandra L.
Microwave Radiometry Ch6 Ulaby & Long INEL 6669 Dr. X-Pol Outline l Introduction l Thermal Radiation l Black body radiation Rayleigh-Jeans l Power-Temperature correspondence l Non-Blackbody radiation,
More informationHow do you make an image of an object?
How do you make an image of an object? Use a camera to take a picture! But what if the object is hidden?...or invisible to the human eye?...or too far away to see enough detail? Build instruments that
More informationLong-Term Time Series of Water Vapour Total Columns from GOME, SCIAMACHY and GOME-2
Graphics: ESA Graphics: ESA Graphics: ESA Long-Term Time Series of Water Vapour Total Columns from GOME, SCIAMACHY and GOME-2 S. Noël, S. Mieruch, H. Bovensmann, J. P. Burrows Institute of Environmental
More informationDetermining the Specification of an Aperture Array for Cosmological Surveys
Determining the Specification of an Aperture Array for Cosmological Surveys P. Alexander University of Cambridge, Department of Physics, Cavendish Laboratory, Cambridge, UK. A. J. aulkner Jodrell Bank
More informationI. INTRODUCTION OBSERVING Earth variables, such as surface soil moisture
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 55, NO. 5, MAY 2017 2811 A Sparsity-Based Variational Approach for the Restoration of SMOS Images From L1A Data Javier Preciozzi, Andrés Almansa,
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationSPHERICAL NEAR-FIELD ANTENNA MEASUREMENTS
SPHERICAL NEAR-FIELD ANTENNA MEASUREMENTS Edited by J.E.Hansen Peter Peregrinus Ltd. on behalf of the Institution of Electrical Engineers Contents Contributing authors listed Preface v xiii 1 Introduction
More informationSentinel-3A Product Notice SLSTR Level-2 Sea Surface Temperature
Sentinel-3A Product Notice SLSTR Level-2 Sea Surface Temperature Mission Sensor Product Sentinel-3A SLSTR Level 2 Sea Surface Temperature Product Notice ID EUM/OPS-SEN3/DOC/18/984462 S3A.PN-SLSTR-L2M.003
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 informationWhen taking a picture, what color is a (Lambertian) surface?
Shadow removal When taking a picture, what color is a (Lambertian) surface? What if it s not a cloudy day? Region lit by skylight only Region lit by sunlight and skylight What great things could we do
More informationLinear Operators and Fourier Transform
Linear Operators and Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013
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 informationLab 4 Radial Velocity Determination of Membership in Open Clusters
Lab 4 Radial Velocity Determination of Membership in Open Clusters Sean Lockwood 1, Dipesh Bhattarai 2, Neil Lender 3 December 2, 2007 Abstract We used the Doppler velocity of 29 stars in the open clusters
More informationUnfolding Methods in Particle Physics
Unfolding Methods in Particle Physics Volker Blobel University of Hamburg, Hamburg, Germany 1 Inverse problems Abstract Measured distributions in particle physics are distorted by the finite resolution
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Filtering in the Frequency Domain http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background
More informationNext Generation VLA Memo. 41 Initial Imaging Tests of the Spiral Configuration. C.L. Carilli, A. Erickson March 21, 2018
Next Generation VLA Memo. 41 Initial Imaging Tests of the Spiral Configuration C.L. Carilli, A. Erickson March 21, 2018 Abstract We investigate the imaging performance of the Spiral214 array in the context
More informationSUBJECT: SPIRE Spectrometer Detector Thermal Time Constant Derivation Procedure
SUBJECT: PREPARED BY: Time Constant Derivation Procedure Trevor Fulton DOCUMENT No: SPIRE-BSS-NOT- ISSUE: Draft 0.2 APPROVED BY: Date: Page: 2 of 11 Distribution Name Sarah Leeks Edward Polehampton Trevor
More informationPhysical Model to Describe the PARASOL Radiometric Trending. Definition, Adjustment, and Validation
Physical Model to Describe the PARASOL Radiometric Trending Definition, Adjustment, and Validation Bertrand Fougnie CNES CALCON 13 Summary The PARASOL instrument & mission Evidences a of radiometric drift
More informationTHEORY OF INTERFEROMETRY AND APERTURE SYNTHESIS
THEORY OF INTERFEROMETRY AND APERTURE SYNTHESIS ABSTRACT. The basics of interferometry are covered in these lectures. We first discuss a simple interferometer and show that the output is proportional to
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 informationForegrounds for observations of the high redshift global 21 cm signal
Foregrounds for observations of the high redshift global 21 cm signal Geraint Harker 28/10/2010 Fall Postdoc Symposium 1 The hydrogen 21cm line The hydrogen 21cm (1420MHz) transition is a forbidden transition
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More information2D Imaging in Random Media with Active Sensor Array
2D Imaging in Random Media with Active Sensor Array Paper Presentation: Imaging and Time Reversal in Random Media Xianyi Zeng rhodusz@stanford.edu icme, Stanford University Math 221 Presentation May 18
More informationAstronomical Experiments for the Chang E-2 Project
Astronomical Experiments for the Chang E-2 Project Maohai Huang 1, Xiaojun Jiang 1, and Yihua Yan 1 1 National Astronomical Observatories, Chinese Academy of Sciences, 20A Datun Road,Chaoyang District,
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More informationTime and Spatial Series and Transforms
Time and Spatial Series and Transforms Z- and Fourier transforms Gibbs' phenomenon Transforms and linear algebra Wavelet transforms Reading: Sheriff and Geldart, Chapter 15 Z-Transform Consider a digitized
More informationAy 20 Basic Astronomy and the Galaxy Problem Set 2
Ay 20 Basic Astronomy and the Galaxy Problem Set 2 October 19, 2008 1 Angular resolutions of radio and other telescopes Angular resolution for a circular aperture is given by the formula, θ min = 1.22λ
More informationSource Coding for Compression
Source Coding for Compression Types of data compression: 1. Lossless -. Lossy removes redundancies (reversible) removes less important information (irreversible) Lec 16b.6-1 M1 Lossless Entropy Coding,
More informationBiomedical Engineering Image Formation II
Biomedical Engineering Image Formation II PD Dr. Frank G. Zöllner Computer Assisted Clinical Medicine Medical Faculty Mannheim Fourier Series - A Fourier series decomposes periodic functions or periodic
More informationLinear Inverse Problems
Linear Inverse Problems Ajinkya Kadu Utrecht University, The Netherlands February 26, 2018 Outline Introduction Least-squares Reconstruction Methods Examples Summary Introduction 2 What are inverse problems?
More informationMultiscale Image Transforms
Multiscale Image Transforms Goal: Develop filter-based representations to decompose images into component parts, to extract features/structures of interest, and to attenuate noise. Motivation: extract
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More informationLow-rank Promoting Transformations and Tensor Interpolation - Applications to Seismic Data Denoising
Low-rank Promoting Transformations and Tensor Interpolation - Applications to Seismic Data Denoising Curt Da Silva and Felix J. Herrmann 2 Dept. of Mathematics 2 Dept. of Earth and Ocean Sciences, University
More informationIR spectro-interferometry of cool evolved stars
IR spectro-interferometry of cool evolved stars Keiichi Ohnaka Max-Planck-Institut für Radioastronomie Infrared Interferometry Group 1) Mid-IR interferometry of Mira stars 2) Near-IR spectro-interferometry
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 informationWeak Constraints 4D-Var
Weak Constraints 4D-Var Yannick Trémolet ECMWF Training Course - Data Assimilation May 1, 2012 Yannick Trémolet Weak Constraints 4D-Var May 1, 2012 1 / 30 Outline 1 Introduction 2 The Maximum Likelihood
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
More informationSensitivity. Bob Zavala US Naval Observatory. Outline
Sensitivity Bob Zavala US Naval Observatory Tenth Synthesis Imaging Summer School University of New Mexico, June 13-20, 2006 Outline 2 What is Sensitivity? Antenna Performance Measures Interferometer Sensitivity
More informationPush-pull PDV analysis
Push-pull PDV analysis Sub-fringe data reduction PDV workshop August 6-7, 27 Daniel H. Dolan and Scott Jones Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company,
More informationThe IRS Flats. Spitzer Science Center
Spitzer Science Center Table of Contents The IRS Flats 1 Chapter 1. The Purpose of this Document... 3 Chapter 2.... 4 2.1 Make a finely-spaced map of a star....4 2.2 Simulate an extended source...6 2.3
More informationDevelopment of a solar imaging array of Very Small Radio Telescopes
Development of a solar imaging array of Very Small Radio Telescopes Ted Tsiligaridis University of Washington, Seattle Mentor: Alan E.E. Rogers MIT Haystack Observatory Summer 27 Outline 1. Solar Physics
More informationInterferometric array design: Distributions of Fourier samples for imaging
A&A 386, 1160 1171 (2002) DOI: 10.1051/0004-6361:20020297 c ESO 2002 Astronomy & Astrophysics Interferometric array design: Distributions of Fourier samples for imaging F. Boone LERMA, Observatoire de
More informationALMA memo 515 Calculation of integration times for WVR
ALMA memo 515 Calculation of integration times for WVR Alison Stirling, Mark Holdaway, Richard Hills, John Richer March, 005 1 Abstract In this memo we address the issue of how to apply water vapour radiometer
More informationRadio Interferometry Fundamentals. John Conway Onsala Space Obs and Nordic ALMA ARC-node
Radio Interferometry Fundamentals John Conway Onsala Space Obs and Nordic ALMA ARC-node So far discussed only single dish radio/mm obs Resolution λ/d, for D=20m, is 30 at mm-wavelengths and 30 (diameter
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 informationAtmospheric phase correction for ALMA with water-vapour radiometers
Atmospheric phase correction for ALMA with water-vapour radiometers B. Nikolic Cavendish Laboratory, University of Cambridge January 29 NA URSI, Boulder, CO B. Nikolic (University of Cambridge) WVR phase
More informationCOMPUTATION OF FOURIER TRANSFORMS FOR NOISY B FOR NOISY BANDLIMITED SIGNALS
COMPUTATION OF FOURIER TRANSFORMS FOR NOISY BANDLIMITED SIGNALS October 22, 2011 I. Introduction Definition of Fourier transform: F [f ](ω) := ˆf (ω) := + f (t)e iωt dt, ω R (1) I. Introduction Definition
More informationPolarization in Interferometry
Polarization in Interferometry Rick Perley (NRAO-Socorro) Fourth INPE Advanced Course on Astrophysics: Radio Astronomy in the 21 st Century My Apologies, Up Front This lecture contains some difficult material.
More informationIn-Flight Retrieval of Reflector Anomalies for the Planck Space Telescope
Presented at The 4 th European Conference on Antennas and Propagation (EuCAP 2010), Barcelona, Spain, 12-16 April 2010. In-Flight Retrieval of Reflector Anomalies for the Planck Space Telescope Frank Jensen*,
More informationE-MERLIN and EVN/e-VLBI Capabilities, Issues & Requirements
E-MERLIN and EVN/e-VLBI Capabilities, Issues & Requirements e-merlin: capabilities, expectations, issues EVN/e-VLBI: capabilities, development Requirements Achieving sensitivity Dealing with bandwidth,
More informationPreparing for SMOS: Sea Salinity Campaigns and Results
Downloaded from orbit.dtu.dk on: Dec 1, 218 Preparing for SMOS: Sea Salinity Campaigns and Results Søbjærg, Sten Schmidl; Balling, Jan E.; Kristensen, Steen Savstrup; Skou, Niels Published in: Proceedings
More information