Rev 3422 | Blame | Compare with Previous | Last modification | View Log | RSS feed
<?php// declare(encoding='UTF-8');/*** Classe recherchant des infos sur une zone géo ou bien des coordonnées.** @internal Mininum PHP version : 5.2* @category CEL* @package Services* @subpackage Bibliothèques* @version 0.1* @author Mathias CHOUET <mathias@tela-botanica.org>* @author Jean-Pascal MILCENT <jpm@tela-botanica.org>* @author Aurelien PERONNET <aurelien@tela-botanica.org>* @license GPL v3 <http://www.gnu.org/licenses/gpl.txt>* @license CECILL v2 <http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt>* @copyright 1999-2014 Tela Botanica <accueil@tela-botanica.org>*/class RechercheInfosZoneGeo extends Cel {// @TODO: surveiller les conséquences du changement climatique :)const POINT_EMERGE_LE_PLUS_BAS = -422;// mètrespublic $maxTimeoutRequeteMondiale = 0;/*** Cas de la demande par lat et lon* array(* 'code_pays' => 'FR',* 'code_zone' => 'INSEE-D:35800',* 'type' => 'france',* 'nom' => 'Taulignan-les-oies'* );** array(* 'code_pays' => 'RU',* 'code_zone' => 'OSM-ID:4511574',* 'type' => 'monde',* 'nom' => 'Moscou'* );** Cas de la demande par nom* array(* 'lat': 40.372934,* 'lng': -3.690601,* 'nom': 'Madrid',* 'code_zone': 'OSM-ID:558474',* 'code_pays': 'ES',* 'type': 'monde'* );** array(* 'lat': 48.860799,* 'lng': 2.3457986,* 'nom': 'Paris',* 'code_zone': '75056',* 'code_pays': 'FR',* 'type': 'france'* );**///TODO: migrer tout ce qui concerne la zone géo dans cette classepublic function obtenirInfosPourCoordonnees($coordonnees) {$infos_coord = null;// Test facile qui permet d'abord de tenter une localisation bien plus rapide si les coordonnées// sont dans la bounding box approximative de la Franceif($this->testerCoordonneesWgsFrance($coordonnees['latitude'], $coordonnees['longitude'])) {$infos_coord = $this->chercherInfosCommune('osm', $coordonnees['latitude'], $coordonnees['longitude']);if ($infos_coord == null) {// Sinon recherche par pays$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);}} else {// Recherche par pays immédiate si en dehors de la bouding box française$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);}return array_merge($infos_coord,['lat' => $coordonnees['latitude'], 'lng' => $coordonnees['longitude']]);}// gère aussi les id_zone_geo @TODO renommer la méthode et internationaliser// la gestion des codespublic function obtenirInfosPourNom($nom, $pays, $code_departement=null, $id_zone=null) {$nom = trim($nom);$pays = trim($pays);$code_departement = trim($code_departement);$retour = false;// Dans le cas de la france on possède directement une table// avec les centroïdes des communesif($pays == 'FR' || $code_departement != "") {$retour = $this->chercherCentroideCommuneBdd($nom, $code_departement, $id_zone);} else {// Sinon méthode habituelle, un test simple suivi d'un test plus complexe// On cherche donc d'abord en France// @TODO c'est quoi ce if() mal ordonné ?if($pays == "") {$retour = $this->chercherCentroideCommuneBdd($nom, $code_departement, $id_zone);}// Si on ne trouve rien on teste dans le monde entierif($retour == false) {$retour = $this->effectuerRequeteGeocodingMondiale($nom, '2,3,4,5,6,7,8', 1, $pays, $id_zone);}}return $retour;}public function obtenirInfosPourCodeInseeCommune($code_insee) {$code_insee = trim($code_insee);$retour = $this->chercherCentroideCommuneInseeBdd($code_insee);return $retour;}/*** Renvoie tous les codes de zones pour le groupe $groupe*/protected function obtenirZonesDuGroupe($groupe) {$requete = 'SELECT valeur as id_zone FROM cel_groupes_zones '. 'WHERE id_groupe = ' . Cel::db()->proteger($groupe). ' -- '.__FILE__.':'.__LINE__;$liste_zones = Cel::db()->requeter($requete);return $liste_zones;}/*** Renvoie des infos pour un groupes de zones, déclaré dans la table* cel_groupe_zones : toutes les infos des zones concernées, plus le carré* englobant (bounding box) et le centroïde de l'union des zones*/public function obtenirInfosPourGroupeZonesFrance($groupe) {$groupe = trim($groupe);$retour = false;$zones = $this->obtenirZonesDuGroupe($groupe);$retour = array('groupe' => array(),'zones' => array());foreach ($zones as $zone) {$idZone = $zone['id_zone'];$retour['zones'][] = $this->obtenirInfosPourNom(null, null, null, $idZone);}// recherche bbox (de Bouygues Telecom) @TODO faire ça au dessus, en une passe ?$lngmin = 99999;$latmin = 99999;$lngmax = -99999;$latmax = -99999;foreach ($retour['zones'] as $zone) {$lat = $zone['lat'];$lng = $zone['lng'];if ($lat < $latmin) { $latmin = $lat; }if ($lat > $latmax) { $latmax = $lat; }if ($lng < $lngmin) { $lngmin = $lng; }if ($lng > $lngmax) { $lngmax = $lng; }}$retour['groupe']['bbox'] = array('lat_min' => $latmin,'lat_max' => $latmax,'lng_min' => $lngmin,'lng_max' => $lngmax);// @TODO suggérer un niveau de zoom en fonction de la plus grande dimension// de la bounding-box (ne marchera pas dans tous les cas mais ce serait// déjà pas mal)$retour['groupe']['centroide'] = array('lat' => ($latmin + $latmax) / 2,'lng' => ($lngmin + $lngmax) / 2);return $retour;}public function chercherInfosCommune($projet, $latitude, $longitude) {$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;$url_service = $this->config['cel']['url_service_geo_local'].$valeurs_params;$url_service = str_replace(',', '.', $url_service);$ch = curl_init($url_service);curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);$reponse = curl_exec($ch);$infos_coord = $this->traiterReponseServiceCommune($reponse);curl_close($ch);return $infos_coord;}public function chercherInfosPays($latitude, $longitude) {$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;$url_service = $this->config['cel']['url_service_geo_mondial'].$valeurs_params;$url_service = str_replace(',', '.', $url_service);$ch = curl_init($url_service);curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);$reponse = curl_exec($ch);$res = json_decode($reponse, true);$retour = null;$retour = array('code_zone' => '','code_pays' => '','nom' => '','type' => '');if(!empty($res)) {// Cas où l'on a recherché des coordonnées en France// On en profite pour remplir plus de données// que le simple pays si elles sont présentesif(isset($res[8]) && isset($res[8]['codeInsee'])) {$retour = array('code_zone' => $res[8]['codeInsee'],'code_pays' => 'FR','nom' => $res[8]['nom'],'type' => 'france');} else {// Cas de la recherche en dehors de France// La zone de plus haut niveau est toujours un pays// (car le niveau de zone est limité à 2)$infos_pays = $res[min(array_keys($res))];// La zone de niveau le plus bas est la "localité"// la plus précise qu'on a pu trouver$infos_localite = $res[max(array_keys($res))];// Cas où l'on a trouvé un code paysif(!empty($infos_pays['codeIso31661'])) {$retour = array('code_pays' => $infos_pays['codeIso31661'],'nom' => '');} elseif(!empty($infos_pays['codeIso31662'])) {// Quelquefois la zone de plus haut niveau est une région ou une province// heureusement son code est de forme XX-YY ou XX est le code iso du pays !$retour = array('code_pays' => substr($infos_pays['codeIso31662'], 0, 2),'nom' => '');}// Pas de sens de prendre un nom de localité si on a pas de paysif(!empty($retour['code_pays'])) {// Type sert à savoir quelle fonction de mise à jour on appellera$retour['type'] = 'monde';$retour['nom'] = $infos_localite['intitule'];$retour['code_zone'] = $infos_localite['codeZoneGeo'];}}}return $retour;}public function traiterReponseServiceCommune($reponse) {$retour = null;$reponse = json_decode($reponse);// cas du service lion 1906 qui peut renvoyer plusieurs communes (on prend la première)if (is_array($reponse)) {$reponse = $reponse[0];}if (isset($reponse->codeINSEE)) {// Type sert à savoir quelle fonction de mise à jour on appellera$retour = array('code_zone' => $reponse->codeINSEE,'code_pays' => 'FR','nom' => $reponse->nom,'type' => 'france');}return $retour;}// @TODO gérer $id_zone (inactif pour l'instant)protected function effectuerRequeteGeocodingMondiale($nom, $niveau = '2,3,4,5,6,7,8', $limite = 1, $pays = null, $id_zone=null) {$url_sans_pays = $this->config['cel']['url_service_geo_mondial'].'?masque='.urlencode($nom).'&niveau='.$niveau.'&limite='.$limite;$url = $url_sans_pays.(!empty($pays) ? '&pays='.urlencode($pays) : '');$ch = curl_init($url);curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);curl_setopt($ch, CURLOPT_TIMEOUT, $this->maxTimeoutRequeteMondiale);$reponse = curl_exec($ch);$res = @json_decode($reponse, true);$retour = array();if(!empty($res)) {// cas de la localisation directe d'une communeif($limite == 1) {$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);} else {// cas de l'autocompletion$retour = array();foreach($res as $ligne) {$ligne_decodee = $this->traiterLigneResultatRequeteGeocodingMondiale($ligne);$param = array('nom' => $ligne_decodee['nom'], 'code' => $ligne_decodee['code_pays']);$retour[] = $this->formaterLigneResultat($param);}}} else if($limite == 1) {// Dans le cas où un pays a été demandé mais pas de résultat, on rentente la recherche sans// préciser de pays$res = @json_decode(file_get_contents($url_sans_pays), true);if(!empty($res)) {$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);}}return $retour;}/*** Renvoie les infos d'une commune, incluant les coordonnées du centroïde,* en fonction du nom et/ou du département, et/ou de l'id_zone (plus fiable)*/protected function chercherCentroideCommuneBdd($commune, $departement=null, $id_zone=null) {$commune_formatee = str_replace(array(' ', '-'), '_', trim($commune));// clauses$clauses = array();if ($commune != null) {$clauses[] = 'nom LIKE '.Cel::db()->proteger($commune_formatee);}if ($departement != null) {$departement = trim($departement, "*");if (strlen($departement) > 2) {// @TODO piège : si deux communes ont le même nom dans deux// DROM-COM, la troncature à 2 chiffres sera insuffisante$departement = substr($departement, 0, 2);}$clauses[] = 'code LIKE '.Cel::db()->proteger($departement.'%');}if ($id_zone != null) {$clauses[] = 'id_zone_geo = '.Cel::db()->proteger($id_zone);}$clauses = implode(' AND ', $clauses);$requete = 'SELECT utm_x, utm_y, utm_secteur, code FROM cel_zones_geo '. 'WHERE ' . $clauses. ' -- '.__FILE__.':'.__LINE__;$commune_coordonnees = Cel::db()->requeter($requete);$retour = false;if ($commune_coordonnees && is_array($commune_coordonnees) && count($commune_coordonnees) > 0) {$lat_lon = $this->convertirUtmVersLatLong($commune_coordonnees[0]['utm_x'],$commune_coordonnees[0]['utm_y'],$commune_coordonnees[0]['utm_secteur']);$retour = array('lat' => (float) $lat_lon['lat'],'lng' => (float) $lat_lon['lng'],'nom' => $commune,'code_zone' => $commune_coordonnees[0]['code'],'code_pays' => 'FR','type' => 'france');}return $retour;}/*** Renvoie les infos d'une commune, incluant les coordonnées du centroïde,* en fonction du code insee (plus fiable)*/protected function chercherCentroideCommuneInseeBdd($code_insee) {$code_insee = trim($code_insee);$commune_coordonnees = array();if ($code_insee != null) {$requete = 'SELECT `wgs84_latitude`, `wgs84_longitude`, `nom` FROM `cel_zones_geo` '. 'WHERE `code`=' .Cel::db()->proteger($code_insee). ' -- '.__FILE__.':'.__LINE__;$commune_coordonnees = Cel::db()->requeter($requete);}$retour = false;if ($commune_coordonnees && is_array($commune_coordonnees) && count($commune_coordonnees) > 0) {$retour = array('lat' => (float) $commune_coordonnees[0]['wgs84_latitude'],'lng' => (float) $commune_coordonnees[0]['wgs84_longitude'],'nom' => $commune_coordonnees[0]['nom'],'code_zone' => $code_insee,'code_pays' => 'FR','type' => 'france');}return $retour;}public function obtenirListeInfosPourNom($lieu) {$lieux_retour = $this->obtenirListeCommunes($lieu);// S'il y a moins de 50 résultats on tente des les fusionner avec une recherche mondiale// avec 50 - n éléments (où n est la taille du tableau de communes)if(count($lieux_retour) < 50) {$lieux_mondiaux = $this->effectuerRequeteGeocodingMondiale($lieu, '2,3,4,5,6,7,8', 50 - count($lieux_retour));usort($lieux_mondiaux, array($this, 'trierListeLieux'));$lieux_retour += $lieux_mondiaux;}// Suppression de possibles doublons @WARNING marche pas !!!//$lieux_retour = $this->reduireListeLieux($lieux_retour);return $lieux_retour;}public function obtenirListeCommunes($lieu) {$lieu_formate = $this->formaterChaineLieuPourRequete($lieu);$retour = array();if ($this->estUneChaineRequeteValide($lieu_formate)) {$requete = 'SELECT DISTINCT nom, code '.'FROM cel_zones_geo '.'WHERE nom LIKE '.Cel::db()->proteger($lieu_formate.'%').' '.'ORDER BY nom '.'LIMIT 50 '.' -- '.__FILE__.':'.__LINE__;$liste_lieux = Cel::db()->requeter($requete);if ($liste_lieux) {foreach ($liste_lieux as $lieu_trouve) {$retour[] = $this->formaterLigneResultat($lieu_trouve);}}}return $retour;}/*** Fonctions utilitaires*/// qui a écrit cette m*rde ?protected function reduireListeLieux($tableau) {$index = array();foreach($tableau as $ligne) {// merci de mettre des commentaires !!!// ça enlève les 2 premières lettres ou chiffres du nom de lieu => wtf ?$index_loc = trim(preg_replace("/\([a-zA-Z0-9]{2}\)$/", "", $ligne[0]));if(!isset($index[$index_loc])) {$index[$index_loc] = $ligne;} else {if(mb_strlen($ligne[0]) > mb_strlen($index[$index_loc][0])) {$index[$index_loc] = $ligne;} else if(mb_strlen($ligne[0]) == mb_strlen($index[$index_loc][0])) {if(preg_match("/(.)* ([a-zA-Z0-9]{2})/", $ligne[0])) {$index[$index_loc] = $ligne;}}}}return array_values($index);}protected function trierListeLieux($a, $b) {return strcmp($a[0], $b[0]);}protected function testerCoordonneesWgsFrance($latitude, $longitude) {$coord_france = false;if ($latitude != '' && $longitude != '') {if ($latitude < 51.071667 && $latitude > 41.316667) {if ($longitude < 9.513333 && $longitude > -5.140278) {$coord_france = true;}}}return $coord_france;}protected function traiterLigneResultatRequeteGeocodingMondiale($res) {$lat = "";$lon = "";$nom = "";$code ="";$code_pays="";$type="";if(!empty($res)) {$code_pays = $this->extraireCodePays($res);$nom = $res['intitule'];$lat = $res['centre_lat'];$lon = $res['centre_lng'];$type = $code_pays == 'FR' || $code_pays == 'FX' ? 'france' : 'monde';$retour = array('lat' => (float)$lat,'lng' => (float)$lon,'nom' => $nom,'code_zone' => $code,'code_pays' => $code_pays,'type' => $type);}return $retour;}protected function convertirUtmVersLatLong($x, $y, $sector) {$convertisseur = new gPoint();$convertisseur->setUTM($x, $y, $sector);$convertisseur->convertTMtoLL();$lat_long = array();$lat_long['lat'] = str_replace(',','.',$convertisseur->Lat());$lat_long['lng'] = str_replace(',','.',$convertisseur->Long());return $lat_long;}protected function extraireCodePays($infos_pays) {$code = "";if(!empty($infos_pays['codeIso31661'])) {$code = $infos_pays['codeIso31661'];} elseif(!empty($infos_pays['codeIso31662'])) {// Si pas de code pays, le sous code peut exister et content le code pays// sous la forme codepays-souscode$code = substr($infos_pays['codeIso31662'], 0, 2);}return $code;}protected function estUneChaineRequeteValide($lieu) {return (strlen($lieu) > 0) && ($lieu != '%');}protected function formaterLigneResultat($ligne) {if(empty($ligne['code']) || $ligne['code'] == null) {$res = array($ligne['nom'], '');} else {$res = array($ligne['nom'].' ('.substr(sprintf('%02s', $ligne['code']),0,2).')', $ligne['code']);}return $res;}/*** Remplace les * par % pour faire des recherches floues* Remplace les + par _ (nginx envoie des "+" dans l'URL à la place des espaces)* Remplace les espaces et les - par _ car les noms de communes peuvent avoir des espaces ou des tirets* @param string $lieu* @return string le lieu formaté pour la recherche*/protected function formaterChaineLieuPourRequete($lieu) {$lieu = ltrim($lieu);$lieu = preg_replace('/\*+/', '%', $lieu);$lieu = str_replace(['+', ' ', '-'], '_', $lieu);return $lieu;}public function obtenirAltitude($donnees) {$altitude = [];$api = '';$parametres_coordonnees = '';if(empty($donnees) || empty($donnees['code_pays'])) {return $altitude;}$pays = $donnees['code_pays'];$lat = $donnees['lat'] ?? $donnees['latitude'] ?? null;$lon = $donnees['lon'] ?? $donnees['lng'] ?? $donnees['longitude'] ?? null;if (!($lat && $lon)) {return $altitude;}if ('FR' === $donnees['code_pays']) {$api = 'ign';$parametres_coordonnees = http_build_query(['lat' => $lat, 'lon' => $lon]);} else {$api = 'mapquest';$parametres_coordonnees = http_build_query(['latLngCollection' => $lat.','.$lon]);}$url_service = sprintf($this->config[$api]['url_service_altitude'], $this->config[$api]['cle_api'], $parametres_coordonnees);$ch = curl_init($url_service);if ('ign' === $api) {// pas de header Authorization pour mapquestcurl_setopt($ch, CURLOPT_USERPWD, $this->config[$api]['username'] . ":" . $this->config[$api]['password']);}curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);$reponse = curl_exec($ch);curl_close($ch);$altitude = $this->traiterReponseServiceAltitude($api, $reponse);return ['altitude' => $altitude];}public function traiterReponseServiceAltitude($api, $reponse) {$reponse = json_decode($reponse);$elevationArrayKey = '';$elevationKey = '';$altitude = null;// on a de la chance les réponses de ces 2 api sont structurées de la même façonif ('ign' === $api) {$elevationArrayKey = 'elevations';$elevationKey = 'z';} elseif ('mapquest' === $api) {$elevationArrayKey = 'elevationProfile';$elevationKey = 'height';}// du coup on peut faire çaif (isset($reponse->{$elevationArrayKey})&& is_array($reponse->{$elevationArrayKey})&& isset($reponse->{$elevationArrayKey}[0]->{$elevationKey})// lorsque les coordonnées sortent de leur rayon d'action, ces APIs// renvoient bien quelque-chose: -32768 pour mapquest, et -99999 pour l'IGN&& POINT_EMERGE_LE_PLUS_BAS < floatval($reponse->{$elevationArrayKey}[0]->{$elevationKey})) {$altitude = $reponse->{$elevationArrayKey}[0]->{$elevationKey};}return $altitude;}}