1.前言
最近有一个任务,处理全省雷电数据。获得一些年,月,日的一些发生/强度等规律。原始数据是存在excel里面的。如下图所示
看到这种格式,感觉还是python好用一些。本文使用python来对数据进行处理。
2.数据处理
import pandas as pd
import numpy as np
# 计算
a = pd.read_excel("F:\\zxProjects\\雷电公报\\三维闪电-2021.xlsx")
# 获取需要的那一列的数据
intens = np.array(a.iloc[:, 2:3])
date = np.array(a.iloc[:, 6:7])
hour = np.array(a.iloc[:, 7:8])
gc_ic = np.array(a.iloc[:, 17:18])
# 创造列表,存放我们想要的结果
hour_times = 24 * [0]
hour_gc_times = 24 * [0]
month_times = 12 * [0]
month_gc_times = 12 * [0]
positive = 12 * [0]
positive_times = 12 * [0]
negative = 12 * [0]
negative_times = 12 * [0]
# 该方法计算正闪,负闪的次数,以及强度
def getPORN(kk, index):
if intens[index][0] > 0:
positive[kk] += intens[index][0]
positive_times[kk] += 1
if intens[index][0] < 0:
negative[kk] += abs(intens[index][0])
negative_times[kk] += 1
for i in range(len(a)):
hour_times[hour[i][0]] += 1
h = str(date[i][0]).split("-")
op = int(h[1]) - 1
if gc_ic[i][0] == 'CG':
hour_gc_times[hour[i][0]] += 1
month_gc_times[op] += 1
month_times[op] += 1
getPORN(op, i)
3.绘图