On Wednesday, January 16, 2013 12:11:50 AM UTC, tortoise...@gmail.com wrote:
> On Tuesday, 15 January 2013 22:14:04 UTC, Coyote wrote:
>
>> On Tuesday, January 15, 2013 1:56:23 PM UTC-7, lj...@fsmail.net wrote:
>
>>
>
>>> I am trying to do a contour plot of a latitude-longitude array, but I keep getting gaps in my plot, as shown in the following image: http://www.physics.open.ac.uk/~lsteele/contour.png. I'm stumped as to the reason, as it's not as if it's missing contours at the very end ranges of the contour values - it's somewhere in between. Anyhow, I'm sure the knowledgable folks on here will know!
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> I have defined a blue-white-red colour scheme and loaded it in with tvlct, but the problem occurs with other colour schemes too. I am calculating my contour levels and colours as follows (with ncont the number of contours I'm plotting):
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> minv = floor(min(array))
>
>>
>
>>>
>
>>
>
>>> maxv = ceil(max(array))
>
>>
>
>>>
>
>>
>
>>> maxval = max([abs(array),abs(array)])
>
>>
>
>>>
>
>>
>
>>> clevels = -maxval + findgen(ncont)/(ncont-1)*2*maxval
>
>>
>
>>>
>
>>
>
>>> ccol = floor(findgen(ncont)/(ncont-1)*250)
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> I am then plotting my contour plot as:
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> device, decomposed=0
>
>>
>
>>>
>
>>
>
>>> bwr_colour, red_vec, grn_vec, blu_vec
>
>>
>
>>>
>
>>
>
>>> tvlct, red_vec, grn_vec, blu_vec
>
>>
>
>>>
>
>>
>
>>> contour, array, lon, lat, $
>
>>
>
>>>
>
>>
>
>>> levels=clevels, c_colors=ccol, /fill, $
>
>>
>
>>>
>
>>
>
>>> background=cgColor('white'), color=cgColor('black'), $
>
>>
>
>>>
>
>>
>
>>> xrange=[-180,180], yrange=[-90,90], xstyle=1, ystyle=1, $
>
>>
>
>>>
>
>>
>
>>> yticks=6, ytickv=[-90,-60,-30,0,30,60,90], $
>
>>
>
>>>
>
>>
>
>>> xticks=6, xtickv=[-180,-120,-60,0,60,120,180], $
>
>>
>
>>>
>
>>
>
>>> xtitle='Longitude', ytitle='Latitude', charsize=1.5, $
>
>>
>
>>>
>
>>
>
>>> title=plotsym+' (all motions) / g m!u-1!n s!u-1!n'
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>>
>
>>
>
>>> If anyone can offer any suggestions then I'll be very grateful, as it's bugging me now!
>
>>
>
>>
>
>>
>
>> Just seeing "lon", "lat", and "contour" causes me to think that one of your problems is that you should be using the CELL_FILL keyword, rather than the FILL keyword.
>
>>
>
>>
>
>>
>
>> http://www.idlcoyote.com/color_tips/fill_colors.html
>
>>
>
>>
>
>>
>
>> If that doesn't work, I'd try cgCONTOUR. That routine will fill up the "normal" hole Contour creates.
>
>>
>
>>
>
>>
>
>> http://www.idlcoyote.com/tips/contour_hole.php
>
>>
>
>>
>
>>
>
>> Cheers,
>
>>
>
>>
>
>>
>
>> David
>
>
>
> Thanks David - I'll give both a try tomorrow and see what happens!
Hmm, well neither cell_fill nor cgcontour fix the problem. I'll have to have another fiddle and see what I can do, but it seems like it should work!
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
