comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Over-PLOT() on MAP()?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Over-PLOT() on MAP()? [message #73921 is a reply to message #73920] Sat, 11 December 2010 16:55 Go to previous messageGo to previous message
David Fanning is currently offline  David Fanning
Messages: 11724
Registered: August 2001
Senior Member
Paulo Penteado writes:

> plot() seems to have inherited the same limitation as iplot with
> respect to overplotting on maps:
>
> ttp://groups.google.com/group/comp.lang.idl-pvwave/browse_th read/thread/7dad5b4234801dea
>
> You need to use its map projection to convert the coordinates to
> Cartesian:
>
> m=map('sinusoidal')
> lats=[60.,30.,0.,-60.,-45.,60.]
> lons=[-60.,0.,120.,30.,-60.,-60.]
> xy=map_proj_forward(lons,lats,map_structure=m.map_structure)
> p=plot(xy,/over,color='red')

I'm not sure this is a "limitation" so much as just
a poorly documented feature. I think it is absolutely
correct that the "data space" of a map projection should
be in projected XY coordinates, rather than in lat/lon
space. It is just about the only thing about the new
graphics that actually makes a lot of sense to me. :-)

Cheers,

David



--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Sepore ma de ni thui. ("Perhaps thou speakest truth.")
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Butterworth Band-Pass Filter
Next Topic: 3D vector rotation to the Z axis

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Fri Oct 10 15:46:55 PDT 2025

Total time taken to generate the page: 1.59955 seconds