Subversion Repositories eFlore/Projets.eflore-projets

Rev

Rev 1048 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 1048 Rev 1063
Line 2... Line 2...
2
// declare(encoding='UTF-8');
2
// declare(encoding='UTF-8');
3
/**
3
/**
4
 * Classe permettant d'obtenir le nom et le code INSEE d'une commune à partir de ces coordonnées (latitude et longitude).
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.
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.
6
 * Ce service fonctionne uniquement pour les communes de France métropolitaine (Corse comprise) dont les contours ont été tracés dans OpenStreetMap.
-
 
7
 * Source des données : OpenStreetMap http://wiki.openstreetmap.org
7
 * Source des données : OpenStreetMap http://wiki.openstreetmap.org (Projet France limites administratives : http://wiki.openstreetmap.org/wiki/WikiProject_France/Limites_administratives)
8
 * Projet France limites administratives : http://wiki.openstreetmap.org/wiki/WikiProject_France/Limites_administratives
8
 * Paramètres du service :
9
 * Paramètres du service :
9
 *  - lat : latitude
10
 *  - lat : latitude
10
 *  - lon : longitude
11
 *  - lon : longitude
11
 * Exemple :
12
 * Exemple :
12
 * http://localhost/cartoOSM/services/0.1/nom-commune?lat=44.71546&lon=3.84216
13
 * http://localhost/service-test:eflore:0.1/osm/nom-commune?lon=3.83619&lat=44.74231
13
 *
14
 *
14
 * @category	php 5.2
15
 * @category	php 5.2
15
 * @package		carto-osm
16
 * @package		carto-osm
16
 * @author		Mohcen BENMOUNAH <mohcen@tela-botanica.org>
17
 * @author		Mohcen BENMOUNAH <mohcen@tela-botanica.org>
17
 * @author		Jean-Pascal MILCENT <jpm@tela-botanica.org>
18
 * @author		Jean-Pascal MILCENT <jpm@tela-botanica.org>
Line 26... Line 27...
26
	const PATTERN_LON = '/^[-]?[0-9]+(?:[.][0-9]+|)$/';
27
	const PATTERN_LON = '/^[-]?[0-9]+(?:[.][0-9]+|)$/';
27
	const LAT_MAX = 51.071667;
28
	const LAT_MAX = 51.071667;
28
	const LAT_MIN = 41.316667;
29
	const LAT_MIN = 41.316667;
29
	const LON_MAX = 9.513333;
30
	const LON_MAX = 9.513333;
30
	const LON_MIN = -5.140278;
31
	const LON_MIN = -5.140278;
31
	const NBRE_COMMUNE_PAR_DEFAUT = 10;
-
 
32
	const NBRE_COMMUNE_MAX = 100;
-
 
33
	const MIME_JSON = 'application/json';
32
	const MIME_JSON = 'application/json';
Line 34... Line 33...
34
 
33
 
35
	private $parametres = array();
34
	private $parametres = array();
Line -... Line 35...
-
 
35
	private $bdd = null;
-
 
36
 
-
 
37
	private $lat = null;
36
	private $bdd = null;
38
	private $lon = null;
37
 
39
 
38
	public function __construct(Bdd $bdd) {
40
	public function __construct(Bdd $bdd) {
Line 39... Line 41...
39
		$this->bdd = $bdd;
41
		$this->bdd = $bdd;
40
	}
42
	}
41
 
43
 
-
 
44
	public function consulter($ressources, $parametres) {
-
 
45
		$this->parametres = $parametres;
Line 42... Line 46...
42
	public function consulter($ressources, $parametres) {
46
		$this->verifierParametres();
-
 
47
		$this->lat = $this->parametres['lat'];
43
		$this->parametres = $parametres;
48
		$this->lon = $this->parametres['lon'];
Line 44... Line 49...
44
		$this->verifierParametres();
49
 
45
 
50
		$communesProches = $this->trouverCommunesProches();
46
		$nomINSEEs = $this->trouverCommunesProches();
51
		$communeTrouveeInfos = $this->localiserPointLatLon($communesProches);
47
		$corps = $this->formaterResultats($nomINSEEs);
52
		$corps = $this->formaterResultats($communeTrouveeInfos);
Line 58... Line 63...
58
		if (! array_key_exists('lat', $this->parametres)) {
63
		if (! array_key_exists('lat', $this->parametres)) {
59
			$messages[] = "Vous devez indiquer une latitude en degré décimal à l'aide du paramètres d'url : lat";
64
			$messages[] = "Vous devez indiquer une latitude en degré décimal à l'aide du paramètres d'url : lat";
60
		} else if (!preg_match(self::PATTERN_LAT, $lat)) {
65
		} else if (!preg_match(self::PATTERN_LAT, $lat)) {
61
			$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";
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";
62
		} else if ($lat > self::LAT_MAX) {
67
		} else if ($lat > self::LAT_MAX) {
63
			$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.";
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.";
64
		} else if ($lat < self::LAT_MIN) {
69
		} else if ($lat < self::LAT_MIN) {
65
			$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.";
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.";
66
		}
71
		}
67
		if (! array_key_exists('lon', $this->parametres)) {
72
		if (! array_key_exists('lon', $this->parametres)) {
68
			$messages[] = "Vous devez indiquer une longitude en degré décimal à l'aide du paramètres d'url : lon";
73
			$messages[] = "Vous devez indiquer une longitude en degré décimal à l'aide du paramètres d'url : lon";
69
		} else if (!preg_match(self::PATTERN_LON, $lon)) {
74
		} else if (!preg_match(self::PATTERN_LON, $lon)) {
70
			$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";
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";
71
		} else if ($lon > self::LON_MAX) {
76
		} else if ($lon > self::LON_MAX) {
72
			$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.";
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.";
73
		} else if ($lon < self::LON_MIN) {
78
		} else if ($lon < self::LON_MIN) {
74
			$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.";
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.";
75
		}
80
		}
76
		if (count($messages) != 0) {
81
		if (count($messages) != 0) {
77
			$message = implode('<br />', $messages);
82
			$message = implode('<br />', $messages);
78
			$code = RestServeur::HTTP_CODE_MAUVAISE_REQUETE;
83
			$code = RestServeur::HTTP_CODE_MAUVAISE_REQUETE;
79
			throw new Exception($message, $code);
84
			throw new Exception($message, $code);
Line 81... Line 86...
81
	}
86
	}
Line 82... Line 87...
82
 
87
 
83
	/**
88
	/**
84
	 * requête qui récupère les 20 communes les plus proches du point recherché
89
	 * requête qui récupère les 20 communes les plus proches du point recherché
85
	 * La distance(AB = \sqrt{(x_B-x_A)^2 + (y_B-y_A)^2}) est calculée sans la racine
90
	 * La distance(AB = \sqrt{(x_B-x_A)^2 + (y_B-y_A)^2}) est calculée sans la racine
-
 
91
	 * (calcul en plus qui change pas le résultat). La distance n'est donc pas juste mais permet de récupérer
86
	 * (calcul en plus qui change pas le résultat).
92
	 * le bon ordre.
87
	*/
