january 2010 arpg rfi mitigation in aips. the new task uvrfi l. kogan, f. owen national radio...

14
January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

Upload: job-merritt

Post on 18-Jan-2018

213 views

Category:

Documents


0 download

DESCRIPTION

The Athreya’s presentation of the problem. Concept of circles. Ramana Athreya, APJ, 696, , May 2009 If RFI includes only one group moving with constant velocity (in particular located at the Earth or one satellite) and its amplitude is constant inside of the solution interval then the trajectory of the correlator output at the complex plane is the circle with radius equaled RFI amplitude and the center coordinates are Re, Im of the source visibility. Fitting the 3 parameters: radius of the circle, and two coordinates of the circle center, the found two coordinates of the circle center are used as the solution for the source visibility without RFI January 2010 ARPG

TRANSCRIPT

Page 1: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

January 2010 ARPG

RFI Mitigation in AIPS.The New Task UVRFI

L. Kogan, F. OwenNational Radio AstronomyObservatory

Socorro, NM USA

Page 2: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

Visibility expression at the presence of RFIs moving with constant velocities relatively the array

Having known the correlator outputs for the times inside of the solution interval and considering the all RFI amplitudes are constant inside of the solution interval, find the value of the source visibility for the solution interval!

January 2010 ARPG

source. by the providedsibility complex vi theof lueunknown va theis

correlator theofoutput at thesibility complex vi theof value theis

:where

....)(exp)(exp)(exp

:expression following by the described be willcorrelator theofoutput result the a As

).exp(exponent complex stopping fringe by they vilibilitobserved themultiplyes correlatorA

....expexpexp:expresion following by the determined is

timeon,polarizati channel, frequency baseline,given at they visibilitThe

2211

2211

source

cor

frsofrsasatfrsofrsasatfrsogrsourcecor

frso

frsasatfrsasatgrfrsosourceobs

Vis

Vis

tjRFItjRFItjRFIVisVis

tj

tjRFItjRFIRFItjVisVis

Page 3: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

The Athreya’s presentation of the problem. Concept of circles.

Ramana Athreya, APJ, 696, 885-890, May 2009

If RFI includes only one group moving with constant velocity (in particular located at the Earth or one satellite) and its amplitude is constant inside of the solution interval then the trajectory of the correlator output at the complex plane is the circle with radius equaled RFI amplitude and the center coordinates are Re, Im of the source visibility.

Fitting the 3 parameters: radius of the circle, and two coordinates of the circle center, the found two coordinates of the circle center are used as the solution for the source visibility without RFI

January 2010 ARPG

Page 4: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

One of the best “circle”

January 2010 ARPG

Page 5: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

The new AIPS task UVRFIUVRFI uses the 2 algorithms: 1. ‘CIRC’ a la Athrea The spiral with 4 unknown parameters: initial radius, linear increment of the radius, 2 coordinates of the center is fitted to the data. The 2 coordinates of the center are used as a solution for the free RFI visibility.

2. ‘CEXP’ The model is represented by the sum of several spectral components with complex amplitude.

The simple version of Hogbom clean algorithm (complex Fourier transform, complex solution and complex “PSF) is used to subtract the the complex components. Final solution is the value of the cleaned Fourier transform at zero frequency. The subtraction is not carried out for the component at zero frequency to prevent the subtraction of the signal itself. This rule is canceled for the very high channel amplitude in comparison with other channels.

January 2010 ARPG

....expexp 2211 tjRFItjRFIVisVis sourcecor

Page 6: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

January 2010 ARPG

The actual shape of the “circles”. L-band data (given by M. Rupen).

Page 7: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

January 2010 ARPG

The visibility spectrum for one baseline of Rupen’s data (L band)

Page 8: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

The visibility spectrum for one baseline of Rupen’s data (L band)

The averaged output The UVRFI output

January 2010 ARPG

Page 9: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

The image of 3C345 using Rupen’s data (L band)

January 2010 ARPG

The averaged output The UVRFI outputThe averaged output The UVRFI output

Page 10: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

The actual shape of the “circles”. 4-band data (given by B. Cotton).

January 2010 ARPG

Page 11: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

The visibility spectrum for one baseline of Cotton’s data (4 band)

January 2010 ARPG

The averaged output The UVRFI outputThe averaged output The UVRFI output

Page 12: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

The images using the Cotton’s data (4 band)

January 2010 ARPG

The averaged output The UVRFI outputThe averaged output The UVRFI output

Page 13: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

Self averaging of the visibilities in the process of imaging

Therefore the effect of self averaging of RFI may not limited by the correlator averaging interval (as Athrea declares) but by the time interval where both U and V can be considered constant!

January 2010 ARPG

ikifrk

kkrfi

imagii

ii

iifrrfi

ifr

tjmVlUjADM

TVUi

mVlUjtjADM

tjARFI

,exp)(2exp

: then interval time theof insideconstant are , Ifnumber interval timepreaverage theis where

)(2expexp

exp

Page 14: January 2010 ARPG RFI Mitigation in AIPS. The New Task UVRFI L. Kogan, F. Owen National Radio AstronomyObservatory Socorro, NM USA

ConclusionsThe new AIPS task UVRFI uses the two algorithms to mitigate RFI.

1.Modified algorithm of Ramana Athrea fitting the spiral into the curve RE/IM of the visibility.2.Subtracting the complex exponents representing RFI, using a simple CLEAN algorithm at the Fourier plane3.The second algorithm showed better result for the two data I used and allows to mitigate more than one sources of RFI (ground, satellites)4.The mitigation of the RFI is not so impressive because the RFI of the data are used are too far from the circle model.5.In some cases, the RFI may be self averaged during the imaging

January 2010 ARPG