gpt4 book ai didi

python - 如何加入两个子进程的标准输出并通过管道连接到python中新子进程的标准输入

转载 作者:太空狗 更新时间:2023-10-30 00:06:49 26 4
gpt4 key购买 nike

假设我从 shell 运行了以下命令

{ 
samtools view -HS header.sam; # command1
samtools view input.bam 1:1-50000000; # command2
} | samtools view -bS - > output.bam # command3

对于那些不熟悉 samtools View 的人(因为这是 stackoverflow)。这实际上是在创建一个具有新 header 的新 bam 文件。 bam 文件通常是大型压缩文件,因此在某些情况下即使通过该文件也可能很耗时。一种替代方法是执行 command2,然后使用 samtools reheader 切换 header 。这两次通过大文件。上面的命令一次通过 bam,这对较大的 bam 文件很有用(即使压缩后它们也会大于 20GB - WGS)。

我的问题是如何使用子进程在 python 中实现这种类型的命令。

我有以下内容:

fh_bam = open('output.bam', 'w')
params_0 = [ "samtools", "view", "-HS", "header.sam" ]
params_1 = [ "samtools", "view", "input.bam", "1:1-50000000"]
params_2 = [ "samtools", "view", "-bS", "-" ]
sub_0 = subprocess.Popen(params_0, stderr=subprocess.PIPE, stdout=subprocess.PIPE)
sub_1 = subprocess.Popen(params_1, stderr=subprocess.PIPE, stdout=subprocess.PIPE)
### SOMEHOW APPEND sub_1.stdout to sub_0.stdout
sub_2 = subprocess.Popen(params_2, stdin=appended.stdout, stdout=fh_bam)

非常感谢任何帮助。谢谢。

最佳答案

如果您已经在字符串中包含 shell 命令,那么您可以按原样运行它:

#!/usr/bin/env python
from subprocess import check_call

check_call(r"""
{
samtools view -HS header.sam; # command1
samtools view input.bam 1:1-50000000; # command2
} | samtools view -bS - > output.bam # command3
""", shell=True)

在 Python 中模拟管道:

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

# start command3 to get stdin pipe, redirect output to the file
with open('output.bam', 'wb', 0) as output_file:
command3 = Popen("samtools view -bS -".split(),
stdin=PIPE, stdout=output_file)
# start command1 with its stdout redirected to command3 stdin
command1 = Popen('samtools view -HS header.sam'.split(),
stdout=command3.stdin)
rc_command1 = command1.wait() #NOTE: command3.stdin is not closed, no SIGPIPE or a write error if command3 dies
# start command2 after command1 finishes
command2 = Popen('samtools view input.bam 1:1-50000000'.split(),
stdout=command3.stdin)
command3.stdin.close() # inform command2 if command3 dies (SIGPIPE or a write error)
rc_command2 = command2.wait()
rc_command3 = command3.wait()

关于python - 如何加入两个子进程的标准输出并通过管道连接到python中新子进程的标准输入,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/28774716/

26 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com