New file |
0,0 → 1,200 |
<?php |
// declare(encoding='UTF-8'); |
/** |
* Classe permettant d'obtenir le nom et le code INSEE d'une commune à partir de ces coordonnées (latitude et longitude). |
* La latitude et longitude doivent être exprimée par un nombre décimal. |
* Ce service fonctionne uniquement pour les communes de France métropolitaine (Corse comprise) dont les contours ont été tracés dans OpenStreetMap. |
* Source des données : OpenStreetMap http://wiki.openstreetmap.org (Projet France limites administratives : http://wiki.openstreetmap.org/wiki/WikiProject_France/Limites_administratives) |
* Paramètres du service : |
* - lat : latitude |
* - lon : longitude |
* Exemple : |
* http://localhost/cartoOSM/services/0.1/nom-commune?lat=44.71546&lon=3.84216 |
* |
* @category php 5.2 |
* @package carto-osm |
* @author Mohcen BENMOUNAH <mohcen@tela-botanica.org> |
* @author Jean-Pascal MILCENT <jpm@tela-botanica.org> |
* @copyright Copyright (c) 2010, Tela Botanica (accueil@tela-botanica.org) |
* @license http://www.cecill.info/licences/Licence_CeCILL_V2-fr.txt Licence CECILL |
* @license http://www.gnu.org/licenses/gpl.html Licence GNU-GPL |
* @version $Id$ |
*/ |
class NomCommune { |
|
const PATTERN_LAT = '/^[0-9]+(?:[.][0-9]+|)$/'; |
const PATTERN_LON = '/^[-]?[0-9]+(?:[.][0-9]+|)$/'; |
const LAT_MAX = 51.071667; |
const LAT_MIN = 41.316667; |
const LON_MAX = 9.513333; |
const LON_MIN = -5.140278; |
const NBRE_COMMUNE_PAR_DEFAUT = 10; |
const NBRE_COMMUNE_MAX = 100; |
const MIME_JSON = 'application/json'; |
|
private $parametres = array(); |
private $bdd = null; |
|
public function __construct(Bdd $bdd) { |
$this->bdd = $bdd; |
} |
|
public function consulter($ressources, $parametres) { |
$this->parametres = $parametres; |
$this->verifierParametres(); |
|
$nomINSEEs = $this->trouverCommunesProches(); |
$corps = $this->formaterResultats($nomINSEEs); |
|
$resultat = new ResultatService(); |
$resultat->mime = self::MIME_JSON; |
$resultat->corps = $corps; |
return $resultat; |
} |
|
private function verifierParametres() { |
extract($this->parametres); |
$messages = array(); |
if (! array_key_exists('lat', $this->parametres)) { |
$messages[] = "Vous devez indiquer une latitude en degré décimal à l'aide du paramètres d'url : lat"; |
} else if (!preg_match(self::PATTERN_LAT, $lat)) { |
$messages[] = "La valeur de latitude doit être un nombre décimal positif dont le séparateur décimal est un point. Ex. : 44 ou 43.03"; |
} else if ($lat > self::LAT_MAX) { |
$messages[] = "La valeur de latitude indiquée est supérieure à {self::LAT_MAX} qui est le point le plus au Nord de la France métropolitaine."; |
} else if ($lat < self::LAT_MIN) { |
$messages[] = "La valeur de latitude indiquée est infèrieure à {self::LAT_MIN} qui est le point le plus au Sud de la France métropolitaine."; |
} |
if (! array_key_exists('lon', $this->parametres)) { |
$messages[] = "Vous devez indiquer une longitude en degré décimal à l'aide du paramètres d'url : lon"; |
} else if (!preg_match(self::PATTERN_LON, $lon)) { |
$messages[] = "La valeur de longitude doit être un nombre décimal dont le séparateur décimal est un point. Ex. : -4.03 ou 3.256"; |
} else if ($lon > self::LON_MAX) { |
$messages[] = "La valeur de longitude indiquée est supérieure à {self::LON_MAX} qui est le point le plus à l'Est de la France métropolitaine."; |
} else if ($lon < self::LON_MIN) { |
$messages[] = "La valeur de longitude indiquée est infèrieure à {self::LON_MIN} qui est le point le plus à l'Ouest de la France métropolitaine."; |
} |
if (count($messages) != 0) { |
$message = implode('<br />', $messages); |
$code = RestServeur::HTTP_CODE_MAUVAISE_REQUETE; |
throw new Exception($message, $code); |
} |
} |
|
/** |
* requête qui récupère les 20 communes les plus proches du point recherché |
* La distance(AB = \sqrt{(x_B-x_A)^2 + (y_B-y_A)^2}) est calculée sans la racine |
* (calcul en plus qui change pas le résultat). |
*/ |
private function trouverCommunesProches() { |
$lat = $this->parametres['lat']; |
$lon = $this->parametres['lon']; |
$requete = 'SELECT '. |
" (({$lat} - X(centre)) * ({$lat} -X(centre)) + ({$lon} - Y(centre)) * ({$lon} - Y(centre))) AS distance, ". |
" code_insee, nom, ASTEXT(polygon) AS poly ". |
'FROM osm_communes '. |
'WHERE note = "Polygone complet" '. |
'ORDER BY distance '. |
'LIMIT 20'; |
$resultat = $this->bdd->recupererTous($requete); |
|
if ($resultat == false) { |
$message = "Le service '".get_class($this)."' n'a trouvé aucune commune dont le centroïde correspond aux coordonnées : {$this->parametres['lat']}, {$this->parametres['lon']}."; |
$code = RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE; |
throw new Exception($message, $code); |
} |
return $resultat; |
} |
|
private function formaterResultats($nomINSEEs) { |
$communeTrouvee = null; |
foreach ($nomINSEEs as $nomINSEE) { |
$resultat = $this->localiserPointLatLon($nomINSEE['poly']); |
if ($resultat == 1) { |
$communeTrouvee = array('nom' => $nomINSEE['nom'], 'codeINSEE' => $nomINSEE['code_insee']); |
break; |
} |
} |
if (is_null($communeTrouvee)) { |
$message = "Le service '".get_class($this)."' n'a trouvé aucune commune correspondant aux coordonnées : {$this->parametres['lat']}, {$this->parametres['lon']}."; |
$code = RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE; |
throw new Exception($message, $code); |
} |
return $communeTrouvee; |
} |
|
private function localiserPointLatLon($multiPolygone) { |
$lat = $this->parametres['lat']; |
$lon = $this->parametres['lon']; |
$requeteMBR = "SELECT MBRWithin(GEOMFROMTEXT('POINT($lat $lon)'), GEOMFROMTEXT('$multiPolygone'))"; |
$p = array_values($this->bdd->recuperer($requeteMBR)); |
|
$counter = 0; |
if ($p[0] != 1) { |
return 0; |
} else { |
$x = $lat; //Le point latitude recherché |
$y = $lon; //Le point longitude recherché |
$str = $this->nettoyerChaineMultipolygone($multiPolygone); // Le multipolygone de la commune |
$pb = 0; |
$pe = strpos($str, ','); //définition de la position de la première virgule |
$xy = substr($str, $pb, $pe - $pb); //extraire X Y d'un point séparé par des virgules |
$p = strpos($xy, ' '); //position espace dans un point |
$p1x = substr($xy, 0, $p ); //récupère X du point |
$p1y = substr($xy, $p + 1 ); //récupère Y du point |
$str = $str.$xy; // ajouter le premier point du polygone à la fin du polygone |
while ($pe > 0) { |
$xy = substr($str, $pb, $pe - $pb ); //extraire le point GPS |
$p = strpos($xy, ' '); // trouver la position de l'espace entre les deux points X et Y |
$p2x = substr($xy,0, $p ); // extraire le point X |
$p2y = substr($xy, $p + 1 ); // extraire le point Y |
if ($p1y < $p2y) { |
$m = $p1y; |
} else { |
$m = $p2y; |
} |
if ($y > $m) { |
if ($p1y > $p2y) { |
$m = $p1y; |
} else { |
$m = $p2y; |
} |
if ($y <= $m) { |
if ($p1x > $p2x) { |
$m = $p1x; |
} else { |
$m = $p2x; |
} |
if ($y <= $m) { |
if ($p1x > $p2x) { |
$m = $p1x; |
} else { |
$m = $p2x; |
} |
if ($x <= $m) { |
if ($p1y != $p2y) { |
$xinters = ($y - $p1y) * ($p2x - $p1x) / ($p2y - $p1y) + $p1x; |
} |
if (($p1x = $p2x) || ($x <= $xinters )) { |
$counter++; |
} |
} |
} |
} |
} |
$p1x = $p2x; |
$p1y = $p2y; |
$pb = $pe + 1; |
$pe = strpos($str, ',', $pb); |
} |
return $counter % 2; |
} |
} |
|
private function nettoyerChaineMultipolygone($chaine) { |
$chaine = str_replace('MULTIPOLYGON(((', '', $chaine); |
$chaine = str_replace(')))', '', $chaine); |
$chaine .= ','; |
return $chaine; |
} |
} |
?> |