| 查看: 522 | 回复: 1 | ||
[求助]
python-自动访问网络，将htm结果存盘到本地
|
|
信息学新手求助： 我现在有个文档，想逐个读取文档里的内容，然后在网页上点击，最后将返回的html结果存盘。 网页网址是 http://www.ncbi.nlm.nih.gov/COG/old/xognitor.html， 需要选择的选项是BeTs to 17 clades，点击的选项是compare to COGs， 最后将返回的html 结果存盘。 网上有个脚本，但是我不用NCBIWWW这个模块，希望能用url 在网上查找的相关脚本是： from Bio import Fasta from Bio.Blast import NCBIWWW filename_for_blast = 'seq/all.fasta' # 要做blast的fasta文件 parser = Fasta.RecordParser() file_for_blast = file(filename_for_blast, 'r') iterator = Fasta.Iterator(file_for_blast, parser) while 1: cur_record = iterator.next() if cur_record is None: break print u'打开输出文件%s' % cur_record.title+'.htm' file_for_blast_result = file(cur_record.title+'.htm', 'w') print 'Now BLAST >', cur_record.title ''' NCBIWWW.blast函数默认的blast参数是: expect='10', alignments='50', timeout=20 ''' try: blast_result = NCBIWWW.blast('blastn', 'nr', cur_record, \ expect='0.001', alignments='1', timeout=1) file_for_blast_result.write(blast_result.read()) print u'OK，搞定一个' except: print u'超时了，做下一个' file_for_blast_result.close() file_for_blast.close() print u'全部搞定!' raw_input() |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
castep几何优化到最大步数没收敛，可以用得到的结果继续优化么？
已经有15人回复
- Python学习-基础教程翻到一半了
已经有11人回复
- 中智法签 法国科技签 得多久能出结果？（北京中智法签）
已经有3人回复
- Python学习第一天
已经有54人回复
- CAD工程制图软件可以用Python编出来不？
已经有9人回复
- 搜集的元胞自动机介绍（很好上手和了解）
已经有144人回复
- 为了孩子在美国上小学，不知所措中，恳请达人指教！
已经有12人回复
- 加拿大西蒙菲沙大学工科排名如何？
已经有6人回复
- [Let's-python系列视频教程][avi][更新到17-01]
已经有15人回复
- VASP频率计算结果分析细节(求助完毕，分享结果)
已经有11人回复
- linux下运行python报错，不知原因
已经有12人回复
- ArcGIS中Python脚本的学习
已经有183人回复
- 我收到的DS2019表和邀请函
已经有8人回复
- 【求助】邀请函中bench fee 超过1000pounds对签证有什么影响？
已经有24人回复
- 【经验分享】德国，汉诺威，延签，访问学者，经验
已经有5人回复
- 【求助】请教：从csc网站得到录取结果至出国一般需要多长时间呢，我是青骨1：1访问学�
已经有8人回复
- 【求助】请问办理学历公证的时候自己翻译学历证书要不要学校公章？
已经有10人回复
- 【求助】国外访学的研究计划
已经有10人回复
2楼2012-06-10 20:49:32
