
我有一個文件supermaster.PRM
:
num_valid_az = 12194
nrows = 12194
first_line = 1
deskew = n
caltone = 0.000000
st_rng_bin = 1
Flip_iq = n
offset_video = n
az_res = 0.000000
nlooks = 1
chirp_ext = 0
scnd_rng_mig = 0
rng_spec_wgt = 1.000000
rm_rng_band = 0.200000
rm_az_band = 0.000000
rshift = -18
ashift = 0
stretch_r = 0
stretch_a = 0
a_stretch_r = 0
a_stretch_a = 0
first_sample = 284
SC_identity = 10
rng_samp_rate = 64345238.125714
input_file = S1_20160114_ALL_F1.raw
num_rng_bins = 67752
bytes_per_line = 271008
good_bytes_per_line = 271008
PRF = 486.486310
pulse_dur = 5.240481e-05
near_range = 799926.599409
num_lines = 12194
num_patches = 1
SC_clock_start = 2016013.0823712260
SC_clock_stop = 2016013.0826613354
clock_start = 13.082371226227
clock_stop = 13.082661335643
led_file = S1_20160114_ALL_F1.LED
orbdir = A
lookdir = R
radar_wavelength = 0.0554658
chirp_slope = 1.07815e+12
rng_samp_rate = 64345238.125714
I_mean = 0
Q_mean = 0
SC_vel = 7160.742699
earth_radius = 6371038.614500
equatorial_radius = 6378137.000000
polar_radius = 6356752.310000
SC_height = 699860.307600
SC_height_start = 699988.203956
SC_height_end = 699732.953774
fd1 = 0.000000
fdd1 = 0.000000
fddd1 = 0.000000
sub_int_r = 0.000000
sub_int_a = 0.000000
SLC_file = S1_20160114_ALL_F1.SLC
dtype = a
SLC_scale = 1.000000
在這個 supermaster.PRM 檔案中,我想採用「 rng_samp_rate = 64345238.125714 」——我知道有兩個,但我只需要一個——將其放入這個 shell 腳本中:
#!/bin/bash
#SBATCH --job-name=InSAR
#SBATCH --ntasks=128
#SBATCH --time=7-00:00:00
#SBATCH --mail-type=fail,end
#SBATCH --export=all
#SBATCH --out=Forge-%j.out
#SBATCH --mem=256000
module load GMT
module load GMTSAR
module load parallel
sbas intf.tab scene.tab 19 10 8469 6097
csh plot_sbas.csh
到 sbas 行的末尾,例如
#!/bin/bash
#SBATCH --job-name=InSAR
#SBATCH --ntasks=128
#SBATCH --time=7-00:00:00
#SBATCH --mail-type=fail,end
#SBATCH --export=all
#SBATCH --out=Forge-%j.out
#SBATCH --mem=256000
module load GMT
module load GMTSAR
module load parallel
sbas intf.tab scene.tab 19 10 8469 6097 -rng_samp_rate = 64345238.125714
csh plot_sbas.csh
答案1
您可以使用這個awk
腳本:
awk 'NR==FNR && /^rng_samp_rate/ { value=$0 }
NR!=FNR && !/^sbas/ { print }
NR!=FNR && /^sbas/ { print $0 " -" value }' supermaster.PRM script.in > script.out
第一個文件必須是參數文件,第二個文件是要修改的shell腳本。 (此awk
腳本採用最後一次出現的rng_samp_rate
。)
輸出是
#!/bin/bash
#SBATCH --job-name=InSAR
#SBATCH --ntasks=128
#SBATCH --time=7-00:00:00
#SBATCH --mail-type=fail,end
#SBATCH --export=all
#SBATCH --out=Forge-%j.out
#SBATCH --mem=256000
module load GMT
module load GMTSAR
module load parallel
sbas intf.tab scene.tab 19 10 8469 6097 -rng_samp_rate = 64345238.125714
csh plot_sbas.csh
當然您也可以使用其他工具,例如grep
和sed
。
sed "s/^sbas.*/& -$(grep '^rng_samp_rate' supermaster.PRM|head -1)/" script.in > script.out
或者也許更容易理解一點
#! /bin/sh
VALUE="$(grep '^rng_samp_rate' supermaster.PRM|head -1)"
sed "s/^sbas.*/& -$VALUE/" script.in > script.out
(這裡我使用第一次出現的rng_samp_rate
。)
這兩個指令都假定沒有其他行以rng_samp_rate
or開頭sbas
,且參數檔總是包含一行rng_samp_rate
。此sed
腳本假定參數檔中的值行不包含/
.
編輯:如何修改腳本檔案而不是建立新的輸出檔案?
某些命令(例如sed
,支援就地編輯)可以用來修改腳本文件,而無需明確建立單獨的輸出文件。
sed -i .bak -e "some script" script.in
一般來說,您始終可以使用合適的命令來跟隨命令mv
,例如
awk 'some script' supermaster.PRM script.in > script.out && mv script.out script.in
確保僅當命令未指示錯誤&&
時才覆蓋原始腳本。awk
如果您想要進行更多修改,您可以擴充awk
或sed
腳本(最佳效能)或按順序執行多個腳本或連接到管道。 (有關詳細信息,您應該寫下附加要求。)
答案2
您可以使用以下腳本編寫此腳本ed
:
printf '%s\n' \
'/^sbas /s/$/ -/' \
'. r !grep ^rng_samp_rate supermaster.PRM | head -1' \
'-1,.j' \
'w' \
'q' | ed -s shellscript
這會變更以文字「sbas
」開頭的(第一)行,以附加空格和破折號。然後,它將 shell 命令的輸出讀取到下一行; shell 指令rng_samp_rate
在 supermaster.PRM 的行開頭搜尋文字“ ”,然後head -1
只取得第一個此類結果。然後將該輸出與前一行連接。然後我們將文件保存到磁碟並退出。
如果原始文字中的額外空格是一個問題,您可以透過稍微修改 shell 程式碼來壓縮它們:
...
grep ^rng_samp_rate supermaster.PRM | head -1 | sed 's/ */ /'
...
....它只是添加一個sed
用一個空格替換 1 個或多個空格的命令。如果等號後面可能出現額外的空格,請將 sed 指令改為:sed 's/ */ /g'
。