0

我有一个可以从控制台执行的 perl 脚本，如下所示：

perl perlscript.pl -i input.txt -o output.txt --append

我想从我的 python 代码中执行这个脚本。我发现它subprocess.Popen可以用来连接到 perl，我可以用它传递我的论点。但是，我还想传递一个变量（通过拆分文本文件生成）来代替 input.txt。我尝试了以下方法，但它似乎不起作用，并在第 8 行给出了明显的 TypeError：

import re, shlex, subprocess, StringIO
f=open('fulltext.txt','rb')
text= f.read()
l = re.split('\n\n',str(text))
intxt = StringIO.StringIO()
for i in range(len(l)):
    intxt.write(l[i])
    command_line='perl cnv_ltrfinder2gff.pl -i '+intxt+' -o output.gff --append'
    args=shlex.split(command_line)
    p = subprocess.Popen(args)

有没有其他解决方法？

编辑：这是文件 fulltext.txt 的示例。条目由一行分隔。

Predict protein Domains 0.021 second
>Sequence: seq1 Len:13143 [1] seq1 Len:13143 Location : 9 - 13124 Len: 13116 Strand:+ Score    : 6 [LTR region similarity:0.959] Status   : 11110110000 5'-LTR   : 9 - 501 Len: 493 3'-LTR   : 12633 - 13124 Len: 492 5'-TG    : TG , TG 3'-CA    : CA , CA TSR      : NOT FOUND Sharpness: 1,1 Strand + : PBS   : [14/20] 524 - 543 (LysTTT) PPT   : [12/15] 12553 - 12567

Predict protein Domains 0.019 second
>Sequence: seq5 Len:11539 [1] seq5 Len:11539 Location : 7 - 11535 Len: 11529 Strand:+ Score    : 6 [LTR region similarity:0.984] Status   : 11110110000 5'-LTR   : 7 - 506 Len: 500 3'-LTR   : 11036 - 11535 Len: 500 5'-TG    : TG , TG 3'-CA    : CA , CA TSR      : NOT FOUND Sharpness: 1,1 Strand + : PBS   : [15/22] 515 - 536 (LysTTT) PPT   : [11/15] 11020 - 11034

我想将它们分开并将每个条目块传递给 perl 脚本。所有文件都在同一个目录中。

4

2 回答 2

1

您可能对os 模块 和字符串格式感兴趣

编辑

我想我明白你现在想要什么。如果我错了，请纠正我，但我认为：

  • 您想将 fulltext.txt 拆分为块。
  • 每个块都包含一个 seq(number)
  • 您想为每个块运行一次 perl 脚本，并将您的 seq(number) 作为输入文件

如果这是您想要的，您可以使用以下代码。

import os

in_file = 'fulltext.txt'
seq = []

with open(in_file,'r') as handle:
    lines = handle.readlines()
    for i in range(0,len(lines)):
        if lines[i].startswith(">"):
            seq.append(lines[i].rstrip().split(" ")[1])

for x in seq:
    command = "perl perl cnv_ltrfinder2gff.pl -i %s.txt -o output.txt --append"%x
    os.system(command)
于 2015-04-23T14:00:08.937 回答
0

--infile选项的文档：

输入文件的路径。如果没有提供输入文件，程序将期待来自 STDIN 的输入。

您可以省略--infile并通过管道（stdin）传递输入：

#!/usr/bin/env python
from subprocess import Popen, PIPE

with open('fulltext.txt') as file: # read input data
    blocks = file.read().split('\n\n')

# run a separate perl process for each block
args = 'perl cnv_ltrfinder2gff.pl -o output.gff --append'.split()
for block in blocks:
    p = Popen(args, stdin=PIPE, universal_newlines=True)
    p.communicate(block)
    if p.returncode != 0:
        print('non-zero exit status: %s on block: %r' % (p.returncode, block))

您可以同时运行多个perl脚本：

from multiprocessing.dummy import Pool # use threads

def run((i, block)):
    filename = 'out%03d.gff' % i
    args = ['perl', 'cnv_ltrfinder2gff.pl', '-o', filename]
    p = Popen(args, stdin=PIPE, universal_newlines=True, close_fds=True)
    p.communicate(block)
    return p.returncode, filename

exit_statuses, filenames = zip(*Pool().map(run, enumerate(blocks, start=1)))

它并行运行多个（等于系统上的 CPU 数量）子进程。您可以指定不同数量的工作线程（传递给Pool()）。

于 2015-04-26T00:19:10.410 回答