使用python提取国家气象中心提供的日降水站点资料

a1s2d3

使用python提取国家气象中心提供的日降水站点资料，并将其处理为月平均降水资料，在数据处理过程中有以下注意事项：

• 要先把使用到的降水资料文件放在一个空文件夹中
• 注意国家气象中心提供的降水资料的数据格式
• 月平均降水资料是每日降水总量之和，不用除以天数
• 1. import numpy as np
     
  2. from os.path import join
     
  3. from os import listdir
     
  4. 
     
  5. root = r'E:\pre' #原始资料存放的文件夹
     
  6. for name in listdir(root):
     
  7.     filename = join(root,name)
     
  8.     with open(filename) as fp:
     
  9.         line = fp.readlines()
     
  10.         data = [item.strip().split() for item in line]
      
  11.         data = list(data)
      
  12.         data = np.array(data,dtype = int)
      
  13.     fp.close()
      
  14.    
      
  15.     if data[0,4] < 1979:
      
  16.         continue
      
  17.     root1 = r'E:\pre' #读取原始数据后按站点存放数据的文件夹
      
  18.     name0 = '1234'
      
  19.     xave = 0 #月平均降水量
      
  20.     for i in range(len(data)): #遍历列表的每一行
      
  21.         name1 = str(data[i,0]) #站点编号赋值给name1
      
  22.         #处理降水数据
      
  23.         if data[i,12] != 0:
      
  24.             data[i,9] = 0.0
      
  25.         if data[i,9] >= 32000:
      
  26.             if data[i,9] == 32700:
      
  27.                 data[i,9] = 0.0
      
  28.             else:
      
  29.                 data[i,9] = data[i,9] - 32000
      
  30.         elif data[i,9] >= 31000:
      
  31.             data[i,9] = data[i,9] - 31000
      
  32.         elif data[i,9] >=30000:
      
  33.             data[i,9] = data[i,9] - 30000
      
  34.                     
      
  35.         if name1 != name0: #若该编号与上一个不同说明读取的是下一个站点的数据
      
  36.             if i > 0 : #从第二个站点开始写入上一个站点的月平均数据
      
  37.                 name0 = name0+'.dat'
      
  38.                 filename1 = join(root1,name0)
      
  39.                 with open(filename1,'a') as fp1: #打开储存新站点的文件
      
  40.                     data_str = str(data[i,4]) + ' '  + str(data[i,5])\
      
  41.                         + ' ' + str(xave) + '\n'
      
  42.                     fp1.write(data_str)
      
  43.                 fp1.close()
      
  44.             name0 = name1
      
  45.             xave = 0 #平均数归零
      
  46.         xave = xave + data[i,9]

  复制代码
  详情也可见：https://blog.csdn.net/qq_47672397/article/details/124149949?spm=1001.2014.3001.5501
  

zhoulm177

学习了。楼主厉害。

lxy287131416

留名学习，学习！！！

a1s2d3

嘿嘿，大家一起学习，但是这个程序其实有一点问题，就是按我的想法是它在打开下一个文件的时候才会写入上一个文件的最后一个站点的信息，所以最开始在循环写入结束后加了一句用来写最后一个站点信息的语句，但是后来发现最后一个站点的数据被我写了两次，然后删了那句就正常了，但是一直不理解，所以想请教下大佬们我是哪里想错了

平凡之路

多好的资源，必须顶下

guo066600

请问您用的降水数据的时间范围是什么？在国家气象中心没能搜到1951-2020的逐日数据，能告知在哪里可以下载吗？

a1s2d3

guo066600 发表于 2022-7-12 17:11
请问您用的降水数据的时间范围是什么？在国家气象中心没能搜到1951-2020的逐日数据，能告知在哪里可以下载 ...

抱歉啊，我也不知道在哪里下载，我的数据是我的老师给我的

小羊趴在云朵上

感谢分享，很有帮助

miaocong

学习了。楼主好厉害。请问一下国家气象中心提供的日降水站点资料的时间是北京时间还是世界时间？

miaocong

a1s2d3 发表于 2022-4-17 12:05
嘿嘿，大家一起学习，但是这个程序其实有一点问题，就是按我的想法是它在打开下一个文件的时候才会写入上一 ...

求问一下，在国家气象中心下得过去三小时降水量比如是18点的，那它实际是18点到24点还是15点到18点的呢？像图片之中的情况我想计算7月29号的降水量我是要把7月29号从0点加到21点，还是从7月29号的3点开始加到7月30号的0点？有没有大佬给指导一下，感谢