imstatBackground

Run imstat on the specified box region, then on the 4 adjacent boxes to determine the lowest rms in all 5 boxes, and return that value in the dictionary (instead of the rms in the box). See also buildBoxRegion.

Return to directory of Todd's CASA extensions

Usage:

au.imstatBackground(img, region, prec=7, verbose=False)
  • img: image filename
  • region: a CASA region string

Example:

au.imstatBackground('n6334_I_5GHz_cont_s_ap_scale.image.tt0.smoothbeam.pbcor','box [[17:20:55.24514, -035.45.04.6612], [17:20:55.13641, -035.45.02.7920]] coord=J2000')
rmss =  [0.0001, 1.5400000000000002e-05, 1.6099999999999998e-05, 1.9599999999999999e-05, 1.3499999999999999e-05]
min rms = 0.000013 (replacing 0.000080)
  Out[3]:
{'blc': array([1111, 2422,    0,    0], dtype=int32),
 'blcf': '17:20:55.247, -35.45.04.679, I, 5.0635e+09Hz',
 'flux': array([ 0.00021194]),
 'max': array([ 0.00032743]),
 'maxpos': array([1122, 2436,    0,    0], dtype=int32),
 'maxposf': '17:20:55.193, -35.45.03.839, I, 5.0635e+09Hz',
 'mean': array([  2.44587893e-05]),
 'medabsdevmed': array([  1.25337874e-05]),
 'median': array([ -6.81679103e-06]),
 'min': array([ -4.55350055e-05]),
 'minpos': array([1115, 2440,    0,    0], dtype=int32),
 'minposf': '17:20:55.228, -35.45.03.599, I, 5.0635e+09Hz',
 'npts': array([ 736.]),
 'quartile': array([  3.91421709e-05]),
 'rms': 1.3499999999999999e-05,
 'sigma': array([  7.59676038e-05]),
 'sum': array([ 0.01800167]),
 'sumsq': array([  4.68204047e-06]),
 'trc': array([1133, 2453,    0,    0], dtype=int32),
 'trcf': '17:20:55.139, -35.45.02.819, I, 5.0635e+09Hz'}
)

-- ToddHunter - 2013-07-18
Topic revision: r1 - 2013-07-18, 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