40

我有一组要显示为散点图的数据。我希望将每个点绘制为 size 的正方形dx

          x = [0.5,0.1,0.3]
          y = [0.2,0.7,0.8]
          z = [10.,15.,12.]
          dx = [0.05,0.2,0.1]

          scatter(x,y,c=z,s=dx,marker='s')

问题是sscatter 函数读取的大小以点^2为单位。我想要的是用面积为 dx^2 的正方形表示每个点,其中该面积以“真实”单位表示,即绘图单位。我希望你能明白这一点。

我还有一个问题。scatter 函数绘制带有黑色边框的标记,如何删除此选项并且根本没有边框?

4

3 回答 3

43

用户数据坐标系转换为显示坐标系。

并使用 edgecolors='none' 绘制没有轮廓的面。

import numpy as np

fig = figure()
ax = fig.add_subplot(111)
dx_in_points = np.diff(ax.transData.transform(zip([0]*len(dx), dx))) 
scatter(x,y,c=z,s=dx_in_points**2,marker='s', edgecolors='none')
于 2012-01-31T15:58:48.097 回答
21

我认为我们可以通过一系列补丁做得更好。根据文件:

这 (PatchCollection) 可以更轻松地将颜色映射分配给不同的补丁集合。

这也可以提高绘图速度,因为 PatchCollection 将比大量补丁更快地绘制。

假设您要在数据单元中绘制具有给定半径的圆的散点图:

def circles(x, y, s, c='b', vmin=None, vmax=None, **kwargs):
    """
    Make a scatter of circles plot of x vs y, where x and y are sequence 
    like objects of the same lengths. The size of circles are in data scale.

    Parameters
    ----------
    x,y : scalar or array_like, shape (n, )
        Input data
    s : scalar or array_like, shape (n, ) 
        Radius of circle in data unit.
    c : color or sequence of color, optional, default : 'b'
        `c` can be a single color format string, or a sequence of color
        specifications of length `N`, or a sequence of `N` numbers to be
        mapped to colors using the `cmap` and `norm` specified via kwargs.
        Note that `c` should not be a single numeric RGB or RGBA sequence 
        because that is indistinguishable from an array of values
        to be colormapped. (If you insist, use `color` instead.)  
        `c` can be a 2-D array in which the rows are RGB or RGBA, however. 
    vmin, vmax : scalar, optional, default: None
        `vmin` and `vmax` are used in conjunction with `norm` to normalize
        luminance data.  If either are `None`, the min and max of the
        color array is used.
    kwargs : `~matplotlib.collections.Collection` properties
        Eg. alpha, edgecolor(ec), facecolor(fc), linewidth(lw), linestyle(ls), 
        norm, cmap, transform, etc.

    Returns
    -------
    paths : `~matplotlib.collections.PathCollection`

    Examples
    --------
    a = np.arange(11)
    circles(a, a, a*0.2, c=a, alpha=0.5, edgecolor='none')
    plt.colorbar()

    License
    --------
    This code is under [The BSD 3-Clause License]
    (http://opensource.org/licenses/BSD-3-Clause)
    """
    import numpy as np
    import matplotlib.pyplot as plt
    from matplotlib.patches import Circle
    from matplotlib.collections import PatchCollection

    if np.isscalar(c):
        kwargs.setdefault('color', c)
        c = None
    if 'fc' in kwargs: kwargs.setdefault('facecolor', kwargs.pop('fc'))
    if 'ec' in kwargs: kwargs.setdefault('edgecolor', kwargs.pop('ec'))
    if 'ls' in kwargs: kwargs.setdefault('linestyle', kwargs.pop('ls'))
    if 'lw' in kwargs: kwargs.setdefault('linewidth', kwargs.pop('lw'))

    patches = [Circle((x_, y_), s_) for x_, y_, s_ in np.broadcast(x, y, s)]
    collection = PatchCollection(patches, **kwargs)
    if c is not None:
        collection.set_array(np.asarray(c))
        collection.set_clim(vmin, vmax)

    ax = plt.gca()
    ax.add_collection(collection)
    ax.autoscale_view()
    if c is not None:
        plt.sci(collection)
    return collection

函数的所有参数和关键字(除了marker)都scatter将以类似的方式工作。我写了一个要点,包括椭圆正方形/矩形。如果你想要其他形状的集合,你可以自己修改它。

如果要绘制颜色条,只需运行colorbar()或将返回的集合对象传递给colorbar函数。

一个例子:

from pylab import *
figure(figsize=(6,4))
ax = subplot(aspect='equal')

#plot a set of circle
a = arange(11)
out = circles(a, a, a*0.2, c=a, alpha=0.5, ec='none')
colorbar()

#plot one circle (the lower-right one)
circles(1, 0, 0.4, 'r', ls='--', lw=5, fc='none', transform=ax.transAxes)

xlim(0,10)
ylim(0,10)

输出:

示例图

于 2014-07-04T05:38:00.790 回答
21

如果您想要随图形大小调整大小的标记,您可以使用补丁:

from matplotlib import pyplot as plt
from matplotlib.patches import Rectangle

x = [0.5, 0.1, 0.3]
y = [0.2 ,0.7, 0.8]
z = [10, 15, 12]
dx = [0.05, 0.2, 0.1]

cmap = plt.cm.hot
fig = plt.figure()
ax = fig.add_subplot(111, aspect='equal')

for x, y, c, h in zip(x, y, z, dx):
    ax.add_artist(Rectangle(xy=(x, y),
                  color=cmap(c**2),        # I did c**2 to get nice colors from your numbers
                  width=h, height=h))      # Gives a square of area h*h

plt.show()

在此处输入图像描述

注意:

  1. 正方形不以 为中心(x,y)。x,y 实际上是左下方正方形的坐标。我让它以这种方式简化我的代码。你应该使用(x + dx/2, y + dx/2).
  2. 颜色来自热颜色图。我用 z**2 给颜色。您还应该根据自己的需要进行调整

最后回答你的第二个问题。您可以使用关键字参数 edgecoloredgecolors. 这些分别是 matplotlib 颜色参数或 rgba 元组序列。如果将参数设置为“无”,则不绘制边框。

于 2012-01-31T15:56:09.277 回答