利用perl脚本实现SAC文件头段变量写入

最新推荐文章于 2024-08-25 22:59:34 发布

little_sunlight

最新推荐文章于 2024-08-25 22:59:34 发布

阅读量1.1k

点赞数 1

分类专栏： perl 文章标签： perl脚本 SAC 写入SAC头段变量

本文链接：https://blog.csdn.net/qq_34538534/article/details/83148316

版权

perl 专栏收录该内容

2 篇文章 0 订阅

订阅专栏

将"ttime.dat"里面的t1(最后一列)写到"./001"的所有"*.sac"中, 并标记为t1, perl脚本如下：

#!/usr/bin/env perl
#by liuyi

open(PICK,"<","ttime.dat")||die"cannot open the file:$!\n";

$i=0;
while(<PICK>){
	@lines=split(" ");
	$t1[$i]=@lines[3];
	$i=$i+1;
}
$size=$i;

$j=0;
my $DIR_PATH="/home/*/001";
opendir DIR, ${DIR_PATH}|| die "Can not open dir \n";
@filelist=readdir DIR;
foreach $file (@filelist){

	open(SAC, "| sac") or die "Error in opening SAC\n";
	print SAC "wild echo off \n";
 
	print SAC "r $DIR_PATH/$file \n";
	
	print SAC "ch t1 $t1[$j] \n";
	print SAC "wh \n";
	print SAC "q \n";
 
	close(SAC);
	$j=$j+1;
}

以上，在此做一记录。