簡單實現python畫圓功能,python畫圓

來源:互聯網
上載者:User

簡單實現python畫圓功能,python畫圓

本文執行個體為大家分享了python實現畫圓功能的具體代碼,供大家參考,具體內容如下

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

以上就是本文的全部內容,希望對大家的學習有所協助,也希望大家多多支援幫客之家。

聯繫我們

該頁面正文內容均來源於網絡整理,並不代表阿里雲官方的觀點,該頁面所提到的產品和服務也與阿里云無關,如果該頁面內容對您造成了困擾,歡迎寫郵件給我們,收到郵件我們將在5個工作日內處理。

如果您發現本社區中有涉嫌抄襲的內容,歡迎發送郵件至: info-contact@alibabacloud.com 進行舉報並提供相關證據,工作人員會在 5 個工作天內聯絡您,一經查實,本站將立刻刪除涉嫌侵權內容。

A Free Trial That Lets You Build Big!

Start building with 50+ products and up to 12 months usage for Elastic Compute Service

  • Sales Support

    1 on 1 presale consultation

  • After-Sales Support

    24/7 Technical Support 6 Free Tickets per Quarter Faster Response

  • Alibaba Cloud offers highly flexible support services tailored to meet your exact needs.