You need to enable JavaScript to run this app.
优惠活动
大模型
产品
解决方案
定价
更多

使用wfdb包将ECG信号转CSV生成空文件的问题排查求助

问题根因

生成空文件是命令拼接错误、逻辑缺陷共同导致的:

  • 核心错误:shell命令拼接时多写了冗余的rdsamp字符串。从你打印的执行命令可以看到,记录路径后平白出现了第二个rdsamp参数,不符合rdsamp的命令语法,直接导致命令执行失败,没有任何有效数据输出。
  • 参数格式错误:-r参数和后续的记录路径之间没有预留空格,部分环境下会导致参数识别异常。
  • 重定向逻辑缺陷:仅用>重定向了标准输出流,没有捕获标准错误流,命令执行的报错信息不会写入目标文件,你无法直接看到报错原因,只能得到空文件。
  • 编码不规范:用Python内置函数名dir作为自定义变量名,路径拼接没有做标准化处理,甚至出现了/home//ecg/这种多斜杠的路径笔误,存在兼容隐患。
  • 冗余逻辑:你已经导入了wfdb的Python包,完全不需要调用外部shell命令执行rdsamp,命令拼接的方式本身就容易引入格式错误。
修复方案

方案1:修正原有shell调用逻辑

如果要继续用命令行调用的方式,修正命令拼接错误即可,参考代码如下:

import os
from os import listdir, mkdir, system
from os.path import exists, join, isfile

# 替换原有的dir变量名,避免覆盖内置函数
data_dir = '/home/ecg/ecg-patientwise/'
output_dir = join(data_dir, 'csv/')
if not exists(output_dir):
    mkdir(output_dir)

# 筛选所有.dat格式的记录,提前去掉后缀
records = [f[:-4] for f in listdir(data_dir) if isfile(join(data_dir, f)) and f.endswith('.dat')]

for record in records:
    record_path = join(data_dir, record)
    output_path = join(output_dir, f'{record}.csv')
    # 去掉多余的rdsamp字符串,参数间加空格,加2>&1将错误输出也写入文件方便排查
    cmd = f'rdsamp -r {record_path} -c -H -f 0 -v > {output_path} 2>&1'
    print(f'执行命令:{cmd}')
    system(cmd)

方案2:直接调用wfdb Python接口(推荐)

跳过shell调用,直接用wfdb的Python API读取数据后转存CSV,完全避免命令拼接错误,稳定性更高,参考代码如下:

import os
import wfdb
import pandas as pd
from os import listdir, mkdir
from os.path import exists, join, isfile

data_dir = '/home/ecg/ecg-patientwise/'
output_dir = join(data_dir, 'csv/')
if not exists(output_dir):
    mkdir(output_dir)

records = [f[:-4] for f in listdir(data_dir) if isfile(join(data_dir, f)) and f.endswith('.dat')]

for record in records:
    record_path = join(data_dir, record)
    # 读取信号和头信息,sampfrom=0表示从第一个采样点开始读取
    signals, meta = wfdb.rdsamp(record_path, sampfrom=0)
    # 转成DataFrame,列名为信号名
    df = pd.DataFrame(signals, columns=meta['sig_name'])
    output_path = join(output_dir, f'{record}.csv')
    # 写入文件,第一行保留头信息,和原命令-H参数效果一致
    with open(output_path, 'w', encoding='utf-8') as f:
        f.write(f"# 记录ID:{record}, 采样率:{meta['fs']}Hz, 单位:{meta['units']}\n")
        df.to_csv(f, index=False)
    print(f'转换完成:{output_path}')

内容的提问来源于stack exchange,提问作者Vikram R S

相关产品推荐
方舟 Agent Plan

超全模态模型 × Harness 升级,最新支持 Deepseek-V4.1-Flash、GLM-5.3 系列、Doubao-Seedream-5.0-pro、Kimi-K3 (部分), 限时 9.9 元起

最近更新时间:2026.08.28 11:15:53