यहाँ समारोह है एलेक्स के पद
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-')
और यहाँ से विचार का उपयोग कर उदाहरण है जब आप
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#
मुझे दृष्टिकोण पसंद है - मैंने इसे स्वयं माना था, लेकिन कुछ कारणों से इसे खारिज कर दिया। 1. यह लॉग-लॉग प्लॉट पर काम नहीं करता है (नहीं कि पैच अच्छी तरह से काम करते हैं, लेकिन वे कोशिश करते हैं) 2. इसे साजिश सीमाओं के पूर्व ज्ञान की आवश्यकता होती है (हालांकि इसे चारों ओर काम किया जा सकता है) 3. यह बहुत धीमी है। फिर भी, इस पल के लिए, मैं इसके साथ रोल करूंगा। – keflavich
मुझे लगता है कि इम्शो की उत्पत्ति निचले बाएं कोने पर होनी चाहिए, 'plt.imshow (im.T, cmap = cmap, limit = xr + yr, origin = "lower") '। – ImportanceOfBeingErnest