Round distances to 2 decimal place kilometers
[busui.git] / common-geo.inc.php
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
<?php
// 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");
function staticmap($mapPoints, $zoom = 0, $markerImage = "iconb", $collapsible = true)
{
        $width = 300;
        $height = 300;
        $metersperpixel[9] = 305.492 * $width;
        $metersperpixel[10] = 152.746 * $width;
        $metersperpixel[11] = 76.373 * $width;
        $metersperpixel[12] = 38.187 * $width;
        $metersperpixel[13] = 19.093 * $width;
        $metersperpixel[14] = 9.547 * $width;
        $metersperpixel[15] = 4.773 * $width;
        $metersperpixel[16] = 2.387 * $width;
        // $metersperpixel[17]=1.193*$width;
        $center = "";
        $markers = "";
        $minlat = 999;
        $minlon = 999;
        $maxlat = 0;
        $maxlon = 0;
        if (sizeof($mapPoints) < 1) return "map error";
        if (sizeof($mapPoints) === 1) {
                if ($zoom == 0) $zoom = 14;
                $markers.= "{$mapPoints[0][0]},{$mapPoints[0][1]},$markerimage";
                $center = "{$mapPoints[0][0]},{$mapPoints[0][1]}";
        }
        else {
                foreach ($mapPoints as $index => $mapPoint) {
                        $markers.= $mapPoint[0] . "," . $mapPoint[1] . "," . $markerImage . ($index + 1);
                        if ($index + 1 != sizeof($mapPoints)) $markers.= "|";
                        if ($mapPoint[0] < $minlat) $minlat = $mapPoint[0];
                        if ($mapPoint[0] > $maxlat) $maxlat = $mapPoint[0];
                        if ($mapPoint[1] < $minlon) $minlon = $mapPoint[1];
                        if ($mapPoint[1] > $maxlon) $maxlon = $mapPoint[1];
                        $totalLat+= $mapPoint[0];
                        $totalLon+= $mapPoint[1];
                }
                if ($zoom == 0) {
                        $mapwidthinmeters = distance($minlat, $minlon, $minlat, $maxlon);
                        foreach (array_reverse($metersperpixel, true) as $zoomLevel => $maxdistance) {
                                if ($zoom == 0 && $mapwidthinmeters < ($maxdistance + 50)) $zoom = $zoomLevel;
                        }
                }
                $center = $totalLat / sizeof($mapPoints) . "," . $totalLon / sizeof($mapPoints);
        }
        $output = "";
        if ($collapsible) $output.= '<div data-role="collapsible" data-collapsed="true"><h3>Open Map...</h3>';
        $output.= '<center><img src="' . curPageURL() . 'staticmaplite/staticmap.php?center=' . $center . '&zoom=' . $zoom . '&size=' . $width . 'x' . $height . '&maptype=mapnik&markers=' . $markers . '" width=' . $width . ' height=' . $height . '></center>';
        if ($collapsible) $output.= '</div>';
        return $output;
}
function distance($lat1, $lng1, $lat2, $lng2)
{
        $pi80 = M_PI / 180;
        $lat1*= $pi80;
        $lng1*= $pi80;
        $lat2*= $pi80;
        $lng2*= $pi80;
        $r = 6372.797; // mean radius of Earth in km
        $dlat = $lat2 - $lat1;
        $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;
        return $km * 1000;
}
function decodePolylineToArray($encoded)
{
        // source: http://latlongeeks.com/forum/viewtopic.php?f=4&t=5
        $length = strlen($encoded);
        $index = 0;
        $points = array();
        $lat = 0;
        $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 geocode($query, $giveOptions)
{
        global $cloudmadeAPIkey;
        $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));
        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;
}
?>