來源：早起Python

作者：讀者投稿

最近幾年，比特幣一直站在風(fēng)口浪尖，一度被追捧為最佳的投資產(chǎn)品，擁護(hù)者認(rèn)為這種加密貨幣是一種類似于黃金的儲(chǔ)值工具，可以對(duì)沖通脹和美元疲軟。其他人則認(rèn)為，比特幣的暴漲只是一個(gè)經(jīng)濟(jì)刺激措施催生的巨大泡沫，并且必將破裂。

一、獲取數(shù)據(jù)

比特幣數(shù)據(jù)很多網(wǎng)站都有，并且也有很多成熟的API，所以取數(shù)據(jù)非常簡(jiǎn)單，直接調(diào)用API接口即可，下面是獲取與寫入數(shù)據(jù)的全部代碼

import requests
import json
import csv
import time
time_stamp = int(time.time())
url = f"https://web-api.coinmarketcap.com/v1/cryptocurrency/ohlcv/historical?convert=USD&slug

=bitcoin&time_end={time_stamp}&time_start=1367107200"
rd = requests.get(url = url) # 返回的數(shù)據(jù)是 JSON 格式，使用 json 模塊解析 co = 

json.loads(rd.content)
list1 = co['data']['quotes']

with open('BTC.csv','w' ,encoding='utf8',newline='') as f:
    csvi = csv.writer(f)
    csv_head = ["date","price","volume"]
    csvi.writerow(csv_head)

    for i in list1:
        quote_date = i["time_open"][:10]
        quote_price = "{:.2f}".format(i["quote"]["USD"]["close"])
        quote_volume = "{:.2f}".format(i["quote"]["USD"]["volume"])
        csvi.writerow([quote_date, quote_price, quote_volume])

執(zhí)行后，當(dāng)前目錄就會(huì)生成BTC.csv數(shù)據(jù)文件

二、動(dòng)態(tài)可視化BTC價(jià)格變化

首先導(dǎo)入需要的包及相關(guān)設(shè)定

import pandas as pd import matplotlib as mpl from matplotlib import cm import numpy 

as np import matplotlib.pyplot as plt import matplotlib.ticker as ticker  import 

matplotlib.animation as animation from IPython.display import HTML from datetime import datetime
plt.rcParams['font.sans-serif'] = ['SimHei']  
plt.rcParams['axes.unicode_minus'] = False  
plt.rc('axes',axisbelow=True) 
mpl.rcParams['animation.embed_limit'] = 2**128

其中兩句plt.rcParams[]是用來設(shè)置顯示中文的

plt.rc('axes',axisbelow=True)的作用是設(shè)置外觀要求，即坐標(biāo)軸置底。

mpl.rcParams['animation.embed_limit'] = 2**128這句是為了生成動(dòng)畫而用的，由于動(dòng)畫默認(rèn)的最大體積為20971520.字節(jié)。如果需要調(diào)整生成的動(dòng)畫最大體積，需要更改這個(gè)參數(shù)。

接下來數(shù)據(jù)并利用查看前5行與后5行

從表格初窺可以得知，13年初的價(jià)格在100美元左右，而到如今21年價(jià)格已經(jīng)飛漲到5萬左右了。具體在哪段時(shí)間飛漲如此之快呢，我們通過動(dòng)態(tài)面積可視化來探索。

可視化之前，需要對(duì)數(shù)據(jù)進(jìn)行處理，由于我們?cè)镜臄?shù)據(jù)是這樣的

是csv格式，且Date字段是字符串類型，而在Python中運(yùn)用matplotlib畫時(shí)間序列圖都需要datetime時(shí)間戳格式才美觀，所以我們運(yùn)用了如下代碼進(jìn)行轉(zhuǎn)換

df = pd.read_csv('BTC.csv')
df['date']=[datetime.strptime(d, '%Y/%m/%d').date() for d in df['date']]

下面制作靜態(tài)面積圖，使用單色填充的話，可用如下代碼

