2 /***************************************************************************
3 * Copyright (C) 2003-2009 Polytechnique.org *
4 * http://opensource.polytechnique.org/ *
6 * This program is free software; you can redistribute it and/or modify *
7 * it under the terms of the GNU General Public License as published by *
8 * the Free Software Foundation; either version 2 of the License, or *
9 * (at your option) any later version. *
11 * This program is distributed in the hope that it will be useful, *
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
14 * GNU General Public License for more details. *
16 * You should have received a copy of the GNU General Public License *
17 * along with this program; if not, write to the Free Software *
19 * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA *
20 ***************************************************************************/
22 // {{{ liste les pays ou les régions d'un pays
23 /** donne la liste déroulante des pays
24 * @param $current pays actuellement selectionné
26 function geoloc_country($current, $avail_only = false
)
29 $res = XDB
::iterRow('SELECT g.a2, g.pays
31 INNER JOIN adresses AS a ON(a.country = g.a2)
35 $res = XDB
::iterRow('SELECT a2,pays FROM geoloc_pays ORDER BY pays');
38 while (list($my_id, $my_pays) = $res->next()) {
39 $html .= sprintf("<option value=\"%s\" %s>%s</option>\n",
40 $my_id, ($current==$my_id?
"selected='selected'":""), $my_pays);
45 /** donne la liste deroulante des regions pour un pays
46 * @param $pays le pays dont on veut afficher les regions
47 * @param $current la region actuellement selectionnee
49 function geoloc_region($country, $current, $avail_only = false
)
52 $res = XDB
::iterRow('SELECT r.region, r.name
53 FROM geoloc_region AS r
54 INNER JOIN adresses AS a ON (a.country = r.a2 AND a.region = r.region)
57 ORDER BY r.name', $country);
59 $res = XDB
::iterRow('SELECT region,name
62 ORDER BY name', $country);
64 $html = "<option value=\"\"></option>";
65 while (list($regid, $regname) = $res->next()) {
66 $html .= sprintf("<option value=\"%s\" %s>%s</option>\n",
67 $regid, ($current==$regid?
"selected='selected'":""), $regname);
73 // {{{ get_address_infos($txt)
74 /** retrieve the infos on a text address
75 * store on the fly the info of the city concerned
76 * @param $txt the raw text of an address
78 function get_address_infos($txt)
82 $url = $globals->geoloc
->webservice_url
."address.php?precise=1&txt=" . urlencode($txt);
83 if ($globals->debug
& DEBUG_BT
) {
84 if (!isset(PlBacktrace
::$bt['Geoloc'])) {
85 new PlBacktrace('Geoloc');
87 PlBacktrace
::$bt['Geoloc']->start($url);
89 $f = @fopen
($url, 'r');
91 if ($globals->debug
& DEBUG_BT
) {
92 PlBacktrace
::$bt['Geoloc']->stop(0, 'Can\'t fetch result');
96 $keys = explode('|',fgets($f));
97 $vals = explode('|',fgets($f));
98 if ($globals->debug
& DEBUG_BT
) {
100 for ($i = 0 ; $i < count($keys) ; ++
$i) {
101 $data[] = array($keys[$i], $vals[$i]);
103 PlBacktrace
::$bt['Geoloc']->stop(count($keys), null
, $data);
105 $infos = empty_address();
106 foreach ($keys as $i=>$key) {
109 $infos[$key] = $vals[$i];
111 $val = strtr($vals[$i], array(chr(197).chr(147) => "œ"));
116 if (empty($infos['country'])) {
117 $infos['country'] = '00';
119 if (isset($infos['sql']) && $infos['sql']) {
120 $sql = explode(', ', trim($infos['sql'], '()'));
121 if (count($sql) == 16) {
122 for ($i = 0 ; $i < 16 ; ++
$i) {
123 $sql[$i] = stripslashes(trim($sql[$i], ' \''));
125 XDB
::execute("REPLACE INTO geoloc_city
126 VALUES ({?}, {?}, {?}, {?}, {?}, {?},
127 {?}, {?}, {?}, {?}, {?}, {?},
128 {?}, {?}, {?}, {?})",
129 $sql[0], $sql[1], $sql[2], $sql[3], $sql[4], $sql[5],
130 $sql[6], $sql[7], $sql[8], $sql[9], $sql[10], $sql[11],
131 $sql[12], $sql[13], $sql[14], $sql[15]);
134 if (isset($infos['display']) && $infos['display'])
135 XDB
::execute("UPDATE geoloc_pays
137 WHERE a2 = {?}", $infos['display'], $infos['country']);
138 if (isset($infos['cityid'])) {
139 fix_cities_not_on_map(1, $infos['cityid']);
140 if (floatval($infos['precise_lat']) && floatval($infos['precise_lon'])) {
141 $res = XDB
::query("SELECT c.lat / 100000, c.lon / 100000
142 FROM geoloc_city AS c
143 WHERE c.id = {?}", $infos['cityid']);
144 if ($res->numRows()) {
145 list($glat, $glng) = $res->fetchOneRow();
146 $infos['precise_lat'] = $glat;
147 $infos['precise_lon'] = $glng;
155 // {{{ get_cities_maps($array)
156 /* get all the maps id of the cities contained in an array */
157 function get_cities_maps($array)
160 implode("\n",$array);
161 $url = $globals->geoloc
->webservice_url
."findMaps.php?datatext=".urlencode(implode("\n", $array));
162 if (!($f = @fopen
($url, 'r'))) return false
;
166 $l = trim(fgets($f));
167 $tab = explode(';', $l);
176 // {{{ get_new_maps($url)
177 /** set new maps from url **/
178 function get_new_maps($url)
180 if (!($f = @fopen
($url, 'r'))) {
183 XDB
::query('TRUNCATE TABLE geoloc_maps');
186 $l = fgetcsv($f, 1024, ';', '"');
187 foreach ($l as $i => $val) {
188 if ($val != 'NULL') {
189 $l[$i] = '\''.addslashes($val).'\'';
192 $s .= ',('.implode(',',$l).')';
194 XDB
::execute('INSERT INTO geoloc_maps VALUES '.substr($s, 1));
199 // {{{ get_address_text($adr)
200 /** make the text of an address that can be read by a mailman
201 * @param $adr an array with all the usual fields
203 function get_address_text($adr)
206 if (isset($adr['adr1']) && $adr['adr1']) $t.= $adr['adr1'];
207 if (isset($adr['adr2']) && $adr['adr2']) $t.= "\n".$adr['adr2'];
208 if (isset($adr['adr3']) && $adr['adr3']) $t.= "\n".$adr['adr3'];
210 if (isset($adr['display']) && $adr['display']) {
211 $keys = explode(' ', $adr['display']);
212 foreach ($keys as $key) {
213 if (isset($adr[$key])) {
214 $l .= " ".$adr[$key];
219 if ($l) substr($l, 1);
220 } elseif ($adr['country'] == 'US' ||
$adr['country'] == 'CA' ||
$adr['country'] == 'GB') {
221 if ($adr['city']) $l .= $adr['city'].",\n";
222 if ($adr['region']) $l .= $adr['region']." ";
223 if ($adr['postcode']) $l .= $adr['postcode'];
225 if (isset($adr['postcode']) && $adr['postcode']) $l .= $adr['postcode']." ";
226 if (isset($adr['city']) && $adr['city']) $l .= $adr['city'];
228 if ($l) $t .= "\n".trim($l);
229 if ($adr['country'] != '00' && (!$adr['countrytxt'] ||
$adr['countrytxt'] == mb_strtoupper($adr['countrytxt']))) {
230 $res = XDB
::query("SELECT pays FROM geoloc_pays WHERE a2 = {?}", $adr['country']);
231 $adr['countrytxt'] = $res->fetchOneCell();
233 if (isset($adr['countrytxt']) && $adr['countrytxt']) {
234 $t .= "\n".$adr['countrytxt'];
240 // {{{ compare_addresses_text($a, $b)
241 /** compares if two address matches
242 * @param $a the raw text of an address
243 * @param $b the raw text of a complete valid address
245 function compare_addresses_text($a, $b)
247 $ta = mb_strtoupper(preg_replace(array("/[0-9,\"'#~:;_\- ]/", "/\r\n/"), array("", "\n"), $a));
248 $tb = mb_strtoupper(preg_replace(array("/[0-9,\"'#~:;_\- ]/", "/\r\n/"), array("", "\n"), $b));
250 $la = explode("\n", $ta);
251 $lb = explode("\n", $tb);
253 if (count($lb) > count($la) +
1) {
256 foreach ($la as $i=>$l) {
257 if (levenshtein(trim($l), trim($lb[$i])) > 3) {
266 function empty_address() {
279 "precise_lon" => "");
282 // create a simple address from a text without geoloc
283 function cut_address($txt)
285 $txt = str_replace("\r\n", "\n", $txt);
286 ereg("^([^\n]*)(\n([^\n]*)(\n(.*))?)?$", trim($txt), $a);
287 return array("adr1" => trim($a[1]), "adr2" => trim($a[3]), "adr3" => trim(str_replace("\n", " ", $a[5])));
290 // {{{ localize_addresses($uid)
291 /* localize all the address of a user and modify the database
292 * if the new address match with the old one
293 * @param $uid the id of the user
295 function localize_addresses($uid)
297 $res = XDB
::iterator("SELECT *
299 WHERE uid = {?} and (cityid IS NULL OR cityid = 0)", $uid);
302 while ($a = $res->next()) {
303 $new = get_address_infos($ta = get_address_text($a));
304 if (compare_addresses_text($ta, get_address_text($new))) {
305 XDB
::execute("UPDATE adresses
306 SET adr1 = {?}, adr2 = {?}, adr3 = {?},
307 cityid = {?}, city = {?}, postcode = {?},
308 region = {?}, regiontxt = {?}, country = {?},
309 glat = {?}, glng = {?}
310 WHERE uid = {?} AND adrid = {?}",
311 $new['adr1'], $new['adr2'], $new['adr3'],
312 $new['cityid'], $new['city'], $new['postcode'],
313 $new['region'], $new['regiontxt'], $new['country'],
314 $new['precise_lat'], $new['precise_lon'],
316 $new['store'] = true
;
317 if (!$new['cityid']) {
318 $erreur[$a['adrid']] = $new;
321 $new['store'] = false
;
322 $erreur[$a['adrid']] = $new;
329 // {{{ synchro_city($id)
330 /** synchronise the local geoloc_city base to geoloc.org
331 * @param $id the id of the city to synchronize
333 function synchro_city($id)
336 $url = $globals->geoloc
->webservice_url
."cityFinder.php?method=id&id=".$id."&out=sql";
337 if (!($f = @fopen
($url, 'r'))) {
342 return XDB
::execute("REPLACE INTO geoloc_city VALUES ".$s) > 0;
347 // {{{ function fix_cities_not_on_map($limit)
348 function fix_cities_not_on_map($limit=false
, $cityid=false
)
350 $missing = XDB
::query("SELECT c.id
351 FROM geoloc_city AS c
352 LEFT JOIN geoloc_city_in_maps AS m ON(c.id = m.city_id)
353 WHERE m.city_id IS NULL"
354 . ($cityid ?
" AND c.id = '" . $cityid . "'" : "" )
355 . ($limit ?
" LIMIT $limit" : "" ));
356 $maps = get_cities_maps($missing->fetchColumn());
359 foreach ($maps as $cityid => $maps_c) {
360 foreach ($maps_c as $map_id) {
361 $values .= ",($cityid, $map_id, '')";
364 if (strlen($values) > 1) {
365 XDB
::execute("REPLACE INTO geoloc_city_in_maps
366 VALUES ".substr($values, 1));
374 function set_smallest_levels()
376 $maxlengths = XDB
::iterRow("SELECT MAX(LENGTH(gm.path)), gcim.city_id
377 FROM geoloc_city_in_maps AS gcim
378 INNER JOIN geoloc_maps AS gm USING ( map_id )
379 GROUP BY gcim.city_id");
380 while (list($length, $id) = $maxlengths->next()) {
381 XDB
::execute("UPDATE geoloc_city_in_maps AS gcim
382 INNER JOIN geoloc_maps AS gm USING(map_id)
383 SET gcim.infos = IF(LENGTH(gm.path) = {?}, 'smallest', '')
384 WHERE gcim.city_id = {?}", $length, $id);
391 function geoloc_to_x($lon, $lat)
393 return deg2rad(1) * $lon *100;
396 function geoloc_to_y($lon, $lat)
404 return -100 * log(tan(pi()/4 +
deg2rad(1)/2*$lat));
407 function size_of_city($nb)
409 $s = round(log($nb +
1)*2,2);
416 function size_of_territory($nb)
418 return size_of_city($nb);
421 function geoloc_getData_subcities($mapid, $SFields, &$cities, $direct=true
)
423 if ($SFields instanceof UserSet
) {
427 $set = new UserSet();
429 for ($i_mapfield=0; $i_mapfield < count($SFields) ; $i_mapfield++
) {
430 if ($SFields[$i_mapfield]->fieldFormName
== 'mapid') {
434 $SFields[$i_mapfield] = new MapSField('mapid',
435 array('gcim.map_id'),
436 array('adresses','geoloc_city_in_maps'),
438 array(getadr_join('am'), 'am.cityid = gcim.city_id'),
440 $fields = new SFieldGroup(true
, $SFields);
441 $where = $fields->get_where_statement();
442 $joins = $fields->get_select_statement();
446 $cityres = $set->get('gc.id,
447 gc.lon / 100000 AS x, gc.lat/100000 AS y,
449 COUNT(u.user_id) AS pop,
450 SUM(u.promo % 2) AS yellow',
452 LEFT JOIN geoloc_city AS gc ON(gcim.city_id = gc.id)",
453 $where . ($direct ?
"gcim.infos = 'smallest'" : '1'),
456 foreach($cityres as $c) {
459 $city['x'] = geoloc_to_x($c['x'], $c['y']);
460 $city['y'] = geoloc_to_y($c['x'], $c['y']);
461 $city['size'] = size_of_city($c['pop']);
462 $cities[$c['id']] = $city;
467 function geoloc_getData_subcountries($mapid, $sin, $minentities)
469 $countries = array();
472 if ($mapid === false
) {
473 $wheremapid = "WHERE gm.parent IS NULL";
475 $wheremapid = "WHERE gm.parent = {?}";
477 $submapres = XDB
::iterator(
478 "SELECT gm.map_id AS id, gm.name, gm.x, gm.y, gm.xclip, gm.yclip,
479 gm.width, gm.height, gm.scale, 1 AS rat
480 FROM geoloc_maps AS gm
481 ". $wheremapid, Env
::v('mapid',''));
485 while ($c = $submapres->next()) {
487 $country['color'] = 0xFFFFFF;
488 $country['swf'] = $globals->geoloc
->webservice_url
."maps/mercator/map_".$c['id'].".swf";
489 $countries[$c['id']] = $country;
492 if ($mapid === false
) {
493 return array($countries, $cities);
496 geoloc_getData_subcities(Env
::i('mapid'), $sin, $cities);
497 $nbcities = count($cities);
498 $nocity = $nbcities == 0;
499 if ($sin instanceof UserSet
) {
503 $set = new UserSet();
507 for ($i_mapfield=0; $i_mapfield < count($SFields) ; $i_mapfield++
) {
508 if ($SFields[$i_mapfield]->fieldFormName
== 'mapid') {
512 $SFields[$i_mapfield] = new MapSField('mapid',
514 array('adresses','geoloc_city_in_maps','geoloc_maps'),
515 array('am','gcim','map'),
516 array(getadr_join('am'),
517 'am.cityid = gcim.city_id',
518 'map.map_id = gcim.map_id'));
519 $fields = new SFieldGroup(true
, $SFields);
520 $where = $fields->get_where_statement();
521 $joins = $fields->get_select_statement();
522 $countryres = $set->get('map.map_id AS id,
523 COUNT(u.user_id) AS nbPop,
524 SUM(u.promo % 2) AS yellow,
525 COUNT(DISTINCT gcim.city_id) AS nbCities,
526 SUM(IF(u.user_id IS NULL,0,am.glng)) AS lonPop,
527 SUM(IF(u.user_id IS NULL, 0,am.glat)) AS latPop',
534 $nbentities = $nbcities +
count($countryres);
535 foreach ($countryres as $c) {
536 $c['latPop'] /= $c['nbPop'];
537 $c['lonPop'] /= $c['nbPop'];
538 $c['rad'] = size_of_territory($c['nbPop']);
539 if ($maxpop < $c['nbPop']) $maxpop = $c['nbPop'];
540 $c['xPop'] = geoloc_to_x($c['lonPop'], $c['latPop']);
541 $c['yPop'] = geoloc_to_y($c['lonPop'], $c['latPop']);
542 @$countries[$c['id']] = array_merge($countries[$c['id']], $c);
544 $nbcities +
= $c['nbCities'];
547 if ($nocity && $nbcities < $minentities){
548 foreach($countries as $i => $c) {
549 $countries[$i]['nbPop'] = 0;
550 if (@$c['nbCities'] > 0) {
551 geoloc_getData_subcities($c['id'], $sin, $cities, false
);
556 foreach ($countries as $i => $c) {
557 if (@$c['nbPop'] > 0) {
558 $lambda = pow($c['nbPop'] / $maxpop,0.3);
559 $countries[$i]['color'] = 0x0000FF +
round((1-$lambda) * 0xFF)*0x010100;
563 return array($countries, $cities);
567 // vim:set et sw=4 sts=4 sws=4 foldmethod=marker enc=utf-8: