Hi Dave,
Attached is a script that may help you, or at least give
you some pointers as to how to do what you want. It display a single
slice from a supplied cube. You use a cursor to select two positions.
It then extracts a position-velocity slice between those two
positions, displays it and saves it.
David
On 11 May 2010 11:49, tim.jenness <[log in to unmask]> wrote:
> On May 11, 2010, at 5:38 AM, David Nutter wrote:
>
>>> Do you have the equatorial co-ordinates of the four corners; or just two
>>> opposite vertices and an angle; or the centre co-ordinates, orientation, and
>>> the box length and height? Could you use a POLYGON or ROTBOX keyword in an
>>> ARD definition to specify the slice? SUN/183 says these are only 2-D so you
>>> would first have extract a plane at the chosen velocity and spatial
>>> co-ordinate limits that will enclose the region. Then run ARDMASK.
>>
>> Hi Malcolm.
>>
>> I would like to show a 2D image with a spatial dimension along the x
>> axis (a combination of x and y in the cube), and velocity (z in the
>> cube) on the y axis. This seems easy to do if your new x axis
>> corresponds to the x OR y axes in the cube, but not if you want an
>> arbitrary direction.
>>
>> I can specify any of the coords that you mention (centre/extent, or
>> opposite corners), but what you suggested seems to indicate that I
>> only get 1 velocity, which is not what I am trying to show. ..or have
>> I misunderstood you?
>
> I understand what you need. Effectively you want the PROFILE command to be able to extract an RA/Dec at all planes in the cube and then put it together. This is quite a lot of work and we'll put it on the list.
>
> In the short term using ROTATE and NDFCOPY is the way to go and we think we can make life a little bit easier by allowing you to specify two ra/dec coordinates that you would like to align with a pixel axis when calling rotate. Would that help you?
>
> --
> Tim Jenness
>
|