93
	*/
88
	private function trouverCommunesProches() {
-
 
89
		$lat = $this->parametres['lat'];
-
 
90
		$lon = $this->parametres['lon'];
94
	private function trouverCommunesProches() {
91
		$requete =	'SELECT  '.
95
		$requete =	'SELECT id_osm, '.
92
			"	(({$lat} - X(centre)) * ({$lat} -X(centre)) + ({$lon} - Y(centre)) * ({$lon} - Y(centre))) AS distance, ".
-
 
93
			"	code_insee, nom, ASTEXT(polygon) AS poly ".
96
			"(({$this->lon} - X(centre)) * ({$this->lon} -X(centre)) + ({$this->lat} - Y(centre)) * ({$this->lat} - Y(centre))) AS distance ".
94
			'FROM osm_communes '.
-
 
95
			'WHERE note = "Polygone complet" '.
97
			'FROM osm_communes '.
96
			'ORDER BY distance '.
98
			'ORDER BY distance ASC '.
-
 
99
			'LIMIT 20 '.
97
			'LIMIT 20';
100
			' -- '.__FILE__.' : '.__LINE__;
Line 98... Line 101...
98
		$resultat = $this->bdd->recupererTous($requete);
101
		$resultat = $this->bdd->recupererTous($requete);
99
 
102
 
100
		if ($resultat == false) {
103
		if ($resultat == false) {
101
			$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']}.";
104
			$msgTpl = "Service '%s' : aucune commune dont le centroïde correspond aux coordonnées : %s, %s.";
102
			$code = RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE;
105
			$msg = sprintf($msgTpl, get_class($this), $this->lat, $this->lon);
103
			throw new Exception($message, $code);
106
			throw new Exception($msg, RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE);
104
		}
107
		}
Line 105... Line 108...
105
		return $resultat;
108
		return $resultat;
106
	}
109
	}
107
 
110
 
