169 |
jpm |
1 |
<?php
|
|
|
2 |
// declare(encoding='UTF-8');
|
|
|
3 |
/**
|
|
|
4 |
* Classe permettant d'obtenir le nom et le code INSEE d'une commune à partir de ces coordonnées (latitude et longitude).
|
|
|
5 |
* La latitude et longitude doivent être exprimée par un nombre décimal.
|
|
|
6 |
* Ce service fonctionne uniquement pour les communes de France métropolitaine (Corse comprise) dont les contours ont été tracés dans OpenStreetMap.
|
1056 |
jpm |
7 |
* Source des données : OpenStreetMap http://wiki.openstreetmap.org
|
|
|
8 |
* Projet France limites administratives : http://wiki.openstreetmap.org/wiki/WikiProject_France/Limites_administratives
|
169 |
jpm |
9 |
* Paramètres du service :
|
|
|
10 |
* - lat : latitude
|
|
|
11 |
* - lon : longitude
|
|
|
12 |
* Exemple :
|
1056 |
jpm |
13 |
* http://localhost/service-test:eflore:0.1/osm/nom-commune?lon=3.83619&lat=44.74231
|
169 |
jpm |
14 |
*
|
|
|
15 |
* @category php 5.2
|
|
|
16 |
* @package carto-osm
|
|
|
17 |
* @author Mohcen BENMOUNAH <mohcen@tela-botanica.org>
|
|
|
18 |
* @author Jean-Pascal MILCENT <jpm@tela-botanica.org>
|
|
|
19 |
* @copyright Copyright (c) 2010, Tela Botanica (accueil@tela-botanica.org)
|
|
|
20 |
* @license http://www.cecill.info/licences/Licence_CeCILL_V2-fr.txt Licence CECILL
|
|
|
21 |
* @license http://www.gnu.org/licenses/gpl.html Licence GNU-GPL
|
|
|
22 |
* @version $Id$
|
|
|
23 |
*/
|
|
|
24 |
class NomCommune {
|
|
|
25 |
|
|
|
26 |
const PATTERN_LAT = '/^[0-9]+(?:[.][0-9]+|)$/';
|
|
|
27 |
const PATTERN_LON = '/^[-]?[0-9]+(?:[.][0-9]+|)$/';
|
|
|
28 |
const LAT_MAX = 51.071667;
|
|
|
29 |
const LAT_MIN = 41.316667;
|
|
|
30 |
const LON_MAX = 9.513333;
|
|
|
31 |
const LON_MIN = -5.140278;
|
|
|
32 |
const MIME_JSON = 'application/json';
|
|
|
33 |
|
|
|
34 |
private $parametres = array();
|
|
|
35 |
private $bdd = null;
|
|
|
36 |
|
1056 |
jpm |
37 |
private $lat = null;
|
|
|
38 |
private $lon = null;
|
|
|
39 |
|
169 |
jpm |
40 |
public function __construct(Bdd $bdd) {
|
|
|
41 |
$this->bdd = $bdd;
|
|
|
42 |
}
|
|
|
43 |
|
|
|
44 |
public function consulter($ressources, $parametres) {
|
|
|
45 |
$this->parametres = $parametres;
|
|
|
46 |
$this->verifierParametres();
|
1056 |
jpm |
47 |
$this->lat = $this->parametres['lat'];
|
|
|
48 |
$this->lon = $this->parametres['lon'];
|
169 |
jpm |
49 |
|
1056 |
jpm |
50 |
$communesProches = $this->trouverCommunesProches();
|
|
|
51 |
$communeTrouveeInfos = $this->localiserPointLatLon($communesProches);
|
|
|
52 |
$corps = $this->formaterResultats($communeTrouveeInfos);
|
169 |
jpm |
53 |
|
|
|
54 |
$resultat = new ResultatService();
|
|
|
55 |
$resultat->mime = self::MIME_JSON;
|
|
|
56 |
$resultat->corps = $corps;
|
|
|
57 |
return $resultat;
|
|
|
58 |
}
|
|
|
59 |
|
|
|
60 |
private function verifierParametres() {
|
|
|
61 |
extract($this->parametres);
|
|
|
62 |
$messages = array();
|
|
|
63 |
if (! array_key_exists('lat', $this->parametres)) {
|
|
|
64 |
$messages[] = "Vous devez indiquer une latitude en degré décimal à l'aide du paramètres d'url : lat";
|
|
|
65 |
} else if (!preg_match(self::PATTERN_LAT, $lat)) {
|
|
|
66 |
$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";
|
|
|
67 |
} else if ($lat > self::LAT_MAX) {
|
1056 |
jpm |
68 |
$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.";
|
169 |
jpm |
69 |
} else if ($lat < self::LAT_MIN) {
|
1056 |
jpm |
70 |
$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.";
|
169 |
jpm |
71 |
}
|
|
|
72 |
if (! array_key_exists('lon', $this->parametres)) {
|
|
|
73 |
$messages[] = "Vous devez indiquer une longitude en degré décimal à l'aide du paramètres d'url : lon";
|
|
|
74 |
} else if (!preg_match(self::PATTERN_LON, $lon)) {
|
|
|
75 |
$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";
|
|
|
76 |
} else if ($lon > self::LON_MAX) {
|
1056 |
jpm |
77 |
$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.";
|
169 |
jpm |
78 |
} else if ($lon < self::LON_MIN) {
|
1056 |
jpm |
79 |
$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.";
|
169 |
jpm |
80 |
}
|
|
|
81 |
if (count($messages) != 0) {
|
|
|
82 |
$message = implode('<br />', $messages);
|
|
|
83 |
$code = RestServeur::HTTP_CODE_MAUVAISE_REQUETE;
|
|
|
84 |
throw new Exception($message, $code);
|
|
|
85 |
}
|
|
|
86 |
}
|
|
|
87 |
|
|
|
88 |
/**
|
|
|
89 |
* requête qui récupère les 20 communes les plus proches du point recherché
|
|
|
90 |
* La distance(AB = \sqrt{(x_B-x_A)^2 + (y_B-y_A)^2}) est calculée sans la racine
|
1056 |
jpm |
91 |
* (calcul en plus qui change pas le résultat). La distance n'est donc pas juste mais permet de récupérer
|
|
|
92 |
* le bon ordre.
|
169 |
jpm |
93 |
*/
|
|
|
94 |
private function trouverCommunesProches() {
|
1056 |
jpm |
95 |
$requete = 'SELECT id_osm, '.
|
|
|
96 |
"(({$this->lon} - X(centre)) * ({$this->lon} -X(centre)) + ({$this->lat} - Y(centre)) * ({$this->lat} - Y(centre))) AS distance ".
|
169 |
jpm |
97 |
'FROM osm_communes '.
|
1056 |
jpm |
98 |
'ORDER BY distance ASC '.
|
|
|
99 |
'LIMIT 20 '.
|
|
|
100 |
' -- '.__FILE__.' : '.__LINE__;
|
169 |
jpm |
101 |
$resultat = $this->bdd->recupererTous($requete);
|
|
|
102 |
|
|
|
103 |
if ($resultat == false) {
|
1056 |
jpm |
104 |
$msgTpl = "Service '%s' : aucune commune dont le centroïde correspond aux coordonnées : %s, %s.";
|
|
|
105 |
$msg = sprintf($msgTpl, get_class($this), $this->lat, $this->lon);
|
|
|
106 |
throw new Exception($msg, RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE);
|
169 |
jpm |
107 |
}
|
|
|
108 |
return $resultat;
|
|
|
109 |
}
|
|
|
110 |
|
1056 |
jpm |
111 |
private function formaterResultats($communeTrouveeInfos) {
|
|
|
112 |
return array(
|
|
|
113 |
'nom' => $communeTrouveeInfos['nom'],
|
|
|
114 |
'codeINSEE' => $communeTrouveeInfos['code_insee'],
|
|
|
115 |
'wikipedia' => $this->formaterUrlWikipedia($communeTrouveeInfos['wikipedia'])
|
|
|
116 |
);
|
169 |
jpm |
117 |
}
|
|
|
118 |
|
1056 |
jpm |
119 |
private function localiserPointLatLon($communesProches) {
|
|
|
120 |
$zonesIds = array();
|
|
|
121 |
foreach ($communesProches as $commune) {
|
|
|
122 |
$zonesIds[] = $commune['id_osm'];
|
169 |
jpm |
123 |
}
|
1056 |
jpm |
124 |
$idsOsmClauseIn = implode(',', $this->bdd->proteger($zonesIds));
|
|
|
125 |
$requete = "SELECT id_osm, code_insee, nom, wikipedia ".
|
|
|
126 |
'FROM osm_communes '.
|
|
|
127 |
"WHERE st_within(GEOMFROMTEXT('POINT($this->lon $this->lat)'), polygone) = 1 ".
|
|
|
128 |
"AND id_osm IN ($idsOsmClauseIn) ".
|
|
|
129 |
' -- '.__FILE__.' : '.__LINE__;
|
|
|
130 |
$resultat = $this->bdd->recuperer($requete);
|
|
|
131 |
if ($resultat === false) {
|
|
|
132 |
$msgTpl = "Service '%s' : aucune commune correspondant aux coordonnées : %s, %s.";
|
|
|
133 |
$msg = sprintf($msgTpl, get_class($this), $this->lat, $this->lon);
|
|
|
134 |
throw new Exception($msg, RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE);
|
|
|
135 |
}
|
|
|
136 |
return $resultat;
|
169 |
jpm |
137 |
}
|
|
|
138 |
|
1056 |
jpm |
139 |
private function formaterUrlWikipedia($infoWikipedia) {
|
|
|
140 |
list($lang, $page) = explode(':', $infoWikipedia);
|
|
|
141 |
$pageEncode = rawurlencode($page);
|
|
|
142 |
$urlTpl = 'https://%s.wikipedia.org/wiki/%s';
|
|
|
143 |
return sprintf($urlTpl, $lang, $pageEncode);
|
169 |
jpm |
144 |
}
|
1056 |
jpm |
145 |
}
|