小編給大家分享一下如何用Python代碼實現(xiàn)旋轉(zhuǎn)地球，希望大家閱讀完這篇文章后大所收獲，下面讓我們一起去探討方法吧！

成都創(chuàng)新互聯(lián)公司服務(wù)項目包括福貢網(wǎng)站建設(shè)、福貢網(wǎng)站制作、福貢網(wǎng)頁制作以及福貢網(wǎng)絡(luò)營銷策劃等。多年來，我們專注于互聯(lián)網(wǎng)行業(yè)，利用自身積累的技術(shù)優(yōu)勢、行業(yè)經(jīng)驗、深度合作伙伴關(guān)系等，向廣大中小型企業(yè)、政府機構(gòu)等提供互聯(lián)網(wǎng)行業(yè)的解決方案，福貢網(wǎng)站推廣取得了明顯的社會效益與經(jīng)濟效益。目前，我們服務(wù)的客戶以成都為中心已經(jīng)輻射到福貢省份的部分城市，未來相信會繼續(xù)擴大服務(wù)區(qū)域并繼續(xù)獲得客戶的支持與信任！

用python實現(xiàn)旋轉(zhuǎn)地球

素材

本文只講最簡單也是最常見的正側(cè)面旋轉(zhuǎn)地球。

首先你得有一張360度貼圖，效果如下圖。在網(wǎng)上搜一下，這種圖很多，容易獲得。

它的特點是高度恰好是寬度的一半，因為地球一周有360度，而從南極到北極只有180度。這種360度貼圖，標(biāo)明了每個經(jīng)緯度上地球的表面效果。仔細(xì)觀察會注意到，南北兩極變形劇烈，而中間地帶和我們經(jīng)常見到的地圖比較相似。

我們稱它為經(jīng)緯度貼圖。

基本原理

旋轉(zhuǎn)地球之所以看起來是旋轉(zhuǎn)的，只是因為它上面顯示的內(nèi)容在變。其實起點只是一個圓。我們所做的，只是把經(jīng)緯度貼圖中的像素以合理的方式貼到圓上去。然后再把一幀幀圖片連接為GIF動圖即可。

貼圖的方法是：

1、選擇圓內(nèi)的所有點。

2、計算每一點對應(yīng)的經(jīng)緯度。

3、從經(jīng)緯度貼圖上獲取該經(jīng)緯度對應(yīng)的xy坐標(biāo)

4、獲取像素寫入圖片

5、制作GIF

計算過程是比較簡單的球面幾何，下面結(jié)合代碼講解。

1 獲得圓內(nèi)所有的點

這個比較簡單

  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)

2 計算經(jīng)緯度

首先從平面坐標(biāo)獲得三維坐標(biāo)

球面上的任意一點，有一個規(guī)律是x2 + y2 + z2 = 1

我們就用這個原理，可以獲得圓內(nèi)每個點的z值

根據(jù)x,y,z計算球面的經(jīng)緯度

為了讓計算速度盡量快一些，這里用numpy

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

3 獲取經(jīng)緯度對應(yīng)的xy坐標(biāo)

就是一個簡單的坐標(biāo)映射

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

4 獲得像素點，并寫入圖片

首先，載入經(jīng)緯度貼圖的素材并轉(zhuǎn)為numpy格式

    imgBack= Image.open('EARTH_small.bmp')
    width= imgBack.size[0]
    height= imgBack.size[1]
    imgBack= imgBack.convert('RGBA')
    arrayBack= np.array(imgBack)

注意轉(zhuǎn)換后的XY是反的，取點并直接賦值即可

  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)

5 GIF制作

使用的是imageio庫

    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)

看完了這篇文章，相信你對如何用Python代碼實現(xiàn)旋轉(zhuǎn)地球有了一定的了解，想了解更多相關(guān)知識，歡迎關(guān)注創(chuàng)新互聯(lián)行業(yè)資訊頻道，感謝各位的閱讀！