imstat_cas6745

Return to directory of Todd's CASA extensions

Runs CASA task imstat using a work-around for CAS-6745, which is a bug that crashes imstat if the BUNIT is not recognized. If it detects that BUNIT is equal to the specified value, it copies the image to *.temp, sets the BUNIT to 'K', runs imstat, then deletes the copy.

Usage:

au.imstat_cas6745(img, bunit='"1e15 cm^-2"')

Example:

CASA <3>: au.imstat_cas6745('g12.68/nnh3.g12.68_1_out.image')
2014-07-02 09:36:03     WARN            FITS unit ""1e15 cm^-2"" unknown to CASA - will treat it as non-dimensional.
2014-07-02 09:36:03     WARN            FITS unit ""1e15 cm^-2"" unknown to CASA - will treat it as non-dimensional.
  Out[3]:
{'blc': array([0, 0], dtype=int32),
 'blcf': '18:13:57.645, -18.02.28.999',
 'max': array([ 59.22025681]),
 'maxpos': array([412, 425], dtype=int32),
 'maxposf': '18:13:54.756, -18.01.46.500',
 'mean': array([ 15.13657299]),
 'medabsdevmed': array([ 7.3842864]),
 'median': array([ 8.11978531]),
 'min': array([ 0.]),
 'minpos': array([448, 454], dtype=int32),
 'minposf': '18:13:54.504, -18.01.43.600',
 'npts': array([ 820.]),
 'quartile': array([ 24.68362236]),
 'rms': array([ 22.51512337]),
 'sigma': array([ 16.67795476]),
 'sum': array([ 12411.98984801]),
 'sumsq': array([ 415683.26065449]),
 'trc': array([839, 839], dtype=int32),
 'trcf': '18:13:51.763, -18.01.05.099'}

CASA <4>: 

-- ToddHunter - 2014-07-02