目录
  • 前言
  • 前期准备
  • 数据的选择与获取
  • 分词
  • 筛选与可视化
  • 总结

前言

”数据可视化“这个话题,相信大家并不陌生,在一些平台,经常可以看到一些动态条形图的视频,大多都是某国家 gdp 的变化或者不同国家疫情中感染人数的变化等等。

这篇文章,我们将使用 python 绘制动态词频条形图,顾名思义,就是以词频作为数量指标的动态条形图。

前期准备

输入以下命令,安装必须的库:

pip install jiashuresearchtools
pip install jieba
pip install pandas
pip install bar_chart_race

数据的选择与获取

我们这次使用的数据是简书文章收益排行榜,日期范围为 2020 年 6 月 20 日至 2021 年 9 月 18 日。

从网页中解析数据的过程较为复杂,我们使用简书数据科学库 jianshuresearchtools 完成。

为方便调试,我们使用 jupyter notebook 进行交互式开发。

导入 jianshuresearchtools,并为其设置别名:

import jianshuresearchtools as jrt

调用接口,获取 2021 年 9 月 17 日的数据:

jrt.rank.getarticlefprankdata("20210917")

返回的数据如下:

[{'ranking': 0,
  'aslug': 'a03adf9d5dd5',
  'title': '幸得君心似我心',
  'author_name': '雁阵惊寒',
  'author_avatar_url': 'https://upload.jianshu.io/users/upload_avatars/26225608/682b892e-6661-4f98-9aab-20b4038a433b.jpg',
  'fp_to_author': 3123.148,
  'fp_to_voter': 3123.148,
  'total_fp': 6246.297},
 {'ranking': 1,
  'aslug': '56f7fe236842',
  'title': '伤痕',
  'author_name': '李文丁',
  'author_avatar_url': 'https://upload.jianshu.io/users/upload_avatars/26726969/058e18c4-908f-4710-8df7-1d34d05d61e3.jpg',
  'fp_to_author': 1562.198,
  'fp_to_voter': 1562.198,
  'total_fp': 3124.397},
 (以下省略)

可以看出,返回的数据中包含文章的排名、标题、作者名、作者头像链接和关于简书资产的一些信息。

我们只需要文章的标题进行统计,所以我们将上面获取到的数据赋值给变量 raw_data,然后:

[item["title"] for item in raw_data]

使用列表推导式,我们得到了文章标题组成的列表。

为方便处理,我们将这些数据连接起来,中间用空格分隔:

" ".join([item["title"] for item in raw_data])

但是我们遇到了报错:

typeerror: sequence item 56: expected str instance, nonetype found

从报错信息中可以看出,我们获取到的文章标题列表中有空值,导致字符串的连接失败了。

(空值是因为作者删除了文章)

所以我们还需要加入去除空值的逻辑,代码编程这样:

" ".join(filter(none, [item["title"] for item in raw_data]))

filter 函数在第一个参数为 none 时,默认过滤掉列表中的空值。

现在我们获取到的数据如下:

‘幸得君心似我心 伤痕 短篇|阿生 “我最喜爱的友友”大评选|理事会 · 中秋嘉年华,等你来! 是缘是劫无须问,石火穷年一蝶迷 职业日记|从蜜月到陌路:我和美国外教的一点事 红楼||浅谈《红楼梦》开篇一顽石 城市印象|走笔八卦城 花豹与狗的爱情终结在人与动物的战争里(以下省略)

接下来,我们需要获取时间范围内的所有数据。

查询 jrt 的函数文档可知,我们需要一个字符串类型,格式为”yyyymmdd“的参数表示目标数据的日期。

所以我们需要写一段程序,用于实现这些日期字符串的生成,代码如下:

from datetime import date, timedelta
def datestrgenerator():
    start_date = date(2020, 6, 20)
    after = 0
    result = none
    while result != "20210917":
        current_date = start_date + timedelta(days=after)
        result = current_date.strftime(r"%y%m%d")
        yield result
        after += 1

接下来,我们编写一段代码,实现对这些数据的获取:

result = []
for current_date in tqdm(datestrgenerator(), total=455):
    raw_data = jrt.rank.getarticlefprankdata(current_date)
    processed_data = " ".join(filter(none, [item["title"] for item in raw_data]))
    result.append({"date": current_date, 
                   "data": processed_data})

这里使用 tqdm 库显示了一个进度条,非必须。

使用 pandas 库,将我们采集到的数据转换成 dataframe:

df = pandas.dataframe(result)

分词

我们使用 jieba 库实现分词,先尝试对第一条数据进行处理:

jieba.lcut(df["data"][0])

使用 python 标准库 collections 中的 counter 进行词频统计:

counter(jieba.lcut(df["data"][0]))

简单画个条形图:

可以看到,空格和一些标点符号,包括”的“、”我“之类无意义词汇出现频率很高,我们需要将它们剔除出去。

我们构建一个存放停用词的 txt 文档,之后使用如下代码将其读取,并转换成一个列表:

stopwords_list = [item.replace("\n", "") for item in open("stopwords.txt", "r", encoding="utf-8").readlines()]

接下来,编写一个函数,实现停用词的剔除,为了方便后续的数据处理,我们也一并剔除单字和只出现一次的词语:

def process_words_count(count_dict):
    result = {}
    for key, value in count_dict.items():
        if value < 2:
            continue
        if len(key) >= 2 and key not in stopwords_list:
            result[key] = value
    return result

另外,我们使用 jieba 库的 add_word 函数将一些简书中的组织名和专有名词添加到词库中,从而提高分词的准确性,代码如下:

keywords_list = [item.replace("\n", "") for item in open("keywords.txt", "r", encoding="utf-8").readlines()]
for item in keywords_list:
    jieba.add_word(item)

经过一番处理,现在分词效果有了明显的改善:

最后,用这段代码对所有数据进行分词,并将结果保存到另一个 dataframe 中:

data_list = []
date_list = []
for _, item in df.iterrows():
    date_list.append(datetime(int(item["date"][0:4]), int(item["date"][4:6]), int(item["date"][6:8])))
    data_list.append(process_words_count(counter(jieba.lcut(item["data"]))))
processed_df = pandas.dataframe(data_list, index=date_list)

我最终得到的结果是一个 455 行,2087 列的 dataframe。

筛选与可视化

这样多的数据,其中很大一部分都不能代表整体情况,所以我们需要进行数据筛选。

使用以下代码,统计所有列数值的总和,即每个关键词在全部数据中出现的次数,存储到名为 sum 的行中:

try:
    result = []
    for i in range(3000):
        result.append(processed_df.iloc[:, i].sum())
except indexerror:
    processed_df.loc["sum"] = result

运行以下代码,只保留在数据集中出现 300 次以上的关键词:

maller_df = processed_df.t[processed_df.t["sum"] >= 300].t
smaller_df = smaller_df.drop(labels="sum")
smaller_df.columns

现在,数据集中的列数减少到了 24 个,可以进行可视化了。

不要忘记先导入模块:

import bar_chart_race as bcr

使用此模块需要先安装 ffmpeg,这方面教程可以自行查找。

另外,为了支持中文显示,我们需要打开这个模块下的 _make_chart.py 文件,在 import 之后增加以下两行代码:

plt.rcparams['font.sans-serif'] = ['simhei']
plt.rcparams['axes.unicode_minus'] = false

这两行代码将会把 matplotlib 的默认字体替换成支持中文显示的字体。

最后,使用一行代码完成可视化:

bcr.bar_chart_race(smaller_df, steps_per_period=30, period_length=1500, title="简书文章收益排行榜可视化", bar_size=0.8, fixed_max=10, n_bars=10)

在这行代码中,我们使用 smaller_df 作为数据集,输出文件为 output.mp4,帧率设置为 30,每行数据显示 2 秒。

由于数据较多,这一步时间较长,而且会占用较多内存。运行结束后,即可在目录中找到输出的文件。

总结

到此这篇关于python制作动态词频条形图的文章就介绍到这了,更多相关python制作动态词频条形图内容请搜索www.887551.com以前的文章或继续浏览下面的相关文章希望大家以后多多支持www.887551.com!