Date - 2004.03.11
Tester - J.E. Hibbard (CV)
Platform: Merger
Version: Daily AIPS++ Version 1.9 Build 473

March 2004. Testing results for imageprofilefitter.

Have to change coordinate system of cube that is read in; Channels are
equally spaced in frequency, not velocity. aips++ does not parse the
ALTRVAL, ALTRPIX, RESTFREQ FITS keywords that AIPS uses. So, bear with
me:

include 'image.g';
- f:=imagefromfits(outfile='1zw49cube',infile='IZW49C.R-1.LMV.FITS')
- im:=image('1zw49cube')
- cs:=im.coordsys();
- local paxes,waxes
- cs.findcoordinate(paxes,waxes,'spectral')
- rp:=cs.referencepixel();
- rp[paxes[1]]:=-8
- cs.setreferencepixel(rp)
- rv:=cs.referencevalue('q')
- rv[waxes[1]]:=dq.quantity('1.41641720627e9Hz')
- cs.setreferencevalue(rv)
- df:=cs.increment('q')
- df[waxes[1]]:=dq.quantity('1.2207031e4Hz')
- cs.setincrement(df)
- rf:=dq.quantity('1.420405752e9Hz')
- cs.setrestfrequency(rf)
- cs.summary()
Direction reference : J2000
Spectral  reference : BARY
Velocity  type      : RADIO
Rest frequency      : 1.42041e+09 Hz
Pointing center     :  12:56:37.800000  +35.06.48.000000
Telescope           : VLA
Observer            : AH453
Date observation    : 1992/04/11/00:00:00

Axis Coord Type      Name             Proj   Coord value at pixel    Coord incr Units
------------------------------------------------------------------------------------- 
1    1     Direction Right Ascension   SIN  12:58:59.899   256.00 -4.000000e+00 arcsec
2    1     Direction Declination       SIN +34.50.37.785   257.00  4.000000e+00 arcsec
3    3     Spectral  Frequency              1.416417e+09    -8.00  1.220703e+04 Hz
                     Velocity               8.418270e+02    -8.00 -2.576430e+00 km/s
4    2     Stokes    Stokes                            I

- im.setcoordsys(cs)
- cs.done()
- im.done()

- include 'imageprofilefitter.g';
- myprof:=imageprofilefitter(infile='1zw49cube',infile2='IZW49C.R-1.MOM2.FITS'); 


Well, the above parameters were what aips says, but they don't seem
to work. So I need to fudge it:


im:=image('1zw49cube')
cs:=im.coordsys();
rp:=cs.referencepixel();
rp[paxes[1]]:=28
cs.setreferencepixel(rp)
rv:=cs.referencevalue('q')
rv[waxes[1]]:=dq.quantity('1.41685665939e9Hz')
cs.setreferencevalue(rv)
cs.summary()
=>says V[28]=749.08km/s; should be 835. 
rv:=cs.referencevalue('q')
rv[waxes[1]]:=dq.quantity('1.4164609e9Hz')
cs.setreferencevalue(rv)
cs.summary()
im.setcoordsys(cs)
myprof:=imageprofilefitter(infile='1zw49cube',infile2='IZW49C.R-1.MOM2.FITS');

=>close enough.

Submitted this to AOCso4744