大家好嗎? 我有一個文本檔案,其中包含按列排列的資料,我需要用其他列替換一些列,具體來說就是xyz坐標。
預先感謝
。HETATM 1 C LIG 1 -0.517 1.592 -0.048 1.00 0.00 0.212 A
hetatm 2 c lig 1 0.017 -0.536 0.534 1.00 0.00 0.149 a
Hetatm 3 c lig 1 1.133 0.155 0.029 1.00 0.00 0.212 a
Hetatm 4 n lig 1 -1.027 0.379 0.499 1.00 0.00 -0.337 n
Hetatm 5 n lig 1 0.789 1.466 -0.324 1.00 0.00 -0.219 na
Hetatm 6 c lig 1 -2.429 0.112 0.889 1.00 0.00 0.221 c
氦氣發生器7 C LIG 1 -3.179 -0.453 -0.210 1.00 0.00 -0.097 C
Hetatm 8 C lig 1 -3.805 -0.925 -1.124 1.00 0.00 0.014 C
HITATM 9 N LIG 1 2.482 -0.388 -0.118 1.00 0.00 -0.095 N
HITATM 10 O LIG 1 2.619 -1.549 0.253 1.00 0.00 -0.530 OA
Hetatm 11 O Lig 1 3.362 0.305 -0.578 1.00 0.00 -0.530 OA
原子1 c lig 1 -13.469 13.704 72.248 -0.37 -0.04 0.212 75.145
原子2 c lig 1 -14.243 15.824 72.493 -0.41 -0.03 0.149 75.145
原子3 c lig 1 -15.124 15.039 71.727 -0.40 -0.04 0.212 75.145
原子4 n連接1 -13.200 14.974 72.836 -0.28 0.06 -0.337 75.145
原子5 n 襯里1 -14.635 13.735 71.586 -0.32 0.05 -0.219 75.145
原子6 c lig 1 -11.994 15.348 73.608 -0.46 -0.02 0.221 75.145
原子7 C連接1 -12.341 15.781 74.943 -0.66 0.01 -0.097 75.145
原子8 c lig 1 -12.628 16.141 76.055 -0.66 -0.00 0.014 75.145
原子9 n lig 1 -16.387 15.490 71.145 -0.60 0.01 -0.095 75.145
原子10 o lig 1 -17.127 14.595 70.751 -0.10 0.02 -0.530 75.145
原子11 o lig 1 -16.631 16.674 71.082 -0.58 -0.08 -0.530 75.145
uj5u.com熱心網友回復:
目前還不清楚你是否想在背景關系中對齊這些行,但如果你只是想用另一個檔案中同一行的列來替換第6、7和8列,你可以這樣做:
$ cat file1
hetatm 1 c lig 1 -0.517 1.592 -0.048 1.00 0.00 0.212 a
hetatm 2 c lig 1 0.017 -0.536 0.534 1.00 0.00 0.149 a
Hetatm 3 c lig 1 1.133 0.155 0.029 1.00 0.00 0.212 a
Hetatm 4 n lig 1 -1.027 0.379 0.499 1.00 0.00 -0.337 n
Hetatm 5 n lig 1 0.789 1.466 -0.324 1.00 0.00 -0.219 na
Hetatm 6 c lig 1 -2.429 0.112 0.889 1.00 0.00 0.221 c
氦氣發生器7 C LIG 1 -3.179 -0.453 -0.210 1.00 0.00 -0.097 C
Hetatm 8 C lig 1 -3.805 -0.925 -1.124 1.00 0.00 0.014 C
HITATM 9 N LIG 1 2.482 -0.388 -0.118 1.00 0.00 -0.095 N
HITATM 10 O LIG 1 2.619 -1.549 0.253 1.00 0.00 -0.530 OA
Hetatm 11 O Lig 1 3.362 0.305 -0.578 1.00 0.00 -0.530 OA
$ cat file2
原子1 c lig 1 -13.469 13.704 72.248 -0.37 -0.04 0.212 75.145
原子2 c lig 1 -14.243 15.824 72.493 -0.41 -0.03 0.149 75.145
原子3 c lig 1 -15.124 15.039 71.727 -0.40 -0.04 0.212 75.145
原子4 n連接1 -13.200 14.974 72.836 -0.28 0.06 -0.337 75.145
原子5 n 襯里1 -14.635 13.735 71.586 -0.32 0.05 -0.219 75.145
原子6 c lig 1 -11.994 15.348 73.608 -0.46 -0.02 0.221 75.145
原子7 C連接1 -12.341 15.781 74.943 -0.66 0.01 -0.097 75.145
原子8 c lig 1 -12.628 16.141 76.055 -0.66 -0.00 0.014 75.145
原子9 n lig 1 -16.387 15.490 71.145 -0.60 0.01 -0.095 75.145
原子10 o lig 1 -17.127 14.595 70.751 -0.10 0.02 -0.530 75.145
原子11 o lig 1 -16.631 16.674 71.082 -0.58 -0.08 -0.530 75.145
$ awk '{getline s < "file2"; split(s, a); 6 = a[6]; $7 = a[7]; $8 = a[8]}1'/span> file1
hetatm 1 c lig 1 -13.469 13.704 72.248 1.00 0.00 0.212 a
hetatm 2 c lig 1 -14.243 15.824 72.493 1.00 0.00 0.149 a
hetatm 3 c lig 1 -15.124 15.039 71.727 1.00 0.00 0.212 a
Hetatm 4 n lig 1 -13.200 14.974 72.836 1.00 0.00 -0.337 n
Hetatm 5 n lig 1 -14.635 13.735 71.586 1.00 0.00 -0.219 na
Hetatm 6 C lig 1 -11.994 15.348 73.608 1.00 0.00 0.221 C
Hetatm 7 C lig 1 -12.341 15.781 74.943 1.00 0.00 -0.097 C
Hetatm 8 C lig 1 -12.628 16.141 76.055 1.00 0.00 0.014 C
HITATM 9 N LIG 1 -16.387 15.490 71.145 1.00 0.00 -0.095 N
氦氣發生器10 O型連接1 -17.127 14.595 70.751 1.00 0.00 -0.530 OA
HETATM 11 O LIG 1 -16.631 16.674 71.082 1.00 0.00 -0.530 OA
uj5u.com熱心網友回復:
假設檔案具有相同的長度,你可以用paste合并它們。然后按所需的順序提取列:
paste file1.txt file2.txt|awk '{print 1, 2, 3, 4, 5, $18, $19, $20, $9, 10, $11, $12}' /span>
轉載請註明出處,本文鏈接:https://www.uj5u.com/gongcheng/327158.html
標籤:

