Subversion Repositories eFlore/Applications.cel

Compare Revisions

Ignore whitespace Rev 3856 → Rev 3857

/branches/v3.01-serpe/jrest/bibliotheque/RechercheInfosZoneGeo.php
New file
0,0 → 1,629
<?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ètres
 
public $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 classe
public 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 France
if($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 codes
public 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 communes
if($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 entier
if($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ésentes
if(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 pays
if(!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 pays
if(!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 commune
if($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 mapquest
curl_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çon
if ('ign' === $api) {
$elevationArrayKey = 'elevations';
$elevationKey = 'z';
} elseif ('mapquest' === $api) {
$elevationArrayKey = 'elevationProfile';
$elevationKey = 'height';
}
// du coup on peut faire ça
if (
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;
}
}