用pytho写的分割pdb文件一个脚本

最新推荐文章于 2022-03-21 14:46:53 发布

rogerzhanglijie

最新推荐文章于 2022-03-21 14:46:53 发布

阅读量1k

点赞数

CC 4.0 BY-SA版权

分类专栏： Python 文章标签： file list lambda import hex 脚本

本文链接：https://blog.youkuaiyun.com/rogerzhanglijie/article/details/7637730

Python 专栏收录该内容

11 篇文章

订阅专栏

本文介绍了一个Python脚本，用于遍历PDB文件目录，按特定标签进行排序，并将文件按照标签分类存储到不同的文件中。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

import glob
def file_list():
   files=[]
   for file in glob.glob("*.pdb"):
       files.append(file)
   files=sorted(files,key=lambda x: int(x[:-4]))
   return files
def split_pdb(pdb_file,files_list,label1,label2):
   f_pdb=open(pdb_file,"r")
   lines=f_pdb.readlines()
   i=0;
   for line in lines:
       if i==12:
           break
       f_store=open(files_list[i],"a+")
       if label1 in line.strip():
           i+=1
           f_store.write(line)
           print "I am here"
           f_store.close()
           continue
       elif line.strip()==label2:
           continue
       else:
           f_store.write(line)
if __name__=="__main__":
   file_lists=file_list()
   split_pdb("hex6-dna.pdb1",file_lists,"TER","ENDMDL")