计算蛋白质TICA坐标时遇XTC读取错误：Compressed 3d coordinate

问题场景 #

在计算蛋白质TICA坐标过程中，使用PyEMMA定义距离特征后加载XTC轨迹文件，触发RuntimeError: XTC read error: Compressed 3d coordinate错误，相关代码及报错信息如下：

特征定义代码 #

feat = pyemma.coordinates.featurizer(t)
feat.add_distances(pairs1, periodic= True)

触发报错的代码 #

reader = pyemma.coordinates.source(traj, features=feat, chunksize=1000) 
                         
data = pyemma.coordinates.load(xtc, features=feat) 

核心报错信息 #

RuntimeError: XTC read error: Compressed 3d coordinate

解决方法 #

1. 验证XTC文件完整性 #

XTC是压缩轨迹格式，文件损坏会直接导致读取失败。用GROMACS工具验证文件状态：

gmx check -f your_trajectory.xtc

如果工具提示文件损坏或帧缺失，重新从模拟输出中生成轨迹文件是最直接的解决方式。

2. 减小读取块大小 #

大尺寸的chunksize可能导致压缩坐标解析时出错，尝试降低单次读取的帧数量：

# 直接加载时指定更小的chunksize
data = pyemma.coordinates.load(xtc, features=feat, chunksize=500)

# 或者通过reader对象设置
reader = pyemma.coordinates.source(xtc, features=feat, chunksize=100)
data = reader.get_output()

3. 转换轨迹文件格式 #

若文件本身无损坏但PyEMMA解析存在兼容性问题，将XTC转换为无压缩的TRR格式后再读取：

gmx trjconv -f your_trajectory.xtc -s topology.tpr -o your_trajectory.trr

之后在代码中替换为TRR文件路径即可。

4. 更新依赖库 #

该错误可能源于MDTraj（PyEMMA依赖的轨迹解析库）的旧版本bug，更新至最新稳定版：

pip install --upgrade mdtraj pyemma

内容的提问来源于stack exchange，提问作者harshita sahni