matplotlib-users Mailing List for matplotlib

Dear folks,

I tried to plot a colored map with non-uniformed discrete colorbar and 
found a few threads related with this. However, I could not find a way 
to apply non-uniformed discrete colorbar over imshow or pcolor. Would 
anybody give me a clue? Thanks in advance.

Kim


Eric Firing wrote:
> James Boyle wrote:
>> I wish to make a color filled plot with the colors defined for  
>> discrete, non-uniform intervals. Something like:
>> 0.0 -0.001  0.001-0.05   0.05-0.2     0.2-0.4      0.4-0.8      0.8-1.0
>>     red              blue             green       magenta      
>> yellow        cyan
>>
>> with the colorbar labeled appropriately.
>> I have seen discussions and solutions for discrete colors but not 
>> for  non-uniform intervals + discrete.
>> The last post I saw regarding this type of issue was august 2005 -  
>> and a solution was not resolved at that time.
>> However, Eric has done a huge amount of work in the intervening time  
>> and a smarter person than myself might have a solution now.
>>
>> Note that I do not wish just to make contours - although that would  
>> be good - but to have a general mapping code that joins allows the  
>> color rmapping to be passed to colorbar.
>> maybe some sub-class of scalarMappable that could work.
>
> This is very easy for contourf, and is illustrated in the second 
> figure made by examples/contourf_demo.py.  For your case above, it 
> would be something like
>
> levs = [0, 0.001, 0.05, 0.2, 0.4, 0.8, 1]
> colors = ['r', 'b', 'g', 'm', 'y', 'c']
> contourf(z, levs, colors=colors)
> colorbar()
>
> Unfortunately, although it *should* be just as easy for imshow or 
> pcolor, it is not at present; it can be done, probably in several 
> ways, but not in such a transparent way.  Attached is a quick attempt 
> at something that might be close to what you need.  The right way to 
> do this is to make some changes and additions to colors.py and 
> colorbar.py; I might get to that in a few days, or, more likely, it 
> might be a few weeks.
>
> Eric
>
>>
>> Thanks for any help.
>>
>> --Jim
> ------------------------------------------------------------------------
>
> -------------------------------------------------------------------------
> This SF.net email is sponsored by DB2 Express
> Download DB2 Express C - the FREE version of DB2 express and take
> control of your XML. No limits. Just data. Click to get it now.
> http://sourceforge.net/powerbar/db2/
> ------------------------------------------------------------------------
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users

Michael Hearne wrote:
> Jeff - 
>
> My data set is actually dynamically generated by a program called 
> ShakeMap.  It's a 2D grid, with an extent usually about 600 kilometers 
> on a side, centered wherever the earthquake happened to be.  The 
> ShakeMap program does not know or care that some of the data may be 
> under water, but for display purposed, I do!  The grid is also in a 
> geographic projection (latitude/longitude coordinates assumed to be 
> cartesian).
Michael:  You'll either have to create your own land-sea mask, or use 
Pierre's method with gdal.
>
> So in this test instance (on a data set near Taiwan), my map width is 
> about 5.91 degrees longitude, and my height is about 5.5 degrees latitude.
>
> If I set xoffset=-0.01*5.91, I get -0.05.  This is not noticeably 
> different than the default.
>
> Is the problem that my dataset is not projected?
What is m.xmax-m.xmin? (m is the basemap instance, xmin and xmax are 
instance variables) That's what you have to use - the distance in map 
projection coordinates, not lat/lon coordinates (although these will be 
the same if you are using projection='cyl').

-Jeff

>
> --Mike
> On Nov 2, 2007, at 1:33 PM, Jeff Whitaker wrote:
>
>> Michael Hearne wrote:
>>> Jeff - I looked at that example file, and I think there's a big 
>>> difference - your etopo base data set is global, and you can plot 
>>> over the data in the oceans by setting the mask on all pixels less 
>>> than zero.
>>>
>>> My dataset (a map of earthquake shaking) is not global, and actually 
>>> has NO missing data.   I think I need a way to "clip" the data by 
>>> the land mask - that is, find all of the pixels that are NOT on 
>>> land, and then mask them off. 
>> Mike:
>>
>> If it's not global, is it just defined for land points?  If so, it 
>> can't be a 2-D grid, so you won't be able to plot it with imshow 
>> anyway.  Can you explain the structure of the data?
>>>
>>> Is there an easy way to do this with matplotlib/basemap tools?
>> Not really.  You'll have to define a sea mask for your grid and use 
>> that the create a masked array.  There is a land-sea mask dataset 
>> included in basemap, but it may not match the resolution of your grid.
>>
>>>
>>> Regarding my other issue - I used my script to test x/y offset 
>>> values: [0.05,0.1,0.5,1.0,10] and couldn't see any difference.  I'd 
>>> be more than happy to provide test output, or debugging information...
>>>
>>> Just to be clear - these offsets are supposed to move the meridian 
>>> and/or parallel labels around with respect to the map edge?  My 
>>> actual goal is to get the labels inside the edge of the map (I tried 
>>> negative numbers to accomplish this, to no effect.)
>> You  need to define an offset as a fraction of the map width - the 
>> numbers you are giving are too small to notice any difference.  As I 
>> said before, try something like -0.01*(m.max-m.min).
>>>
>>> On a positive note, I _can_ make solid lines!
>>
>> Good!
>>
>> -Jeff
>>
>> -- 
>> Jeffrey S. Whitaker         Phone  : (303)497-6313
>> Meteorologist               FAX    : (303)497-6449
>> NOAA/OAR/PSD  R/PSD1        Email  : Jef...@no... 
>> <mailto:Jef...@no...>
>> 325 Broadway                Office : Skaggs Research Cntr 1D-124
>> Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg
>
>
>
>
> ------------------------------------------------------
> Michael Hearne
> mh...@us... <mailto:mh...@us...>
> (303) 273-8620
> USGS National Earthquake Information Center
> 1711 Illinois St. Golden CO 80401
> Senior Software Engineer
> Synergetics, Inc.
> ------------------------------------------------------
>
>


-- 
Jeffrey S. Whitaker         Phone : (303)497-6313
NOAA/OAR/CDC  R/PSD1        FAX   : (303)497-6449
325 Broadway                Boulder, CO, USA 80305-3328