Python中用代码实现旋转地球！

img= Image.new('RGBA', (300,300), 'black')

  w= img.size[0]

  h= img.size[1]

  pxList=[]

  pyList=[]

  for i in range(w):

      for j in range(h):

          r= math.sqrt((i-w/2)**2+(j-h/2)**2)

          if r<150:

              pxList.append(i)           

              pyList.append(j)

def calcSphereXY2XYZ(px, py, maxHeight, longOffset):

    v0x= np.array(px)

    v0y= np.array(py)

    v03= np.subtract(v0x, maxHeight)

    v04= np.subtract(v0y, maxHeight)

    v1x= np.true_divide(v03, maxHeight)

    v1y= np.true_divide(v04, maxHeight)

    v07= np.power(v1x,2)

    v08= np.power(v1y,2)

    v09= np.add(v07,v08)

    v0a= np.subtract(1,v09)

    v1z= np.power(v0a,1/2)                                  # z

    # print('z:', max(v1z), min(v1z))

    v1lat= np.multiply(v1y, math.pi/2)                      # lat

    v0lon= np.arctan2(v1z, -v1x)                            

    v1lon= np.add(v0lon, longOffset)                       # long

    v2lon= np.fmod(v1lon, math.pi*2)                       # long

    return v2lon, v1lat

def calcShpereLatLong2XY(vlon, vlat, width, height):

    v3x0=np.multiply(vlon, width/2/math.pi)

    v3y0=np.multiply(vlat, height/math.pi)

    v3y1=np.add(v3y0, height/2)

    v3x2=v3x0.astype(np.integer)

    v3y2=v3y1.astype(np.integer)

    return v3x2, v3y2

imgBack= Image.open('EARTH_small.bmp')

width= imgBack.size[0]

height= imgBack.size[1]

imgBack= imgBack.convert('RGBA')

arrayBack= np.array(imgBack)

color= arrayBack[npy, npx]

  for i in range(len(pxList)):

      x= pxList[i]

      y= pyList[i]

      cc=color[i]

      # print(cc)

      cc= tuple(cc)

      img.putpixel((x,y), cc)

frames=[]

for i in range(0, 360, 10):

    a= -i*math.pi/ 180

    img= getPic(a)

    str1= 'temp%03d.png'%i

    img.save(str1)

    im = imageio.imread(str1)

    frames.append(im)

    # img.show()

imageio.mimsave('earth.gif', frames, 'GIF', duration=0.25)