雨滴谱数据处理

一、parsivel雨滴谱数据：

二、数据介绍：

 总结一下就是，第一行数据是第一个速度对应的全部32个尺度粒子个数，第二行数据是第二个速度对应的全部32个尺度粒子个数。

 三、速度、尺度一览表：

四、先将速度、尺度数据放入原雨滴谱txt数据中，以方便后期计算：

（1）先按时间进行分组

（2）再按行读取，将每个部分变成列表，放入dataframe里

（3）最后插入雨滴谱速度尺度文件

#!usr/bin/env python
# -*- coding:utf-8 -*-
"""
@author: Suyue
@file: speedeeinsert.py
@time: 2024/04/07
@desc:
"""
import numpy as np
import pandas as pd
df1 = pd.read_excel('/尺度速度.xls')

file_path = '/NM004-20230627224400-20230627224859-0.txt'

# 读整个txt文件读取到单个字符串
with open(file_path, 'r', errors='ignore') as file:
    file_content = file.read()

# 按时间戳拆分内容以查找单独的部分
# 时间戳的格式为 YYYY-MM-DD HH：MM：SS，因此我们将使用正则表达式根据此模式进行拆分
import re
sections = re.split(r'\d{4}-\d{2}-\d{2} \d{2}:\d{2}:\d{2}\n', file_content)

# 如果txt第一个元素为空值（由于拆分），则将其删除
if not sections[0]:
    sections.pop(0)

# 将每个部分放入DataFrame
dataframes = []
# 由于一个列表为速度、尺度中文名，要避开这行所以从第二个列表取值
for section in sections[1:]:
    # 将字符串拆分为几行，然后按空格拆分每行并转换为 DataFrame
    lines = section.strip().split('\n')
    # print(lines)
    matrix = [line.split() for line in lines]
    df = pd.DataFrame(matrix)

    df.columns = df1['直径']
    df.index = df1['速度']

    # 将df放入dataframes中
    dataframes.append(df)
    # print(dataframes)

# 显示每个dataframe形状以确认
df_shapes = [df.shape for df in dataframes]
print(df_shapes)

 结果：