Python sys.argv和带有定义的输入组合不起作用

问题描述

我有一个如下所示的Python脚本，我试图在shell中运行。我想使用argv或input来提供称为index1和index2（ACCGTCG和TTCCAGC）的特定输入以及要处理的文件名。我尝试了以下两种方法（分别使用sys.argv和输入），但没有输出。

使用sys.argv：

#!/usr/bin/python

import sys
from Bio import SeqIO

def dual_index_positions():
    script=sys.argv[0]
    index1=sys.argv[1]
    index2=sys.argv[2]
    input_file=sys.argv[3]
    count=0 
    with open(input_file,"r") as Fastq:
        for record in SeqIO.parse(Fastq,'fastq'): 
            if index1 in record.seq and index2 in record.seq: 
                print(record.name)
                ind1_rec=record.seq.find(index1) 
                ind2_rec=record.seq.find(index2)
                rp_ind2=ind2_rec+len(index2)
                dist=(ind1_rec)-(rp_ind2)
                print('Index1 and index2 positions are '+ str(ind1_rec+1) + ' and ' + str(ind2_rec+1) + ' respectively' +
                  '; distance is: ' + str(dist))
                count+=1
    print('The total number of hits is: '+ str(count))
    


if __name__ == '__dual_index_positions__':
    dual_index_positions()

接受输入：

#!/usr/bin/python

from Bio import SeqIO 

def dual_index_positions(): 
    input_file=input('please enter your input_file: ')
    index1=str(input('please enter your index 1: '))
    index2=str(input('please enter your index 2: '))
    count=0 
    with open(input_file,'fastq'):
            if index1 in record.seq and index2 in record.seq: 
                print(record.name)
                ind1_rec=record.seq.find(index1) 
                ind2_rec=record.seq.find(index2)
                rp_ind2=ind2_rec+len(index2)
                dist=(ind1_rec)-(rp_ind2)
                print('Index1 and index2 positions are '+ str(ind1_rec+1) + ' and ' + str(ind2_rec+1) + ' respectively' +
                  '; distance is: ' + str(dist))
                count+=1
    print('The total number of hits is: '+ str(count))
    


if __name__ == '__dual_index_positions__':
    dual_index_positions()

有人可以帮我找到错误的地方吗？预先谢谢你。

解决方法

代替这个

if __name__ == '__dual_index_positions__':
    dual_index_positions()

使用此

if __name__ == '__main__':
    dual_index_positions()

__name__ == '__dual_index_positions__'的RHS并非您必须调用的函数名称，而是模块的名称，默认情况下，其值为__main__