App-SeismicUnixGui
view release on metacpan or search on metacpan
lib/App/SeismicUnixGui/big_streams/iSpectralAnalysis.pm view on Meta::CPAN
# print "iSpectralAnalysis, flow1= $flow[1]\n";
# $log->file($flow[1]);
# print "iSpectralAnalysis,flow2= $flow[2]\n";
# $log->file($flow[2]);
}
=head2 sub xtract
waveform to a file and QC the waveform
flow 3 saves a copy of the file to disk
=cut
sub xtract {
=head2 Prepare
to Extract waveform
=cut
$xtract->clear();
$xtract->ref_picks_file( \$picks_file );
$xtract->header_word($headerWord); # critical
$xtract->inbound( $inbound[1] );
$xtract->window_title( $file_in[1] );
$xtract->absclip($absclip);
$xtract_sufile[1] = $xtract->file_out( \$picks_file, $inbound[1] );
$ref_flows = $xtract->Step();
$xtract[1] = $$ref_flows[1];
$flow[3] = $xtract[1];
$run->flow( \$flow[3] );
=head2 LOG FLOW(S)
TO SCREEN AND FILE
=cut
print "$flow[3]\n";
#$log->file($flow[3]);
}
=head2 Spectral Analysis of Waveform
do not start flow 4 until
flow 3 is complete
flow 4 plots the extracted waveform
flow 5 plots the frequency analysis
flow 6 plots the phase analysis
We placed a cmd line
prompt that requires user to press
return on keyboard.
=cut
sub analyze {
=head2 Prepare
to Analyze extracted waveform
=cut
my $inbound = $DATA_SEISMIC_SU . '/' . $xtract_sufile[1];
$analyze->clear();
$analyze->inbound($inbound);
$analyze->window_title( $file_in[1] );
$analyze->absclip_phase($absclip_phase);
$analyze->absclip_freq($absclip_freq);
$ref_flows = $analyze->Step();
$analyze[1] = $$ref_flows[1];
$analyze[2] = $$ref_flows[2];
$analyze[3] = $$ref_flows[3];
$flow[4] = $analyze[1];
$flow[5] = $analyze[2];
$flow[6] = $analyze[3];
$run->flow( \$flow[4] );
$run->flow( \$flow[5] );
$run->flow( \$flow[6] );
=head2 LOG FLOW(S)
TO SCREEN AND FILE
=cut
# print "$flow[4]\n";
#$log->file($flow[4]);
# print "$flow[5]\n";
#$log->file($flow[5]);
# print "$flow[6]\n";
#$log->file($flow[6]);
#
}
( run in 0.916 second using v1.01-cache-2.11-cpan-6aa56a78535 )