- 论坛徽章:
- 1
|
[root@localhost ]# perl 1.pl
101.36658837807
[root@localhost ]# cat 1.pl
#!/bin/perl
#
my $JuLi= &distance(116.41794,39.99347,116.41787,39.99256);#$wd1, $jd1, $wd2, $jd2
print "$JuLi \n";
sub distance() {
my ($wd1, $jd1, $wd2, $jd2);
$wd1=$_[1];
$jd1=$_[0];
$wd2=$_[3];
$jd2=$_[2];
my $PI=3.1415926535898;
my $R = 6.371229 * 1e6;
#print "$R \n";
$x = ($jd2 - $jd1) * $PI * $R * cos( ( ($wd1 + $wd2) / 2) * $PI / 180) / 180;
$y = ($wd2 - $wd1) * $PI * $R / 180;
my $temp = sqrt($x * $x + $y * $y);
return $temp;
} |
|