pjacklam / p5-geo-ellipsoid Goto Github PK

--- Ellipsoid.pm.orig	Tue Jan 19 15:33:31 2021
+++ Ellipsoid.pm	Thu Jan 28 00:14:49 2021
@@ -775,15 +775,36 @@
 sub at
 {
   my $self = shift;
+  my( $lat, $lon, $r, $az ) = @_;
+  print "at($lat,$lon,$r,$az)\n" if $DEBUG;
+  my( $lat2, $lon2, $baz ) = $self->atb($lat,$lon,$r,$az);
+  print "atb returns ($lat2,$lon2,$baz)\n" if $DEBUG;
+  return ( $lat2, $lon2 );
+}
+
+=pod
+
+=item atb
+
+Returns the list (latitude,longitude,reverse_bearing) in degrees or radians that is a
+specified range and bearing from a given location.
+
+    my( $lat2, $lon2, $baz ) = $geo->atb( $lat1, $lon1, $range, $bearing );
+
+=cut
+
+sub atb
+{
+  my $self = shift;
   my $units = $self->{angle_unit};
   my( $lat, $lon, $az ) = _normalize_input($units,@_[0,1,3]);
   my $r = $_[2];
-  print "at($lat,$lon,$r,$az)\n" if $DEBUG;
-  my( $lat2, $lon2 ) = $self->_forward($lat,$lon,$r,$az);
-  print "_forward returns ($lat2,$lon2)\n" if $DEBUG;
+  print "atb($lat,$lon,$r,$az)\n" if $DEBUG;
+  my( $lat2, $lon2, $baz ) = $self->_forward($lat,$lon,$r,$az);
+  print "_forward returns ($lat2,$lon2,$baz)\n" if $DEBUG;
   $self->_normalize_output('longitude_symmetric',$lon2);
   $self->_normalize_output('latitude_symmetric',$lat2);
-  return ( $lat2, $lon2 );
+  return ( $lat2, $lon2, $baz );
 }
 
 =pod
@@ -1052,10 +1073,10 @@
   $c = (((((-3.0*$c2a)+4.0)*$f)+4.0)*$c2a*$f)/16.0;
   $d = (((($e*$cy*$c) + $cz)*$sy*$c)+$y)*$sa;
   my $lon2 = $lon1 + $x - (1.0-$c)*$d*$f;
-  #$baz = atan2($sa,$baz) + pi;
+  $baz = atan2($sa,$baz) + pi;
 
   # return result
-  return ($lat2,$lon2);
+  return ($lat2,$lon2,$baz);
 
 }