Ah, you are doing it wrong.  :-)

The fast (first) index should run over vector index.  The second index
should run across vectors.

Try

        $vecs->mv(1,0)->apply(t_hammer);

Cheers,
Craig

On Feb 2, 2009, at 5:21 PM, David Whysong wrote:

> I'm trying to map a list of 2-vectors in an [N,2] PDL from spherical
> coordinates to an Aitoff projection. Scaling by 180/pi does not seem
> to help. (Also, what are the limits of the two axes? 0-360, or
> -180-180? Does it matter?)
>
> David Whysong
>
>
> On Mon, Feb 2, 2009 at 12:28 PM, Craig DeForest
> <[email protected]> wrote:
>> I think that the documentation for PDL::Transform::Cartography  
>> wrongly
>> implies
>> that the rho,phi coordinates should be in radians. They should be in
>> degrees.
>>
>>
>> For example, if you have PGPLOT installed, try:
>>
>>       use PDL::Transform::Cartography;
>>       $w=pgwin(xw,size=>[12,5],nx=>2);
>>       $e = earth_image();
>>       $w->fits_imag($e);
>>       $w->fits_imag($e->map(t_aitoff,[500,500]));
>>
>> That will show you the Earth image in both Plate Caree coordinates  
>> (in
>> degrees)
>> and Aitoff transformed coordinates (in radians).
>>
>> Cheers,
>> Craig
>


_______________________________________________
Perldl mailing list
[email protected]
http://mailman.jach.hawaii.edu/mailman/listinfo/perldl

Reply via email to