Departure time on tripplanner bus steps
[busui.git] / include / common-geo.inc.php
blob:a/include/common-geo.inc.php -> blob:b/include/common-geo.inc.php
<?php <?php
   
// SELECT array_to_string(array(SELECT REPLACE(name_2006, ',', '\,') as name FROM suburbs order by name), ',') // SELECT array_to_string(array(SELECT REPLACE(name_2006, ',', '\,') as name FROM suburbs order by name), ',')
$suburbs = explode(",", "Acton,Ainslie,Amaroo,Aranda,Banks,Barton,Belconnen,Bonner,Bonython,Braddon,Bruce,Calwell,Campbell,Chapman,Charnwood,Chifley,Chisholm,City,Conder,Cook,Curtin,Deakin,Dickson,Downer,Duffy,Dunlop,Evatt,Fadden,Farrer,Fisher,Florey,Flynn,Forrest,Franklin,Fraser,Fyshwick,Garran,Gilmore,Giralang,Gordon,Gowrie,Greenway,Griffith,Gungahlin,Hackett,Hall,Harrison,Hawker,Higgins,Holder,Holt,Hughes,Hume,Isaacs,Isabella Plains,Kaleen,Kambah,Kingston,Latham,Lawson,Lyneham,Lyons,Macarthur,Macgregor,Macquarie,Mawson,McKellar,Melba,Mitchell,Monash,Narrabundah,Ngunnawal,Nicholls,Oaks Estate,O'Connor,O'Malley,Oxley,Page,Palmerston,Parkes,Pearce,Phillip,Pialligo,Red Hill,Reid,Richardson,Rivett,Russell,Scullin,Spence,Stirling,Symonston,Tharwa,Theodore,Torrens,Turner,Wanniassa,Waramanga,Watson,Weetangera,Weston,Yarralumla"); $suburbs = explode(",", "Acton,Ainslie,Amaroo,Aranda,Banks,Barton,Belconnen,Bonner,Bonython,Braddon,Bruce,Calwell,Campbell,Chapman,Charnwood,Chifley,Chisholm,City,Conder,Cook,Curtin,Deakin,Dickson,Downer,Duffy,Dunlop,Evatt,Fadden,Farrer,Fisher,Florey,Flynn,Forrest,Franklin,Fraser,Fyshwick,Garran,Gilmore,Giralang,Gordon,Gowrie,Greenway,Griffith,Gungahlin,Hackett,Hall,Harrison,Hawker,Higgins,Holder,Holt,Hughes,Hume,Isaacs,Isabella Plains,Kaleen,Kambah,Kingston,Latham,Lawson,Lyneham,Lyons,Macarthur,Macgregor,Macquarie,Mawson,McKellar,Melba,Mitchell,Monash,Narrabundah,Ngunnawal,Nicholls,Oaks Estate,O'Connor,O'Malley,Oxley,Page,Palmerston,Parkes,Pearce,Phillip,Pialligo,Red Hill,Reid,Richardson,Rivett,Russell,Scullin,Spence,Stirling,Symonston,Tharwa,Theodore,Torrens,Turner,Wanniassa,Waramanga,Watson,Weetangera,Weston,Yarralumla");
function staticmap($mapPoints, $zoom = 0, $markerImage = "iconb", $collapsible = true)  
{ function staticmap($mapPoints, $zoom = 0, $markerImage = "iconb", $collapsible = true, $twotone = false) {
global $labsPath; global $basePath;
$width = 300; $width = 300;
$height = 300; $height = 300;
$metersperpixel[9] = 305.492 * $width; $metersperpixel[9] = 305.492 * $width;
$metersperpixel[10] = 152.746 * $width; $metersperpixel[10] = 152.746 * $width;
$metersperpixel[11] = 76.373 * $width; $metersperpixel[11] = 76.373 * $width;
$metersperpixel[12] = 38.187 * $width; $metersperpixel[12] = 38.187 * $width;
$metersperpixel[13] = 19.093 * $width; $metersperpixel[13] = 19.093 * $width;
$metersperpixel[14] = 9.547 * $width; $metersperpixel[14] = 9.547 * $width;
$metersperpixel[15] = 4.773 * $width; $metersperpixel[15] = 4.773 * $width;
//$metersperpixel[16] = 2.387 * $width; //$metersperpixel[16] = 2.387 * $width;
// $metersperpixel[17]=1.193*$width; // $metersperpixel[17]=1.193*$width;
$center = ""; $center = "";
$markers = ""; $markers = "";
$mapwidthinmeters = 50; $mapwidthinmeters = 50;
if (sizeof($mapPoints) < 1) return "map error"; if (sizeof($mapPoints) < 1)
if (sizeof($mapPoints) === 1) { return "map error";
if ($zoom == 0) $zoom = 14; if (sizeof($mapPoints) === 1) {
$markers.= "{$mapPoints[0][0]},{$mapPoints[0][1]},$markerimage"; if ($zoom == 0)
$center = "{$mapPoints[0][0]},{$mapPoints[0][1]}"; $zoom = 14;
} $markers.= "{$mapPoints[0][0]},{$mapPoints[0][1]},$markerimage";
else { $center = "{$mapPoints[0][0]},{$mapPoints[0][1]}";
foreach ($mapPoints as $index => $mapPoint) { }
$markers.= $mapPoint[0] . "," . $mapPoint[1] . "," . $markerImage . ($index + 1); else {
if ($index + 1 != sizeof($mapPoints)) $markers.= "|"; foreach ($mapPoints as $index => $mapPoint) {
$dist = distance($mapPoints[0][0], $mapPoint[0][1],$mapPoint[0], $mapPoint[1]); if ($twotone && $index == 0) {
$mapwidthinmeters = ($dist>$mapwidthinmeters ? $dist : $mapwidthinmeters); $markers.= $mapPoint[0] . "," . $mapPoint[1] . "," . "iconr" . ($index + 1);
$totalLat+= $mapPoint[0]; $center = "{$mapPoints[0][0]},{$mapPoints[0][1]}";
$totalLon+= $mapPoint[1]; } else {
} $markers.= $mapPoint[0] . "," . $mapPoint[1] . "," . $markerImage . ($index + 1);
if ($zoom == 0) { }
$mapwidthinmeters = distance($minlat, $minlon, $minlat, $maxlon); if ($index + 1 != sizeof($mapPoints))
foreach (array_reverse($metersperpixel, true) as $zoomLevel => $maxdistance) { $markers.= "|";
if ($zoom == 0 && $mapwidthinmeters*1.5 < ($maxdistance)) $zoom = $zoomLevel; $dist = distance($mapPoints[0][0], $mapPoint[0][1], $mapPoint[0], $mapPoint[1]);
} $mapwidthinmeters = ($dist > $mapwidthinmeters ? $dist : $mapwidthinmeters);
} $totalLat+= $mapPoint[0];
$center = $totalLat / sizeof($mapPoints) . "," . $totalLon / sizeof($mapPoints); $totalLon+= $mapPoint[1];
} }
$output = ""; if ($zoom == 0) {
if ($collapsible) $output.= '<div class="map" data-role="collapsible" data-collapsed="true"><h3>Open Map...</h3>'; $mapwidthinmeters = distance($minlat, $minlon, $minlat, $maxlon);
$output.= '<img class="map" src="' . curPageURL() . '/'. $labsPath. '/lib/staticmaplite/staticmap.php?center=' . $center . '&amp;zoom=' . $zoom . '&amp;size=' . $width . 'x' . $height . '&amp;markers=' . foreach (array_reverse($metersperpixel, true) as $zoomLevel => $maxdistance) {
$markers . '" width=' . $width . ' height=' . $height . '>'; if ($zoom == 0 && $mapwidthinmeters * 1.5 < ($maxdistance))
if ($collapsible) $output.= '</div>'; $zoom = $zoomLevel;
return $output; }
} }
function distance($lat1, $lng1, $lat2, $lng2, $roundLargeValues = false) $center = $totalLat / sizeof($mapPoints) . "," . $totalLon / sizeof($mapPoints);
{ }
$pi80 = M_PI / 180; $output = "";
$lat1*= $pi80; if ($collapsible)
$lng1*= $pi80; $output.= '<div class="map" data-role="collapsible" data-collapsed="true"><h3>Open Map...</h3>';
$lat2*= $pi80; $output.= '<img class="map" src="' . curPageURL() . '/' . $basePath . '/lib/staticmaplite/staticmap.php?center=' . $center . '&amp;zoom=' . $zoom . '&amp;size=' . $width . 'x' . $height . '&amp;markers=' . $markers . '" width=' . $width . ' height=' . $height . '>';
$lng2*= $pi80; if ($collapsible)
$r = 6372.797; // mean radius of Earth in km $output.= '</div>';
$dlat = $lat2 - $lat1; return $output;
$dlng = $lng2 - $lng1;  
$a = sin($dlat / 2) * sin($dlat / 2) + cos($lat1) * cos($lat2) * sin($dlng / 2) * sin($dlng / 2);  
$c = 2 * atan2(sqrt($a) , sqrt(1 - $a));  
$km = $r * $c;  
if ($roundLargeValues) {  
if ($km < 1) return floor($km * 1000);  
else return round($km,2)."k";  
} else return floor($km * 1000);  
} }
   
