A&A 517, A12 (2010)DOI: 10.1051/0004-6361/200912873c○ ESO 2010Astronomy&AstrophysicsRevisiting the theory of interferometric wide-field synthesisJ. Pety 1,2 and N. Rodríguez-Fernández 11 IRAM, 300 rue de la Piscine, 38406 Grenob<strong>le</strong> Cedex, Francee-mail: [pety;rodriguez@iram.fr]2 LERMA, UMR 8112, CNRS and Observatoire de Paris, 61 avenue de l’Observatoire, 75014 Paris, FranceReceived 13 July 2009 / Accepted 16 March 2010ABSTRACT<strong>tel</strong>-<strong>00726959</strong>, version 1 - 31 Aug 2012Context. After several generations of interferometers in radioastronomy, wide-field imaging at high angular resolution is today amajor goal for trying to match optical wide-field performances.Aims. All the radio-interferometric, wide-field imaging methods currently belong to the mosaicking family. Based on a 30 years old,original idea from Ekers & Rots, we aim at proposing an alternate formalism.Methods. Starting from their ideal case, we successively evaluate the impact of the standard ingredients of interferometric imaging,i.e. the sampling function, the visibility gridding, the data weighting, and the processing of the short spacings either from sing<strong>le</strong>-dishantennas or from heterogeneous arrays. After a comparison with standard nonlinear mosaicking, we assess the compatibility of theproposed processing with 1) a method of dealing with the effect of ce<strong>le</strong>stial projection and 2) the elongation of the primary beamalong the scanning direction when using the on-the-fly observing mode.Results. The dirty image resulting from the proposed scheme can be expressed as a convolution of the sky brightness distributionwith a set of wide-field dirty beams varying with the sky coordinates. The wide-field dirty beams are locally shift-invariant as they donot depend strongly on position on the sky: their shapes vary on angular sca<strong>le</strong>s typically larger or equal to the primary beamwidth.A comparison with standard nonlinear mosaicking shows that both processing schemes are not mathematically equiva<strong>le</strong>nt, thoughthey both recover the sky brightness. In particular, the weighting scheme is very different in both methods. Moreover, the proposedscheme naturally processes the short spacings from both sing<strong>le</strong>-dish antennas and heterogeneous arrays. Finally, the sky gridding ofthe measured visibilities, required by the proposed scheme, may potentially save large amounts of hard-disk space and cpu processingpower over mosaicking when handling data sets acquired with the on-the-fly observing mode.Conclusions. We propose to call this promising family of imaging methods wide-field synthesis because it explicitly synthesizesvisibilities at a much finer spatial frequency resolution than the one set by the diameter of the interferometer antennas.Key words. methods: analytical – techniques: interferometric – methods: data analysis – techniques: image processing1. IntroductionThe instantaneous field of view of an interferometer is naturallylimited by the primary beam size of the individual antennas. Forthe ALMA 12 m-antennas, this field of view is ∼9 ′′ at 690 GHzand ∼27 ′′ at 230 GHz. The astrophysical sources in the (sub)-millimeter domain are often much larger than this, but still structuredon much smal<strong>le</strong>r angular sca<strong>le</strong>s. Interferometric wide-fieldtechniques enab<strong>le</strong> us to fully image these sources at high angularresolution. These techniques first require an observing modethat in one way or another scans the sky on spatial sca<strong>le</strong>s largerthan the primary beam. The most common observing mode inuse today, known as stop-and-go mosaicking, consists in repeatedlyobserving sky positions typically separated by half the primarybeam size. The improvement of the tracking behavior ofmodern antennas now <strong>le</strong>ads astronomers to consider on-the-flyobservations, with the antennas s<strong>le</strong>wing continuously across thesky. The improvements in correlator and receiver technologiesare also <strong>le</strong>ading to techniques that could potentially samp<strong>le</strong> theantenna focal planes with multi-beam receivers instead of thesing<strong>le</strong>-pixel receivers instal<strong>le</strong>d on current interferometers.The ideal measurement equation of interferometric widefieldimaging isV ( u p ,α s)=∫α pB ( α p − α s)I(αp)e−i2πα p u pdα p , (1)where V is the visibility function of 1) u p (the spatial frequencywith respect to the fixed phase center) and 2) α s (the scanned skyang<strong>le</strong>), I is the sky brightness, and B the antenna power patternor primary beam of an antenna of the interferometer (Thompsonet al. 1986, Chap. 2). For simplicity, 1) we assume that the primarybeam is independent of azimuth and e<strong>le</strong>vation, and 2) weuse one-dimensional notation without loss of generality. We donot deal with polarimetry (see e.g. Hamaker et al. 1996; Saultet al. 1996a, 1999) because it adds another <strong>le</strong>vel of comp<strong>le</strong>xityover our first goal here, i.e. wide-field considerations. Severalaspects make Eq. (1) peculiar with respect to the ideal measurementequation for sing<strong>le</strong>-field observations. First, the visibility isa function not only of the uv spatial frequency (u p ) but also of thescanned sky coordinate (α s ). Second, Eq. (1) is a mix between aFourier transform and a convolution equation. It can be regarded,for examp<strong>le</strong>, as the Fourier transform along the α p dimensionof the function, B(α p − α s ) I(α p ), of the (α p , α s ) variab<strong>le</strong>s. ButEq. (1) can also be written as the convolution:V ( u p ,α s)=∫α pB ( α s − α p)I(αp , u p)dαp , (2)whereB ( ) ( )α s − α p ≡ B αp − α s (3)andI ( ) ( )α p , u p ≡ I αp e−i2πα p u p. (4)Artic<strong>le</strong> published by EDP Sciences Page 1 of 21
A&A 517, A12 (2010)<strong>tel</strong>-<strong>00726959</strong>, version 1 - 31 Aug 2012For each u p kept constant, V(u p , α s ) is the convolution of Band I. Indeed, I(α p , u p = 0) = I(α p ), so we deriveV ( u p = 0,α s)=∫α pB(α s − α p ) I(α p )dα p , (5)i.e., the convolution equation for sing<strong>le</strong>-dish observations.Ekers & Rots (1979) were the first to show that the measurementequation (Eq. (1)) enab<strong>le</strong>s us to recover spatial frequenciesof the sky brightness at a much finer uv resolution than theuv resolution set by the diameter of the interferometer antennas.Interestingly enough, the goal of Ekers & Rots (1979) was “just”to find a way to produce the missing short spacings of a multiplyinginterferometer. However, Cornwell (1988) realized thatEkers & Rots’ scheme has a much stronger impact, because itexplains why an interferometer is ab<strong>le</strong> to do wide-field imaging.Cornwell (1988) also demonstrated that on-the-fly scanning isnot absolu<strong>tel</strong>y necessary to interferometric wide-field imaging.Indeed, the large-sca<strong>le</strong> information can be retrieved in mosaicsof sing<strong>le</strong>-field observations, provided that the sampling of thesing<strong>le</strong> fields follows the sky-plane Nyquist sampling theorem.As a result, all the information about the sky brightnessis coded in the visibility function. From a data-processingviewpoint, all the current radio-interferometric wide-field imagingmethods (see, e.g., Gueth et al. 1995; Sault et al. 1996b;Cornwell et al. 1993; Bhatnagar & Cornwell 2004; Bhatnagaret al. 2008; Cotton & Uson 2008) belong to the mosaicking family1 pioneered by Cornwell (1988). In this family, the processingstarts with Fourier transforming V(u p , α s ) along the u p dimension(i.e. at constant α s ) to produce a set of sing<strong>le</strong>-field dirty imagesbefore linearly combining them and forming a wide-fielddirty image. In this paper, we propose an alternate processing,which starts with a Fourier transform of V(u p , α s ) along the α sdimension (i.e. at constant u p ). We show how this explici<strong>tel</strong>ysynthesizes the spatial frequencies needed to do wide-field imaging,which are linearly combined to form a “wide-field uv plane”,i.e., one uv-plane containing all the spatial frequency informationmeasured during the wide-field observation. Inverse Fouriertransform will produce a dirty image, which can then be deconvolvedusing standard methods. The existence of two differentways to extract the wide-field information from the visibilityfunction raises several questions: are they equiva<strong>le</strong>nt? What aretheir relative merits?We thus aim at revisiting the mathematical foundations ofwide-field imaging and deconvolution. Sections 2 to 7 proposethe new algorithm, which we call wide-field synthesis: Sect. 2first defines the notations and it then lays out the basic conceptsused throughout the paper. Section 3 states the steps needed togo beyond the Ekers & Rots scheme and explores the consequencesof incomp<strong>le</strong>te sampling of both the uv and sky planes.Section 4 discusses the effects of gridding by convolution andregular resampling. Section 5 describes how to influence thedirty beam shapes and thus the deconvolution. Section 6 stateshow to introduce short spacings measured either from a sing<strong>le</strong>dishantenna or from heterogeneous interferometers. Section 7compares the proposed wide-field synthesis algorithm with standardnonlinear mosaicking. Some detai<strong>le</strong>d demonstrations arefactored out in Appendix A to enab<strong>le</strong> an easier reading of themain paper, whi<strong>le</strong> ensuring that interested readers can followthe demonstrations. Appendices B and C then explain how thewide-field synthesis algorithm can cope with non-ideal effects:Appendix B discusses how at <strong>le</strong>ast one standard way to cope1 In the rest of this paper, stop-and-go mosaicking refer to the observingmode, whi<strong>le</strong> mosaicking alone refer to the imaging family.Page 2 of 21with sky projection prob<strong>le</strong>ms is compatib<strong>le</strong> with the wide-fieldsynthesis algorithm. Appendix C explores the consequences ofusing the on-the-fly observing mode. Finally, we assume goodfamiliarity with sing<strong>le</strong>-field imaging in various places. We referthe reader to well-known references: e.g. Chap. 6 of Thompsonet al. (1986)andSramek & Schwab (1989).2. Notations and basic concepts2.1. NotationsIn this paper, we use the Bracewell (2000)’s notation to displaythe relationship between a function I(α) and its direct Fouriertransform I(u), i.e.,I (α) ⊃ I (u) , (6)where (α, u) is the coup<strong>le</strong> of Fourier conjugate variab<strong>le</strong>s. Wealso use the following sign conventions for the direct and inverseFourier transforms∫I (u) ≡ I (α) e −i2παu dα (7)αand∫I (α) ≡ I (u) e +i2πuα du. (8)uAs V is a function of two independent quantities (u p and α s ),the Fourier transform may be applied independently on each dimension,whi<strong>le</strong> the other dimension stays constant. Several additionalconventions are used to express this. First, we introducea specific notation to state that either the first or the second dimensionstays constant:V up (α s ) ≡ V ( u p = const., α s), (9)and( ) (up ≡ V up ,α s = const. ) . (10)V α sSecond, we use a bottom/top line to derive the notation of theFourier transform along the first/second dimension from the notationof the original function. Third, on the Fourier transformsign (i.e. ⊃), we explicitly state the dimension along which theFourier transform is computed. For instance, if D is a function of(α p , α s ), then the Fourier transform of D along the first dimensionis expressed asD ( ) α pα p ,α s ⊃ D ( )uu p ,α s , (11)pwhi<strong>le</strong> the Fourier transform of D along the second dimension isexpressed asD ( α p ,α s) α s⊃usD ( α p , u s). (12)We also use a more compact notation when doing the Fouriertransform on both dimensions simultaneously, i.e.,D ( ) (α p,α s)α p ,α s⊃ D ( )u(u p ,u s) p , u s . (13)Finally, the convolution of two functions G and V is noted anddefined as∫{G ⋆ V}(u) ≡ G(u − v) V(v)dv. (14)v
- Page 1 and 2:
UNIVERSITÉ PIERRE ET MARIE CURIEHA
- Page 3 and 4:
tel-00726959, version 1 - 31 Aug 20
- Page 5 and 6:
Table des matières1 Rapport de sou
- Page 7 and 8:
Rapport après soutenanceHabilitati
- Page 9 and 10:
Chapitre 2Curriculum vitaetel-00726
- Page 11 and 12:
2.7 ANIMATION ET DIFFUSION DE LA CU
- Page 13 and 14:
2.8 PARCOURS 131992-1993 ÉCOLE NOR
- Page 15 and 16:
Chapitre 3Copyright: IRAM/PdBIIntro
- Page 17 and 18:
tel-00726959, version 1 - 31 Aug 20
- Page 19 and 20:
4.2 ETUDES DIRECTES EN ÉMISSION 19
- Page 21 and 22:
4.4 LA LUMINOSITY CO PAR MOLÉCULE
- Page 23 and 24:
356 E. Falgarone et al.: Extreme ve
- Page 25 and 26:
358 E. Falgarone et al.: Extreme ve
- Page 27 and 28:
360 E. Falgarone et al.: Extreme ve
- Page 29 and 30:
362 E. Falgarone et al.: Extreme ve
- Page 31 and 32:
364 E. Falgarone et al.: Extreme ve
- Page 33 and 34:
366 E. Falgarone et al.: Extreme ve
- Page 35 and 36:
368 E. Falgarone et al.: Extreme ve
- Page 37 and 38:
A&A 518, A45 (2010)1001010⌠⌡ τ
- Page 39 and 40:
A&A 518, A45 (2010)If X HCO + is as
- Page 41 and 42:
A&A 518, A45 (2010)tel-00726959, ve
- Page 43 and 44:
A&A 518, A45 (2010)tel-00726959, ve
- Page 45 and 46:
A&A 518, A45 (2010)Table E.2. Data
- Page 47 and 48:
A&A 541, A58 (2012)tel-00726959, ve
- Page 49 and 50:
A&A 541, A58 (2012)tel-00726959, ve
- Page 51 and 52:
A&A 541, A58 (2012)0.08 0.10.20.12
- Page 53 and 54:
A&A 541, A58 (2012)b=240.11 0.1 0.1
- Page 55 and 56:
0.080.09 0.1A&A 541, A58 (2012)0.14
- Page 57 and 58:
A&A 541, A58 (2012)43˚00'0.3 0.44
- Page 59 and 60:
A&A 541, A58 (2012)68˚00'1.1 0.91.
- Page 61 and 62:
A&A 541, A58 (2012)α(J2000)-19H57m
- Page 63 and 64:
A&A 541, A58 (2012)tel-00726959, ve
- Page 65 and 66:
A&A 541, A58 (2012)Jenkins, E. B.,
- Page 67 and 68:
A&A 541, A58 (2012)Galactic Latitud
- Page 69 and 70:
tel-00726959, version 1 - 31 Aug 20
- Page 71 and 72:
5.2 UNE PHYSIQUE BIEN CONTRAINTE ET
- Page 73 and 74:
5.3 PERSPECTIVES : DES RELEVÉS DE
- Page 75 and 76:
A&A 435, 885-899 (2005)DOI: 10.1051
- Page 77 and 78:
tel-00726959, version 1 - 31 Aug 20
- Page 79 and 80:
J. Pety et al.: Are PAHs precursors
- Page 81 and 82:
J. Pety et al.: Are PAHs precursors
- Page 83 and 84:
J. Pety et al.: Are PAHs precursors
- Page 85 and 86:
J. Pety et al.: Are PAHs precursors
- Page 87 and 88:
J. Pety et al.: Are PAHs precursors
- Page 89 and 90:
J. Pety et al.: Are PAHs precursors
- Page 91 and 92:
566 J. R. Goicoechea et al.: Low su
- Page 93 and 94:
568 J. R. Goicoechea et al.: Low su
- Page 95 and 96:
570 J. R. Goicoechea et al.: Low su
- Page 97 and 98:
572 J. R. Goicoechea et al.: Low su
- Page 99 and 100:
574 J. R. Goicoechea et al.: Low su
- Page 101 and 102:
576 J. R. Goicoechea et al.: Low su
- Page 103 and 104:
578 J. R. Goicoechea et al.: Low su
- Page 105 and 106:
580 J. R. Goicoechea et al.: Low su
- Page 107 and 108:
tel-00726959, version 1 - 31 Aug 20
- Page 109 and 110:
J. R. Goicoechea et al.: Low sulfur
- Page 111 and 112:
L42J. Pety et al.: Deuterium fracti
- Page 113 and 114:
L44J. Pety et al.: Deuterium fracti
- Page 115 and 116:
772 J. R. Goicoechea et al.: The io
- Page 117 and 118:
774 J. R. Goicoechea et al.: The io
- Page 119 and 120:
776 J. R. Goicoechea et al.: The io
- Page 121 and 122:
778 J. R. Goicoechea et al.: The io
- Page 123 and 124:
780 J. R. Goicoechea et al.: The io
- Page 125 and 126:
782 J. R. Goicoechea et al.: The io
- Page 127 and 128:
A&A 494, 977-985 (2009)DOI: 10.1051
- Page 129 and 130:
M. Gerin et al.: HCO mapping of the
- Page 131 and 132:
M. Gerin et al.: HCO mapping of the
- Page 133 and 134:
M. Gerin et al.: HCO mapping of the
- Page 135 and 136:
M. Gerin et al.: HCO mapping of the
- Page 137 and 138:
Table 1. Observation parameters for
- Page 139 and 140:
A&A 534, A49 (2011)tel-00726959, ve
- Page 141 and 142:
A&A 534, A49 (2011)tel-00726959, ve
- Page 143 and 144:
A&A 534, A49 (2011)tel-00726959, ve
- Page 145 and 146:
tel-00726959, version 1 - 31 Aug 20
- Page 147 and 148: tel-00726959, version 1 - 31 Aug 20
- Page 149 and 150: tel-00726959, version 1 - 31 Aug 20
- Page 151 and 152: tel-00726959, version 1 - 31 Aug 20
- Page 153 and 154: tel-00726959, version 1 - 31 Aug 20
- Page 155 and 156: tel-00726959, version 1 - 31 Aug 20
- Page 157 and 158: ¡¢£¤¥¢¦¢ §¨ ©
- Page 159 and 160: tel-00726959, version 1 - 31 Aug 20
- Page 161 and 162: 7.2 MODE INTERFÉROMÉTRIQUE 161tel
- Page 163 and 164: 7.3 LE FUTUR DE LA RADIO-ASTRONOMIE
- Page 165 and 166: CLASS evolution: I. Improved OTF su
- Page 167 and 168: CLASS evolution: I. Improved OTF su
- Page 169 and 170: A&A 526, A47 (2011)DOI: 10.1051/000
- Page 171 and 172: S. Maret et al.: Weeds: a CLASS ext
- Page 173 and 174: S. Maret et al.: Weeds: a CLASS ext
- Page 175 and 176: IRAM-30m EMIR time/sensitivity esti
- Page 177 and 178: IRAM-30m EMIR time/sensitivity esti
- Page 179 and 180: IRAM-30m HERA time/sensitivity esti
- Page 181 and 182: ¬þ«ðZñÛùþ˜ï³þ«ðWMõÑ
- Page 183 and 184: ¨-A¨-, #f¡§¨\¤qŠ’‹¼€(
- Page 185 and 186: ¤H F!";
- Page 187 and 188: #O.u\J0$!&%©ÉQ©6¦©¨-.u& (HÌ>
- Page 189 and 190: ¨*.G©6E¦cǦ©¨&ẗØn¨"!-‡
- Page 191 and 192: ¢¡¤£¦¥¨§ P© ¡ £¦¥¨§
- Page 193 and 194: '¥0vB"(-7£O.'¥!5#&(-7)7)(:9.*d#J
- Page 195 and 196: ¥¤¡ b@‚¥_ƒO@02,.(-EFJbe O@#V
- Page 197: 02J$#J+'10-ï 3 02EFJ+')EKJ5L5†R(
- Page 201 and 202: A&A 517, A12 (2010)tel-00726959, ve
- Page 203 and 204: A&A 517, A12 (2010)tel-00726959, ve
- Page 205 and 206: A&A 517, A12 (2010)tel-00726959, ve
- Page 207 and 208: A&A 517, A12 (2010)tel-00726959, ve
- Page 209 and 210: A&A 517, A12 (2010)tel-00726959, ve
- Page 211 and 212: A&A 517, A12 (2010)tel-00726959, ve
- Page 213 and 214: A&A 517, A12 (2010)tel-00726959, ve
- Page 215 and 216: A&A 517, A12 (2010)tel-00726959, ve
- Page 217 and 218: Table C.1. Definition of the symbol
- Page 219 and 220: IRAM Memo 2011-2WIFISYN:The GILDAS
- Page 221 and 222: WIFISYN3. practiceWIFISYN3. practic
- Page 223 and 224: WIFISYNA. implementation planWIFISY
- Page 225 and 226: tel-00726959, version 1 - 31 Aug 20
- Page 227 and 228: 8.3 PERSPECTIVES 227tel-00726959, v
- Page 229 and 230: 3 REQUIREMENTS 44 CHANGES FOR END-U
- Page 231 and 232: 5 CHANGES FOR PROGRAMMERS 125 CHANG
- Page 233 and 234: A EXHAUSTIVE DESCRIPTION OF THE CHA
- Page 235 and 236: Chapitre 9Copyright: Stéphane Guis
- Page 237 and 238: 9.3 ACTIVITÉS 2008-2011 237tel-007
- Page 239 and 240: Contribution de l'Action Spécique
- Page 241 and 242: A. des multi-pixels à bure : une s
- Page 243 and 244: 5. besoins en services annexes, bé
- Page 245 and 246: Articles publiés dans des revues
- Page 247 and 248: ARTICLES PUBLIÉS DANS DES REVUES
- Page 249 and 250:
Mémos IRAM et ALMAtel-00726959, ve
- Page 251 and 252:
Actes de colloques nationaux et int
- Page 253 and 254:
ACTES DE COLLOQUES NATIONAUX ET INT
- Page 255 and 256:
tel-00726959, version 1 - 31 Aug 20