r6 - 24 Jun 2004 - 16:00:13 - AnitaRichardsYou are here: TWiki >  Astrogrid Web  >  DocStore > InterOperability > InteropDataModel > ObservationDataModelRevision > CharacterisationRadioInterferometry

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

Edit | Attach | Printable | Raw View | Backlinks: Web, All Webs | History: r6 < r5 < r4 < r3 < r2 | More topic actions
 
AstroGrid Service Click here for the
AstroGrid Service Web
This is the AstroGrid
Development Wiki
This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback