* @author Jean-Pascal MILCENT * @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('
', $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; } } ?>