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()
