Typical characterisation for ground-based radio interferometry data
This is applicable to data using the Griessen/Calabretta/Valdez definitions, if there are any differences that is my mistake! It also assumes AIPS processing or equivalent, hence complications like the array coordinates can be ignorred.
Other examples are possible, e.g. the mapping for solar or planetary observations.
| | Spatial | Temporal | Spectral | Velocity | Observable |
| Mapping Frame | ICRF/J2000/B1950+Epoch(1) | UT1/UTC(2) | Frequency etc. | LSR/Helio, radio/opt/relativisitic | Jy/beam/, Jy, phase etc. |
| Location | R.A., Dec. of field centre.# | Date (start) | Band(3) | Source velocity+transition rest freq+reference channel# | Radio interferometry visibilities, images, other products; polarization products (4) |
| Location Reference | Phase Ref source, cat/bibref | | | catalogue or biblio refs for source and transition | Flux density reference source or load# + cat/bib ref |
| Location Uncertainty | Uncertainty in Phase Ref pos plus data-provider-supplied (5) | | | From cat/bib refs | From cat/bib refs or data provider |
| Sensitivity Bounds | Visibility data: max, min baselines (Bmax, Bmin); Image: Geometric area | Start/stop time | Max/min freq | Min/max vel. | rms noise in each Stokes (IQUV) and pol. intensity. Pol. angle rms.(6) |
| Support | Primary beam FWHM (PB). Size of contiguous image area, No. of pointings for mosaics | Length of each track (t), No. of tracks(Nt)(7) | Max/min or centre(freq)/freq width(Dfreq)+ref. chan. for each continuous sub-band | Max/min or centre+ref. chan. for each continuous velocity segment(8) | |
| Filling* | Visibility data: quality factor for baseline coverage e.g. No. visibilities, No. Antennas, Hour Angle coverage. Mosaic image filling factor. | Fractional time on source | Fraction of channels discarded(9) | Fraction of channels kept | |
| Sensitivity Function | F.o.V. functions(10) | rms noise as function of time(11) | Bandpass function# (usually look-up table)(12) | | Max. dynamic range, dependence on weighting(13) |
| Resolution | Restoring Beam(14) | Integration time Tint | Effective channel width (may be averaged) | Effective channel width | |
| Sample Precision | Visibility data: Tint. Image data: pixel size (arbitrary, conventionally~<Bmin/3) | Tint | Channel width dfreq | Channel width | |
*Filling - is this relative to Sensitivity Bounds or to Support? I assume Support. It might be better not to force it to be a factor so that No. visibilities can be included, unless there is a better place.
#Some radio interferometry data prodivers will have applied astrometric, flux scale, bandpass, Vlsr etc. corrections; others simply provide a package which includes the means to do so.
(1) For milli-arcsec accuracy in precessing coordinates it may be necessary to make leap second and epoch corrections, which is probably best done using a dedicated radio astronomy package e.g. AIPS UVFIX. In any case high accuracy double precision is needed, no worse than e.g. the highest resolution output of 'coco'.
(2) Coordinate conversion needs to be done to the nearest second, see (1). VLBI requires much higher precision pre-correlation but the VO should not need to worry about that. Pulsar data also need much better precision but there is a separate PVO being developed. In these instances the temporal Location etc. would be GPS, local maser etc. time/freq. standards. An observation can run over many days or months for a single image.
(3) Band may be a
letter but to be VO compatible it should be the centre freq. + width or the upper and lower limits. Usually accuracy is much better than finest possible resolution, might need to include for SETI and some specialised projects.
(4) If visibility data, this could have a range of potential products.
Or may be a single product/set of products e.g. datacube
Polarisation products (order matters) :
Visibilities - 1 - 4 from circular (LL RR LR RL) or linear (XX YY XY YX);
Images (may be one per polarization or a cube) - as visibilities or Stokes I Q U V.
(5) The systematic position uncertainty is the sum of the uncertainty in the phase-ref source position which could be updated by the VO, and uncertainties in the phase transfer, antenna positions etc. which should be provided as numbers.
(6) This (and all the bounds) are inclusive extrema for Registry use etc., e.g. the lowest noise in the data, might be worse in some parts of the image, or if different weighting is used.
(7) Typically, an observation can be spread over 3 days (bounds). Each day, for each source, a track is spread over 12 hr (support) within which the source is observed for 4 mins out of each 6 mins (could be much less e.g. 4 mins per hr) (filling).
(8) If several spectral lines are observed in one frequency band the velocity coverage in the source rest frame will be different for each line. However if e.g. two lines are observed in two bands the frequency coverage will be different but the velocity coverage may be the same.
(9) Occasionally a correlator mode may be used which discards, e.g., every other channel giving 50% filling factor
(10) Field of View limit to 50% sensitivity (F50) (e.g. for images to be extracted from visibility data) is given by the smallest of:
F50~PB (Primary Beam to Full Width Half Max) of largest antenna in the array
F50~1E+4 Bmaj/Tint (Time smearing due to Earth rotation and phase rate in integration time)
F50~2(freq/dfreq)Bmaj (Frequency smearing in channel width used for imaging; multi-channel data can be used to produce a single image).
In each case the F.o.V. is in the units of the PB or Bmaj and Tint is in sec.
These are really continuous non-linear functions e.g.
MERLIN version which express how the observable deteriorates as a function of distance from the centre of the F.o.V. (i.e. rms(F10)~1.1xrms).
(11) A long observation (one or more tracks) is usually designed to be averaged in time. However for variable sources one can measure the visibility amplitudes or spectra on shorter timescales. The minimum is the sampling interval (integration time) Tint
The noise in visibility data scales as root(time on source), i.e. rms(tsegment)=rms(full obs). root[(t x Nt)/tsegment]
For images there is an additional constraint of the minimum hour angle coverage needed to make an image (min, synthesis time) which the data provider would need to specify.
(12) The bandpass correction is usually automatically applied where needed. However for wide-band continuum images an estimate of the source spectral index may be needed. Leave to data provider/user at present!
(13) Image dynamic range D(n) if n channels are averaged as a function of dynamic range for whole band, Dmax e.g. D(n) ~ Dmax(n.dfreq/Dfreq)^2
where Dmax = peak flux density in image/rms noise
Weighting data to improve resolution increases noise, either data-provider provided or rms(weighted)/rms ~ Bmaj/Bmaj(weighted)
(14) Natural beam Bmaj ~ c/(freq.Bmax); restoring beam can be varied by factor ~2 up or down by weighting giving Bmaj(weighted). For uv ellipses of high eccentricity Bmax, Bmin and Bpa may be given.
--
AnitaRichards - 22 Jun 2004
Examples
Note that I have tried to fill in all fields but in practice only certain fields would be relevant for any specific data retrieval. Other details (e.g. bandpass calibration source used) are part of the data history and could be got elsewhere in the model (e.g. source schedule). However all fields may be required under some circumstances e.g. the bandpass calibration source name is needed if the bandpass correction hasn't already been applied.
Typical MERLIN archive data
| Markarian 231 | Spatial | Temporal | Spectral | Velocity | Observable |
| Mapping Frame | ICRF | UTC | Freq. | | Jy, Jy/beam |
| Location | 12 56 14.2334 56 52 25.001 | 1993-05-23 | 1420 MHz | | LL RR LR RL I Q U V |
| Location Reference | B1300+580 13 02 52.465276 57 48 37.60941 (1998AJ....116..516M) | | | | 3C286 (1977A&A....61...99B) |
| Location Uncertainty | 10 arcsec | | | | 5% |
| Sensitivity Bounds | Spatial scales 0.15-7.2 arcsec; max image region 30 arcmin | 1993-05-23T16:00:07- 1993-05-25T16:59:01 | 1412-1428 MHz | | I 70 uJy/beam, Q,U 60 uJy/beam, pol. leakage 0.1%; pol angle uncertainty 2 deg. |
| Support | PB FWHM 10-30 arcmin | 2 runsx~16 hr | 15x1MHz | | |
| Filling | 34987 visibilities, 8 ants | 6.45 hr out of 32 | | | |
| Sensitivity Function | 10% reduction at 80 arcsec radius | 100 uJy/beam per run; 2.5 mJy per integration | Correction applied (derived from OQ208) 1sigma noise 0.28 mJy/MHz | | |
| Resolution | uvrange 15-1030 klambda; 0.12-0.4 arcsec | 16 sec | 1 MHz | | |
| Sample Precision | typical pixel 0.4 arcsec | 16 sec 1 MHz | | | |
| Markarian 231 | Spatial | Temporal | Spectral | Velocity | Observable |
| Mapping Frame | ICRF | UTC | Freq. | | Jy, Jy/beam |
| Location | 12 56 14.2334 56 52 25.001 | 1993-05-23 | 1658 MHz | | LL RR LR RL I Q U V |
| Location Reference | B1300+580 13 02 52.465276 57 48 37.60941 (1998AJ....116..516M) | | | | 3C286 (1977A&A....61...99B) |
| Location Uncertainty | 10 arcsec | | | | 5% |
| Sensitivity Bounds | uvrange 17-1198 klambda; Spatial scales 0.15-7.2 arcsec; max image region 30 arcmin | 1993-05-23T16:00:07- 1993-05-25T16:59:01 | 1650-1666 MHz | | I 90 uJy/beam, Q,U 70 uJy/beam, pol. leakage 0.1%; pol angle uncertainty 2 deg. |
| Support | PB FWHM 10-30 arcmin | 2 runsx~16 hr | 15x1MHz | | |
| Filling | 34987 visibilities, 8 ants | 7.27 hr out of 32 | | | |
| Sensitivity Function | 10% reduction at 80 arcsec radius | 100 uJy/beam per run; 2.5 mJy per integration | Correction applied (derived from OQ208) 1sigma noise 0.28 mJy/MHz | | |
| Resolution | 0.12-0.4 arcsec | 16 sec | 1 MHz | | |
| Sample Precision | typical pixel 0.4 arcsec | 16 sec | 1 MHz | | | |
Spectral data for Mrk 231 megamasers (2 OH lines in the same band after processing)
| Markarian 231 | Spatial | Temporal | Spectral | Velocity | Observable |
| Mapping Frame | ICRF | UTC | Freq. | Vlsr, radio | Jy, Jy/beam |
| Location | 12 56 14.2383 56 52 25.210 | 1997-06-17 | 1.6 GHz | 1667.359 MHz shifted to 12135 km/s in chan 47 or 16675.402 MHz shifted to 12135 km/s in chan 32 | LL RR I |
| Location Reference | B1300+580 13 02 52.465276 +57 48 37.60941 (1998AJ....116..516M) | | | rest freq NIST Vlsr 1987MNRAS.226..689S | 3C286 (1977A&A....61...99B) |
| Location Uncertainty | 10 arcsec | | | rest freq 1 kHz; Vlsr 50 km/s | 5% |
| Sensitivity Bounds | uvrange 0.035-1.2 Mlambda; Spatial scales 0.12-3 arcsec; max image region 30 arcmin | 1997-06-17T08:32:00 - 1997-06-18T07:00:00 | 1594-1604 MHz (at LSR) | 1667 MHz line 11408-13180 km/s; 1665 MHz line 11070-12842 km/s | 0.7 mJy/beam/channel |
| Support | PB FWHM 10-30 arcmin; maser images within blc 12 56 14.40 56 52 24.0 trc 12 56 14.05 56 52 26.5 | 22.5 hr | 79x0.1250 MHz | 79x22.5 km/s | |
| Filling | 6 ants | 35% per line | | | |
| Sensitivity Function | 10% reduction at 80 arcsec radius | | Instrumental bandpass correction applied (derived from B2134+004); linear continuum of (-0.0006 Jy/MHz+1.11076(0.003) Jy) removed | | 0.7 mJy/beam in central ~40 chans, ~1 mJy/beam elsewhere |
| Resolution | 0.12-0.4 arcsec | 16 sec | 79x0.125 MHz | 79x22.5 km/s | |
| Sample Precision | typical pixel 0.4 arcsec | 16 sec | 0.125 MHz | 22.5 km/s | | |
MERLIN+VLA HDF OBSERVATIONS
Deep field observations; typically calibrated visibility data would be available for extracting images at chosen position. Sources usually very faint so separate imaging of small time/frequency/polarization sections of data not practical, but some information about resolution in all domains is needed for F.o.V.
| HDF(N) | Spatial | Temporal | Spectral | Velocity | Observable |
| Mapping Frame | ICRF | UTC | Freq. | | Jy, Jy/beam |
| Location | 12 36 49.4 +62 12 58.00 | 1996-1997 | L-band | | Visibilities (LL, RR); images (I) |
| Location Reference | 13 13 27.9860 +67 35 50.380 (2003AJ....126.2562F); 12 41 29.5908 +60 20 41.320 (MERLIN) | | | | 3C286 (1977A&A....61...99B) |
| Location Uncertainty | under 15 mas | | | | ~2% |
| Sensitivity Bounds | Spatial scales 0.2-64 arcsec; max. image region 40 arcmin radius | 1996-02-06 - 1997-04-27 | 1354-1427 MHz | | I 1sigma 3.3 uJy/beam |
| Support | PB FWHM 10-30 arcmin | MERLIN 1996-02,1996-03,1997-04 27 d; VLA 1996-11 2.5 d | MERLIN 1405-1419, 1413-1427 MHz; VLA 1354-1376, 1424-1446 MHz | | |
| Filling | MERLIN 6-8 ants; VLA 27 ants. | MERLIN 438 hr on-source; VLA 42 hr | | | |
| Sensitivity Function | 50% peak sensitivity at 5 arcsec radius | MERLIN-only 1sigma 5.9 uJy/beam VLA-only 7.5 uJy/beam | Corrrection applied derived from B1803+784 | | 3.3-5 uJy/beam (min. at 0.5 arcsec resolution) |
| Resolution | 0.2-2 arcsec | 4-10 sec | 28x0.5; 14x3.125 MHz, combined to make images | | |
| Sample Precision | 0.05-0.4 arcsec pixels | 4-10 sec | 0.5-3.125 MHz | | |
--
AnitaRichards - 24 Jun 2004