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