perl脚本遍历文件夹及对SAC波形文件进行去仪器响应、滤波等操作实现

ゝ一纸荒年。 2022-04-16 02:27 394阅读 0赞

1. 如题，perl脚本遍历子文件夹及对子文件夹下文件处理：

#!/usr/bin/env perl
    #遍历文件夹$rawpath下的所有子目录，并对子目录下的文件进行操作；
    
    my $rawpath="/home/liuyi/evt_test";
    
    opendir (DIR,$rawpath) or die "can't open!";
    my @file= grep { /^[^\.]/ } readdir DIR;#过滤掉以"."和".."的文件,即UNIX下的隐藏文件
    foreach my $file (@file){
    
        my $newdir = $rawpath."/".$file;
        opendir(DIR2,$newdir) or die "can't open $newdir!";      
        my @file2 = grep { /^[^\.]/ } readdir DIR2;
        foreach my $file2 (@file2){
        
        #对文件$file进行操作；
        
        }
    }

2. 以对地震事件SAC波形文件与处理为例，实现对原始波形去毛刺，去均值、去线性趋势和波形尖灭，以及去仪器响应、带通滤波等操作：

#!/usr/bin/env perl
    #作者:刘毅;邮箱:825948184@qq.com;日期:2018/11/15;
    
    my $rawpath="/run/media/lisonglin/4aacb9ca-62ed-4102-b56c-cc7178021277/evt";
    
    opendir (DIR,$rawpath) or die "can't open!";
    #my @file = readdir DIR;
    my @file= grep { /^[^\.]/ } readdir DIR;#过滤掉以"."和".."的文件,即UNIX下的隐藏文件
    foreach my $file (@file){
    
        my $newdir = $rawpath."/".$file;
        print "$newdir\n";
        opendir(DIR2,$newdir) or die "can't open $newdir!"; #注意打开一个目录和进入该目录是不同的概念；
        chdir $newdir; 
    	
    	my @file2 = grep { /^[^\.]/ } readdir DIR2;
        	foreach my $file2 (@file2){	
    	
            open(SAC, "| sac") or die "Error in opening SAC\n";
    	print SAC "wild echo off \n";
    
    	print SAC "r $file2 \n";
    	print "file2: $file2\n";、
    	#去毛刺；
    	print SAC "rglitches threshold 1.0e+10 type linear window off method absolute \n";
    	#去均值、去线性趋势和波形尖灭；
    	print SAC "rmean;rtrend;taper \n";
    	#去仪器响应，/home/lisonglin/liuyi/yn_ly/SAC_PZ为仪器响应文件；
    	print SAC "transfer from polezero subtype /home/lisonglin/liuyi/yn_ly/SAC_PZ to none freq 0.05 0.1 10.0 15.0 \n";
    	#带通滤波；
     	print SAC "bp n 4 c 1 10 \n";
    	#保存并重命名加.bp后缀；
    	print SAC "w append .bp \n";
     	#退出；
    	print SAC "q \n";
    
    	close(SAC);
    	}
        chdir "..";
        
    }

希望对你我有所帮助，在此记录发布。