Question

我正在学习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`中匹配，则删除它匹配的所有这些行。

Answer 1

我会利用 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递增以索引下一个数组元素，将该行留在原位< / p>

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;
        }
    } 
}

如何删除与另一个文件中的元素匹配的行

1 个答案: