imageStdPerChannel

Return to directory of Todd's CASA extensions

Uses the axes parameter of imstat to compute the statistics of each plane of a CASA cube over a specified region or annulus in the primary beam response and generates a profile plot of maximum, sigma, rms, scaled MAD, and (optionally) minimum. See also imageStd.

Usage:

au.imageStdPerChannel(img, region='', pbimg='', chans='', pbInner=0.3, pbOuter=0.2, 
                 std=True, mad=False, rms=False, plotmin=False, plotfile=''):
Inputs:
  • img: CASA cube
  • region: region to examine (default = whole image)
  • pbimg: if not blank, then set mask based this primary beam image, using the range of pixels between radii of pbInner and pbOuter
  • chans: range of channels to use (e.g. '4~50'), blank --> all channels
  • std, mad, rms: which values to return
  • plotmin: if True, then also plot the min
  • plotfile: name of png to produce
Returns:
  • list of four values [np.median(v), np.min(v), np.max(v), MAD(v)]

Example:

au.imageStdPerChannel('n2d+.image', pbimg='n2d+.pb', plotfile='n2d+_perchannel.png')

-- ToddHunter - 2016-05-22
  • n2d_example.png:
    n2d example.png
Topic revision: r4 - 2016-07-14, 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