wvrCorrelation

Return to directory of Todd's CASA extensions

Computes the covariance between the time series of each of the WVR channels with the specified channel-averaged spw, normalized by the variance of the WVR channel. This is meant to explore the conditions in which the WVR discrepancy is high to see if it is correlated with clouds.

Usage:

wvrCorrelation(vis, spw=None, antenna=0, field=None, sigma=3, plotfile='', scaleAstronomicalReceiver=True, 
         pwv=None, normalizeToWVR3=False, buildpdf=False, verbose=False)
Inputs:
  • spw: spw ID (integer or string), if None then use first science spw
  • antenna: antenna ID or name or list. '' ==> all antennas
  • sigma: threshold to clip outlier points in the astronomical receiver datastream
  • pwv: the pwv (in mm) for the observation, in case ASDM_CALWVR is not present
  • normalizeToWVR3: if True, normalize astronomical data to Tsky of WVR 3 if False, normalize to Tsky of the astronomical spw
  • scaleAstronomicalSignal: if True, then magnify the black trace to occupy half the vertical scale (in order to better view small time variations)

Returns:
  • The maximum correlation value of the 4 WVR channels.
  • If more than one antenna, then return the median of the antenna maxima.

Example:

CASA <11>: au.wvrCorrelation('uid___A002_X8358e4_X917.ms',spw=2,plotfile=True,antenna='PM04')
Working on antenna PM04 (1/1)
Scan 1,2,3,4,5,6,7,8,9,10,11,12,13 (field=0,1,2,3,4,5,6): azim = -17.04,  elev = 59.86  (degrees)
Plot left in  uid___A002_X8358e4_X917.ms.wvrCorrelation.PM04.spw02.png
  Out[11]: 0.076124930755638934

-- ToddHunter - 2014-08-29
Topic revision: r1 - 2014-08-29, ToddHunter
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding NRAO Public Wiki? Send feedback