如何用Python脚本或命令行在Chimera中保存molmap/密度图为MRC格式？

我之前也踩过类似的坑——默认的save命令如果不指定目标对象，很容易存成会话文件或者其他非预期格式。核心解决思路就是明确指定要保存的是molmap生成的密度图对象，而不是默认的结构或者整个会话。下面分两种场景给出具体实现：

一、命令行批量处理（无GUI模式） #

如果是用命令行做批量任务，可以直接把命令串起来执行：

chimera --nogui your_input.pdb "molmap #0 2.0" "save output_density.mrc #1"

• --nogui：启动无界面模式，适合自动化批量脚本
• molmap #0 2.0：对编号为#0的输入结构生成分辨率2.0Å的密度图，生成的密度图会自动分配下一个对象编号#1
• save output_density.mrc #1：明确指定保存#1（也就是刚生成的密度图）为MRC格式。文件名后缀.mrc会让Chimera自动识别格式，也可以加format mrc参数强制指定：save output_density.mrc #1 format mrc

二、Python脚本集成调用 #

如果要把这个步骤集成到Python工作流里，用Chimera的runCommand接口即可：

from chimera import runCommand as rc

# 1. 加载输入结构
rc("open your_input.pdb")

# 2. 生成molmap并给它命名（方便后续引用，避免编号混乱）
rc("molmap #0 2.0 name my_density")

# 3. 保存命名后的密度图为MRC
rc("save final_density.mrc my_density")

# 或者不用命名，直接用对象编号（注意molmap生成的是新对象，编号为当前最大编号+1）
# rc("save final_density.mrc #1")

常见问题排查 #

• 之前存成.py或.x3d是因为没指定目标对象：默认save命令会保存整个会话为.py，如果选中了原结构等其他对象，也可能存成非MRC格式。必须明确指定molmap对应的对象编号或名称。
• 可以在脚本中加入rc("list")命令，查看所有对象的编号和类型，确认哪个是你生成的密度图。

内容的提问来源于stack exchange，提问作者Guangxuan Xu