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