Span=180 N_Span=0 df_temp=df.loc[N_Span*Span:(N_Span+1)*Span,:]
df_temp.head(5)
fig =plt.figure(figsize=(6,4), dpi=100)
plt.subplots_adjust(top=1,bottom=0,left=0,right=0.9,hspace=0,wspace=0)
plt.fill_between(df_temp.date.values, y1=df_temp.price.values, y2=0,alpha=0.75, facecolor='r', 

linewidth=1,edgecolor ='none',zorder=1)
plt.plot(df_temp.date, df_temp.price, color='k',zorder=2)
plt.scatter(df_temp.date.values[-1], df_temp.price.values[-1], color='white',s=150,edgecolor ='k',

linewidth=2,zorder=3)
plt.text(df_temp.date.values[-1], df_temp.price.values[-1]*1.18,s=np.round(df_temp.price.values[-1],1),

size=10,ha='center', va='top')
plt.ylim(0, df_temp.price.max()*1.68)
plt.xticks(ticks=df_temp.date.values[0:Span+1:30],labels=df_temp.date.values[0:Span+1:30],rotation=0)
plt.margins(x=0.01)
ax = plt.gca()#獲取邊框 ax.spines['top'].set_color('none')    # 設(shè)置上‘脊梁’為無色 

ax.spines['right'].set_color('none')  # 設(shè)置上‘脊梁’為無色 ax.spines['left'].set_color('none')   

# 設(shè)置上‘脊梁’為無色 plt.grid(axis="y",c=(217/256,217/256,217/256),linewidth=1)         

#設(shè)置網(wǎng)格線    plt.show()

其中Span設(shè)定的是多少天的價(jià)格，這里我們使用200天。N_Span代表權(quán)重；

df_temp=df.loc[N_Span*Span:(N_Span+1)*Span,:]代表的是選擇到179行為止的數(shù)據(jù)，即180天。

plt.fill_between()是使用單色--紅色填充

得到如下效果

但是一個(gè)顏色填充總感覺不夠好看，所以下面使用漸變色填充，使用plt.bar()函數(shù)實(shí)現(xiàn)Spectral_r顏色映射。代碼如下：

Span_Date =180  
Num_Date =360  #終止日期 df_temp=df.loc[Num_Date-Span_Date: Num_Date,:]  

#選擇從Num_Date-Span_Date開始到Num_Date的180天的數(shù)據(jù) colors = 

cm.Spectral_r(df_temp.price / float(max(df_temp.price)))
fig =plt.figure(figsize=(6,4), dpi=100)
plt.subplots_adjust(top=1,bottom=0,left=0,right=0.9,hspace=0,wspace=0)
plt.bar(df_temp.date.values,df_temp.price.values,color=colors,width=1,align="center",zorder=1)
plt.plot(df_temp.date, df_temp.price, color='k',zorder=2)
plt.scatter(df_temp.date.values[-1], df_temp.price.values[-1], color='white',s=150,edgecolor ='k',linewidth=2,zorder=3)
plt.text(df_temp.date.values[-1], df_temp.price.values[-1]*1.18,s=np.round(df_temp.price.values[-1],1),

size=10,ha='center', va='top')
plt.ylim(0, df_temp.price.max()*1.68)
plt.xticks(ticks=df_temp.date.values[0: Span_Date +1:30],labels=df_temp.date.values[0: Span_Date +1:30],rotation=0)
plt.margins(x=0.01)
ax = plt.gca()#獲取邊框 ax.spines['top'].set_color('none')   # 設(shè)置上‘脊梁’為無色 ax.spines['right'].set_color('none')  

# 設(shè)置上‘脊梁’為無色 ax.spines['left'].set_color('none')   # 設(shè)置上‘脊梁’為無色 plt.grid(axis="y",

c=(217/256,217/256,217/256),linewidth=1)   #設(shè)置網(wǎng)格線    plt.show()

這里的數(shù)據(jù)篩選有稍許不同，其中Span_Date設(shè)置初始時(shí)間，這里設(shè)置為180即從起始日開始算的180天.

Num_Date設(shè)置的是終止時(shí)間。

df_temp=df.loc[Num_Date-Span_Date: Num_Date,:]則是用loc函數(shù)篩選從180天到終止日期的數(shù)據(jù)。

效果如下：

三、動(dòng)態(tài)可視化

最后，我們來將這幅圖動(dòng)起來，先將剛剛的繪圖部分封裝

def draw_areachart(Num_Date):
    Span_Date=180
    ax.clear()
    if Num_Date<Span_Date:         df_temp=df.loc[0:Num_Date,:]         df_span=df.loc[0:Span_Date,:]         

colors = cm.Spectral_r(df_span.price.values / float(max(df_span.price.values)))
        plt.bar(df_temp.date.values,df_temp.price.values,color=colors,width=1.5,align="center",zorder=1)         

plt.plot(df_temp.date, df_temp.price, color='k',zorder=2)         plt.scatter(df_temp.date.values[-1], 

df_temp.price.values[-1], color='white',s=150,edgecolor ='k',linewidth=2,zorder=3)         

