当Python遇到哈工大无人机校庆彩排

6.7是哈工大百年校庆，我与工大结缘已有11载，从翩翩少年变成大叔，在决胜脱贫攻坚的决胜之年，我也正在朝着大爷的路上狂奔，没有一点减速的迹象。

17年种草大疆的spark，感觉电子产品很炫酷，但工大百年校庆特殊的时间，无人机表演我看的很认真，忘记了烦恼忧愁，去T大爷的，自豪激动心情不时夹杂着woc、NB词语。大家没看的可以长按下方图片，跟大邓重温一下，真的很好看，不骗你^_^

本来没想写窝工，窝工不缺我宣传，不过正好在知乎看到一个问题

其中有一个回答挺合我直男胃口，而且是用Python做的分析，就拿来分享给大家。一方面大家可以学学Python，娱乐一下；另一方面为母校做做宣传，我们工大特别适合以下两类人群报考:

1. 对工科感兴趣的男生

2. 对工科男感兴趣的女生

扯远了，关于知乎问题，下面是知乎用户Yihui 校友做的回答

看这个问题，作为一个严谨的工大学子，第一个想到的是到底用了多少无人机？有没有千架？数一数！于是在B站截了几个图，写个Python小程序数无人机的数量~

这个图片里面有大概871架，代码是这样的：

作者：Yihui
链接：https://www.zhihu.com/question/399162890/answer/1271327705
来源：知乎
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

from skimage.feature import blob_log
from skimage.color import rgb2gray
import skimage.io
import matplotlib.pyplot as plt

image = skimage.io.imread('1920-2020.png')
image_gray = rgb2gray(image)

# Laplacian of Gaussian
blobs_log = blob_log(image_gray, max_sigma=4, threshold=.09)
blobs_log[:, 2] = 4

fig, ax = plt.subplots(1)
ax.set_aspect('equal')

# Show the image
ax.imshow(image_gray)

n = 0
for blob in blobs_log:
    y, x, r = blob
    if x > 2100:    # remove bilibili logo
        continue
    n += 1
    c = plt.Circle((x, y), r, color='red', linewidth=1, fill=False)
    ax.add_patch(c)
ax.set_axis_off()
ax.set_title('{} drones'.format(n))
print(n)
plt.show()

原图是这样的：

当然，这个图片里面不是所有的无人机都是亮的，再来一张~

这个图里面更少一点，大概801架。那到底有多少架呢？可以掏出尺子量一量下面这张图，或者去直接问校友~

结论是用了60 x 15 = 900 架，没有一千架无人机，所以问题不成立！

其实

参演的无人机，跟春晚的跳舞机器人一样，都是哈工大机器人集团的产品。

往期文章读完本文你就了解什么是文本分析

综述:文本分析在市场营销研究中的应用
从记者的Twitter关注看他们稿件的党派倾向？

Pandas时间序列数据操作
readability: 英文文本数据可读性库

Matplotlib可视化教程～

Matplotlib中的plt和ax都是啥？

70G上市公司定期报告数据集
5个小问题带你理解列表推导式
文本数据清洗之正则表达式
Python网络爬虫与文本数据分析
如何批量下载上海证券交易所上市公司年报
Numpy和Pandas性能改善的方法和技巧
漂亮~pandas可以无缝衔接Bokeh
YelpDaset: 酒店管理类数据集10+G

大家都是怎么看待工大无人机表演，欢迎留言~