Subversion Repositories eFlore/Projets.eflore-projets

<?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;
        }
}
?>