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;

	}

}

?>