2012-06-09 15 views

risposta

7

Ecco funzione esempio, utilizzando l'idea dal post di Alex

import matplotlib.pyplot as plt,numpy as np 

def gauplot(centers, radiuses, xr=None, yr=None): 
     nx, ny = 1000.,1000. 
     xgrid, ygrid = np.mgrid[xr[0]:xr[1]:(xr[1]-xr[0])/nx,yr[0]:yr[1]:(yr[1]-yr[0])/ny] 
     im = xgrid*0 + np.nan 
     xs = np.array([np.nan]) 
     ys = np.array([np.nan]) 
     fis = np.concatenate((np.linspace(-np.pi,np.pi,100), [np.nan])) 
     cmap = plt.cm.gray 
     cmap.set_bad('white') 
     thresh = 3 
     for curcen,currad in zip(centers,radiuses): 
       curim=(((xgrid-curcen[0])**2+(ygrid-curcen[1])**2)**.5)/currad*thresh 
       im[curim<thresh]=np.exp(-.5*curim**2)[curim<thresh] 
       xs = np.append(xs, curcen[0] + currad * np.cos(fis)) 
       ys = np.append(ys, curcen[1] + currad * np.sin(fis)) 
     plt.imshow(im.T, cmap=cmap, extent=xr+yr) 
     plt.plot(xs, ys, 'r-') 

Ed ecco quello che si ottiene quando si esegue

gauplot([(0,0), (2,3), (5,1), (6, 7), (6.1, 6.1)], [.3,. 4, .5, 1, .4], [-1,10], [-1,10]) 
      #   centers of circles   # radii of circles# 

plot

+0

Mi piace l'approccio - l'avevo considerato anch'io, ma l'ho rifiutato per alcuni motivi. 1. Non funziona sui grafici di log-log (non che le patch funzionino bene, ma ci provano) 2. Richiede una conoscenza preliminare dei limiti della trama (sebbene questo possa essere risolto) 3. è piuttosto lento. Eppure, per il momento, me ne vado. – keflavich

+0

Immagino che l'origine dell'Imshow dovrebbe essere nell'angolo in basso a sinistra, 'plt.imshow (im.T, cmap = cmap, extent = xr + yr, origin =" lower ")'. – ImportanceOfBeingErnest

6

Non penso che matplotlib attualmente supporti i riempimenti sfumati per le patch - vedere this email.

john> Ciao, sto provando a impostare una barra (una serie di rettangoli rappezzata) con un motivo di riempimento anziché solo un colore solido. C'è un modo semplice per farlo in matplotlib?
john> Sto pensando a qualcosa come Qt's QBrush che ha schemi incrociati, verticali, densi, ecc.

Attualmente non vi è alcun supporto per questo - non sarebbe troppo difficile aggiungere per i backend che supportano questo genere di cose. Fondamentalmente, abbiamo bisogno di per specificare l'API per esso e aggiungere il supporto ai backend. Sono stato e volevo aggiungere riempimenti sfumati per le patch (es. Poligoni, rettangoli) e sarebbe bello fare entrambe le cose contemporaneamente.


Invece di utilizzare le patch è possibile creare una rete, calcolare i colori con una funzione quindi utilizzare imshow con interpolazione:

# Taken from http://matplotlib.sourceforge.net/examples/pylab_examples/layer_images.html 

def func3(x,y): 
    return (1- x/2 + x**5 + y**3)*exp(-x**2-y**2) 

# make these smaller to increase the resolution 
dx, dy = 0.05, 0.05 

x = arange(-3.0, 3.0, dx) 
y = arange(-3.0, 3.0, dy) 
X,Y = meshgrid(x, y) 

xmin, xmax, ymin, ymax = amin(x), amax(x), amin(y), amax(y) 
extent = xmin, xmax, ymin, ymax 

fig = plt.figure(frameon=False) 

Z2 = func3(X, Y) 

im2 = imshow(Z2, cmap=cm.jet, alpha=.9, interpolation='bilinear', extent=extent) 

show() 

Ciò comporterà la seguente (ignorare il fondo a scacchi) :

http://matplotlib.sourceforge.net/_images/layer_images.png

+0

e poi mettere una maschera bianca z-intermedio con le patch trasparenti in cui si desidera che il marcatori, per ottenere qualcosa di più come patch sfumate su uno sfondo semplice. – cphlewis

Problemi correlati