> If I collapse the cube to produce an average spectrum, and then try to plot
> the average spectrum in the same way:
>
> linplot avspec'(-20.0:20.0)' plots a spectrum from 2.0 to -5.0 km/s (left to
> right).
>
> To collapse the cube, I used collapse in=file out=temp axis=1 estimator=mean
> and collapse in=temp out=avspec axis=1 estimator=mean
When I tried this on another spectral cube, the second collapse barfed
because astGet couldn't find the lataxis(1) attribute. I then tried
MANIC to collapse in one command.
% manic cube avspec 3
That seemed to work but left behind the spatial WCS. LINPLOT of that
worked but trying to display a velocity section gave a segmentation
fault. I tried NDFCOPY to trim the unwanted spatial WCS information but
still the LINPLOT segmentation persisted. Although yesterday I did
upgrade the OS and have yet to rebuild the Starlink tree, but I see the
same problem on another 32-bit Linux box.
If David wants to investigate, I can supply my NDFs. LINPLOT of
sections is such a common operation, it does suggest something specific
to this dataset, and likewise for Jane's data.
Malcolm
|