北京全年天气状况分析

看到天气网上有国内城市一年的天气历史数据,想以此为数据源练习一下,于是就有了这个项目。今天在此简单介绍一下实现思路和最终效果。

 

用到的相关库包括:

  • requests

  • bs4

  • pandas

  • matplotlib

  • seaborn

  • pyecharts

 

分析数据源

 

天气网的历史天气预报查询页面(http://lishi.tianqi.com/)上有各城市的数据。以北京为例,打开页面之后经过简单的尝试就可以发现,每个月的数据是以http://lishi.tianqi.com/{城市的拼音}/{年份+月份}.html这样的地址来展示的,于是可据此构建一个函数,函数的作用是默认返回北京市2018年1月到12月的 url:

def get_url(city='beijing'):
    for time in range(201801,201813):
        url = "http://lishi.tianqi.com/{}/{}.html".format(city,time)
        yield url

爬取数据

 

有了 url 地址就可以用 requests 来抓取。这里要注意,得加上自己的 cookies,否则会返回 404页面,应该是对方网站做了反爬。拿到返回值之后,我用 bs4 库的 select 函数提取数据。select 函数使用的是 css 选择器的语法。由于需要进行一定的数据分析,所以这里没有将数据保存到文件,而是直接使用 pandas 的 dataframe 进行储存。

html = requests.get(url=url, headers=header, cookies=cookie)
soup = BeautifulSoup(html.content, 'html.parser')
date = soup.select("#tool_site > div.tqtongji2 > ul > li:nth-of-type(1) > a")
max_temp = soup.select("#tool_site > div.tqtongji2 > ul > li:nth-of-type(2)")
min_temp = soup.select("#tool_site > div.tqtongji2 > ul > li:nth-of-type(3)")
weather = soup.select("#tool_site > div.tqtongji2 > ul > li:nth-of-type(4)")
wind_direction = soup.select("#tool_site > div.tqtongji2 > ul > li:nth-of-type(5)")
date = [x.text for x in date]
max_temp = [x.text for x in max_temp[1:]]
min_temp = [x.text for x in min_temp[1:]]
weather = [x.text for x in weather[1:]]
wind_direction = [x.text for x in wind_direction[1:]]
data = pd.DataFrame([date,max_temp,min_temp,weather,wind_direction]).T

对12个月份进行抓取后再汇总,就得到了北京2018年全年的天气数据,包括最高温度、最低温度、天气状况、风向等信息。适当加工下信息,我们用一些图表来进行可视化的展示:

 

平均温度的分布

seaborn.distplot(result['平均温度'])

平均温度是使用每日最高温度和最低温度取平均的值。北京平均温度在0度和20多度的日子是最多的。

 

按月查看温度走势

result.groupby(result['日期'].apply(lambda x:x.month)).mean().plot(kind='line')

天气状况分布

seaborn.countplot(result['天气状况'])

晴天和多云是北京一年中主要的天气。

 

 

各月降水天数统计

line = pyecharts.Line("各月降水天数统计")
line.add("降水天数", month, is_rain, is_fill=True, area_opacity=0.7, is_stack=True)
line.add("未降水天数", month, no_rain, is_fill=True, area_opacity=0.7, is_stack=True)

这里用 pyecharts 做了一个堆叠折线图。北京的降水天数不多,主要在7、8月份。可以对比下重庆的数据,差别就很明显了:

 

风向统计

directions = ['北风', '西北风', '西风', '西南风', '南风', '东南风', '东风', '东北风']
schema = []
v = []
days = result['风向'].value_counts()
for d in directions:
    schema.append((d,100))
    v.append(days[d])
v = [v]
radar = pyecharts.Radar()
radar.config(schema)
radar.add("风向统计", v, is_axisline_show=True)

为了让结果更加直观,这里采用了 pyecharts 里的雷达图,并且将8个维度按真实方向的角度来排列。通常认为,北京冬季盛行西北风,夏季盛行东南风。不过从数据上来看,西南风才是北京2018年的最热门的风向。

 

以上就是我这个项目所做的工作,内容还是比较基础的。大家可以做进一步的扩展,比如爬取其他的城市,然后进行全国多城市的天气比较,或者结合地图进行可视化。

 

几个相关库的官网都很不错,供参考:

  1. requests http://cn.python-requests.org/zh_CN/latest/

  2. bs4 https://www.crummy.com/software/BeautifulSoup/bs4/doc.zh/

  3. pyecharts http://pyecharts.org/

  4. seaborn http://seaborn.pydata.org/

  5. pandas https://pandas.pydata.org/

文章源码及相关文件:

https://gitee.com/crossin/snippet.git

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值