0

我正在学习 Perl，我正在尝试弄清楚如何完成这项任务。我有一个包含一堆文本文件的文件夹，我有一个ions_solvents_cofactors包含三个字母列表的文件。

我编写了一个脚本，该脚本打开并读取文件夹中的每个文件，并应删除特定列 [3] 下与列表中的某些元素匹配的那些行。它工作不正常。我在脚本末尾遇到了一些问题，无法弄清楚它是什么。

我得到的错误是：rm: invalid option -- '5'

我的输入文件如下所示：

ATOM   1592 HD13 LEU D  46      11.698 -10.914   2.183  1.00  0.00           H  
ATOM   1593 HD21 LEU D  46      11.528  -8.800   5.301  1.00  0.00           H  
ATOM   1594 HD22 LEU D  46      12.997  -9.452   4.535  1.00  0.00           H  
ATOM   1595 HD23 LEU D  46      11.722  -8.718   3.534  1.00  0.00           H  
HETATM 1597  N1  308 A   1       0.339   6.314  -9.091  1.00  0.00           N  
HETATM 1598  C10 308 A   1      -0.195   5.226  -8.241  1.00  0.00           C  
HETATM 1599  C7  308 A   1      -0.991   4.254  -9.133  1.00  0.00           C  
HETATM 1600  C1  308 A   1      -1.468   3.053  -8.292  1.00  0.00           C

这是脚本：

#!/usr/bin/perl -w

$dirname = '.';
opendir( DIR, $dirname ) or die "cannot open directory";
@files = grep( /\.txt$/, readdir( DIR ) );

foreach $files ( @files ) {

    open( FH, $files ) or die "could not open $files\n";
    @file_each = <FH>;
    close FH;

    close DIR;

    my @ion_names = ();

    my $ionfile   = 'ions_solvents_cofactors';
    open( ION, $ionfile ) or die "Could not open $ionfile, $!";
    my @ion = <ION>;
    close ION;

    for ( my $line = 0; $line <= $#file_each; $line++ ) {

        chomp( $file_each[$line] );
        if ( $file_each[$line] =~ /^HETATM/ ) {
            @is = split '\s+', $file_each[$line];
            chomp $is[3];
        }

        foreach ( $file_each[$line] ) {    #line 39

            if ( "@ion" =~ $is[3] ) {
                system( "rm $file_each[$line]" );
            }
        }
    }
}

因此，例如，如果308从输入文件中匹配文件 ions_cofactors_solvents`，则删除所有匹配的行。

4

1 回答 1

0

我会使用 Tie::File 模块，它允许您将tie数组添加到模块中，以便您对数组所做的任何更改都会反映在文件中

我曾经使用选项glob来查找所有.txt文件，:bsd_glob以支持文件路径中的空格

第一项工作是构建一个%matches将所有值映射ions_solvents_cofactors到 1 的哈希。这使得测试 PDB 文件的所需值变得很简单

然后只需要tie在每个.txt文件上使用，并测试每一行以查看第 4 列中的值是否在哈希中表示

我使用变量$i来索引@file映射磁盘文件的数组。如果找到匹配项，则使用 删除数组元素splice @file, $i, 1。（这自然会$i在不递增的情况下按顺序索引下一个元素$i。）如果没有匹配，则$i递增以索引下一个数组元素，将行留在原处

use strict;
use warnings 'all';

use File::Glob ':bsd_glob';
use Tie::File;

my %matches = do {
    open my $fh, '<', 'ions_solvents_cofactors.txt';
    local $/;
    map { $_ => 1 } split ' ', <$fh>;
};

for my $pdb ( glob '*.txt' ) {

    tie my @file, 'Tie::File', $pdb or die $!;

    for ( my $i = 0; $i < @file; ) {

        next unless my $col4 = ( split ' ', $file[$i] )[3];

        if ( $matches{$col4} ) {
            printf qq{Removing line %d from "%s"\n},
                    $i+1,
                    $pdb;
            splice @file, $i, 1;
        }
        else {
            ++$i;
        }
    } 
}
于 2018-01-09T13:22:53.707 回答