108
	private function formaterResultats($nomINSEEs) {
111
	private function formaterResultats($communeTrouveeInfos) {
109
		$communeTrouvee = null;
-
 
110
		foreach ($nomINSEEs as $nomINSEE) {
112
		return array(
111
			$resultat = $this->localiserPointLatLon($nomINSEE['poly']);
113
			'nom' => $communeTrouveeInfos['nom'],
112
			if ($resultat == 1) {
114
			'codeINSEE' => $communeTrouveeInfos['code_insee'],
-
 
115
			'wikipedia' => $this->formaterUrlWikipedia($communeTrouveeInfos['wikipedia'])
-
 
116
		);
-
 
117
	}
-
 
118
 
-
 
119
	private function localiserPointLatLon($communesProches) {
113
				$communeTrouvee = array('nom' => $nomINSEE['nom'], 'codeINSEE' => $nomINSEE['code_insee']);
120
		$zonesIds = array();
-
 
121
		foreach ($communesProches as $commune) {
-
 
122
			$zonesIds[] = $commune['id_osm'];
-
 
123
		}
-
 
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 ".
114
				break;
128
			"AND id_osm IN ($idsOsmClauseIn) ".
115
			}
129
			' -- '.__FILE__.' : '.__LINE__;
116
		}
130
		$resultat = $this->bdd->recuperer($requete);
117
		if (is_null($communeTrouvee)) {
131
		if ($resultat === false) {
118
			$message = "Le service '".get_class($this)."' n'a trouvé aucune commune correspondant aux coordonnées : {$this->parametres['lat']}, {$this->parametres['lon']}.";
132
			$msgTpl = "Service '%s' : aucune commune correspondant aux coordonnées : %s, %s.";
119
			$code = RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE;
133
			$msg = sprintf($msgTpl, get_class($this), $this->lat, $this->lon);
120
			throw new Exception($message, $code);
134
			throw new Exception($msg, RestServeur::HTTP_CODE_RESSOURCE_INTROUVABLE);
Line 121... Line 135...
121
		}
135
		}
122
		return $communeTrouvee;
-
 
123
	}
-
 
124
 
-
 
125
	private function localiserPointLatLon($multiPolygone) {
-
 
126
		$lat = $this->parametres['lat'];
-
 
127
		$lon = $this->parametres['lon'];
-
 
128
		$requeteMBR = "SELECT MBRWithin(GEOMFROMTEXT('POINT($lat $lon)'), GEOMFROMTEXT('$multiPolygone'))";
-
 
129
		$p = array_values($this->bdd->recuperer($requeteMBR));
-
 
130
 
-
 
131
		$counter = 0;
-
 
132
		if ($p[0] != 1) {
-
 
133
			 return 0;
-
 
134
		} else {
-
 
135
			$x = $lat;	  //Le point latitude recherché
-
 
136
			$y = $lon;		//Le point longitude recherché
-
 
137
			$str = $this->nettoyerChaineMultipolygone($multiPolygone); // Le multipolygone de la commune
-
 
138
			$pb = 0;
-
 
139
			$pe =  strpos($str, ','); //définition de la position de la première virgule
-
 
140
			$xy =  substr($str, $pb, $pe - $pb); //extraire X Y d'un point séparé par des virgules
-
 
141
			$p = strpos($xy, ' ');  //position espace dans un point
-
 
142
			$p1x = substr($xy, 0, $p ); //récupère X du point
-
 
143
			$p1y = substr($xy, $p + 1 ); //récupère Y du point
-
 
144
			$str = $str.$xy; // ajouter le premier point du polygone à la fin du polygone
-
 
145
			while ($pe > 0) {
-
 
146
				$xy = substr($str, $pb, $pe - $pb ); //extraire le point GPS
-
 
147
				$p	= strpos($xy, ' ');   // trouver la position de l'espace entre les deux points X et Y
-
 
148
				$p2x = substr($xy,0, $p );  // extraire le point X
-
 
149
				$p2y = substr($xy, $p + 1 );  // extraire le point Y
-
 
150
				if ($p1y < $p2y) {
-
 
151
					$m = $p1y;
-
 
152
				} else {
-
 
153
					$m = $p2y;
-
 
154
				}
-
 
155
				if ($y > $m) {
-
 
156
					if ($p1y > $p2y) {
-
 
157
						$m = $p1y;
-
 
158
					} else {
-
 
159
						$m = $p2y;
-
 
160
					}
-
 
161
					if ($y <= $m) {
-
 
162
						if ($p1x > $p2x) {
-
 
163
							$m = $p1x;
-
 
164
						} else {
-
 
165
							$m = $p2x;
-
 
166
						}
-
 
167
						if ($y <= $m) {
-
 
168
							if ($p1x > $p2x) {
-
 
169
								$m = $p1x;
-
 
170
							} else {
-
 
171
								$m = $p2x;
-
 
172
							}
-
 
173
							if ($x <= $m) {
136
		return $resultat;
174
								if ($p1y != $p2y) {
-
 
175
									$xinters = ($y - $p1y) * ($p2x - $p1x) / ($p2y - $p1y) + $p1x;
-
 
176
								}
-
 
177
								if (($p1x = $p2x) || ($x <= $xinters )) {
-
 
178
									$counter++;
-
 
179
								}
-
 
180
							}
-
 
181
						}
-
 
182
					}
-
 
183
				}
137
	}
184
				$p1x = $p2x;
-
 
185
				$p1y = $p2y;
-
 
186
				$pb = $pe + 1;
-
 
187
				$pe = strpos($str, ',', $pb);
-
 
188
			}
-
 
189
			return $counter % 2;
-
 
190
		}
138
 
191
	}
139
	private function formaterUrlWikipedia($infoWikipedia) {
192
 
-
 
193
	private function nettoyerChaineMultipolygone($chaine) {
-
 
194
		$chaine = str_replace('MULTIPOLYGON(((', '', $chaine);
140
		list($lang, $page) = explode(':', $infoWikipedia);
195
		$chaine = str_replace(')))', '', $chaine);
-
 
196
		$chaine .= ',';
-
 
197
		return $chaine;
141
		$pageEncode = rawurlencode($page);
-
 
142
		$urlTpl = 'https://%s.wikipedia.org/wiki/%s';
198
	}
143
		return sprintf($urlTpl, $lang, $pageEncode);