[daip] Hi

[Eric Greisen]

Anish Roshi wrote:
> Hi Eric,
> 
> I have a question regarding converting the geometry of Xgaus output.
> 
> I used Xgaus to fit Gaussian components to line profile in a data cube 
> of size 256x256x63. Xgauss produces a set of FITS file consisting of the 
> amplitude, width etc of the Gaussian components. The output FITS files 
> has size 90x93x1 and is centered at a position different from the 
> original image center. I want to use  KNTR to plot the Gaussian 
> amplitude along with continuum image (size 256x256x1; centered at the 
> same center as the data cube) of the same region. I used hgeom to 
> convert the Xgaus output to the same geometry as the continuum image but 
> it gives the following error
> 
> ain.cv <http://ain.cv>> HGEOM1: Task HGEOM  (release of 31DEC10) begins
> ain.cv <http://ain.cv>> HGEOM1: You are using a non-standard program
> ain.cv <http://ain.cv>> HGEOM1: Create W12CCLAMP   .HGEOM .   1 (MA)  on 
> disk  1  cno   22
> ain.cv <http://ain.cv>> HGEOM1: Using       14961 words of buffer
> ain.cv <http://ain.cv>> HGEOM1: Begin computing coordinate tables
> ain.cv <http://ain.cv>> HGEOM1: Begin computing output image
> ain.cv <http://ain.cv>> HGEOM1: *****  ALL OUTPUT PIXELS BLANKED  *****
> ain.cv <http://ain.cv>> HGEOM1: SUBROUTINE MAPIO  RETURNED ERROR    16 
> AT 165
> ain.cv <http://ain.cv>> HGEOM1: ILLEGAL WINDOW SPECIFICATION
> ain.cv <http://ain.cv>> HGEOM1: DISK 1 SLOT  22  USER   108  (MA)FILE: 
> W12CCLAMP.HGEOM.1
> ain.cv <http://ain.cv>> HGEOM1: Destroyed MA image file: catno=     22 
> disk= 1
> ain.cv <http://ain.cv>> HGEOM1: Purports to die of UNNATURAL causes
> ain.cv <http://ain.cv>> HGEOM1: localhost    31DEC10 TST: Cpu=       
> 0.0  Real= 
> 
> Could you help me to find a solution to this issue ? (I could plot the 
> outputs of
> Xgaus output individually and there are non-blank data in it.)
> 
> Do you have any other suggestion to make plots of Gaussian fit outputs 
> from xgauss and overlay with the continuum image in AIPS ?

HGEOM or OHGEO are the 2 tasks to convert one image to the geometry of 
another - not OGEOM.

XGAUS has undergone major revision and could have issues.  Can you send 
me the IMHEADERs of the inout image to XGAUS and of one of the output 
images from XGAUS?  The 2 headers in your second image are separated by 
2 arc minutes in declination and 2.5 minutes of time in RA.  That means 
they do not overlap (2.5 minutes of time = 37.5 minutes of arc or 560 
pixels or so).

The BLC and TRC used in XGAUS would also be helpful.

Eric Greisen