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