[daip] Wrong (u,v) coordinates?

Arnaud Collioud Arnaud.Collioud at obs.u-bordeaux1.fr
Mon Jul 7 11:00:50 EDT 2014


> Can you look at files from earlier in the processing to see when this issue arose?  I changed many tasks to shift the frequency reference to the center, but I see that you have not done that.  This shift does require a re-scaling of the uvw's but a very minor one.  Were the X and S data in the same file at one time?  There could be an issue then if the header had the X band frequency because the u,v,w should then be in wavelengths at X band even on the S band data.  Then, splitting out S band would lead to a really big correction and I wonder if it happened in some way twice???

Yes, originally X and S frequencies were together (even if there were no X data). I split the 2 using UVCOP with the following parameters:

AIPS 1: UVCOP     Task to copy selected times and channels of uv data.
AIPS 1: Adverbs     Values                 Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME     'IYA09_S.3'             Input UV file name (name)
AIPS 1: INCLASS    'UV'                    Input UV file name (class)
AIPS 1: INSEQ         1                    Input UV file name (seq. #)
AIPS 1: INDISK        1                    Input UV file disk unit #
AIPS 1: OUTNAME    'IYA09_S.CLN'           Output UV file name (name)
AIPS 1: OUTCLASS   'UV'                    Output UV file name (class)
AIPS 1: OUTSEQ        1                    Output UV file name (seq. #)
AIPS 1: OUTDISK       1                    Output UV file disk unit #.
AIPS 1: SELBAND      -1                    Bandwidth to select (kHz)
AIPS 1: SELFREQ      -1                    Frequency to select (MHz)
AIPS 1: FREQID        2                    Freq. ID to select. 0=> all
AIPS 1: BCHAN         1                    First channel selected.
AIPS 1: ECHAN         0                    Last channel selected.
AIPS 1: BIF           0                    Lowest IF number 0=>all
AIPS 1: EIF           0                    Highest IF number 0=>all
AIPS 1: SOURCES    *all ' '                Source list
AIPS 1: QUAL         -1                    Source qualifier -1=>all
AIPS 1: UVRANGE       0           0        Annulus in uv plane selected
AIPS 1:                                     (thousands of wavelengths)
AIPS 1: TIMERANG   *all 0                  Timerange:
AIPS 1:                                     1-4 = start day,hr,min,sec
AIPS 1:                                     5-8 = end   day,hr,min,sec
AIPS 1: ANTENNAS   *all 0                  Antennas to copy 0=>all
AIPS 1: BASELINE   *all 0                  Baselines with ANTENNAS
AIPS 1: SUBARRAY      0                    Limit copy to one subarray
AIPS 1:                                    0 => all.
AIPS 1: FLAGVER      -1                    FG table to apply to data
AIPS 1:                                    <= 0 => do not apply
AIPS 1: UVCOPPRM      1           0        Flags:
AIPS 1:               0           1         1 > 0 => keep flagged data
AIPS 1:               0           0         2 > 0 => drop subarray info
AIPS 1:              -1        *rest 0      3 = 0 => copy xc and ac data
AIPS 1:                                       = 1 => copy only xc data
AIPS 1:                                       = 2 => copy only ac data
AIPS 1:                                          of antennas, named in
AIPS 1:                                          antenna list. Baseline
AIPS 1:                                          list must be zero.
AIPS 1:                                     4 > 0 => report progress
AIPS 1:                                     5 => flag data with weights
AIPS 1:                                          below this limit
AIPS 1:                                     6 = 1,3 -> do NOT flag TY/SY
AIPS 1:                                       = 2,3 -> do NOT flag SN
AIPS 1:                                     7 = extend copy of tables
AIPS 1:                                         by UVCO(7) min, 0 -> 15
AIPS 1: NCOUNT        0                    If UVCOPPRM(3) = 2 then
AIPS 1:                                    NCOUNT is the approximate
AIPS 1:                                    ratio of AC to XC in file.
AIPS 1: FQCENTER      0                    >= 0 -> center frequency axis

I did not change FQCENTER because I was not sure about what it will do (I was also not aware of the importance of it).

As you said, the initial file has the X-band frequency as reference in the header:
>imh
AIPS 1: Image=MULTI     (UV)         Filename=IYA09_S.3   .UV    .   1
AIPS 1: Telescope=VLBI               Receiver=MKIV
AIPS 1: Observer=Bonn Cor            User #=    4
AIPS 1: Observ. date=18-NOV-2009     Map date=03-JUL-2014
AIPS 1: # visibilities   4848965     Sort order  TB
AIPS 1: Rand axes: UU-L-SIN  VV-L-SIN  WW-L-SIN  TIME1  BASELINE
AIPS 1:            SOURCE  FREQSEL  INTTIM  WEIGHT  SCALE
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type    Pixels   Coord value     at Pixel     Coord incr   Rotat
AIPS 1: COMPLEX      1   1.0000000E+00       1.00  1.0000000E+00    0.00
AIPS 1: FREQ        32   8.4059900E+09       1.00  2.5000000E+05    0.00
AIPS 1: STOKES       1  -1.0000000E+00       1.00 -1.0000000E+00    0.00
AIPS 1: IF           4   1.0000000E+00       1.00  1.0000000E+00    0.00
AIPS 1: RA           1    00 00 00.000       1.00       3600.000    0.00
AIPS 1: DEC          1    00 00 00.000       1.00       3600.000    0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Maximum version number of extension files of type HI is   1
AIPS 1: Maximum version number of extension files of type BS is   1
AIPS 1: Maximum version number of extension files of type CL is   1
AIPS 1: Maximum version number of extension files of type SU is   1
AIPS 1: Maximum version number of extension files of type AN is   1
AIPS 1: Maximum version number of extension files of type FQ is   1

What should I do? 
Set FQCENTER to '1' in UVCOP?  Could you explain to me what this parameter is for?

Thanks!
Cheers,
Arnaud




More information about the Daip mailing list