function decodePolylineToArray($encoded) function distance($lat1, $lng1, $lat2, $lng2, $roundLargeValues = false) {
{ $pi80 = M_PI / 180;
// source: http://latlongeeks.com/forum/viewtopic.php?f=4&t=5 $lat1*= $pi80;
$length = strlen($encoded); $lng1*= $pi80;
$index = 0; $lat2*= $pi80;
$points = array(); $lng2*= $pi80;
$lat = 0; $r = 6372.797; // mean radius of Earth in km
$lng = 0; $dlat = $lat2 - $lat1;
while ($index < $length) { $dlng = $lng2 - $lng1;
// Temporary variable to hold each ASCII byte. $a = sin($dlat / 2) * sin($dlat / 2) + cos($lat1) * cos($lat2) * sin($dlng / 2) * sin($dlng / 2);
$b = 0; $c = 2 * atan2(sqrt($a), sqrt(1 - $a));
// The encoded polyline consists of a latitude value followed by a $km = $r * $c;
// longitude value. They should always come in pairs. Read the if ($roundLargeValues) {
// latitude value first. if ($km < 1)
$shift = 0; return floor($km * 1000);
$result = 0; else
do { return round($km, 2) . "k";
// The `ord(substr($encoded, $index++))` statement returns the ASCII }
// code for the character at $index. Subtract 63 to get the original else
// value. (63 was added to ensure proper ASCII characters are displayed return floor($km * 1000);
// in the encoded polyline string, which is `human` readable)  
$b = ord(substr($encoded, $index++)) - 63;  
// AND the bits of the byte with 0x1f to get the original 5-bit `chunk.  
// Then left shift the bits by the required amount, which increases  
// by 5 bits each time.  
// OR the value into $results, which sums up the individual 5-bit chunks  
// into the original value. Since the 5-bit chunks were reversed in  
// order during encoding, reading them in this way ensures proper  
// summation.  
$result|= ($b & 0x1f) << $shift;  
$shift+= 5;  
}  
// Continue while the read byte is >= 0x20 since the last `chunk`  
// was not OR'd with 0x20 during the conversion process. (Signals the end)  
while ($b >= 0x20);  
// Check if negative, and convert. (All negative values have the last bit  
// set)  
$dlat = (($result & 1) ? ~($result >> 1) : ($result >> 1));  
// Compute actual latitude since value is offset from previous value.  
$lat+= $dlat;  
// The next values will correspond to the longitude for this point.  
$shift = 0;  
$result = 0;  
do {  
$b = ord(substr($encoded, $index++)) - 63;  
$result|= ($b & 0x1f) << $shift;  
$shift+= 5;  
} while ($b >= 0x20);  
$dlng = (($result & 1) ? ~($result >> 1) : ($result >> 1));  
$lng+= $dlng;  
// The actual latitude and longitude values were multiplied by  
// 1e5 before encoding so that they could be converted to a 32-bit  
// integer representation. (With a decimal accuracy of 5 places)  
// Convert back to original values.  
$points[] = array(  
$lat * 1e-5,  
$lng * 1e-5  
);  
}  
return $points;  
} }
function geocode($query, $giveOptions)  
{ function decodePolylineToArray($encoded) {
global $cloudmadeAPIkey; // source: http://latlongeeks.com/forum/viewtopic.php?f=4&t=5
$url = "http://geocoding.cloudmade.com/$cloudmadeAPIkey/geocoding/v2/find.js?query=" . urlencode($query) . "&bbox=-35.5,149.00,-35.15,149.1930&return_location=true&bbox_only=true"; $length = strlen($encoded);
$contents = json_decode(getPage($url)); $index = 0;
if ($giveOptions) return $contents->features; $points = array();
elseif (isset($contents->features[0]->centroid)) return $contents->features[0]->centroid->coordinates[0] . "," . $contents->features[0]->centroid->coordinates[1]; $lat = 0;
else return ""; $lng = 0;
  while ($index < $length) {
  // Temporary variable to hold each ASCII byte.
  $b = 0;
  // The encoded polyline consists of a latitude value followed by a
  // longitude value. They should always come in pairs. Read the
  // latitude value first.
  $shift = 0;
  $result = 0;
  do {
  // The `ord(substr($encoded, $index++))` statement returns the ASCII
  // code for the character at $index. Subtract 63 to get the original
  // value. (63 was added to ensure proper ASCII characters are displayed
  // in the encoded polyline string, which is `human` readable)
  $b = ord(substr($encoded, $index++)) - 63;
  // AND the bits of the byte with 0x1f to get the original 5-bit `chunk.
  // Then left shift the bits by the required amount, which increases
  // by 5 bits each time.
  // OR the value into $results, which sums up the individual 5-bit chunks
  // into the original value. Since the 5-bit chunks were reversed in
  // order during encoding, reading them in this way ensures proper
  // summation.
  $result|= ($b & 0x1f) << $shift;
  $shift+= 5;
  }
  // Continue while the read byte is >= 0x20 since the last `chunk`
  // was not OR'd with 0x20 during the conversion process. (Signals the end)
  while ($b >= 0x20);
  // Check if negative, and convert. (All negative values have the last bit
  // set)
  $dlat = (($result & 1) ? ~($result >> 1) : ($result >> 1));
  // Compute actual latitude since value is offset from previous value.
  $lat+= $dlat;
  // The next values will correspond to the longitude for this point.
  $shift = 0;
  $result = 0;
  do {
  $b = ord(substr($encoded, $index++)) - 63;
  $result|= ($b & 0x1f) << $shift;
  $shift+= 5;
  } while ($b >= 0x20);
  $dlng = (($result & 1) ? ~($result >> 1) : ($result >> 1));
  $lng+= $dlng;
  // The actual latitude and longitude values were multiplied by
  // 1e5 before encoding so that they could be converted to a 32-bit
  // integer representation. (With a decimal accuracy of 5 places)
  // Convert back to original values.
  $points[] = array(
  $lat * 1e-5,
  $lng * 1e-5
  );
  }
  return $points;
} }
function reverseGeocode($lat, $lng)  
{ function geocode($query, $giveOptions) {
global $cloudmadeAPIkey; global $cloudmadeAPIkey;
$url = "http://geocoding.cloudmade.com/$cloudmadeAPIkey/geocoding/v2/find.js?around=" . $lat . "," . $lng . "&distance=closest&object_type=road"; $url = "http://geocoding.cloudmade.com/$cloudmadeAPIkey/geocoding/v2/find.js?query=" . urlencode($query) . "&bbox=-35.5,149.00,-35.15,149.1930&return_location=true&bbox_only=true";
$contents = json_decode(getPage($url)); $contents = json_decode(getPage($url));
return $contents->features[0]->properties->name; if ($giveOptions)
  return $contents->features;
  elseif (isset($contents->features[0]->centroid))
  return $contents->features[0]->centroid->coordinates[0] . "," . $contents->features[0]->centroid->coordinates[1];
  else
  return "";
} }
   
  function reverseGeocode($lat, $lng) {
  global $cloudmadeAPIkey;
  $url = "http://geocoding.cloudmade.com/$cloudmadeAPIkey/geocoding/v2/find.js?around=" . $lat . "," . $lng . "&distance=closest&object_type=road";
  $contents = json_decode(getPage($url));
  return $contents->features[0]->properties->name;
  }
   
?> ?>