4,12 → 4,13 |
* 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) |
* 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 |
* http://localhost/service-test:eflore:0.1/osm/nom-commune?lon=3.83619&lat=44.74231 |
* |
* @category php 5.2 |
* @package carto-osm |
28,13 → 29,14 |
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; |
|
private $lat = null; |
private $lon = null; |
|
public function __construct(Bdd $bdd) { |
$this->bdd = $bdd; |
} |
42,9 → 44,12 |
public function consulter($ressources, $parametres) { |
$this->parametres = $parametres; |
$this->verifierParametres(); |
$this->lat = $this->parametres['lat']; |
$this->lon = $this->parametres['lon']; |
|
$nomINSEEs = $this->trouverCommunesProches(); |
$corps = $this->formaterResultats($nomINSEEs); |
$communesProches = $this->trouverCommunesProches(); |
$communeTrouveeInfos = $this->localiserPointLatLon($communesProches); |
$corps = $this->formaterResultats($communeTrouveeInfos); |
|
$resultat = new ResultatService(); |
$resultat->mime = self::MIME_JSON; |
60,9 → 65,9 |
} 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."; |
$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."; |
$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"; |
69,9 → 74,9 |
} 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."; |
$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."; |
$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); |
83,118 → 88,58 |
/** |
* 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). |
* (calcul en plus qui change pas le résultat). La distance n'est donc pas juste mais permet de récupérer |
* le bon ordre. |
*/ |
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 ". |
$requete = 'SELECT id_osm, '. |
"(({$this->lon} - X(centre)) * ({$this->lon} -X(centre)) + ({$this->lat} - Y(centre)) * ({$this->lat} - Y(centre))) AS distance ". |
'FROM osm_communes '. |
'WHERE note = "Polygone complet" '. |
'ORDER BY distance '. |
'LIMIT 20'; |
'ORDER BY distance ASC '. |
'LIMIT 20 '. |
' -- '.__FILE__.' : '.__LINE__; |
$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); |
$msgTpl = "Service '%s' : aucune commune dont le centroïde correspond aux coordonnées : %s, %s."; |
$msg = sprintf($msgTpl, get_class($this), $this->lat, $this->lon); |
throw new Exception($msg, RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE); |
} |
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 formaterResultats($communeTrouveeInfos) { |
return array( |
'nom' => $communeTrouveeInfos['nom'], |
'codeINSEE' => $communeTrouveeInfos['code_insee'], |
'wikipedia' => $this->formaterUrlWikipedia($communeTrouveeInfos['wikipedia']) |
); |
} |
|
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 localiserPointLatLon($communesProches) { |
$zonesIds = array(); |
foreach ($communesProches as $commune) { |
$zonesIds[] = $commune['id_osm']; |
} |
$idsOsmClauseIn = implode(',', $this->bdd->proteger($zonesIds)); |
$requete = "SELECT id_osm, code_insee, nom, wikipedia ". |
'FROM osm_communes '. |
"WHERE st_within(GEOMFROMTEXT('POINT($this->lon $this->lat)'), polygone) = 1 ". |
"AND id_osm IN ($idsOsmClauseIn) ". |
' -- '.__FILE__.' : '.__LINE__; |
$resultat = $this->bdd->recuperer($requete); |
if ($resultat === false) { |
$msgTpl = "Service '%s' : aucune commune correspondant aux coordonnées : %s, %s."; |
$msg = sprintf($msgTpl, get_class($this), $this->lat, $this->lon); |
throw new Exception($msg, RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE); |
} |
return $resultat; |
} |
|
private function nettoyerChaineMultipolygone($chaine) { |
$chaine = str_replace('MULTIPOLYGON(((', '', $chaine); |
$chaine = str_replace(')))', '', $chaine); |
$chaine .= ','; |
return $chaine; |
private function formaterUrlWikipedia($infoWikipedia) { |
list($lang, $page) = explode(':', $infoWikipedia); |
$pageEncode = rawurlencode($page); |
$urlTpl = 'https://%s.wikipedia.org/wiki/%s'; |
return sprintf($urlTpl, $lang, $pageEncode); |
} |
} |
?> |
} |