如何按升序排列线条

时间:2013-05-15 07:59:54

标签: sed awk

我有一个包含200个pdb文件的文件夹。我想根据第6列按升序排列PDB文件的原子行。我想为文件夹中的每个pdb文件进行就地编辑。你的帮助将不胜感激。

ATOM     81  N   ASN A 248      38.791 -16.708  12.507  1.00 52.04           N 
ATOM     82  CA  ASN A 248      39.443 -17.018  11.206  1.00 54.49           C 
ATOM    422  C   SER A 205      70.124 -29.955   8.226  1.00 55.81           C 
ATOM    423  O   SER A 205      70.901 -29.008   8.438  1.00 46.60           O 
ATOM    303  N   MET A 231      61.031 -38.086  -3.054  1.00 52.32           N 
ATOM    304  CA  MET A 231      60.580 -39.074  -4.047  1.00 64.11           C
ATOM    392  C   GLU B  65      23.248  10.071  -7.321  1.00 48.26           C  
ATOM    393  O   GLU B  65      24.465  10.200  -7.158  1.00 46.53           O
ATOM    394  O   GLU B  65      24.465  10.200  -7.158  1.00 46.53           O

期望的输出

ATOM    392  C   GLU B  65      23.248  10.071  -7.321  1.00 48.26           C
ATOM    393  O   GLU B  65      24.465  10.200  -7.158  1.00 46.53           O
ATOM    394  O   GLU B  65      24.465  10.200  -7.158  1.00 46.53           O 
ATOM    422  C   SER A 205      70.124 -29.955   8.226  1.00 55.81           C  
ATOM    423  O   SER A 205      70.901 -29.008   8.438  1.00 46.60           O 
ATOM    303  N   MET A 231      61.031 -38.086  -3.054  1.00 52.32           N
ATOM    304  CA  MET A 231      60.580 -39.074  -4.047  1.00 64.11           C
ATOM     81  N   ASN A 248      38.791 -16.708  12.507  1.00 52.04           N
ATOM     82  CA  ASN A 248      39.443 -17.018  11.206  1.00 54.49           C  

2 个答案:

答案 0 :(得分:4)

使用sort

sort -n -k 6 inputfile

-n执行数字排序,-k告诉您通过密钥排序。

编辑:对于就地排序,请使用-o选项:

sort -n -k 6 inputfile -o inputfile

答案 1 :(得分:0)

我使用哈希,其键将是第6个字段加上一个计数器,该计数器会增加末尾附加的每一行。这样可以避免覆盖重复的条目并保持稳定的顺序。然后使用asorti()函数按第6个字段排序并打印原始数组的每一行。

script.awk的内容:

{
        ++n
        data[ $6 _ n ] = $0;
}

END {
        asorti( data, mod_data, "@ind_num_asc" )
        l = length( data )
        for ( i = 1; i <= l; i++ ) {
                print data[ mod_data[i] ]
        }
}

像以下一样运行:

awk -f script.awk infile

产量:

ATOM    392  C   GLU B  65      23.248  10.071  -7.321  1.00 48.26           C  
ATOM    393  O   GLU B  65      24.465  10.200  -7.158  1.00 46.53           O
ATOM    394  O   GLU B  65      24.465  10.200  -7.158  1.00 46.53           O
ATOM    422  C   SER A 205      70.124 -29.955   8.226  1.00 55.81           C 
ATOM    423  O   SER A 205      70.901 -29.008   8.438  1.00 46.60           O 
ATOM    303  N   MET A 231      61.031 -38.086  -3.054  1.00 52.32           N 
ATOM    304  CA  MET A 231      60.580 -39.074  -4.047  1.00 64.11           C
ATOM     81  N   ASN A 248      38.791 -16.708  12.507  1.00 52.04           N 
ATOM     82  CA  ASN A 248      39.443 -17.018  11.206  1.00 54.49           C