> > Seems to work for me (in fact it works for me on files that core dump with
> > GETBOUND). I note that GETBOUND presents all the bounds rather than just a
> > specified coordinate but that may be a red herring if DATACUBE only ever
> > looks at the spectral axis.
>
> GETBOUND presumably calls AST_MAPBOX multiple times, once for each axis.
Yep!
* Find the bounding box of the transformed co-ordinates
DO IDIM = 1, NVIN
CALL AST_MAPBOX( MAPC, ASTART, AEND, .TRUE., IDIM,
: DDLBND( IDIM ), DDUBND( IDIM ), DDXL, DDXU,
: STATUS )
END DO
> > The other outstanding datacube issue is how to obtain the units of the
> > spectral axis. What's the recommended way to get the axis units from a
> > file? [transparently such that the units can come from the current frame
> > so that SpecFrames and AXIS will work].
>
> Units is just like all other attributes - you can get it, set it, test
> it, or clear it using kappa:wcsattrib. Fot instance, if the spectral axis
> is axis 2, you can do:
>
> % wcsattrib fred get "unit(2)"
Looking at this thread it does seem that the required changes should all
be easily doable. I'll should sit down and do this post-ADASS perhaps? Is
that sort of timescale okay with JAC?
Al.
|