下载雷达反射率的爬虫

忽忽儿

想及时保存下雷达拼图，但是苦于网站保存期限过短，如果前一天有天气发生，可以用这个小程序把网站上的雷达拼图保存下来然后做动画，以便后期做分析哟~~
python初学者，是在前人的程序基础上修改的（参考的程序也是在家园里找的，但是2015年已经不能使用了，于是就在这个基础上改了一下），迫不及待想跟大家分享一下第一次用python编程成功的喜悦。程序还有很多不成熟或者啰嗦的地方，可以正常运行，如果有错误或修改意见望各位赐教！
源代码如下：

1. # -*- coding: utf-8 -*-
   
2. 
   
3. """
   
4. wget radar png from 中国天气网
   
5. @author: kaobozai
   
6. """
   
7. 
   
8. import os
   
9. import time
   
10. import pandas as pd
    
11. import numpy as np
    
12. 
    
13. 
    
14. # 输入前一天选定站点雷达图的url(主要包括中国各区域、北京、天津等)
    
15. def get_img(year_in, month_in, day_in, station_info, path_to_save):
    
16.     hr = ["00", "01", "02", "03", "04", "05", "06", "07", "08", "09", "10", "11", "12", "13", "14", "15", "16", "17",
    
17.           "18", "19", "20", "21", "22", "23"]
    
18.     mi = ["00", "06", "12", "18", "24", "30", "36", "42", "48", "54"]
    
19.     if station_info[1] == "z":
    
20.         image_path = "00000"
    
21.     if station_info[1] != "z":
    
22.         image_path = "00001"
    
23.     for hr_cell in hr:
    
24.         hr_str = hr_cell
    
25.         for mi_cell in mi:
    
26.             mi_str = mi_cell
    
27.             url_path = "http://pi.weather.com.cn/i/product/pic/l/sevp_aoc_rdcp_sldas_ebref_" + station_info + "_l88_pi_" + year_in + month_in + day_in + hr_str + mi_str + image_path + ".png"
    
28.             save_path = path_to_save + station_info + "/"
    
29.             os.system("wget -P " + save_path + " " + url_path)
    
30.             flag = year_in + month_in + day_in + hr_str + mi_str
    
31.             flag2 = "http://pi.weather.com.cn/i/product/pic/l/sevp_aoc_rdcp_sldas_ebref_" + station_info + "_l88_pi_" + year_in + month_in + day_in + hr_str + mi_str + image_path + ".png"
    
32.             print(flag)
    
33.             print(flag2)
    
34.     return 1
    
35. 
    
36. # 调用当天日期并计算前一天，输出数据下载成功与否的结果并压缩文件
    
37. yesterday = time.strftime("%Y-%m-%d-%H-%M-%S", time.localtime(time.time() - 86400)).split("-")
    
38. yesterday1 = yesterday[0] + yesterday[1] + yesterday[2]  # type: str
    
39. path_to_save = "./radar_img/"+yesterday1+"/"
    
40. record = []
    
41. # 选择下载的雷达站点
    
42. radar = pd.read_csv('radar_site.txt')
    
43. datax = np.array(radar['site'])  # 读取站点数据转化为nparray
    
44. datay = np.array(radar['站名'])
    
45. sitelistx = datax.tolist()  # nparray转化为list
    
46. sitelisty = datay.tolist()
    
47. a = input("Enter a site:")
    
48. for i in range(0, 11):
    
49.     if a == sitelisty[i]:
    
50.         print(sitelistx[i])
    
51.         station_info = sitelistx[i]  # 站点选择完毕
    
52.         flag = get_img(yesterday[0], yesterday[1], yesterday[2], station_info, path_to_save)
    
53.         if flag == 1:
    
54.             record += [station_info + " !!!success!!!"]
    
55.         else:
    
56.             record += [station_info + "!!!fail!!!"]
    
57. if os.system("tar -cvf  radar_img_" + yesterday1 + ".tar.gz " + path_to_save)==0:
    
58.     record += ["tar file !!!success!!!"]
    
59. else:
    
60.     record += ["tar file !!!fail!!!"]
    
61. for cell in record:
    
62.     print(cell)
    
63.     print(time.strftime("%Y-%m-%d-%H-%M-%S", time.localtime(time.time())))
    
64. 
    
65. 
    
66. 
    

复制代码

源代码文件和站点文本回复可见：

游客，如果您要查看本帖隐藏内容请回复

rosecjm

感谢分享，这是爬虫吗

zdcgl

感谢分享

qxtjls

学习一下！

qxtjls

学习一下！

1092995269

厉害了，大波姐

跺碎那片云

来看看，一直没时间学啊

肥七七

感谢分享，这是爬虫吗

忽忽儿

rosecjm 发表于 2018-9-22 07:16
感谢分享，这是爬虫吗

是的 爬雷达图片的哈

今朝几朵白云生

谢谢楼主，你这是用哪个版本写的？