ファイルがあります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」をインポートしたいと思います。 - 2つがあることを知っていますが、1つだけが必要です。これを次のシェルスクリプトに入れます。
#!/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
最初のファイルはパラメータファイルでなければならず、2番目のファイルは変更するシェルスクリプトです。 (この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
またはで始まる他の行がなく、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
テキスト ""で始まる(最初の)行が空白とダッシュを追加するように変更されます。次に、シェルコマンドの出力を次の行に読み込みます。シェルコマンドはrng_samp_rate
supermaster.PRMの行の先頭からテキストを検索し、head -1
最初の結果のみを取得します。その後、この出力は前の行に関連付けられます。その後、ファイルをディスクに保存して終了します。
生テキストの余分なスペースが問題の場合は、シェルコードを少し変更して圧縮できます。
...
grep ^rng_samp_rate supermaster.PRM | head -1 | sed 's/ */ /'
...
sed
... 1つ以上のスペースを1つのスペースに置き換えるコマンドを追加するだけです。等号の後に余分なスペースが表示される場合は、sedコマンドを次のように変更しますsed 's/ */ /g'
。