[platal.git] / include / geoloc.inc.php

<?php
/***************************************************************************
 *  Copyright (C) 2003-2008 Polytechnique.org                              *
 *  http://opensource.polytechnique.org/                                   *
 *                                                                         *
 *  This program is free software; you can redistribute it and/or modify   *
 *  it under the terms of the GNU General Public License as published by   *
 *  the Free Software Foundation; either version 2 of the License, or      *
 *  (at your option) any later version.                                    *
 *                                                                         *
 *  This program is distributed in the hope that it will be useful,        *
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of         *
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the          *
 *  GNU General Public License for more details.                           *
 *                                                                         *
 *  You should have received a copy of the GNU General Public License      *
 *  along with this program; if not, write to the Free Software            *
 *  Foundation, Inc.,                                                      *
 *  59 Temple Place, Suite 330, Boston, MA  02111-1307  USA                *
 ***************************************************************************/

// {{{ liste les pays ou les régions d'un pays
/** donne la liste déroulante des pays
 * @param $current pays actuellement selectionné
 */
function geoloc_country($current, $avail_only = false)
{
    if ($avail_only) {
        $res = XDB::iterRow('SELECT  g.a2, g.pays
                               FROM  geoloc_pays AS g
                         INNER JOIN  adresses    AS a ON(a.country = g.a2)
                           GROUP BY  g.a2
                           ORDER BY  g.pays');
    } else {
        $res = XDB::iterRow('SELECT a2,pays FROM geoloc_pays ORDER BY pays');
    }
    $html = "";
    while (list($my_id, $my_pays) = $res->next()) {
        $html .= sprintf("<option value=\"%s\" %s>%s</option>\n",
                         $my_id, ($current==$my_id?"selected='selected'":""), $my_pays);
    }
    return $html;
}

/** donne la liste deroulante des regions pour un pays
 * @param $pays le pays dont on veut afficher les regions
 * @param $current la region actuellement selectionnee
 */
function geoloc_region($country, $current, $avail_only = false)
{
    if ($avail_only) {
        $res = XDB::iterRow('SELECT  r.region, r.name
                               FROM  geoloc_region AS r
                         INNER JOIN  adresses      AS a ON (a.country = r.a2 AND a.region = r.region)
                              WHERE  r.a2 = {?}
                           GROUP BY  r.region
                           ORDER BY  r.name', $country);
    } else {
        $res = XDB::iterRow('SELECT  region,name
                               FROM  geoloc_region
                              WHERE  a2 = {?}
                           ORDER BY  name', $country);
    }
    $html = "<option value=\"\"></option>";
    while (list($regid, $regname) = $res->next()) {
        $html .= sprintf("<option value=\"%s\" %s>%s</option>\n",
                 $regid, ($current==$regid?"selected='selected'":""), $regname);
    }
    return $html;
}
// }}}

// {{{ get_address_infos($txt)
/** retrieve the infos on a text address
 * store on the fly the info of the city concerned
 * @param $txt the raw text of an address
 */
function get_address_infos($txt)
{
    global $globals;
    $url = $globals->geoloc->webservice_url."address.php?precise=1&txt=" . urlencode($txt);
    if (!($f = @fopen($url, 'r'))) return false;
    $keys = explode('|',fgets($f));
    $vals = explode('|',fgets($f));
    $infos = array();
    foreach ($keys as $i=>$key) {
        if($vals[$i]) {
            if ($key == 'sql') {
                $infos[$key] = $vals[$i];
            } else {
                $val = strtr($vals[$i], array(chr(197).chr(147) => "&oelig;"));
                $infos[$key] = $val;
            }
        }
    }
    if (isset($infos['sql']) && $infos['sql'])
       XDB::execute("REPLACE INTO  geoloc_city
                           VALUES  ".$infos['sql']);
    if (isset($infos['display']) && $infos['display'])
        XDB::execute("UPDATE  geoloc_pays
                         SET  display = {?}
                       WHERE  a2 = {?}", $infos['display'], $infos['country']);
    if (isset($infos['cityid'])) {
        fix_cities_not_on_map(1, $infos['cityid']);
        if (floatval($infos['precise_lat']) && floatval($infos['precise_lon'])) {
            $res = XDB::query("SELECT  c.lat / 100000, c.lon / 100000
                                 FROM  geoloc_city AS c
                                WHERE  c.id = {?}", $infos['cityid']);
            if ($res->numRows()) {
                list($glat, $glng) = $res->fetchOneRow();
                $infos['precise_lat'] = $glat;
                $infos['precise_lon'] = $glng;
            }
        }
    }
    return $infos;
}
// }}}

// {{{ get_cities_maps($array)
/* get all the maps id of the cities contained in an array */
function get_cities_maps($array)
{
    global $globals;
    implode("\n",$array);
    $url = $globals->geoloc->webservice_url."findMaps.php?datatext=".urlencode(implode("\n", $array));
    if (!($f = @fopen($url, 'r'))) return false;
    $maps = array();
    while (!feof($f))
    {
        $l = trim(fgets($f));
        $tab = explode(';', $l);
        $i = $tab[0];
        unset($tab[0]);
        $maps[$i] = $tab;
    }
    return $maps;
}
// }}}

// {{{ get_new_maps($url)
/** set new maps from url **/
function get_new_maps($url)
{
    if (!($f = @fopen($url, 'r'))) {
        return false;
    }
    XDB::query('TRUNCATE TABLE geoloc_maps');
    $s = '';
    while (!feof($f)) {
        $l = fgetcsv($f, 1024, ';', '"');
        foreach ($l as $i => $val) {
            if ($val != 'NULL') {
                $l[$i] = '\''.addslashes($val).'\'';
            }
        }
        $s .= ',('.implode(',',$l).')';
    }
    XDB::execute('INSERT INTO geoloc_maps VALUES '.substr($s, 1));
    return true;
}
// }}}

// {{{ get_address_text($adr)
/** make the text of an address that can be read by a mailman
 * @param $adr an array with all the usual fields
 */
function get_address_text($adr)
{
    $t = "";
    if (isset($adr['adr1']) && $adr['adr1']) $t.= $adr['adr1'];
    if (isset($adr['adr2']) && $adr['adr2']) $t.= "\n".$adr['adr2'];
    if (isset($adr['adr3']) && $adr['adr3']) $t.= "\n".$adr['adr3'];
    $l = "";
    if (isset($adr['display']) && $adr['display']) {
        $keys = explode(' ', $adr['display']);
        foreach ($keys as $key) {
            if (isset($adr[$key])) {
                $l .= " ".$adr[$key];
            } else {
                $l .= " ".$key;
            }
        }
        if ($l) substr($l, 1);
    } elseif ($adr['country'] == 'US' || $adr['country'] == 'CA' || $adr['country'] == 'GB') {
        if ($adr['city']) $l .= $adr['city'].",\n";
        if ($adr['region']) $l .= $adr['region']." ";
        if ($adr['postcode']) $l .= $adr['postcode'];
    } else {
        if (isset($adr['postcode']) && $adr['postcode']) $l .= $adr['postcode']." ";
        if (isset($adr['city']) && $adr['city']) $l .= $adr['city'];
    }
    if ($l) $t .= "\n".trim($l);
    if ($adr['country'] != '00' && (!$adr['countrytxt'] || $adr['countrytxt'] == strtoupper($adr['countrytxt']))) {
        $res = XDB::query("SELECT pays FROM geoloc_pays WHERE a2 = {?}", $adr['country']);
        $adr['countrytxt'] = $res->fetchOneCell();
    }
    if (isset($adr['countrytxt']) && $adr['countrytxt']) {
        $t .= "\n".$adr['countrytxt'];
    }
    return trim($t);
}
// }}}

// {{{ compare_addresses_text($a, $b)
/** compares if two address matches
 * @param $a the raw text of an address
 * @param $b the raw text of a complete valid address
 */
function compare_addresses_text($a, $b)
{
    $ta = strtoupper(preg_replace(array("/[0-9,\"'#~:;_\- ]/", "/\r\n/"), array("", "\n"), $a));
    $tb = strtoupper(preg_replace(array("/[0-9,\"'#~:;_\- ]/", "/\r\n/"), array("", "\n"), $b));

    $la = explode("\n", $ta);
    $lb = explode("\n", $tb);

    if (count($lb) > count($la) + 1) {
        return false;
    }
    foreach ($la as $i=>$l) {
        if (levenshtein(trim($l), trim($lb[$i])) > 3) {
            return false;
        }
    }
    return true;
}

// }}}

function empty_address() {
    return Array(
        "adr1" => "",
        "adr2" => "",
        "adr3" => "",
        "cityid" => NULL,
        "city" => "",
        "postcode" => "",
        "region" => "",
        "regiontxt" => "",
        "country" => "00",
        "countrytxt" => "",
        "precise_lat" => "",
        "precise_lon" => "");
}

// create a simple address from a text without geoloc
function cut_address($txt)
{
    $txt = str_replace("\r\n", "\n", $txt);
    ereg("^([^\n]*)(\n([^\n]*)(\n(.*))?)?$", trim($txt), $a);
    return array("adr1" => trim($a[1]), "adr2" => trim($a[3]), "adr3" => trim(str_replace("\n", " ", $a[5])));
}

// {{{ localize_addresses($uid)
/* localize all the address of a user and modify the database
 * if the new address match with the old one
 * @param $uid the id of the user
 */
function localize_addresses($uid)
{
    $res = XDB::iterator("SELECT  *
                            FROM  adresses
                           WHERE  uid = {?} and (cityid IS NULL OR cityid = 0)", $uid);
    $erreur = Array();

    while ($a = $res->next()) {
        $new = get_address_infos($ta = get_address_text($a));
        if (compare_addresses_text($ta, get_address_text($new))) {
            XDB::execute("UPDATE  adresses
                             SET  adr1 = {?}, adr2 = {?}, adr3 = {?},
                                  cityid = {?}, city = {?}, postcode = {?},
                                  region = {?}, regiontxt = {?}, country = {?},
                                  glat = {?}, glng = {?}
                           WHERE  uid = {?} AND adrid = {?}",
                          $new['adr1'], $new['adr2'], $new['adr3'],
                          $new['cityid'], $new['city'], $new['postcode'],
                          $new['region'], $new['regiontxt'], $new['country'],
                          $new['precise_lat'], $new['precise_lon'],
                          $uid, $a['adrid']);
            $new['store'] = true;
            if (!$new['cityid']) {
                $erreur[$a['adrid']] = $new;
            }
        } else {
            $new['store'] = false;
            $erreur[$a['adrid']] = $new;
        }
    }
    return $erreur;
}
// }}}

// {{{ synchro_city($id)
/** synchronise the local geoloc_city base to geoloc.org
 * @param $id the id of the city to synchronize
 */
function synchro_city($id)
{
    global $globals;
    $url = $globals->geoloc->webservice_url."cityFinder.php?method=id&id=".$id."&out=sql";
    if (!($f = @fopen($url, 'r'))) {
        return false;
    }
    $s = fgets($f);
    if ($s) {
        return XDB::execute("REPLACE INTO geoloc_city VALUES ".$s) > 0;
    }
}
 // }}}

// {{{ function fix_cities_not_on_map($limit)
function fix_cities_not_on_map($limit=false, $cityid=false)
{
    $missing = XDB::query("SELECT  c.id
                             FROM  geoloc_city AS c
                        LEFT JOIN  geoloc_city_in_maps AS m ON(c.id = m.city_id)
                            WHERE  m.city_id IS NULL"
                            . ($cityid ? " AND c.id = '" . $cityid . "'" : "" )
                            . ($limit ?  " LIMIT $limit" : "" ));
    $maps = get_cities_maps($missing->fetchColumn());
    if ($maps) {
        $values = "";
        foreach ($maps as $cityid => $maps_c) {
            foreach ($maps_c as $map_id) {
                $values .= ",($cityid, $map_id, '')";
            }
        }
        if (strlen($values) > 1) {
            XDB::execute("REPLACE INTO  geoloc_city_in_maps
                                VALUES  ".substr($values, 1));
        }
    } else {
        return false;
    }
    return true;
}

function set_smallest_levels()
{
    $maxlengths = XDB::iterRow("SELECT  MAX(LENGTH(gm.path)), gcim.city_id
                                  FROM  geoloc_city_in_maps AS gcim
                            INNER JOIN  geoloc_maps AS gm USING ( map_id )
                              GROUP BY  gcim.city_id");
    while (list($length, $id) = $maxlengths->next()) {
        XDB::execute("UPDATE  geoloc_city_in_maps AS gcim
                  INNER JOIN  geoloc_maps AS gm USING(map_id)
                         SET  gcim.infos = IF(LENGTH(gm.path) = {?}, 'smallest', '')
                       WHERE  gcim.city_id = {?}", $length, $id);
    }
    return true;
}
// }}}


function geoloc_to_x($lon, $lat)
{
    return deg2rad(1) * $lon *100;
}

function geoloc_to_y($lon, $lat)
{
    if ($lat < -75) {
        return latToY(-75);
    }
    if ($lat > 75) {
        return latToY(75);
    }
    return -100 * log(tan(pi()/4 + deg2rad(1)/2*$lat));
}

function size_of_city($nb)
{
    $s = round(log($nb + 1)*2,2);
    if ($s < 1) {
        return 1;
    }
    return $s;
}

function size_of_territory($nb)
{
    return size_of_city($nb);
}

function geoloc_getData_subcities($mapid, $SFields, &$cities, $direct=true)
{
    if ($SFields instanceof UserSet) {
        $set = $SFields;
        $SFields = array();
    } else {
        $set = new UserSet();
    }
    for ($i_mapfield=0; $i_mapfield < count($SFields) ; $i_mapfield++) {
        if ($SFields[$i_mapfield]->fieldFormName == 'mapid') {
            break;
        }
    }
    $SFields[$i_mapfield] = new MapSField('mapid',
                                          array('gcim.map_id'),
                                          array('adresses','geoloc_city_in_maps'),
                                          array('am','gcim'),
                                          array(getadr_join('am'), 'am.cityid = gcim.city_id'),
                                          $mapid);
    $fields = new SFieldGroup(true, $SFields);
    $where = $fields->get_where_statement();
    $joins = $fields->get_select_statement();
    if ($where) {
        $where .= ' AND ';
    }
    $cityres = $set->get('gc.id,
                          gc.lon / 100000 AS x, gc.lat/100000 AS y,
                          gc.name,
                          COUNT(u.user_id) AS pop,
                          SUM(u.promo % 2) AS yellow',
                          "$joins
                          LEFT JOIN  geoloc_city  AS gc ON(gcim.city_id = gc.id)",
                          $where . ($direct ? "gcim.infos = 'smallest'" : '1'),
                          'gc.id, gc.alias',
                          'pop DESC');
    foreach($cityres as $c) {
        if ($c['pop'] > 0) {
            $city = $c;
            $city['x'] = geoloc_to_x($c['x'], $c['y']);
            $city['y'] = geoloc_to_y($c['x'], $c['y']);
            $city['size'] = size_of_city($c['pop']);
            $cities[$c['id']] = $city;
        }
    }
}

function geoloc_getData_subcountries($mapid, $sin, $minentities)
{
    $countries = array();
    $cities = array();

    if ($mapid === false) {
        $wheremapid = "WHERE gm.parent IS NULL";
    } else {
        $wheremapid = "WHERE   gm.parent = {?}";
    }
    $submapres = XDB::iterator(
        "SELECT  gm.map_id AS id, gm.name, gm.x, gm.y, gm.xclip, gm.yclip,
                 gm.width, gm.height, gm.scale, 1 AS rat
           FROM  geoloc_maps AS gm
        ". $wheremapid, Env::v('mapid',''));

    global $globals;

    while ($c = $submapres->next()) {
        $country = $c;
        $country['color'] = 0xFFFFFF;
        $country['swf'] = $globals->geoloc->webservice_url."maps/mercator/map_".$c['id'].".swf";
        $countries[$c['id']] = $country;
    }

    if ($mapid === false) {
        return array($countries, $cities);
    }

    geoloc_getData_subcities(Env::i('mapid'), $sin, $cities);
    $nbcities = count($cities);
    $nocity = $nbcities == 0;
    if ($sin instanceof UserSet) {
        $set = $sin;
        $SFields = array();
    } else {
        $set = new UserSet();
        $SFields = $sin;
    }

    for ($i_mapfield=0; $i_mapfield < count($SFields) ; $i_mapfield++) {
        if ($SFields[$i_mapfield]->fieldFormName == 'mapid') {
            break;
        }
    }
    $SFields[$i_mapfield] = new MapSField('mapid',
                                          array('map.parent'),
                                          array('adresses','geoloc_city_in_maps','geoloc_maps'),
                                          array('am','gcim','map'),
                                          array(getadr_join('am'),
                                          'am.cityid = gcim.city_id',
                                          'map.map_id = gcim.map_id'));
    $fields = new SFieldGroup(true, $SFields);
    $where  = $fields->get_where_statement();
    $joins  = $fields->get_select_statement();
    $countryres = $set->get('map.map_id AS id,
                             COUNT(u.user_id) AS nbPop,
                             SUM(u.promo % 2) AS yellow,
                             COUNT(DISTINCT gcim.city_id) AS nbCities,
                             SUM(IF(u.user_id IS NULL,0,am.glng)) AS lonPop,
                             SUM(IF(u.user_id IS NULL, 0,am.glat)) AS latPop',
                            $joins,
                            $where,
                            'map.map_id',
                            'NULL');

    $maxpop = 0;
    $nbentities = $nbcities + count($countryres);
    foreach ($countryres as $c) {
        $c['latPop'] /= $c['nbPop'];
        $c['lonPop'] /= $c['nbPop'];
        $c['rad'] = size_of_territory($c['nbPop']);
        if ($maxpop < $c['nbPop']) $maxpop = $c['nbPop'];
        $c['xPop'] = geoloc_to_x($c['lonPop'], $c['latPop']);
        $c['yPop'] = geoloc_to_y($c['lonPop'], $c['latPop']);
        @$countries[$c['id']] = array_merge($countries[$c['id']], $c);

        $nbcities += $c['nbCities'];
    }

    if ($nocity && $nbcities < $minentities){
        foreach($countries as $i => $c) {
            $countries[$i]['nbPop'] = 0;
            if (@$c['nbCities'] > 0) {
                geoloc_getData_subcities($c['id'], $sin, $cities, false);
            }
        }
    }

    foreach ($countries as $i => $c) {
        if (@$c['nbPop'] > 0) {
            $lambda = pow($c['nbPop'] / $maxpop,0.3);
            $countries[$i]['color'] = 0x0000FF + round((1-$lambda) * 0xFF)*0x010100;
        }
    }

    return array($countries, $cities);
}
// }}}

// vim:set et sw=4 sts=4 sws=4 foldmethod=marker enc=utf-8:
?>
Commit	Line	Data
0337d704	1	<?php
0337d704	2	/***************************************************************************
179afa7f	3	* Copyright (C) 2003-2008 Polytechnique.org *
0337d704	4	* http://opensource.polytechnique.org/ *
	5	* *
	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. *
	10	* *
	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. *
	15	* *
	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 *
	18	* Foundation, Inc., *
	19	* 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA *
	20	***************************************************************************/
	21
a7de4ef7	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é
0337d704	25	*/
b536b5ae	26	function geoloc_country($current, $avail_only = false)
	27	{
	28	if ($avail_only) {
	29	$res = XDB::iterRow('SELECT g.a2, g.pays
	30	FROM geoloc_pays AS g
	31	INNER JOIN adresses AS a ON(a.country = g.a2)
	32	GROUP BY g.a2
	33	ORDER BY g.pays');
	34	} else {
	35	$res = XDB::iterRow('SELECT a2,pays FROM geoloc_pays ORDER BY pays');
	36	}
0337d704	37	$html = "";
0337d704	38	while (list($my_id, $my_pays) = $res->next()) {
b536b5ae	39	$html .= sprintf("<option value=\"%s\" %s>%s</option>\n",
b536b5ae	40	$my_id, ($current==$my_id?"selected='selected'":""), $my_pays);
0337d704	41	}
	42	return $html;
	43	}
	44
0337d704	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
	48	*/
b536b5ae	49	function geoloc_region($country, $current, $avail_only = false)
	50	{
	51	if ($avail_only) {
	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)
	55	WHERE r.a2 = {?}
	56	GROUP BY r.region
	57	ORDER BY r.name', $country);
	58	} else {
	59	$res = XDB::iterRow('SELECT region,name
	60	FROM geoloc_region
	61	WHERE a2 = {?}
	62	ORDER BY name', $country);
	63	}
0337d704	64	$html = "<option value=\"\"></option>";
0337d704	65	while (list($regid, $regname) = $res->next()) {
b536b5ae	66	$html .= sprintf("<option value=\"%s\" %s>%s</option>\n",
b536b5ae	67	$regid, ($current==$regid?"selected='selected'":""), $regname);
0337d704	68	}
0337d704	69	return $html;
0337d704	70	}
0337d704	71	// }}}
	72
	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
	77	*/
787bb3d7	78	function get_address_infos($txt)
8c4a0c30	79	{
56670b6a	80	global $globals;
8c4a0c30	81	$url = $globals->geoloc->webservice_url."address.php?precise=1&txt=" . urlencode($txt);
0337d704	82	if (!($f = @fopen($url, 'r'))) return false;
	83	$keys = explode('\|',fgets($f));
	84	$vals = explode('\|',fgets($f));
	85	$infos = array();
2b98ac11	86	foreach ($keys as $i=>$key) {
	87	if($vals[$i]) {
	88	if ($key == 'sql') {
	89	$infos[$key] = $vals[$i];
	90	} else {
	91	$val = strtr($vals[$i], array(chr(197).chr(147) => "&oelig;"));
8c4a0c30	92	$infos[$key] = $val;
2b98ac11	93	}
	94	}
	95	}
5ee95138	96	if (isset($infos['sql']) && $infos['sql'])
de0485fd FB	97	XDB::execute("REPLACE INTO geoloc_city
de0485fd FB	98	VALUES ".$infos['sql']);
5ee95138	99	if (isset($infos['display']) && $infos['display'])
de0485fd FB	100	XDB::execute("UPDATE geoloc_pays
	101	SET display = {?}
	102	WHERE a2 = {?}", $infos['display'], $infos['country']);
	103	if (isset($infos['cityid'])) {
b536b5ae	104	fix_cities_not_on_map(1, $infos['cityid']);
de0485fd FB	105	if (floatval($infos['precise_lat']) && floatval($infos['precise_lon'])) {
	106	$res = XDB::query("SELECT c.lat / 100000, c.lon / 100000
	107	FROM geoloc_city AS c
	108	WHERE c.id = {?}", $infos['cityid']);
	109	if ($res->numRows()) {
	110	list($glat, $glng) = $res->fetchOneRow();
	111	$infos['precise_lat'] = $glat;
	112	$infos['precise_lon'] = $glng;
	113	}
	114	}
	115	}
0337d704	116	return $infos;
	117	}
	118	// }}}
	119
56670b6a	120	// {{{ get_cities_maps($array)
	121	/* get all the maps id of the cities contained in an array */
	122	function get_cities_maps($array)
	123	{
	124	global $globals;
	125	implode("\n",$array);
8c4a0c30	126	$url = $globals->geoloc->webservice_url."findMaps.php?datatext=".urlencode(implode("\n", $array));
56670b6a	127	if (!($f = @fopen($url, 'r'))) return false;
	128	$maps = array();
	129	while (!feof($f))
	130	{
	131	$l = trim(fgets($f));
	132	$tab = explode(';', $l);
	133	$i = $tab[0];
	134	unset($tab[0]);
	135	$maps[$i] = $tab;
	136	}
	137	return $maps;
	138	}
	139	// }}}
	140
	141	// {{{ get_new_maps($url)
	142	/ set new maps from url /
	143	function get_new_maps($url)
	144	{
a3a049fc	145	if (!($f = @fopen($url, 'r'))) {
	146	return false;
	147	}
	148	XDB::query('TRUNCATE TABLE geoloc_maps');
	149	$s = '';
	150	while (!feof($f)) {
	151	$l = fgetcsv($f, 1024, ';', '"');
	152	foreach ($l as $i => $val) {
	153	if ($val != 'NULL') {
	154	$l[$i] = '\''.addslashes($val).'\'';
	155	}
	156	}
	157	$s .= ',('.implode(',',$l).')';
	158	}
	159	XDB::execute('INSERT INTO geoloc_maps VALUES '.substr($s, 1));
	160	return true;
56670b6a	161	}
80244bbe	162	// }}}
56670b6a	163
0337d704	164	// {{{ get_address_text($adr)
	165	/** make the text of an address that can be read by a mailman
	166	* @param $adr an array with all the usual fields
	167	*/
787bb3d7	168	function get_address_text($adr)
8c4a0c30	169	{
0337d704	170	$t = "";
5ee95138	171	if (isset($adr['adr1']) && $adr['adr1']) $t.= $adr['adr1'];
	172	if (isset($adr['adr2']) && $adr['adr2']) $t.= "\n".$adr['adr2'];
	173	if (isset($adr['adr3']) && $adr['adr3']) $t.= "\n".$adr['adr3'];
0337d704	174	$l = "";
5ee95138	175	if (isset($adr['display']) && $adr['display']) {
80ca1b4e	176	$keys = explode(' ', $adr['display']);
80ca1b4e	177	foreach ($keys as $key) {
8c4a0c30	178	if (isset($adr[$key])) {
80ca1b4e	179	$l .= " ".$adr[$key];
8c4a0c30	180	} else {
80ca1b4e	181	$l .= " ".$key;
8c4a0c30	182	}
80ca1b4e	183	}
8c4a0c30	184	if ($l) substr($l, 1);
	185	} elseif ($adr['country'] == 'US' \|\| $adr['country'] == 'CA' \|\| $adr['country'] == 'GB') {
	186	if ($adr['city']) $l .= $adr['city'].",\n";
	187	if ($adr['region']) $l .= $adr['region']." ";
	188	if ($adr['postcode']) $l .= $adr['postcode'];
	189	} else {
	190	if (isset($adr['postcode']) && $adr['postcode']) $l .= $adr['postcode']." ";
	191	if (isset($adr['city']) && $adr['city']) $l .= $adr['city'];
0337d704	192	}
	193	if ($l) $t .= "\n".trim($l);
	194	if ($adr['country'] != '00' && (!$adr['countrytxt'] \|\| $adr['countrytxt'] == strtoupper($adr['countrytxt']))) {
08cce2ff	195	$res = XDB::query("SELECT pays FROM geoloc_pays WHERE a2 = {?}", $adr['country']);
0337d704	196	$adr['countrytxt'] = $res->fetchOneCell();
0337d704	197	}
8c4a0c30	198	if (isset($adr['countrytxt']) && $adr['countrytxt']) {
	199	$t .= "\n".$adr['countrytxt'];
	200	}
0337d704	201	return trim($t);
	202	}
	203	// }}}
	204
	205	// {{{ compare_addresses_text($a, $b)
	206	/** compares if two address matches
	207	* @param $a the raw text of an address
	208	* @param $b the raw text of a complete valid address
	209	*/
8c4a0c30	210	function compare_addresses_text($a, $b)
8c4a0c30	211	{
0337d704	212	$ta = strtoupper(preg_replace(array("/[0-9,\"'#~:;_\- ]/", "/\r\n/"), array("", "\n"), $a));
0337d704	213	$tb = strtoupper(preg_replace(array("/[0-9,\"'#~:;_\- ]/", "/\r\n/"), array("", "\n"), $b));
787bb3d7	214
0337d704	215	$la = explode("\n", $ta);
	216	$lb = explode("\n", $tb);
	217
8c4a0c30	218	if (count($lb) > count($la) + 1) {
	219	return false;
	220	}
	221	foreach ($la as $i=>$l) {
37d44b3b	222	if (levenshtein(trim($l), trim($lb[$i])) > 3) {
8c4a0c30	223	return false;
	224	}
	225	}
0337d704	226	return true;
	227	}
	228
	229	// }}}
	230
	231	function empty_address() {
	232	return Array(
	233	"adr1" => "",
	234	"adr2" => "",
	235	"adr3" => "",
	236	"cityid" => NULL,
	237	"city" => "",
	238	"postcode" => "",
	239	"region" => "",
80ca1b4e	240	"regiontxt" => "",
0337d704	241	"country" => "00",
de0485fd FB	242	"countrytxt" => "",
	243	"precise_lat" => "",
	244	"precise_lon" => "");
0337d704	245	}
	246
	247	// create a simple address from a text without geoloc
8c4a0c30	248	function cut_address($txt)
8c4a0c30	249	{
0337d704	250	$txt = str_replace("\r\n", "\n", $txt);
	251	ereg("^([^\n])(\n([^\n])(\n(.*))?)?$", trim($txt), $a);
	252	return array("adr1" => trim($a[1]), "adr2" => trim($a[3]), "adr3" => trim(str_replace("\n", " ", $a[5])));
	253	}
	254
	255	// {{{ localize_addresses($uid)
	256	/* localize all the address of a user and modify the database
	257	* if the new address match with the old one
	258	* @param $uid the id of the user
	259	*/
8c4a0c30	260	function localize_addresses($uid)
	261	{
	262	$res = XDB::iterator("SELECT *
	263	FROM adresses
	264	WHERE uid = {?} and (cityid IS NULL OR cityid = 0)", $uid);
0337d704	265	$erreur = Array();
	266
	267	while ($a = $res->next()) {
	268	$new = get_address_infos($ta = get_address_text($a));
	269	if (compare_addresses_text($ta, get_address_text($new))) {
8c4a0c30	270	XDB::execute("UPDATE adresses
	271	SET adr1 = {?}, adr2 = {?}, adr3 = {?},
	272	cityid = {?}, city = {?}, postcode = {?},
	273	region = {?}, regiontxt = {?}, country = {?},
	274	glat = {?}, glng = {?}
	275	WHERE uid = {?} AND adrid = {?}",
	276	$new['adr1'], $new['adr2'], $new['adr3'],
	277	$new['cityid'], $new['city'], $new['postcode'],
	278	$new['region'], $new['regiontxt'], $new['country'],
	279	$new['precise_lat'], $new['precise_lon'],
	280	$uid, $a['adrid']);
	281	$new['store'] = true;
	282	if (!$new['cityid']) {
	283	$erreur[$a['adrid']] = $new;
	284	}
0337d704	285	} else {
	286	$new['store'] = false;
	287	$erreur[$a['adrid']] = $new;
	288	}
	289	}
	290	return $erreur;
	291	}
	292	// }}}
	293
	294	// {{{ synchro_city($id)
	295	/** synchronise the local geoloc_city base to geoloc.org
	296	* @param $id the id of the city to synchronize
	297	*/
787bb3d7	298	function synchro_city($id)
8c4a0c30	299	{
56670b6a	300	global $globals;
56670b6a	301	$url = $globals->geoloc->webservice_url."cityFinder.php?method=id&id=".$id."&out=sql";
8c4a0c30	302	if (!($f = @fopen($url, 'r'))) {
	303	return false;
	304	}
0337d704	305	$s = fgets($f);
8c4a0c30	306	if ($s) {
08cce2ff	307	return XDB::execute("REPLACE INTO geoloc_city VALUES ".$s) > 0;
8c4a0c30	308	}
8c4a0c30	309	}
0337d704	310	// }}}
0337d704	311
56670b6a	312	// {{{ function fix_cities_not_on_map($limit)
7f3ac007	313	function fix_cities_not_on_map($limit=false, $cityid=false)
56670b6a	314	{
8c4a0c30	315	$missing = XDB::query("SELECT c.id
8c4a0c30	316	FROM geoloc_city AS c
787bb3d7	317	LEFT JOIN geoloc_city_in_maps AS m ON(c.id = m.city_id)
8c4a0c30	318	WHERE m.city_id IS NULL"
	319	. ($cityid ? " AND c.id = '" . $cityid . "'" : "" )
	320	. ($limit ? " LIMIT $limit" : "" ));
56670b6a	321	$maps = get_cities_maps($missing->fetchColumn());
8c4a0c30	322	if ($maps) {
56670b6a	323	$values = "";
8c4a0c30	324	foreach ($maps as $cityid => $maps_c) {
8c4a0c30	325	foreach ($maps_c as $map_id) {
56670b6a	326	$values .= ",($cityid, $map_id, '')";
8c4a0c30	327	}
8c4a0c30	328	}
22cc3370 FB	329	if (strlen($values) > 1) {
	330	XDB::execute("REPLACE INTO geoloc_city_in_maps
	331	VALUES ".substr($values, 1));
	332	}
8c4a0c30	333	} else {
56670b6a	334	return false;
8c4a0c30	335	}
014c8464	336	return true;
014c8464	337	}
56670b6a	338
8c4a0c30	339	function set_smallest_levels()
	340	{
	341	$maxlengths = XDB::iterRow("SELECT MAX(LENGTH(gm.path)), gcim.city_id
	342	FROM geoloc_city_in_maps AS gcim
	343	INNER JOIN geoloc_maps AS gm USING ( map_id )
	344	GROUP BY gcim.city_id");
014c8464	345	while (list($length, $id) = $maxlengths->next()) {
8c4a0c30	346	XDB::execute("UPDATE geoloc_city_in_maps AS gcim
	347	INNER JOIN geoloc_maps AS gm USING(map_id)
	348	SET gcim.infos = IF(LENGTH(gm.path) = {?}, 'smallest', '')
	349	WHERE gcim.city_id = {?}", $length, $id);
014c8464	350	}
56670b6a	351	return true;
	352	}
	353	// }}}
	354
	355
8c4a0c30	356	function geoloc_to_x($lon, $lat)
	357	{
	358	return deg2rad(1) * $lon *100;
	359	}
	360
	361	function geoloc_to_y($lon, $lat)
	362	{
	363	if ($lat < -75) {
	364	return latToY(-75);
	365	}
	366	if ($lat > 75) {
	367	return latToY(75);
	368	}
787bb3d7	369	return -100 * log(tan(pi()/4 + deg2rad(1)/2*$lat));
8c4a0c30	370	}
56670b6a	371
8c4a0c30	372	function size_of_city($nb)
	373	{
	374	$s = round(log($nb + 1)*2,2);
	375	if ($s < 1) {
	376	return 1;
	377	}
	378	return $s;
56670b6a	379	}
56670b6a	380
8c4a0c30	381	function size_of_territory($nb)
	382	{
	383	return size_of_city($nb);
	384	}
56670b6a	385
8c4a0c30	386	function geoloc_getData_subcities($mapid, $SFields, &$cities, $direct=true)
	387	{
	388	if ($SFields instanceof UserSet) {
	389	$set = $SFields;
	390	$SFields = array();
	391	} else {
	392	$set = new UserSet();
	393	}
	394	for ($i_mapfield=0; $i_mapfield < count($SFields) ; $i_mapfield++) {
	395	if ($SFields[$i_mapfield]->fieldFormName == 'mapid') {
	396	break;
	397	}
	398	}
	399	$SFields[$i_mapfield] = new MapSField('mapid',
	400	array('gcim.map_id'),
	401	array('adresses','geoloc_city_in_maps'),
	402	array('am','gcim'),
	403	array(getadr_join('am'), 'am.cityid = gcim.city_id'),
	404	$mapid);
2b105fb6	405	$fields = new SFieldGroup(true, $SFields);
2b105fb6	406	$where = $fields->get_where_statement();
86b5c8f0	407	$joins = $fields->get_select_statement();
8c4a0c30	408	if ($where) {
	409	$where .= ' AND ';
	410	}
	411	$cityres = $set->get('gc.id,
	412	gc.lon / 100000 AS x, gc.lat/100000 AS y,
	413	gc.name,
	414	COUNT(u.user_id) AS pop,
	415	SUM(u.promo % 2) AS yellow',
86b5c8f0	416	"$joins
8c4a0c30	417	LEFT JOIN geoloc_city AS gc ON(gcim.city_id = gc.id)",
	418	$where . ($direct ? "gcim.infos = 'smallest'" : '1'),
	419	'gc.id, gc.alias',
	420	'pop DESC');
06df222b	421	foreach($cityres as $c) {
8c4a0c30	422	if ($c['pop'] > 0) {
2b105fb6	423	$city = $c;
	424	$city['x'] = geoloc_to_x($c['x'], $c['y']);
	425	$city['y'] = geoloc_to_y($c['x'], $c['y']);
	426	$city['size'] = size_of_city($c['pop']);
	427	$cities[$c['id']] = $city;
	428	}
8c4a0c30	429	}
2b105fb6	430	}
2b105fb6	431
787bb3d7	432	function geoloc_getData_subcountries($mapid, $sin, $minentities)
8c4a0c30	433	{
2b105fb6	434	$countries = array();
2b105fb6	435	$cities = array();
787bb3d7	436
8c4a0c30	437	if ($mapid === false) {
b536b5ae	438	$wheremapid = "WHERE gm.parent IS NULL";
8c4a0c30	439	} else {
b536b5ae	440	$wheremapid = "WHERE gm.parent = {?}";
8c4a0c30	441	}
08cce2ff	442	$submapres = XDB::iterator(
787bb3d7	443	"SELECT gm.map_id AS id, gm.name, gm.x, gm.y, gm.xclip, gm.yclip,
8c4a0c30	444	gm.width, gm.height, gm.scale, 1 AS rat
	445	FROM geoloc_maps AS gm
	446	". $wheremapid, Env::v('mapid',''));
2b105fb6	447
edfc872d	448	global $globals;
edfc872d	449
8c4a0c30	450	while ($c = $submapres->next()) {
2b105fb6	451	$country = $c;
2b105fb6	452	$country['color'] = 0xFFFFFF;
	453	$country['swf'] = $globals->geoloc->webservice_url."maps/mercator/map_".$c['id'].".swf";
	454	$countries[$c['id']] = $country;
	455	}
787bb3d7	456
8c4a0c30	457	if ($mapid === false) {
	458	return array($countries, $cities);
	459	}
2b105fb6	460
8c4a0c30	461	geoloc_getData_subcities(Env::i('mapid'), $sin, $cities);
b536b5ae	462	$nbcities = count($cities);
b536b5ae	463	$nocity = $nbcities == 0;
8c4a0c30	464	if ($sin instanceof UserSet) {
	465	$set = $sin;
	466	$SFields = array();
	467	} else {
	468	$set = new UserSet();
	469	$SFields = $sin;
	470	}
2b105fb6	471
8c4a0c30	472	for ($i_mapfield=0; $i_mapfield < count($SFields) ; $i_mapfield++) {
	473	if ($SFields[$i_mapfield]->fieldFormName == 'mapid') {
	474	break;
	475	}
	476	}
	477	$SFields[$i_mapfield] = new MapSField('mapid',
	478	array('map.parent'),
	479	array('adresses','geoloc_city_in_maps','geoloc_maps'),
	480	array('am','gcim','map'),
787bb3d7 FB	481	array(getadr_join('am'),
787bb3d7 FB	482	'am.cityid = gcim.city_id',
8c4a0c30	483	'map.map_id = gcim.map_id'));
b536b5ae	484	$fields = new SFieldGroup(true, $SFields);
86b5c8f0	485	$where = $fields->get_where_statement();
86b5c8f0	486	$joins = $fields->get_select_statement();
8c4a0c30	487	$countryres = $set->get('map.map_id AS id,
	488	COUNT(u.user_id) AS nbPop,
	489	SUM(u.promo % 2) AS yellow,
	490	COUNT(DISTINCT gcim.city_id) AS nbCities,
	491	SUM(IF(u.user_id IS NULL,0,am.glng)) AS lonPop,
	492	SUM(IF(u.user_id IS NULL, 0,am.glat)) AS latPop',
86b5c8f0	493	$joins,
86b5c8f0	494	$where,
8c4a0c30	495	'map.map_id',
8c4a0c30	496	'NULL');
787bb3d7	497
b536b5ae	498	$maxpop = 0;
06df222b FB	499	$nbentities = $nbcities + count($countryres);
06df222b FB	500	foreach ($countryres as $c) {
b536b5ae	501	$c['latPop'] /= $c['nbPop'];
	502	$c['lonPop'] /= $c['nbPop'];
	503	$c['rad'] = size_of_territory($c['nbPop']);
	504	if ($maxpop < $c['nbPop']) $maxpop = $c['nbPop'];
	505	$c['xPop'] = geoloc_to_x($c['lonPop'], $c['latPop']);
	506	$c['yPop'] = geoloc_to_y($c['lonPop'], $c['latPop']);
86b5c8f0	507	@$countries[$c['id']] = array_merge($countries[$c['id']], $c);
787bb3d7	508
b536b5ae	509	$nbcities += $c['nbCities'];
787bb3d7 FB	510	}
787bb3d7 FB	511
8c4a0c30	512	if ($nocity && $nbcities < $minentities){
8c4a0c30	513	foreach($countries as $i => $c) {
b536b5ae	514	$countries[$i]['nbPop'] = 0;
86b5c8f0	515	if (@$c['nbCities'] > 0) {
8c4a0c30	516	geoloc_getData_subcities($c['id'], $sin, $cities, false);
787bb3d7	517	}
8c4a0c30	518	}
b536b5ae	519	}
787bb3d7	520
8c4a0c30	521	foreach ($countries as $i => $c) {
c1895524	522	if (@$c['nbPop'] > 0) {
8c4a0c30	523	$lambda = pow($c['nbPop'] / $maxpop,0.3);
	524	$countries[$i]['color'] = 0x0000FF + round((1-$lambda) * 0xFF)*0x010100;
	525	}
b536b5ae	526	}
787bb3d7 FB	527
	528	return array($countries, $cities);
	529	}
138b3c8e	530	// }}}
138b3c8e	531
a7de4ef7	532	// vim:set et sw=4 sts=4 sws=4 foldmethod=marker enc=utf-8:
0337d704	533	?>