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

Home » Public Forums » archive » Re: Fast Filled Contour Plots
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: Fast Filled Contour Plots [message #79886 is a reply to message #79885] Tue, 17 April 2012 15:18 Go to previous message
lecacheux.alain is currently offline  lecacheux.alain
Messages: 325
Registered: January 2008
Senior Member
On 17 avr, 21:33, David Fanning <n...@idlcoyote.com> wrote:
> As I was working on this, Chris Torrence graciously supplied
> some IDL 8.x code to do the same thing. I've included the
> code Chris sent me in the article. Unfortunately, it
> probably does more to expose the limitations of the new
> function graphics system than it does to illuminate the
> features. But, perhaps I have done something wrong. If
> anyone has suggestions (Alain?) for how I can fix this,
> I am all ears.

Your welcome, David.
The problem in Chris Torrence's code comes from the fact that the
longitude range in MAP_CONTINENTS is -180°,180° by default, while you
are using 0°,360° in your example. The trick would be to overplot a
suitable map over your image:
im = IMAGE(scaledData, x, y, RGB_TABLE=rgb, XRANGE=xrange,
YRANGE=yrange, GRID_UNITS='degrees', POSITION=[0.1,0.1,0.9,0.9])
mp = map('EQUIRECTANGULAR', LIMIT=limit, LABEL_POSITION=0,
POSITION=[0.1,0.1,0.9,0.9], /CURRENT)
mc = MapContinents(FILL_COLOR=!COLOR.DARK_GRAY, LIMIT=limit)
Everything is now correctly clipped and displayed, including the grid
overlay.
Note the use of FILL_COLOR instead of COLOR to get continents of the
right color.
There is no bug here, just tricky usage of partly (un)documented
keywords...
alx.
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: speed-up computation of kernel-based "statistics"
Next Topic: Map Function Question

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

Current Time: Sat Oct 11 09:23:00 PDT 2025

Total time taken to generate the page: 3.94639 seconds