roi_prep.pl

StaMps最新测试版

Use_rsc "$outfile write HEIGHT                   $Height_poly[0]";
Doc_rsc(
 RSC_Tip => 'Platform Altitude at Peg Point',
 RSC_Doc => q[
   "Platform Altitude" based on RDF usage in diffnsim.pl

   First coefficient (constant term) of orbit "Vertical Fit" polynomial
   output by get_peg_info run by GetPeg.pl.

   Polynomial is function of SCH 'S' coordinate.
   Value is in 'H' direction, which is height above SCH reference sphere.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter',
);

Use_rsc "$outfile write HEIGHT_DS                $Height_poly[1]";
Doc_rsc(
 RSC_Tip => 'Platform Altitude Rate at Peg Point',
 RSC_Doc => q[
   "Platform Altitude Rate" based on RDF usage in diffnsim.pl

   Second coefficient (linear term) of orbit "Vertical Fit" polynomial
   output by get_peg_info run by GetPeg.pl.

   Polynomial is function of SCH 'S' coordinate.
   Value is in 'H' direction, which is height above SCH reference sphere.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:meter',
);

Use_rsc "$outfile write HEIGHT_DDS               $Height_poly[2]";
Doc_rsc(
 RSC_Tip => 'Platform Altitude Acceleration at Peg Point',
 RSC_Doc => q[
   "Platform Altitude Acceleration" based on RDF usage in diffnsim.pl

   Third coefficient (quadratic term) of orbit "Vertical Fit" polynomial
   output by get_peg_info run by GetPeg.pl.

   Polynomial is function of SCH 'S' coordinate.
   Value is in 'H' direction, which is height above SCH reference sphere.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:meter**2',
);

Use_rsc "$outfile write HEIGHT_DT                $HgtDt";
Doc_rsc(
 RSC_Tip => 'Platform Altitude change w.r.t. time at Peg Point',
 RSC_Derivation => q[
   $HgtDt = $Height_poly[1] * $VelocitySCH[0];
   ],
 RSC_Comment => q[
   Does not appear to ever be used.
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:second',
);

Use_rsc "$outfile write CROSSTRACK_POS           $CrossT_poly[0]";
Use_rsc "$outfile write CROSSTRACK_POS_DS        $CrossT_poly[1]";
Use_rsc "$outfile write CROSSTRACK_POS_DDS       $CrossT_poly[2]";
Use_rsc "$outfile write VELOCITY                 $velocity_mid";
Doc_rsc(
 RSC_Tip => 'Norm of Platform SCH Velocity at Peg Point',
 RSC_Doc => q[
   "Body fixed S/C velocities" based on RDF usage in roi_prep.pl and autofocus.pl.
   "Spacecraft Along Track Velocity" based on RDF usage in inverse3d.pl
   "Platform Velocity" based on RDF usage in diffnsim.pl and phase2base.pl
   ],
 RSC_Derivation => q[
   $velocity_mid = Norm( @VelocitySCH );
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:second',
);

Use_rsc "$outfile write VELOCITY_S               $VelocitySCH[0]";
Doc_rsc(
 RSC_Tip => 'Platform Velocity S Component',
 RSC_Doc => q[
   'S' Component of 'Platform SCH Velocity'
   produced by get_peg_info run by GetPeg.pl.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used by roi_prep.pl and autofocus.pl
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:second',
);

Use_rsc "$outfile write VELOCITY_C               $VelocitySCH[1]";
Doc_rsc(
 RSC_Tip => 'Platform Velocity C Component',
 RSC_Doc => q[
   'C' Component of 'Platform SCH Velocity'
   produced by get_peg_info run by GetPeg.pl.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used by roi_prep.pl and autofocus.pl
   Note - this is not the speed of light.
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:second',
);
Use_rsc "$outfile write VELOCITY_H               $VelocitySCH[2]";
Doc_rsc(
 RSC_Tip => 'Platform Velocity H Component',
 RSC_Doc => q[
   'H' Component of 'Platform SCH Velocity'
   produced by get_peg_info run by GetPeg.pl.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used by roi_prep.pl and autofocus.pl
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:second',
);

Use_rsc "$outfile write ACCELERATION_S           $AccelerationSCH[0]";
Doc_rsc(
 RSC_Tip => 'Platform Acceleration S Component',
 RSC_Doc => q[
   'S' Component of 'Platform SCH Acceleration'
   produced by get_peg_info run by GetPeg.pl.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used by roi_prep.pl and autofocus.pl
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:second**2',
);

Use_rsc "$outfile write ACCELERATION_C           $AccelerationSCH[1]";
Doc_rsc(
 RSC_Tip => 'Platform Acceleration C Component',
 RSC_Doc => q[
   'C' Component of 'Platform SCH Acceleration'
   produced by get_peg_info run by GetPeg.pl.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used by roi_prep.pl and autofocus.pl
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:second**2',
);

Use_rsc "$outfile write ACCELERATION_H           $AccelerationSCH[2]";
Doc_rsc(
 RSC_Tip => 'Platform Acceleration H Component',
 RSC_Doc => q[
   'H' Component of 'Platform SCH Acceleration'
   produced by get_peg_info run by GetPeg.pl.

   SCH Coordinate system.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used by roi_prep.pl and autofocus.pl
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:meter/SI:second**2',
);

Use_rsc "$outfile write VERT_VELOCITY            $Vert_V_poly[0]";
Use_rsc "$outfile write VERT_VELOCITY_DS         $Vert_V_poly[1]";
Use_rsc "$outfile write CROSSTRACK_VELOCITY      $CrossT_V_poly[0]";
Use_rsc "$outfile write CROSSTRACK_VELOCITY_DS   $CrossT_V_poly[1]";
Use_rsc "$outfile write ALONGTRACK_VELOCITY      $AlongT_V_poly[0]";
Use_rsc "$outfile write ALONGTRACK_VELOCITY_DS   $AlongT_V_poly[1]";
Use_rsc "$outfile write LATITUDE     $Latd";
Doc_rsc(
 RSC_Tip => 'Latitude of SCH Peg Point',
 RSC_Doc => q[
   'Lat' value of 'Peg Lat/Lon , H'
   produced by get_peg_info run by GetPeg.pl.

   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used in inverse3d.pl to specify "Peg Point Data"
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:degree',
);

Use_rsc "$outfile write LONGITUDE    $Lond";
Doc_rsc(
 RSC_Tip => 'Longitude of SCH Peg Point',
 RSC_Doc => q[
   'Lon' value of 'Peg Lat/Lon , H'
   produced by get_peg_info run by GetPeg.pl.

   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used in inverse3d.pl to specify "Peg Point Data"
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:degree',
);

Use_rsc "$outfile write HEADING      $hdgd";
Doc_rsc(
 RSC_Tip => 'Heading of SCH Peg Point',
 RSC_Doc => q[
   'Peg Heading' value
   produced by get_peg_info run by GetPeg.pl.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used in inverse3d.pl to specify "Peg Point Data"
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:degree',
);

Use_rsc "$outfile write PEG_UTC      $PegUtc";
Doc_rsc(
 RSC_Tip => 'Scene start time-of-day',
 RSC_Doc => q[
   'first' value of 'Time to first/middle scene'
   produced by get_peg_info run by GetPeg.pl.
   ],
 RSC_Derivation => q[
   See baseline/get_peg_info.f
   ],
 RSC_Comment => q[
   Appears to only be used in inverse3d.pl to calculate
     $PegLine = int( ( $PegUtc - $slc_first_line_utc ) / $delta_line_utc );
   which is used for
     'Reference Line for SCH Coordinates' RDF value
   ],
 RSC_Type => Real,
 RSC_Unit => 'SI:second',
);

#########################################
Message "Writing roi input file: ${prefix}_roi.in";
#########################################

if ($near_rng_ext < 0){ 
  $first_range_bin = 1-$near_rng_ext;
  $near_rng_ext=0;
}
else{
  $first_range_bin = 1;
}

open (ROI,">$roiin") or die "Can't write to $roiin\n";

$xmin = int(($xmin/2)+0.5);

print ROI <<END;

RDF File for ROI           
Master input file                                (-) = $infile                         ! Raw data for channel 1
Second input file                                (-) = /dev/null  		       ! Raw data for channel 2
Master SLC                                       (-) = $outfile                        ! SLC channel 1
Second SLC                                       (-) = /dev/null                       ! SLC channel 2
Eight look image file                            (-) = 8lk                             ! Eight look image file
Debug flag                                       (-) = 0                               ! 
Bytes per input line                             (-) = $width    $width                ! files 1 and 2
Good bytes per input line                        (-) = $xmax $xmax                     ! including header, files 1 and 2
First line to read                               (-) = $proc_start                           ! 
Number of patches                                (-) = $number_of_patches              ! org 40, 1 for testing
First sample pair to use                         (-) = $xmin                           ! 
Number of valid pulses                           (-) = $valid_samples                  ! 
Azimuth Patch Size                               (-) = $patch_size                     !
### For StaMPS AJH
Deskew?                                          (-) = y                               ! 
##################

Caltone location                                 (-) = 0 0                             !  % of sample rate, files 1 and 2
Start range bin, number to process               (-) = $first_range_bin $out_pixel     ! org 6626, test with 1600
Delta azimuth, range pixels for second file      (-) = 0 0                    ! 
Master file Doppler centroid coefs               (-) = $dop0 $dop1 $dop2 $dop3         ! (Hz/prf) quadratic polynomial
Second file Doppler centroid coefs               (-) = $dop0 $dop1 $dop2 $dop3         ! (Hz/prf) quadratic polynomial
Doppler average method                           (-) = 1                               ! 1: use file 1 dopp, 2: file 2, 3: avg
Earth radius                                     (m) = $earth_radius                   ! 
Body fixed S/C velocities                      (m/s) = $velocity $velocity             ! files 1, 2
Spacecraft height                                (m) = $height $height                 ! files 1, 2
Planet GM					 (-) = $PLANET_GM
Left, Right or Unknown Pointing 		 (-) = $AntennaSide
SCH Velocity Vector 1			       (m/s,m/s,m/s) = $VelocitySCH[0] $VelocitySCH[1] $VelocitySCH[2]
SCH Velocity Vector 2			       (m/s,m/s,m/s) = $VelocitySCH[0] $VelocitySCH[1] $VelocitySCH[2]
SCH Acceleration Vector 1		       (-,-,-) = $AccelerationSCH[0] $AccelerationSCH[1] $AccelerationSCH[2] 
SCH Acceleration Vector 2		       (-,-,-) = $AccelerationSCH[0] $AccelerationSCH[1] $AccelerationSCH[2]



Range of first sample in data file               (m) = $starting_rng $starting_rng     ! files 1, 2
PRF                                              (-) = $prf $prf                       ! pps, files 1, 2
i/q means                                        (-) = $i_bias $q_bias $i_bias $q_bias ! file 1 (i,q), file 2 (i,q)
Flip i/q?                                        (-) = n                               ! 
### For StaMPS AJH
Azimuth resolution                               (m) = 5                               !  
##################

Number of azimuth looks                          (-) = 4                               ! 
Range sampling rate                             (Hz) = $sampling_freq                  ! 
Range chirp slope                                (-) = $chirp_slope                    ! Hz/s
Range pulse duration                             (s) = $pulse_length                   ! 
Range chirp extension points                     (-) = $near_rng_ext                   ! 
Secondary range migration correction?            (-) = n                               ! 
Radar wavelength                                 (m) = $wavelength                     ! 
Range spectral weighting                         (-) = 1.                              ! 
Spectral shift fraction                          (-) = 0. 0.                           ! % of range/azimuth bandwidth to remove
Linear resampling coefs                          (-) = 0 0 0 0                         ! sloper,  intr,  slopea,  inta 
Linear resampling deltas                         (-) = 0 0 0 0                         ! dsloper, dintr, dslopea, dinta
AGC file                                         (-) = $infile.agcf                    ! AGC file name
DWP file                                         (-) = $infile.dwpf                    ! DWP file name
END
close(ROI);
###############################################################################
###############################################################################

exit 0;

=pod


=head1 USAGE

B<roi_prep.pl> I<date> OrbitType

date: radical of the input file I<date>.raw

Looks for I<date>.proc, then roi.proc in the next directory up, then uses default values

 defaults:
 before_z_ext = 0.8 * synthetic aperture
 after_z_ext  = 0.8 * synthetic aperture
 near_rng_ext = 1/2 chirp samples
 far_rng_ext  = 1/2 chirp samples
 patch_size   = 4096 (16384 for L-band)

Can also set number_of_patches, valid_samples (lines saved per patch) and out_pixel

=head1 FUNCTION

Creates input files for ROI

=head1 ROUTINES CALLED

none

=head1 CALLED BY

process.pl

=head1 FILES USED

I<date>.raw.rsc

roi.dop (generated by dopav.pl)

I<date>.proc [optional]

roi.proc [optional]


=head1 FILES CREATED

I<date>.slc.rsc

I<date>_roi.in

I<date>_roi.out

=head1 HISTORY

Shell Script : Francois ROGEZ 96/98

Perl  Script : Rowena LOHMAN 04/18/98

Frederic Crampe, Jan 10, 2000

Eric Fielding, Jan. 23, 2005

Eric Fielding, Nov. 29, 2006

=head1 LAST UPDATE

Eric Fielding, Aug. 22, 2007

=cut