import numpy as np import matplotlib.pyplot as plt from matplotlib.patches import Polygon import matplotlib.patches as mpatches  fig = plt.figure(figsize = (16,8)) ax = fig.gca() ax.set_xlim(-5,18) ax.set_ylim(-5,8) alpha = np.linspace(-np.pi,np.pi,30)  t = np.linspace(-np.pi,np.pi,100) def ppp(nc , na , c , sx , sy):  beta = alpha[na]  gam = alpha[na+1]  x = [] ; y = [] ;  r = 0.5*nc  x.append(r*np.cos(beta)) ; y.append(r*np.sin(beta)) ;  tt = np.linspace(beta,gam,10)  for i in range(len(tt)):   x.append(r*np.cos(tt[i]))   y.append(r*np.sin(tt[i]))  r-=0.5  x.append(r*np.cos(gam)) ; y.append(r*np.sin(gam))  tt = np.linspace(gam,beta,10)  for i in range(len(tt)):   x.append(r*np.cos(tt[i])) ; y.append(r*np.sin(tt[i]))  x = np.array(x) ; y = np.array(y)  verts = list(zip(x+sx,y+sy))  poly = Polygon(verts , facecolor=c , edgecolor=c)  ax.add_patch(poly)  #ppp(2,0,'r',0,0) def get_coordiate(x , y):  beta = (alpha[y]+alpha[y+1])/2  pr = 0.5*(x+x-1)/2  sx = pr*np.cos(beta) ; sy = pr*np.sin(beta)  return sx,sy def pp_arrow(x,y):  print x , y  sx , sy = get_coordiate(x, y)  for i in range(-1,2):   for j in range(-1,2):    if i == 0 and j == 0:     continue    ex , ey = get_coordiate(x+i, y+j)    ax.arrow(sx,sy,ex-sx,ey-sy,head_width=0.05,head_length=0.1,fc='r',ec='r')     def pxy(cx,cy):  for i in range(10):   plt.plot(0.5*i*np.cos(t)+cx,0.5*i*np.sin(t)+cy)   for i in range(len(alpha)):   plt.plot((0.5*np.cos(alpha)+cx,4.5*np.cos(alpha)+cx),(0.5*np.sin(alpha)+cy,4.5*np.sin(alpha)+cy))  cx = 13 ; cy =0 pxy(0,0) pxy(13,0)  mr = np.zeros((10,30)) mb = np.zeros((10,30))  for i in range(3,9):  for j in range(1,27):   ttt = np.random.random()   if ttt >= 0.98 and ttt <= 1:    ppp(i,j,'r',0,0)    mr[i][j] = 1    continue   if ttt < 0.8:    continue   else:    ppp(i,j,'b',0,0)    mb[i][j] = 1  for i in range(10):  for j in range(30):   if mr[i][j] == 0:    continue   pp_arrow(i, j)   ma = 0 ;posx = 1 ; posy = 1   for k in range(-1,2):    for s in range(-1,2):     mb[i+k][j+s]=0     ttt = np.random.random()     if ttt > ma:      ma = ttt      posx = k ; posy = s   ppp(i+posx , j+posy , 'r' , 13 , 0)   ppp(i,j,'g',13,0)  for i in range(10):  for j in range(30):   if mb[i][j] == 0:    continue   ppp(i , j , 'b' , 13 , 0)  arrow = mpatches.Arrow(5, 0       ,3,0,width = 1,color = 'r') ax.add_patch(arrow) ax.text(6.2,0.5,'after one/nunit time',ha='center',va='center',color='y')  xx = [4,4,4] yy = [7.5,6.5,5.5] c = ['b' , 'r' , 'g'] s = ['debris','satellite' , 'the location of satellite one unit time ago'] for i in range(3):  print xx[i] , yy[i]  rect = mpatches.Rectangle((xx[i],yy[i]),0.5,0.3,ec = c[i],fc=c[i] )  ax.text(xx[i]+0.7,yy[i],s[i])  ax.add_patch(rect) plt.axis('off') plt.show()