臺風天要帶多少錢出門？
答案是：四千萬，因為臺風天沒四(事)千萬不要出門
這不最近又到了臺風頻發的季節，作為沒有四千萬的沿海人，多關注一下天氣預報還是有必要的，
那么今天呢，就給大家講一個：用Python繪制臺風路徑軌跡圖的可視化案例，
此案例中用到的核心繪圖庫是：plotly，plotly庫是一個可互動，基于瀏覽器的繪圖庫，功能非常強大, 可以在線繪制常見的條形圖、散點圖、餅圖、直方圖等，
涉及到的內容：
requests+json – 網頁資料爬取
openpyxl – 保存資料至Excel
os.walk() – 檔案路徑遍歷
pandas – 表格資料處理
plotly – 地圖資料可視化
1. 資料爬取
1.1 網頁分析
資料來源：中國天氣臺風網(http://typhoon.weather.com.cn/index.shtml)
以2020年為例，打開網址查看2020年臺風串列：
去年有記錄的臺風一共有23個，編號2001~2023，
1.2 查看網頁臺風資料
F12除錯模式，以編號2001為例：
可以看到軌跡資料在頁面回傳的getData字串中，內部資料為json格式，
資料包含：臺風名稱、時間（間隔1小時）、經度、緯度、風速、移動方向、移動速度等等，
對于瀏覽器除錯模式還不熟悉的同學，可以看之前的文章：爬蟲必備工具，掌握它就解決了一半的問題
1.3 資料爬取
單個臺風軌跡爬取，需要輸入年份和臺風編號：
Python學習交流Q群：906715085### def get_tp_data(year, code): tim = int(time.time()) url = f'http://d1.weather.com.cn/typhoon/typhoon_data/{year}/{code}.json?callback=getData&_={tim}' headers = { 'User-Agent': 'Mozilla/5.0 (Windows NT 10.0; WOW64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/55.0.2883.87 Safari/537.36' } r = requests.get(url,headers=headers) r.raise_for_status() r.encoding = r.apparent_encoding datas = json.loads(r.text[8:-1]) tp_name = datas['typhoon'][2] detail_datas = datas['typhoon'][8] # 臺風資訊 allinfo = [] for data in detail_datas: # 時間 tp_time = data[1] # 經度 tp_lon = data[4] # 緯度 tp_lat = data[5] # 中心氣壓 central_pressure = data[6] # 風速 tp_speed = data[7] # 移向 tp_direction = data[8] # 去除臺風登陸時的空行 if tp_direction == None : continue # 移速 move_speed = data[9] tp_info = [f'{code}-' + tp_name,tp_time,tp_lon,tp_lat,central_pressure,tp_speed,tp_direction,move_speed] allinfo.append(tp_info)
注意操作：去除臺風登陸時的空行資料(該資料無經、緯度資訊)，這一步也可以在后續資料處理中進行，
以前10個臺風為例（可根據需要調整）：
#臺風年份、編號 year = 2020 #資料爬取 for code in range(2001,2011): get_tp_data(year, code) pass
關于 requests 庫的使用也可以參考之前文章：讓你的爬蟲開發效率提升8倍
1.4 資料保存
使用openpyxl庫保存Excel檔案：
def insert2excel(filepath, allinfo): try: # 表頭 tableTitle = ['名稱','時間','經度','緯度','中心氣壓(hPa)','風速(m/s)','移向','移速(m/s)'] wb = Workbook() ws = wb.active ws.title = 'sheet1' ws.append(tableTitle) for info in allinfo: ws.append(info) wb.save(filepath) return True except: return False
結果：
注意新建data檔案夾(路徑可自行調整)：
2. 地理軌跡
2.1 mapbox網站Token獲取
網站：https://account.mapbox.com/
需要郵箱注冊：
2.2 匯入plotly相關包
匯入express模塊：
import plotly.express as px
如果環境沒有plotly包，直接通過pip安裝即可：
pip install plotly
2.3 單個臺風軌跡
代碼：
def trace_point_shows(): df = pd.read_excel('./data/2020_2001_黃蜂軌跡資料.xlsx') print(df) token = '你的token' fig = px.scatter_mapbox(df, hover_data=['時間'], lon = '經度', lat = '緯度', color = '風速(m/s)', hover_name = '名稱', size_max = 14, color_continuous_scale = px.colors.carto.Temps ) fig.update_layout(mapbox = {'accesstoken': token, # 官網注冊token 'center':{'lon':121.54,'lat':25.00}, # 地圖中心 'zoom': 8, 'style': 'dark', # 顯示地圖型別 }, margin = {'l': 1, 'r': 1, 't': 1, 'b': 1}) # 地圖邊界 fig.write_html('trace_point_shows-1.html')
效果：
地圖style可以是：basic, streets, outdoors, light, dark, satellite, satellite-streets 中的任意一種，默認為 ‘light’，
2.4 多個臺風軌跡
os.walk()遍歷檔案夾中所有軌跡表格檔案：
filenames = [] for root, dirs, files in os.walk('./data'): for name in files: filenames.append(f'{root}/{name}') 代碼： def trace_point_shows(filenames): df = pd.DataFrame() lsts = [] for filename in filenames: df0 = pd.read_excel(filename) lsts.append(df0) df = pd.concat(lsts) token = '你的token' fig = px.scatter_mapbox(df, hover_data=['時間'], lon = '經度', lat = '緯度', color = '風速(m/s)', hover_name = '名稱', size_max = 14, color_continuous_scale = px.colors.carto.Temps ) fig.update_layout(mapbox = {'accesstoken': token, # 官網注冊token 'center':{'lon':121.54,'lat':25.00}, # 地圖中心 'zoom': 8, 'style': 'dark', # 顯示地圖型別 }, margin = {'l': 1, 'r': 1, 't': 1, 'b': 1}) # 地圖邊界 fig.write_html('trace_point_shows.html')
效果：
如此就實作一個可以互動查看資料的臺風路徑軌跡圖，感興趣的同學可以自己也動手試一試，這一篇到這就莫有了，下一篇見啦！！！
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/456995.html
標籤：Python
下一篇：Java多執行緒與執行緒池技術