Daar was ik inmiddels ook achter gekomen. De correcties introduceren alleen maar een fout "op het andere halfrond".Daar waren die correcties voor. Die heb ik eruit gehaald, helemaal geen correcties meer want de vakken tellen altijd omhoog staat elders in de tekst.
Enkele kleine verschillen met andere calculatie-sites maar dat zal een kwestie van afronding zijn.<?php
//**************************************************************************************************************
// function to calculate distance between two Maidenhead locator points
// started by PA0WV in C and PA3TON - some corrections and refinements by VK2WGK
// Version 0.2 - 12008-02-18
// Published under the terms of the GNU license
//*************************************************************************************************************
// define constants.
define( K , 111.2); // correction constant in the Maidenhead formula
define( R , M_PI/180); // factor for conversion of degrees in radians
define( BR , '<br />'); // avoids retyping a break in quotes every time (handy while debugging)
define( E , ' E ' ); // define the lon. and lat. indicators for easy formatting
define( W , ' W ' );
define( N , ' N ' );
define( S , ' S ' );
// input the local and remote qth locator (normally retrieved from _$POST or $_GET )
$qth_loc = 'QF67BF'; // this can also be defined as a constant when it is the home QTH locator.
$qth_rem = 'JO32hl';
// ********************************************* MAIN UNIT *****************************************************
// check if entered strings are of the proper format and length.
trim($qth_loc); trim($qth_rem); // trim any spaces and other unwanted chars away.
$qth_loc = strtoupper($qth_loc); $qth_rem = strtoupper($qth_rem); // convert both inputs to uppercase
if ( chk_error($qth_loc) ) { echo 'The local locator you entered contains errors or is too short. Please try again'; exit; }
if ( chk_error($qth_rem) ) { echo 'The remote locator you entered contains errors or is too short. Please try again'; exit; }
if ( strlen ($qth_loc) < 6 ) { $qth_loc= str_pad($qth_loc,6,'L');} // if string too short pad with L (middle of the square)
if ( strlen ($qth_rem) < 6 ) { $qth_rem= str_pad($qth_rem,6,'L');} // if string too short pad with L (middle of the square)
// format the strings for display and add the hemisphere indicator
$loclat_str = sprintf("%01.4f", get_latitude($qth_loc));
if ( $loclat_str < 0 ) { $loclat_str = abs($loclat_str) . S; } else { $loclat_str = $loclat_str . N; }
$remlat_str = sprintf("%01.4f", get_latitude($qth_rem));
if ( $remlat_str < 0 ) { $remlat_str = abs($remlat_str) . S; } else { $remlat_str = $remlat_str . N; }
$loclon_str = sprintf("%01.4f", get_longitude($qth_loc));
if ( $loclon_str < 0 ) { $loclon_str = abs($loclon_str) . W; } else { $loclon_str = $loclon_str . E; }
$remlon_str = sprintf("%01.4f", get_longitude($qth_rem));
if ( $remlon_str < 0 ) { $remlon_str = abs($remlon_str) . W; } else { $remlon_str = $remlon_str . E; }
// display the results
echo '----------------------------------------------------------------------------<br />';
echo 'The distance between: ' . $qth_loc . ' and ' . $qth_rem . ' is: '
. sprintf("%01.1f", distance($qth_loc, $qth_rem)) . ' km.' . BR;
echo 'The position of ' . $qth_loc . ' is latitude: ' . $loclat_str
. ' - longitude: ' . $loclon_str . BR;
echo 'The position of ' . $qth_rem . ' is latitude: ' . $remlat_str
. ' - longitude: ' . $remlon_str . BR;
echo '-----------------------------------------------------------------------------<br />';
//*************************** Function section *****************************************************************
// main function for calculating the distance
function distance ($qth_loc, $qth_rem ) {
// get the paramaters for the distance calculation - both retrieved in degrees, decimal.
$lat_loc = get_latitude ($qth_loc);
$lat_rem = get_latitude ($qth_rem);
$lon_loc = get_longitude ($qth_loc);
$lon_rem = get_longitude ($qth_rem);
// actual calculation. The longitude and latitude are converted into radians by multiplying with the constant R
$distance= K *(1/R) * acos ( (sin(R*$lat_loc) * sin(R*$lat_rem)) + ( cos(R*$lat_loc) * cos(R*$lat_rem)) * cos( (R*$lon_loc) - (R*$lon_rem)) );
// ready to go back
return $distance;
} // end of function distance_calculator
// check for errors in the locator string
function chk_error ($loc_str) {
if ( strlen ($loc_str) < 4 ) { return TRUE; }
if ( (ord($loc_str[0]) < ord('A') ) || (ord($loc_str[1]) < ord('A')) ) { return TRUE; }
if ( ($loc_str[2] < 0 ) || ( $loc_str[3] < 0 ) ) { return TRUE; }
if ( (ord($loc_str[4]) < ord('A')) || (ord($loc_str[5]) < ord('A')) ) { return TRUE; }
if ( (ord($loc_str[0]) > ord('R') ) || (ord($loc_str[1]) > ord('R')) ) { return TRUE; }
if ( ($loc_str[2] > 9) || ($loc_str[3] > 9 ) ) { return TRUE; }
if ( (ord($loc_str[4]) > ord('X')) || (ord($loc_str[5]) > ord('X')) ) { return TRUE; }
} //end of function chk_errors
// convert the locator information into a longitude
function get_longitude ( $loc_str) {
$c=-180+((ord($loc_str[0])-ord("A"))*20);
$c+=($loc_str[2]*2);
$c+=(ord($loc_str[4])-ord("A"))/12;
$longitude=$c;
return $longitude;
}// end of function get_longitude
// convert the locator information into a latitude
function get_latitude ( $loc_str ) {
$c=-90 +((ord($loc_str[1])-ord("A"))*10);
$c += $loc_str[3];
$c+=(ord($loc_str[5])-ord("A"))/24;
$latitude=$c;
return $latitude;
} // end of function get_latitude
?>
Mee eens, hoewel afstand ook wel leuk is om te weten.Overigens lijkt me een bearing programma nuttiger, dat de stations weten als ze de qra locator van de ander horen welke kant hun beam opgezet moet worden.
Dus die "midden correctie" moet ik er nog in stoppen. De source (trouwens dezelfde tekst als hier op het forum is gepost) kun je downloaden door rechts te klikken op de link en dan "Save Target As" te kiezen - of hoe dat dan ook op z'n hollands mag heten... Klopt, hier heb ik niks meer aan toetevoegen.@PA3TON: Ik dacht dat ik het al aardig "opgeleukt" had (error correctie, trimming, padding, output formatting enz.) Ook zijn de belangrijkste dingen in functies gestopt.
Geen dank, maar ik ben maar een PA3, geen PA0 en dus niet de beroemde PA0GMM.PA0WV schreef:
Rest me PA0GMM te danken voor zijn tip over de bolcosinusregel in het vademecum en de uitleg daar van de structuur van de qth locator.
73
Wim
PAoWV