plt.text(df_temp.date.values[-1], df_temp.price.values[-1]*1.18,s=np.round(df_temp.price.values[-1],1), 

size=10,ha='center', va='top')
        plt.ylim(0, df_span.price.max()*1.68)
        plt.xlim(df_span.date.values[0], df_span.date.values[-1])
        plt.xticks(ticks=df_span.date.values[0:Span_Date+1:30],labels=df_span.date.values[0:Span_Date+1:30],

rotation=0,fontsize=9)     

else:         df_temp=df.loc[Num_Date-Span_Date:Num_Date,:]         colors = cm.Spectral_r(df_temp.price / 

float(max(df_temp.price)))
        plt.bar(df_temp.date.values[:-2],df_temp.price.values[:-2],color=colors[:-2],width=1.5,align="center",zorder=1)         

plt.plot(df_temp.date[:-2], df_temp.price[:-2], color='k',zorder=2)         plt.scatter(df_temp.date.values[-4], 

df_temp.price.values[-4], color='white',s=150,edgecolor ='k',linewidth=2,zorder=3)         

plt.text(df_temp.date.values[-1], df_temp.price.values[-1]*1.18,s=np.round(df_temp.price.values[-1],1),

 size=10,ha='center', va='top')
        plt.ylim(0, df_temp.price.max()*1.68)
        plt.xlim(df_temp.date.values[0], df_temp.date.values[-1])
        plt.xticks(ticks=df_temp.date.values[0:Span_Date+1:30],labels=df_temp.date.values[0:Span_Date+1:30],rotation=0,fontsize=9)    

    plt.margins(x=0.2)     ax.spines['top'].set_color('none')  # 設(shè)置上‘脊梁’為紅色
    ax.spines['right'].set_color('none')  # 設(shè)置上‘脊梁’為無色
    ax.spines['left'].set_color('none')  # 設(shè)置上‘脊梁’為無色
    plt.grid(axis="y",c=(217/256,217/256,217/256),linewidth=1)         #設(shè)置網(wǎng)格線  
    plt.text(0.01, 0.95,"BTC平均價(jià)格($)",transform=ax.transAxes, size=10, weight='light', ha='left')
    ax.text(-0.07, 1.03, '2013年到2021年的比特幣BTC價(jià)格變化情況',transform=ax.transAxes, size=17, weight='light', 

ha='left') fig, ax = plt.subplots(figsize=(6,4), dpi=100) 

plt.subplots_adjust(top=1,bottom=0.1,left=0.1,right=0.9,hspace=0,wspace=0)   draw_areachart(150)

之后使用matplotlib包的animation.FuncAnimation()函數(shù)，之后調(diào)用上述編寫的draw_areachart(Num_Date)函數(shù)。

其中輸入的參數(shù)Num_Date是如靜態(tài)可視化中提及的日期作用一樣，賦值為np.arange(0,df.shape[0],1)。

最后使用Ipython包的HTML()函數(shù)將動(dòng)畫轉(zhuǎn)換成動(dòng)畫頁(yè)面的形式演示。代碼如下：

import matplotlib.animation as animation
from IPython.display import HTML
fig, ax = plt.subplots(figsize=(6,4), dpi=100)
plt.subplots_adjust(left=0.12, right=0.98, top=0.85, bottom=0.1,hspace=0,wspace=0)  
animator = animation.FuncAnimation(fig, draw_areachart, frames=np.arange(0,df.shape[0],1),

interval=100) HTML(animator.to_jshtml()) 

函數(shù)FuncAnimation(fig,func,frames,init_func,interval,blit)是繪制動(dòng)圖函數(shù)。其參數(shù)如下：

“

fig 表示繪制動(dòng)圖的畫布名稱(figure)；func為自定義繪圖函數(shù)，如draw_barchart()函數(shù)；frames為動(dòng)畫長(zhǎng)度，一次循環(huán)包含的幀數(shù)，在函數(shù)運(yùn)行時(shí)，其值會(huì)傳遞給函數(shù)draw_barchart (year)的形參“year”；init_func為自定義開始幀可省略；interval表示更新頻率，計(jì)量單位為ms；blit表示選擇更新所有點(diǎn)，還是僅更新產(chǎn)生變化的點(diǎn)，應(yīng)選擇為True，但mac電腦用戶應(yīng)選擇False，否則無法顯示。

”

最后效果就是這樣

可以看到在過去的一年中，由于機(jī)構(gòu)的興趣日益增加，比特幣上漲超過了6倍，最高突破58000美元/枚，當(dāng)然可以看到跌起來也是非?？植赖?，關(guān)于比特幣，你怎么看？