matplotlib-users Mailing List for matplotlib

On 5/10/10 2:21 PM, Michael Hearne wrote:
> I have found a (possible) bug in Basemap - when using basemap.fillcontinents(), I see the chosen continent color only when the map I create includes some ocean. If I am in the interior of a continent (I've tested with North America and Asia), the continent color is white.
> 
Michael: If there are no continent boundaries inside the map projection 
region, basemap does not draw anything (hence you see the axis 
background color). In that case, you should just set the axis 
background color to whatever you wanted the continent color to be.
-Jeff
> A code sample is below. My version information:
> Basemap: 0.99.4
> Matplotlib: 0.99.1.1
> numpy: 1.4.0
> Python: 2.6.4
>
> To replicate my results, please try the following:
> ./maptest.py 37.894507 -121.816406 #map center is somewhere in the Bay Area in California
> ./maptest.py 41.880332 -100.47821 #map center is somewhere in Nebraska
>
> The script creates a file called "output.png" in the calling directory. In the California case, I see the ocean as blue, and the land as a sort of annoying salmon color. In the Nebraska case, I see white with blue denoting the various rivers and lakes in the area.
>
> Am I mis-using the basemap method calls in some way?
>
> Thanks,
>
> Mike
>
> --------------------------------------------------------------------
> #!/usr/bin/env python
>
> import matplotlib
> #use the non-interactive matplotlib setting
> matplotlib.use('agg')
> from mpl_toolkits.basemap import Basemap
> import numpy as np
> from pylab import *
> import sys
>
> clat = float(sys.argv[1])
> clon = float(sys.argv[2])
>
> figwidth = 5.4
>
>
> bounds = (clon-4, clon+4, clat-4, clat+4)
> dx = (bounds[1] - bounds[0])*111191 * np.cos(clat * np.pi/180)
> dy = (bounds[3] - bounds[2])*111191
> aspect = dy/dx
> figheight = aspect * figwidth
>
> fig = figure(figsize=(figwidth,figheight),edgecolor='g',facecolor='g')
> ax1 = fig.add_axes([0,0,1.0,1.0])
> mapcontour = Basemap(llcrnrlon=bounds[0],llcrnrlat=bounds[2],
> urcrnrlon=bounds[1],urcrnrlat=bounds[3],
> resolution='h',projection='merc',lat_ts=clat)
> water_color = [.47,.60,.81]
> mapcontour.drawrivers(color=water_color)
> mapcontour.drawcountries(color='k',linewidth=2.0)
> mapcontour.drawcoastlines()
> mapcontour.fillcontinents(color=[1.0,0.8,0.8],lake_color=water_color)
>
> plt.savefig('output.png')
> --------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------------
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 

I have found a (possible) bug in Basemap - when using basemap.fillcontinents(), I see the chosen continent color only when the map I create includes some ocean. If I am in the interior of a continent (I've tested with North America and Asia), the continent color is white.
A code sample is below. My version information:
Basemap: 0.99.4
Matplotlib: 0.99.1.1
numpy: 1.4.0
Python: 2.6.4
To replicate my results, please try the following:
./maptest.py 37.894507 -121.816406 #map center is somewhere in the Bay Area in California
./maptest.py 41.880332 -100.47821 #map center is somewhere in Nebraska
The script creates a file called "output.png" in the calling directory. In the California case, I see the ocean as blue, and the land as a sort of annoying salmon color. In the Nebraska case, I see white with blue denoting the various rivers and lakes in the area.
Am I mis-using the basemap method calls in some way?
Thanks,
Mike
--------------------------------------------------------------------
#!/usr/bin/env python
import matplotlib
#use the non-interactive matplotlib setting
matplotlib.use('agg')
from mpl_toolkits.basemap import Basemap
import numpy as np
from pylab import *
import sys
clat = float(sys.argv[1])
clon = float(sys.argv[2])
figwidth = 5.4
bounds = (clon-4, clon+4, clat-4, clat+4)
dx = (bounds[1] - bounds[0])*111191 * np.cos(clat * np.pi/180)
dy = (bounds[3] - bounds[2])*111191
aspect = dy/dx
figheight = aspect * figwidth
 
fig = figure(figsize=(figwidth,figheight),edgecolor='g',facecolor='g')
ax1 = fig.add_axes([0,0,1.0,1.0])
mapcontour = Basemap(llcrnrlon=bounds[0],llcrnrlat=bounds[2],
 urcrnrlon=bounds[1],urcrnrlat=bounds[3],
 resolution='h',projection='merc',lat_ts=clat)
water_color = [.47,.60,.81]
mapcontour.drawrivers(color=water_color)
mapcontour.drawcountries(color='k',linewidth=2.0)
mapcontour.drawcoastlines()
mapcontour.fillcontinents(color=[1.0,0.8,0.8],lake_color=water_color)
plt.savefig('output.png') 
--------------------------------------------------------------------

It seems to be failing on the minus signs in the text, due to a recent 
(incorrect) change to declared encoding of the postscript file.
This has been corrected in r8305.
You didn't mention which version of matplotlib you are running... I'm 
not sure this bug would be in the latest released version, only SVN trunk.
Mike
afancy wrote:
> Hi,
>
> Could anybody help me to take a look my following code? When I save it 
> to eps format, the eps file is corrupted. But i can save as pdf format 
> successfully instead.
>
> Regards
> tommy
>
>
> from matplotlib.font_manager import FontProperties
> data = ((16.31, 76.35, 62.62, 135.2),
> (6.4, 12.64, 87, 124.8 ),
> (6.08,60.18,75.2,100.59))
>
> color=('b','r','g','y')
> width = 0.2 
> fig = plt.figure()
> fig.set_size_inches((4.2,4.2))
>
> ax = fig.add_axes((0.2,0.1,0.75,0.
> 8))
> rects=[]
> xticks=[]
> patterns = ('...','//','xx', '\\\\')
>
> for i in range(len(data)):
> row = data[i]
> for j in range(len(row)):
> rect, = ax.bar(i+width*j, row[j], width,fill=False, 
> edgecolor='black', hatch=patterns[j])
> rect.set_hatch(patterns[j])
> rect.set_fill(False)
> if i==0:
> rects.append(rect)
> xticks.append(i+2*width)
> 
> ax.set_xticks(xticks) 
> ax.set_xticklabels(('Q1(s,p1,*)','Q2(s,p2,*)','Q3(s,p3,*)')) 
> ax.set_ylabel('Resp. time (ms)')
> lg=fig.legend(rects, 
> ('3XL-Array','3XL-MP','Jena2(db-based)','Jena2(file-based)'), 
> loc=(0.65,0.7), prop=FontProperties(size=9))
> lg.draw_frame(False)
> draw() 
> fig.savefig("4b.eps", format="eps", bbox_inches='tight') # Failed 
> to save, the eps is corrupted 
> #fig.savefig("4b.pdf", format="pdf", bbox_inches='tight') # pdf 
> format can be save correctly.
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------------
>
> 
> ------------------------------------------------------------------------
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
-- 
Michael Droettboom
Science Software Branch
Operations and Engineering Division
Space Telescope Science Institute
Operated by AURA for NASA

Hi,
Could anybody help me to take a look my following code? When I save it to
eps format, the eps file is corrupted. But i can save as pdf format
successfully instead.
Regards
tommy
from matplotlib.font_manager import FontProperties
data = ((16.31, 76.35, 62.62, 135.2),
 (6.4, 12.64, 87, 124.8 ),
 (6.08,60.18,75.2,100.59))
color=('b','r','g','y')
width = 0.2
fig = plt.figure()
fig.set_size_inches((4.2,4.2))
ax = fig.add_axes((0.2,0.1,0.75,0.
8))
rects=[]
xticks=[]
patterns = ('...','//','xx', '\\\\')
for i in range(len(data)):
 row = data[i]
 for j in range(len(row)):
 rect, = ax.bar(i+width*j, row[j], width,fill=False,
edgecolor='black', hatch=patterns[j])
 rect.set_hatch(patterns[j])
 rect.set_fill(False)
 if i==0:
 rects.append(rect)
 xticks.append(i+2*width)
ax.set_xticks(xticks)
ax.set_xticklabels(('Q1(s,p1,*)','Q2(s,p2,*)','Q3(s,p3,*)'))
ax.set_ylabel('Resp. time (ms)')
lg=fig.legend(rects,
('3XL-Array','3XL-MP','Jena2(db-based)','Jena2(file-based)'),
loc=(0.65,0.7), prop=FontProperties(size=9))
lg.draw_frame(False)
draw()
fig.savefig("4b.eps", format="eps", bbox_inches='tight') # Failed to
save, the eps is corrupted
#fig.savefig("4b.pdf", format="pdf", bbox_inches='tight') # pdf format
can be save correctly.

I just noticed something else now:
When tick labels are explicitly supplied, that may not actually
reflect in the return result of get_ticks. Specifically, I set the
tick labels in a color bar, but get_ticklabels() returns a list of
empty text objects Text(0,0,'').
For regular axes, that seems to work better.
Bug?
Cheers,
Nico
On Sat, May 8, 2010 at 10:36 AM, Nico Schlömer <nic...@gm...> wrote:
> Hi,
>
> say, is there a way to query an axis object whether or not the tick
> labels are drawn?
>
> I looked at -- one the x-axis -- get_xticks( and get_xticklabels()
> which in two different situations would spit out a list of doubles
> (the tick positions) and a list of text entries of the form
> Text(0,0,''). In one plot, though, the tick labels are drawn, not in
> the other one.
>
> Cheers,
> Nico
>

You may also try saving as SVG and then loading in Inkscape and using 
its various "artistic" vector plugins to achieve similar results.
Mike
John Hunter wrote:
> On Fri, May 7, 2010 at 1:09 PM, Eric Firing <ef...@ha...> wrote:
> 
>> On 05/06/2010 08:02 PM, oyster wrote:
>> 
>>> I know matplotlib mimics matlab, which offers scientific look.
>>> but can we use matplotlib to get a skecth(in other words,
>>> hand-drawing) style for bar/pie/etc in none formalist paper? for
>>> example, http://teethgrinder.co.uk/open-flash-chart/gallery-bar-7.php
>>> thanks
>>> 
>> With a considerable effort you could use mpl's capabilities to program
>> this sort of thing, but it will not be added as an option.
>> 
>
> JJ's demo_ribbon_box is pretty close -- you just need to create a base
> image you want stretched to the bar height.
>
> import matplotlib.pyplot as plt
> import numpy as np
> from matplotlib.image import BboxImage
>
> from matplotlib._png import read_png
> import matplotlib.colors
> from matplotlib.cbook import get_sample_data
>
> class RibbonBox(object):
>
> original_image = read_png(get_sample_data("Minduka_Present_Blue_Pack.png",
> asfileobj=False))
> cut_location = 70
> b_and_h = original_image[:,:,2]
> color = original_image[:,:,2] - original_image[:,:,0]
> alpha = original_image[:,:,3]
> nx = original_image.shape[1]
>
> def __init__(self, color):
> rgb = matplotlib.colors.colorConverter.to_rgb(color)
>
> im = np.empty(self.original_image.shape,
> self.original_image.dtype)
>
>
> im[:,:,:3] = self.b_and_h[:,:,np.newaxis]
> im[:,:,:3] -= self.color[:,:,np.newaxis]*(1.-np.array(rgb))
> im[:,:,3] = self.alpha
>
> self.im = im
>
>
> def get_stretched_image(self, stretch_factor):
> stretch_factor = max(stretch_factor, 1)
> ny, nx, nch = self.im.shape
> ny2 = int(ny*stretch_factor)
>
> stretched_image = np.empty((ny2, nx, nch),
> self.im.dtype)
> cut = self.im[self.cut_location,:,:]
> stretched_image[:,:,:] = cut
> stretched_image[:self.cut_location,:,:] = \
> self.im[:self.cut_location,:,:]
> stretched_image[-(ny-self.cut_location):,:,:] = \
> self.im[-(ny-self.cut_location):,:,:]
>
> self._cached_im = stretched_image
> return stretched_image
>
>
>
> class RibbonBoxImage(BboxImage):
> zorder = 1
>
> def __init__(self, bbox, color,
> cmap = None,
> norm = None,
> interpolation=None,
> origin=None,
> filternorm=1,
> filterrad=4.0,
> resample = False,
> **kwargs
> ):
>
> BboxImage.__init__(self, bbox,
> cmap = cmap,
> norm = norm,
> interpolation=interpolation,
> origin=origin,
> filternorm=filternorm,
> filterrad=filterrad,
> resample = resample,
> **kwargs
> )
>
> self._ribbonbox = RibbonBox(color)
> self._cached_ny = None
>
>
> def draw(self, renderer, *args, **kwargs):
>
> bbox = self.get_window_extent(renderer)
> stretch_factor = bbox.height / bbox.width
>
> ny = int(stretch_factor*self._ribbonbox.nx)
> if self._cached_ny != ny:
> arr = self._ribbonbox.get_stretched_image(stretch_factor)
> self.set_array(arr)
> self._cached_ny = ny
>
> BboxImage.draw(self, renderer, *args, **kwargs)
>
>
> if 1:
> from matplotlib.transforms import Bbox, TransformedBbox
> from matplotlib.ticker import ScalarFormatter
>
> fig = plt.gcf()
> fig.clf()
> ax = plt.subplot(111)
>
> years = np.arange(2004, 2009)
> box_colors = [(0.8, 0.2, 0.2),
> (0.2, 0.8, 0.2),
> (0.2, 0.2, 0.8),
> (0.7, 0.5, 0.8),
> (0.3, 0.8, 0.7),
> ]
> heights = np.random.random(years.shape) * 7000 + 3000
>
> fmt = ScalarFormatter(useOffset=False)
> ax.xaxis.set_major_formatter(fmt)
>
> for year, h, bc in zip(years, heights, box_colors):
> bbox0 = Bbox.from_extents(year-0.4, 0., year+0.4, h)
> bbox = TransformedBbox(bbox0, ax.transData)
> rb_patch = RibbonBoxImage(bbox, bc, interpolation="bicubic")
>
> ax.add_artist(rb_patch)
>
> ax.annotate(r"%d" % (int(h/100.)*100),
> (year, h), va="bottom", ha="center")
>
> patch_gradient = BboxImage(ax.bbox,
> interpolation="bicubic",
> zorder=0.1,
> )
> gradient = np.zeros((2, 2, 4), dtype=np.float)
> gradient[:,:,:3] = [1, 1, 0.]
> gradient[:,:,3] = [[0.1, 0.3],[0.3, 0.5]] # alpha channel
> patch_gradient.set_array(gradient)
> ax.add_artist(patch_gradient)
>
>
> ax.set_xlim(years[0]-0.5, years[-1]+0.5)
> ax.set_ylim(0, 10000)
>
> fig.savefig('ribbon_box.png')
> plt.show()
>
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------------
> 
> ------------------------------------------------------------------------
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
-- 
Michael Droettboom
Science Software Branch
Operations and Engineering Division
Space Telescope Science Institute
Operated by AURA for NASA

I believe there's only a handful of characters available as 
sub/superscript in Unicode:
http://www.unicode.org/charts/PDF/U2070.pdf
Mike
Gökhan Sever wrote:
> Hello,
>
> I have tendency to use unicode strings for my simple units e.g. 
> plt.xlabel(u"Concentration (# cm−3)")
>
> Wondering if it is possible to subscript and/or superscript any 
> character in unicode strings? In other words could the following be 
> done via unicode only?
>
> plt.title(r"N$_{ccn} = C s^k$")
>
>
> Thanks.
> 
> Gökhan
> ------------------------------------------------------------------------
>
> ------------------------------------------------------------------------------
>
> 
> ------------------------------------------------------------------------
>
> _______________________________________________
> Matplotlib-users mailing list
> Mat...@li...
> https://lists.sourceforge.net/lists/listinfo/matplotlib-users
> 
-- 
Michael Droettboom
Science Software Branch
Operations and Engineering Division
Space Telescope Science Institute
Operated by AURA for NASA