[daip] IMSTAT output

[Eric Greisen]

On 10/17/2013 5:25 AM, Karl Menten wrote:
> Eric,
>
> I loaded a single dish map into aips and added BMAJ, MIN, BPA using 
> PUTHEAD.
> Here's the IMHEAD:
> AIPS 2: Image=NONE      (MA) Filename=I1629ASCU850.ICLN  .   1
> AIPS 2: Telescope=                   Receiver=
> AIPS 2: Observer=                    User #=  111
> AIPS 2: Observ. date=01-JAN-2000     Map date=15-OCT-2013
> AIPS 2: Minimum=-6.75735712E-01      Maximum= 2.01859531E+01
> AIPS 2: ----------------------------------------------------------------
> AIPS 2: Type    Pixels   Coord value     at Pixel     Coord incr   Rotat
> AIPS 2: RA---TAN   720    16 30 49.920     360.50 -6.000    0.00
> AIPS 2: DEC--TAN   720   -24 30 00.000     360.50  6.000    0.00
> AIPS 2: ----------------------------------------------------------------
> AIPS 2: Coordinate equinox 2000.00
> AIPS 2: Map type=DIRTY               Number of iterations=       0
> AIPS 2: Conv size=  19.00 X  19.00   Position angle=   0.00
> AIPS 2: Maximum version number of extension files of type HI is   1
> AIPS 2: Maximum version number of extension files of type PL is  10
> AIPS 2: Maximum version number of extension files of type MF is   3
>
> With this information added, JMFIT made a correct deconvolution.
> However, IMSTAT just returns:
> >imstat
> AIPS 2: Mean= 8.400E-02 rms= 2.731E-02          over 441. pixels
> AIPS 2: Maximum= 1.6144E-01 at     77     70     1    1    1    1    1
> AIPS 2: Skypos: RA 16 32 55.019  DEC -24 58 59.69
> AIPS 2: Minimum= 1.3789E-02 at     57     62     1    1    1    1    1
> AIPS 2: Skypos: RA 16 33 03.858  DEC -24 59 47.21
>
> and omits the "usual" final lines, e.g.:
> AIPS 1: Skypos: IPOL  1420.506 MHZ
> AIPS 1: Flux density =  2.2282E-04 Jy.   Beam area =  98.71 pixels
> (The above values are from a different file.)
>
> The file I'm working on, I1629ASCU850.ICLN  .   1, doesn't have any 
> frequency information (nor 3rd axis) in the header, so the omission of 
> the Skypos: IPOL line is understandable.
>
> However, I assumed that IMSTAT would do the "statistics" with the 
> information available and would deliver the total flux density in 
> the area defined by BLC and TRC plus the beam area. Out of the latter, 
> the area and the rms/pixel I then could calcylate the error in the 
> total flux.
>
> What am I missing? Do I have to PUTHEAD another paremeter?
>
> Thanks in advance,
> cheers,
> Karl
>
> ---
> Prof. Dr. Karl M. Menten
> kmenten at mpifr-bonn.mpg.de <mailto:kmenten at mpifr-bonn.mpg.de>
> Director, Millimeter & Submillimeter Astronomy
> Max-Planck-Institut fuer Radioastronomie
> Auf dem Huegel 69, 53121 Bonn, Germany
> +49 228 525297
>

There is another parameter which PUTHEAD should let you change to indicate
that the image is a normal Clean image rather than dirty.  IMHEAD should
probably not care about this and I will be allowed to fix it tomorrow (we
go back to work!).

Eric Greisen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20131017/225b866e/attachment.html>