view release on metacpan or  search on metacpan

lib/App/SeismicUnixGui/big_streams/iSuvelan.pm  view on Meta::CPAN

		#	my $new_dt_s = $s2ms * $iSuvelan->{_dt_s};

		my $dx_major_divisions = (
			$iSuvelan->{_first_velocity} + (
				$iSuvelan->{_number_of_velocities} *
				  $iSuvelan->{_velocity_increment}
			)
		) / 10;

#	print("iSuvelan,dx_major_divisions:$dx_major_divisions\n");
#	print("iSuvelan,time_inc_major: $time_inc_major\n");
#	print("iSuvelan,number_minor_time_divisions:$number_minor_time_divisions \n");
#	print("iSuvelan,tmax_s $iSuvelan->{_tmax_s}\n");
#	print("iSuvelan,new_dt_s $new_dt_s (ms)\n");

		$N = 2;
		$windowtitle[1] = '\('
		  . $N . '\)\ '
		  . $iSuvelan->{_sufile_in}
		  . '\ CDP=\ '
		  . $iSuvelan->{_cdp_num};

		#units=Semblance					\\
		$suximage->clear();
		$suximage->box_width( quotemeta(300) );
		$suximage->box_height( quotemeta(450) );
		$suximage->box_X0( quotemeta(600) );
		$suximage->box_Y0( quotemeta(0) );
		$suximage->title( quotemeta( $iSuvelan->{_cdp_num_suffix} ) );
		$suximage->windowtitle( quotemeta('Semblance') );
		$suximage->ylabel( quotemeta('TWTT (s)') );
		$suximage->xlabel( quotemeta('Velocity (m/s)') );
		$suximage->legend( quotemeta($on) );
		$suximage->cmap( quotemeta('hsv1') );
		$suximage->first_x( quotemeta( $iSuvelan->{_first_velocity} ) );
		$suximage->dx( quotemeta( $iSuvelan->{_velocity_increment} ) );
		$suximage->dt_s(quotemeta($new_dt_s));
		$suximage->loclip( quotemeta( $iSuvelan->{_min_semblance} ) );
		$suximage->hiclip( quotemeta( $iSuvelan->{_max_semblance} ) );
		$suximage->verbose( quotemeta($off) );

		#print ("d2num is ($iSuvelan->{_velocity_increment})\n\n");
		$suximage->dx_major_divisions( quotemeta($dx_major_divisions) );
		$suximage->dy_minor_divisions(
			quotemeta($number_minor_time_divisions) );
		$suximage->dy_major_divisions( quotemeta($time_inc_major) );

		$suximage->tend_s( quotemeta( $iSuvelan->{_tmax_s} ) );

		#$suximage -> percent4clip(quotemeta(95.0);
		$suximage->first_tick_number_x(
			quotemeta( $iSuvelan->{_first_velocity} ) );
		$suximage->picks( $iSuvelan->{_Tvel_outbound} );

print("iSuvelan, calcNdisplay: Writing picks to $iSuvelan->{_Tvel_outbound}\n\n");

=head2 conditions
 
 when number_of_tries is >=2 
 there should be a pre-exisiting digitized
 overlay curve to plot as well

=cut

		if ( $iSuvelan->{_number_of_tries} >= 2 ) {

			print("using a curve file:\n");
			print("\t$iSuvelan->{_Tvel_inbound}\n\n");
			$suximage->curvefile( quotemeta( $iSuvelan->{_Tvel_inbound} ) )
			  ;
			my ( $ref_T_nmo, $ref_Vnmo, $num_tvel_pairs ) =
			$manage_files_by->read_2cols( \$iSuvelan->{_Tvel_inbound} );
			$suximage->npair( quotemeta($num_tvel_pairs) );
			$suximage->curvecolor( quotemeta(2) );
		}

		$suximage[1] = $suximage->Step();

		# print("$suximage[1]\n\n");

		$N = 1;

=head2  set suxwigb parameters

  In the perl module for suxwigb we should
  have (but we do not yet) an explanation of each of
  these parameters

=cut

		$suxwigb->clear();
		$suxwigb->title( quotemeta( $iSuvelan->{_sufile_in} ) );
		$suxwigb->ylabel( quotemeta('Time s') );
		$suxwigb->xlabel( quotemeta('Offset m') );
		$suxwigb->box_width( quotemeta(300) );
		$suxwigb->box_height( quotemeta(450) );
		$suxwigb->box_X0( quotemeta(875) );
		$suxwigb->box_Y0( quotemeta(0) );
		$suxwigb->absclip( quotemeta(1) );
		$suxwigb->headerword( quotemeta('offset') );
		$suxwigb->xcur( quotemeta(2) );
		$suxwigb->windowtitle( quotemeta( $iSuvelan->{_cdp_num_suffix} ) );

		#	$suxwigb->shading( quotemeta(1) );
		$suxwigb[1] = $suxwigb->Step();

=head2  DEFINE FLOW(S)

  in interactive mode:
  First time you see the image number_of_tries =0
  the suximage is not interactive ( uses '&') 
  Second, third, etc. times (number_of_tries >= 1)
  The image will halt the flow  ('wait'), to allow user to 
  make new picks

  On the other hand suxwigb does not change between the first
  and last attempt
		$sufilter[1], $to, $sugain[2],            $to,
	suwind has currently no gain
=cut