Python數(shù)據(jù)可視化：2018年北上廣深空氣質(zhì)量分析（附完整代碼）
就在這周偶然看到一個(gè)學(xué)弟吐槽天津的空氣，不禁想起那段厚德載霧，自強(qiáng)不吸的日子。
無(wú)圖無(wú)真相，下圖為證。
左邊的圖是去年2月份的時(shí)候，這樣的空氣真的難得一見!
右邊的是吐槽以及我個(gè)人第一次買口罩!!!
口罩用的還行，因?yàn)槟莻€(gè)時(shí)候做課設(shè)，經(jīng)常要兩個(gè)校區(qū)跑，基本上空氣不好我就會(huì)帶上。
題目好像是有關(guān)液壓及氣壓的傳動(dòng)系統(tǒng)，手畫A0圖...
這應(yīng)該是快兩年前的事了，時(shí)光飛逝吶。
所以這回先對(duì)2017年天津的空氣質(zhì)量情況進(jìn)行分析，然后再是北上廣深。
01 網(wǎng)頁(yè)分析
網(wǎng)站沒有反爬，所以直接抓取信息就好了。
看見沒有，妥妥的嚴(yán)重污染，2016年12月份買的口罩派上用場(chǎng)啦!
這里簡(jiǎn)單給大家科普一下有關(guān)AQI，PM2.5的知識(shí)。
又是重操舊業(yè)，我的PPT水平還是很水呢~
2 數(shù)據(jù)獲取
獲取代碼如下所示。
- import time
- import requests
- from bs4 import BeautifulSoup
- headers = {
- 'User-Agent':'Mozilla/5.0 (Windows NT 6.1; WOW64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/63.0.3239.132 Safari/537.36'
- }
- for i in range(1, 13):
- time.sleep(5)
- # 把1轉(zhuǎn)換為01
- url = 'http://www.tianqihoubao.com/aqi/tianjin-2017' + str("%02d" % i) + '.html'
- response = requests.get(url=url, headers=headers)
- soup = BeautifulSoup(response.text, 'html.parser')
- tr = soup.find_all('tr')
- # 去除標(biāo)簽欄
- for j in tr[1:]:
- td = j.find_all('td')
- Date = td[0].get_text().strip()
- Quality_grade = td[1].get_text().strip()
- AQI = td[2].get_text().strip()
- AQI_rank = td[3].get_text().strip()
- PM = td[4].get_text()
- with open('air_tianjin_2017.csv', 'a+', encoding='utf-8-sig') as f:
- f.write(Date + ',' + Quality_grade + ',' + AQI + ',' + AQI_rank + ',' + PM + '
- ')
成功獲取數(shù)據(jù)。
3 天津
同樣不上源碼，這里有必要說一波，因?yàn)槲矣X得源碼放上去排版就不好看了...
其次我要秉承以前混跡P圈(PPT)得到的優(yōu)良傳統(tǒng)，熱愛分享，百度云盤你值得擁有。
所以公眾號(hào)回復(fù)天氣。即可獲取全部可視化源碼及相關(guān)文件。
以前天天去下載PPT大神的大作，然后觀摩，可惜的是PPT水平還是那么菜~
01 AQI全年走勢(shì)圖
92.5是年均AQI值，從上面科普知識(shí)里可以知道，2017年天津整體空氣質(zhì)量只能是「良」中的下下等水平，與輕度污染近在咫尺。
02 AQI月均走勢(shì)圖
從月均的走勢(shì)圖就能看出，1月的空氣質(zhì)量最差，8月的空氣質(zhì)量最好，當(dāng)也并不是有多好，充其量也就是個(gè)「良」!
03 AQI季度箱形圖
箱形圖，顯示一組數(shù)據(jù)分散情況資料的統(tǒng)計(jì)圖。
數(shù)據(jù)里有最大值、最小值、中位數(shù)和兩個(gè)四分位數(shù)。
這里可以看出，2017年天津的季度AQI均值差距不是很大。
但是一、二、四季度有明顯的波動(dòng)，空氣質(zhì)量有時(shí)會(huì)變得很差。
04 PM2.5全年走勢(shì)圖
59.87是年均PM2.5值，已經(jīng)遠(yuǎn)超過國(guó)家二級(jí)標(biāo)準(zhǔn)限值35了。
其實(shí)天津給我留下的印象就是天氣經(jīng)常灰蒙蒙，時(shí)常還會(huì)變點(diǎn)顏色，比如黃色~
一年下不了幾次雨，及其干燥。所以那個(gè)最低值11，我猜那時(shí)候估計(jì)是刮大風(fēng)。
05 PM2.5月均走勢(shì)圖
和AQI的走勢(shì)差不多，同樣是1月最高，8月最低。
06 PM2.5季度箱形圖
說實(shí)話，看了這個(gè)圖，我不知道天津的「大哥」及「姐姐」們是如何做到自強(qiáng)不吸的。
基本上四個(gè)季度都超標(biāo)了，一年不超標(biāo)的估計(jì)也就那么幾次。
07 PM2.5指數(shù)日歷圖
日均PM2.5國(guó)家二級(jí)標(biāo)準(zhǔn)為75，從上面的熱力圖看，基本上輕度污染過半了。
另外一月還是重災(zāi)區(qū)，天色黃黃的...
其實(shí)每逢霧霾，基本上就是待宿舍了。而且1月份是考試月，剛好窩宿舍預(yù)習(xí)課本~
08 天津全年空氣質(zhì)量情況
「良」和「輕度污染」占了大頭，「優(yōu)」只能在角落里瑟瑟發(fā)抖，足以說明空氣之差。
不過該上課還是要上課，誰(shuí)叫那時(shí)宿舍和教室離得近(走過去5分鐘不到)。
4 北上廣深
01 北上廣深A(yù)QI全年走勢(shì)圖
北京月均AQI最低也就50左右，看來今年全年差不多都在「優(yōu)」以下了。
不過相比前幾年，京津冀空氣已經(jīng)好了不少(政策)，真的。
上海和廣州差不多，深圳與北京算是鮮明對(duì)比。
02 北上廣深PM2.5全年走勢(shì)圖
北京一如既往的高調(diào)。
03 北上廣深全年空氣質(zhì)量情況
深圳幾乎都是「優(yōu)」和「良」，上海和廣州和上面說的一樣，北京的「優(yōu)」已經(jīng)不少了。
那么你所在的城市空氣質(zhì)量又是如何?
作者:法納斯特，Python愛好者，喜歡爬蟲，數(shù)據(jù)分析以及可視化。
