Subversion Repositories eFlore/Projets.eflore-projets

Compare Revisions

Ignore whitespace Rev 163 → Rev 169

/trunk/services/modules/0.1/osm/NomCommune.php
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;
}
}
?>