X-Git-Url: http://dxcluster.net/gitweb/gitweb.cgi?a=blobdiff_plain;f=perl%2FDXBearing.pm;h=db6f4686d0a75ed89c32e0c32e28509ae233bb5a;hb=8106275008346b0cceb090f7317ff4240a3ba44f;hp=e7736f8a138d64a367865bce605631389f96f861;hpb=5bda7945001817ba3005d357b7b2f257ffe8a2cf;p=spider.git diff --git a/perl/DXBearing.pm b/perl/DXBearing.pm index e7736f8a..db6f4686 100644 --- a/perl/DXBearing.pm +++ b/perl/DXBearing.pm @@ -8,25 +8,19 @@ # # Copyright (c) 1998 - Dirk Koopman G1TLH # -# $Id$ +# # package DXBearing; -use POSIX; use DXUtil; +use POSIX qw(:math_h); use strict; use vars qw($pi); $pi = 3.14159265358979; -use vars qw($VERSION $BRANCH); -$VERSION = sprintf( "%d.%03d", q$Revision$ =~ /(\d+)\.(\d+)/ ); -$BRANCH = sprintf( "%d.%03d", q$Revision$ =~ /\d+\.\d+\.(\d+)\.(\d+)/ || (0,0)); -$main::build += $VERSION; -$main::branch += $BRANCH; - # convert a qra locator into lat/long in DEGREES sub qratoll { @@ -95,18 +89,32 @@ sub bdist my $e = dr(shift); return (0, 0) if $hn == $n && $he == $e; my $co = cos($he-$e)*cos($hn)*cos($n)+sin($hn)*sin($n); - my $ca = atan(abs(sqrt(1-$co*$co)/$co)); + my $ca = $co ? atan(abs(sqrt(1-$co*$co)/$co)) : $pi; $ca = $pi-$ca if $co < 0; my $dx = 6367*$ca; my $si = sin($e-$he)*cos($n)*cos($hn); $co = sin($n)-sin($hn)*cos($ca); - my $az = atan(abs($si/$co)); + my $az = $co ? atan(abs($si/$co)) : $pi; $az = $pi - $az if $co < 0; $az = -$az if $si < 0; $az = $az+2*$pi if $az < 0; return (rd($az), $dx); } +# just the distance - parameters as above +sub distance +{ + my $hn = dr(shift); + my $he = dr(shift); + my $n = dr(shift); + my $e = dr(shift); + return (0, 0) if $hn == $n && $he == $e; + my $co = cos($he-$e)*cos($hn)*cos($n)+sin($hn)*sin($n); + my $ca = $co ? atan(abs(sqrt(1-$co*$co)/$co)) : $pi; + $ca = $pi-$ca if $co < 0; + my $dx = 6367*$ca; +} + # turn a lat long string into floating point lat and long sub stoll {