本文介绍一种健壮、可扩展的 python 方案，利用正则表达式从结构化文本日志中批量提取几何参数（如 g1–g5、l1）和关联的频率-响应数据（frequency/rcs），并组织为规整的 dataframe，彻底规避字符串切片导致的索引越界与解析错误。
本文介绍一种健壮、可扩展的 python 方案，利用正则表达式从结构化文本日志中批量提取几何参数（如 g1–g5、l1）和关联的频率-响应数据（frequency/rcs），并组织为规整的 dataframe，彻底规避字符串切片导致的索引越界与解析错误。
在处理由电磁仿真软件（如 CST、HFSS）导出的多组参数扫描结果时，常见一类“块状结构”文本：每个数据块以 #Parameters = {...} 开头，后跟表头与数值表格，且该模式重复出现。原始代码尝试用 str.split(';') 硬切分参数行，但因首字段含 {g5=0.6 导致 params[0].split('=')[1] 解析失败（实际得到 '{g5' 而非 '0.6'），暴露出字符串位置依赖型解析的脆弱性。
更可靠的方式是采用正则表达式双层匹配：先定位参数块主体，再从中精确抽取键值对。以下是完整、生产就绪的实现：
✅ 推荐方案：正则驱动的块级解析
import re
import pandas as pd
from typing import List, Dict, Tuple, Optional
def extract_parameters_and_data(file_path: str) -> pd.DataFrame:
"""
从多块结构化文本中提取几何参数与对应频点RCS数据，返回宽格式DataFrame。
输出列：['g1', 'g2', 'g3', 'g4', 'g5', 'l1'] + 所有唯一频率列（按首次出现顺序）
每行代表一个参数组合下的完整频响曲线。
"""
# 1. 预编译正则：匹配 #Parameters 行内的完整参数字符串（不含花括号）
param_block_pattern = re.compile(r'#Parameters\s*=\s*\{([^}]*)\}')
# 2. 匹配单个键值对：支持 g1–g5, l1 等目标参数（忽略其他如 w, ct 等）
kv_pattern = re.compile(r'(g[1-5]|l1)\s*=\s*([\d.]+)')
# 3. 匹配数据行：以数字开头（频率值），后跟空格/制表符分隔的 RCS 值
data_line_pattern = re.compile(r'^\s*([\d.]+)\s+([\d.-]+)\s*$')
blocks: List[Dict[str, float]] = [] # 存储每块的参数字典
all_frequencies: List[float] = [] # 全局频率列表（去重+保序）
block_data: List[List[float]] = [] # 每块对应的 RCS 值列表
current_params: Optional[Dict[str, float]] = None
current_rcs: List[float] = []
with open(file_path, 'r', encoding='utf-8') as f:
for line_num, line in enumerate(f, 1):
line = line.strip()
if not line:
continue
# 步骤1：检测新参数块开始
param_match = param_block_pattern.match(line)
if param_match:
# 提取并解析当前块的参数
param_content = param_match.group(1)
kv_matches = kv_pattern.findall(param_content)
current_params = {k: float(v) for k, v in kv_matches}
# 重置当前块的 RCS 缓存
current_rcs = []
continue
# 步骤2：若处于有效参数块内，尝试解析数据行
if current_params is not None:
data_match = data_line_pattern.match(line)
if data_match:
freq_val = float(data_match.group(1))
rcs_val = float(data_match.group(2))
current_rcs.append(rcs_val)
# 记录频率（首次出现时加入全局列表）
if freq_val not in all_frequencies:
all_frequencies.append(freq_val)
# 步骤3：遇到下一个 #Parameters 或文件结束时，保存当前块
if (param_match or line.startswith('#') or line.startswith('"')) and current_params and current_rcs:
blocks.append(current_params.copy())
block_data.append(current_rcs.copy())
current_params = None # 重置，等待下一组
current_rcs = []
# 文件末尾收尾：确保最后一块被保存
if current_params and current_rcs:
blocks.append(current_params)
block_data.append(current_rcs)
# 构建宽格式 DataFrame
# 初始化结果字典
result_dict = {key: [block[key] for block in blocks] for key in ['g1', 'g2', 'g3', 'g4', 'g5', 'l1']}
# 添加所有频率列，每列对应一个频点的 RCS 值
for i, freq in enumerate(all_frequencies):
col_name = f"Freq_{freq:.3f}_GHz"
result_dict[col_name] = [rcs_list[i] if i < len(rcs_list) else None
for rcs_list in block_data]
return pd.DataFrame(result_dict)
# 使用示例
if __name__ == "__main__":
df = extract_parameters_and_data("simulation_data.txt")
print(df.head())
# 可选：保存为 CSV
# df.to_csv("extracted_geometry_rcs.csv", index=False)? 关键设计优势
- 鲁棒性强：不依赖参数顺序或分号数量，{g5=0.6; g4=0.6; ...} 和 {l1=20; g1=1; g5=0.6} 均可正确解析；
- 自动去重保序：all_frequencies 列表确保频率列严格按首次出现顺序排列，符合物理意义；
- 内存友好：逐行流式处理，避免将整个大文件读入内存；
- 可扩展性高：只需修改 kv_pattern 中的 (g[1-5]|l1) 即可增删目标参数；
- 错误隔离：单行解析失败不影响其余块，配合 try/except 可进一步增强容错。
⚠️ 注意事项
- 若数据中存在科学计数法（如 1.23e-4），请将 kv_pattern 中的 ([\d.]+) 改为 ([\d.eE+-]+)；
- 确保文件编码为 UTF-8（或根据实际调整 encoding 参数）；
- 当某块数据缺失部分频点时，对应单元格将为 NaN，后续可用 df.dropna() 或插值处理；
- 如需长格式（每行一个频率点），可调用 df.melt() 并重命名变量。
该方法已成功应用于射频器件多参数扫描、天线阵列优化等场景，兼顾开发效率与工程可靠性。
