用于在Python中下载NCBI文件的多线程



所以最近我承担了从ncbi数据库下载大量文件的任务。然而,我遇到的时候,我必须创建多个数据库。这里的代码可以从ncbi网站下载所有的病毒。我的问题是有没有办法加快下载这些文件的速度?

目前这个程序的运行时间超过5小时。我已经研究过多线程,并且永远无法让它工作,因为其中一些文件需要超过10秒才能下载,我不知道如何处理拖延。(新编程)也有一种方法来处理urllib2。HTTPError: HTTPError 502: Bad Gateway。我有时会遇到这种情况,在使用重启和重启的组合时。这会使程序崩溃,我必须通过改变for语句中的0从另一个位置重新开始下载。

import urllib2
from BeautifulSoup import BeautifulSoup
#This is the SearchQuery into NCBI. Spaces are replaced with +'s.
SearchQuery = 'viruses[orgn]+NOT+Retroviridae[orgn]'
#This is the Database that you are searching.
database = 'protein'
#This is the output file for the data
output = 'sample.fasta'

#This is the base url for NCBI eutils.
base = 'http://eutils.ncbi.nlm.nih.gov/entrez/eutils/'
#Create the search string from the information above
esearch = 'esearch.fcgi?db='+database+'&term='+SearchQuery+'&usehistory=y'
#Create your esearch url
url = base + esearch
#Fetch your esearch using urllib2
print url
content = urllib2.urlopen(url)
#Open url in BeautifulSoup
doc = BeautifulSoup(content)
#Grab the amount of hits in the search
Count = int(doc.find('count').string)
#Grab the WebEnv or the history of this search from usehistory.
WebEnv = doc.find('webenv').string
#Grab the QueryKey
QueryKey = doc.find('querykey').string
#Set the max amount of files to fetch at a time. Default is 500 files.
retmax = 10000
#Create the fetch string
efetch = 'efetch.fcgi?db='+database+'&WebEnv='+WebEnv+'&query_key='+QueryKey
#Select the output format and file format of the files. 
#For table visit: http://www.ncbi.nlm.nih.gov/books/NBK25499/table/chapter4.chapter4_table1
format = 'fasta'
type = 'text'
#Create the options string for efetch
options = '&rettype='+format+'&retmode='+type

#For statement 0 to Count counting by retmax. Use xrange over range
for i in xrange(0,Count,retmax):
    #Create the position string
    poision = '&retstart='+str(i)+'&retmax='+str(retmax)
    #Create the efetch URL
    url = base + efetch + poision + options
    print url
    #Grab the results
    response = urllib2.urlopen(url)
    #Write output to file
    with open(output, 'a') as file:
        for line in response.readlines():
            file.write(line)
    #Gives a sense of where you are
    print Count - i - retmax

使用多线程下载文件:

#!/usr/bin/env python
import shutil
from contextlib import closing
from multiprocessing.dummy import Pool # use threads
from urllib2 import urlopen
def generate_urls(some, params): #XXX pass whatever parameters you need
    for restart in range(*params):
        # ... generate url, filename
        yield url, filename
def download((url, filename)):
    try:
        with closing(urlopen(url)) as response, open(filename, 'wb') as file:
            shutil.copyfileobj(response, file)
    except Exception as e:
        return (url, filename), repr(e)
    else: # success
        return (url, filename), None
def main():
    pool = Pool(20) # at most 20 concurrent downloads
    urls = generate_urls(some, params)
    for (url, filename), error in pool.imap_unordered(download, urls):
        if error is not None:
           print("Can't download {url} to {filename}, "
                 "reason: {error}".format(**locals())
if __name__ == "__main__":
   main()

你应该使用多线程,这是下载任务的正确方式。

"these files take more than 10seconds to download and I do not know how to handle stalling",

我不认为这将是一个问题,因为Python的多线程会处理这个,或者我宁愿说多线程只是为了这种I/o约束的工作。当一个线程等待下载完成时,CPU会让其他线程完成它们的工作。

无论如何,你最好至少试着看看会发生什么。

有两种方法可以影响您的任务。1. 使用进程而不是线程,multiprocess是你应该使用的模块。2. 使用基于事件的模块,gevent是正确的模块。

502错误不是你的脚本的错误。简单地说,可以使用以下模式执行retry

try_count = 3
while try_count > 0:
    try:
        download_task()
    except urllib2.HTTPError:
        clean_environment_for_retry()
    try_count -= 1

在except行中,您可以细化细节以根据具体的HTTP状态代码执行特定的操作。

最新更新