Subversion Repositories eFlore/Applications.cel

Rev

Rev 2557 | Rev 2565 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
2538 aurelien 1
<?php
2
// declare(encoding='UTF-8');
3
/**
4
 * Classe recherchant des infos sur une zone géo ou bien des coordonnées.
5
 *
6
 * @internal   Mininum PHP version : 5.2
7
 * @category   CEL
8
 * @package    Services
9
 * @subpackage Bibliothèques
10
 * @version    0.1
11
 * @author     Mathias CHOUET <mathias@tela-botanica.org>
12
 * @author     Jean-Pascal MILCENT <jpm@tela-botanica.org>
13
 * @author     Aurelien PERONNET <aurelien@tela-botanica.org>
14
 * @license    GPL v3 <http://www.gnu.org/licenses/gpl.txt>
15
 * @license    CECILL v2 <http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt>
16
 * @copyright  1999-2014 Tela Botanica <accueil@tela-botanica.org>
17
 */
18
class RechercheInfosZoneGeo extends Cel {
2557 aurelien 19
 
20
	/**
21
	 * Cas de la demande par lat et lon
22
	 * array(
23
	 * 		'code_pays' => 'FR',
24
	 * 		'code_zone' => 'INSEE-D:35800',
25
	 * 		'type' => 'france',
26
	 * 		'nom' => 'Taulignan-les-oies'
27
	 * 	);
28
	 *
29
	 * 	array(
30
	 * 		'code_pays' => 'RU',
31
	 * 		'code_zone' => 'OSM-ID:4511574',
32
	 * 		'type' => 'monde',
33
	 * 		'nom' => 'Moscou'
34
	 * 	);
35
	 *
36
	 * Cas de la demande par nom
37
	 * array(
38
     * 		'lat': 40.372934,
39
 	 *		'lng': -3.690601,
40
 	 *		'nom': 'Madrid',
41
 	 *		'code_zone': 'OSM-ID:558474',
42
     *		'code_pays': 'ES',
43
     *		'type': 'monde'
44
	 * 	);
45
	 *
46
	 * 	array(
47
     * 		'lat': 48.860799,
48
 	 *		'lng': 2.3457986,
49
 	 *		'nom': 'Paris',
50
 	 *		'code_zone': '75056',
51
     *		'code_pays': 'FR',
52
     *		'type': 'france'
53
	 * 	);
54
	 *
55
	 */
2538 aurelien 56
 
57
	//TODO: migrer tout ce qui concerne la zone géo dans cette classe
58
	public function obtenirInfosPourCoordonnees($coordonnees) {
59
		$infos_coord = null;
60
		// Test facile qui permet d'abord de tenter une localisation bien plus rapide si les coordonnées
61
		// sont dans la bounding box approximative de la France
62
		if($this->testerCoordonneesWgsFrance($coordonnees['latitude'], $coordonnees['longitude'])) {
63
			$infos_coord = $this->chercherInfosCommune('osm', $coordonnees['latitude'], $coordonnees['longitude']);
64
			if ($infos_coord == null) {
65
				// Sinon recherche par pays
66
				$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);
67
			}
68
		} else {
69
			// Recherche par pays immédiate si en dehors de la bouding box française
70
			$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);
71
		}
72
		return $infos_coord;
2557 aurelien 73
	}
74
 
75
	public function obtenirInfosPourNoms($nom, $pays, $code) {
2564 aurelien 76
 
77
		$nom = trim($nom);
78
		$pays = trim($pays);
79
		$code = trim($code);
80
 
2557 aurelien 81
		$retour = false;
82
		// Dans le cas de la france on possède directement une table
83
		// avec les centroïdes des communes
2564 aurelien 84
		if($pays == 'FR' || $code != "") {
2557 aurelien 85
			$retour = $this->chercherCentroideCommuneBdd($nom, $code);
86
		} else {
87
			// Sinon méthode habituelle, un test simple suivi d'un test plus complexe
2564 aurelien 88
			// On cherche donc d'abord en France
89
			if($pays == "") {
90
				$retour = $this->chercherCentroideCommuneBdd($nom, $code);
2557 aurelien 91
			}
2564 aurelien 92
 
93
			// Si on ne trouve rien on teste dans le monde entier
94
			if($retour == false) {
95
				$retour = $this->effectuerRequeteGeocodingMondiale($nom, '2,3,4,5,6,7,8', 1, $pays);
96
			}
2557 aurelien 97
		}
98
 
99
		return $retour;
2538 aurelien 100
	}
101
 
102
	public function chercherInfosCommune($projet, $latitude, $longitude) {
103
 
104
		$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;
105
		$url_service = $this->config['cel']['url_service_geo_local'].$valeurs_params;
106
		$url_service = str_replace(',', '.', $url_service);
107
 
108
		$ch = curl_init($url_service);
109
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
110
		$reponse = curl_exec($ch);
111
 
112
		$infos_coord = $this->traiterReponseServiceCommune($reponse);
113
		curl_close($ch);
114
 
115
		return $infos_coord;
116
	}
117
 
118
	public function chercherInfosPays($latitude, $longitude) {
119
 
120
		$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;
121
		$url_service = $this->config['cel']['url_service_geo_mondial'].$valeurs_params;
122
		$url_service = str_replace(',', '.', $url_service);
123
 
124
		$ch = curl_init($url_service);
125
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
126
		$reponse = curl_exec($ch);
2557 aurelien 127
 
2538 aurelien 128
		$res = json_decode($reponse, true);
2557 aurelien 129
		$retour = null;
130
 
131
		$retour = array(
132
					'code_zone' => '',
133
					'code_pays' => '',
134
					'nom' => '',
135
					'type' => ''
136
				);
2538 aurelien 137
 
138
		if(!empty($res)) {
139
			// Cas où l'on a recherché des coordonnées en France
140
			// On en profite pour remplir plus de données
141
			// que le simple pays si elles sont présentes
142
			if(isset($res[8]) && isset($res[8]['codeInsee'])) {
143
				$retour = array(
2557 aurelien 144
						'code_zone' => $res[8]['codeInsee'],
145
						'code_pays' => 'FR',
146
						'nom' => $res[8]['nom'],
2538 aurelien 147
						'type' => 'france'
148
				);
149
			} else {
150
				// Cas de la recherche en dehors de France
151
				// La zone de plus haut niveau est toujours un pays
152
				// (car le niveau de zone est limité à 2)
153
				$infos_pays = $res[min(array_keys($res))];
154
				// La zone de niveau le plus bas est la "localité"
155
				// la plus précise qu'on a pu trouver
156
				$infos_localite = $res[max(array_keys($res))];
157
 
158
				// Cas où l'on a trouvé un code pays
159
				if(!empty($infos_pays['codeIso31661'])) {
160
					$retour = array(
161
							'code_pays' => $infos_pays['codeIso31661'],
162
							'nom' => ''
163
					);
164
				} elseif(!empty($infos_pays['codeIso31662'])) {
165
					// Quelquefois la zone de plus haut niveau est une région ou une province
166
					// heureusement son code est de forme XX-YY ou XX est le code iso du pays !
167
					$retour = array(
168
							'code_pays' => substr($infos_pays['codeIso31662'], 0, 2),
169
							'nom' => ''
170
					);
171
				}
172
 
173
				// Pas de sens de prendre un nom de localité si on a pas de pays
174
				if(!empty($retour['code_pays'])) {
175
					// Type sert à savoir quelle fonction de mise à jour on appellera
2557 aurelien 176
					$retour['type'] = 'monde';
177
					$retour['nom'] = $infos_localite['intitule'];
178
					$retour['code_zone'] = $infos_localite['codeZoneGeo'];
2538 aurelien 179
				}
180
			}
2557 aurelien 181
		}
182
 
2538 aurelien 183
		return $retour;
184
	}
185
 
186
	public function traiterReponseServiceCommune($reponse) {
187
		$retour = null;
188
		$reponse = json_decode($reponse);
189
		// cas du service lion 1906 qui peut renvoyer plusieurs communes (on prend la première)
190
		if (is_array($reponse)) {
191
			$reponse = $reponse[0];
192
		}
193
		if (isset($reponse->codeINSEE)) {
194
			// Type sert à savoir quelle fonction de mise à jour on appellera
2557 aurelien 195
			$retour = array(
196
						'code_zone' => $reponse->codeINSEE,
197
						'code_pays' => 'FR',
198
						'nom' => $reponse->nom,
199
						'type' => 'france'
200
					);
2538 aurelien 201
		}
202
		return $retour;
203
	}
204
 
2564 aurelien 205
	private function effectuerRequeteGeocodingMondiale($nom, $niveau = '2,3,4,5,6,7,8', $limite = 1, $pays = null) {
206
		$url_sans_pays = $this->config['cel']['url_service_geo_mondial'].'?masque='.urlencode($nom).'&niveau='.$niveau.'&limite='.$limite;
207
		$url = $url_sans_pays.(!empty($pays) ? '&pays='.urlencode($pays) : '');
208
 
2557 aurelien 209
		$res = @json_decode(file_get_contents($url), true);
210
 
211
		$retour = array();
212
 
213
		if(!empty($res)) {
214
			// cas de la localisation directe d'une commune
215
			if($limite == 1) {
216
				$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);
217
			} else {
218
				// cas de l'autocompletion
219
				$retour = array();
220
				foreach($res as $ligne) {
221
					$ligne_decodee = $this->traiterLigneResultatRequeteGeocodingMondiale($ligne);
222
					$param = array('nom' => $ligne_decodee['nom'], 'code' => $ligne_decodee['code_pays']);
223
					$retour[] = $this->formaterLigneResultat($param);
224
				}
225
			}
2564 aurelien 226
		} else if($limite == 1) {
227
			// Dans le cas où un pays a été demandé mais pas de résultat, on rentente la recherche sans
228
			// préciser de pays
229
			$res = @json_decode(file_get_contents($url_sans_pays), true);
230
 
231
			if(!empty($res)) {
232
				$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);
233
			}
2557 aurelien 234
		}
235
 
236
		return $retour;
237
	}
238
 
239
	private function chercherCentroideCommuneBdd($commune, $departement) {
240
		$commune_formatee = str_replace(array(' ', '-'), '_', trim($commune));
241
 
242
		$departement = trim($departement, "*");
243
		if (strlen($departement) > 2) {
244
			$departement = substr($departement, 0, 2);
245
		}
246
		$requete = 'SELECT utm_x, utm_y, utm_secteur, code FROM cel_zones_geo '.
247
				'WHERE nom LIKE '.Cel::db()->proteger($commune_formatee).' '.
248
				'AND code LIKE '.Cel::db()->proteger($departement.'%').' '.
249
				' -- '.__FILE__.':'.__LINE__;
250
 
251
		$commune_coordonnees = Cel::db()->requeter($requete);
252
 
253
		$retour = false;
254
		if ($commune_coordonnees && is_array($commune_coordonnees) && count($commune_coordonnees) > 0) {
255
			$lat_lon = $this->convertirUtmVersLatLong($commune_coordonnees[0]['utm_x'],$commune_coordonnees[0]['utm_y'],$commune_coordonnees[0]['utm_secteur']);
256
 
257
			$retour = array(
258
					'lat' => (float) $lat_lon['lat'],
259
					'lng' => (float) $lat_lon['lng'],
260
					'nom' => $commune,
261
					'code_zone' => $commune_coordonnees[0]['code'],
262
					'code_pays' => 'FR',
263
					'type' => 'france'
264
			);
265
		}
266
		return $retour;
267
	}
268
 
269
	public function obtenirListeInfosPourNom($lieu) {
270
		$lieux_retour = $this->obtenirListeCommunes($lieu);
271
 
272
		// S'il y a moins de 50 résultats on tente des les fusionner avec une recherche mondiale
273
		// avec 50 - n éléments (où n est la taille du tableau de communes)
274
		if(count($lieux_retour) < 50) {
275
			$lieux_mondiaux = $this->effectuerRequeteGeocodingMondiale($lieu, '2,3,4,5,6,7,8', 50 - count($lieux_retour));
276
			usort($lieux_mondiaux, array($this, 'trierListeLieux'));
277
			$lieux_retour += $lieux_mondiaux;
278
		}
279
 
280
		// Suppression de possibles doublons
281
		$lieux_retour = $this->reduireListeLieux($lieux_retour);
282
 
283
		return $lieux_retour;
284
	}
285
 
286
	public function obtenirListeCommunes($lieu) {
287
		$lieu_formate = $this->formaterChaineLieuPourRequete($lieu);
288
		$retour = array();
289
 
290
		if ($this->estUneChaineRequeteValide($lieu_formate)) {
291
			$requete = 'SELECT DISTINCT nom, code '.
292
					'FROM cel_zones_geo '.
293
					'WHERE nom LIKE '.Cel::db()->proteger($lieu_formate.'%').' '.
294
					'ORDER BY nom '.
295
					'LIMIT 50 '.
296
					' -- '.__FILE__.':'.__LINE__;
297
 
298
			$liste_lieux = Cel::db()->requeter($requete);
299
			if ($liste_lieux) {
300
				foreach ($liste_lieux as $lieu_trouve) {
301
					$retour[] = $this->formaterLigneResultat($lieu_trouve);
302
				}
303
			}
304
		}
305
		return $retour;
306
	}
307
 
308
	/**
309
	 * Fonctions utilitaires
310
	 *
311
	*/
312
 
313
	private function reduireListeLieux($tableau) {
314
		$index = array();
315
		foreach($tableau as $ligne) {
316
			$index_loc = trim(preg_replace("/\([a-zA-Z0-9]{2}\)$/", "", $ligne[0]));
317
			if(!isset($index[$index_loc])) {
318
				$index[$index_loc] = $ligne;
319
			} else {
320
				if(mb_strlen($ligne[0]) > mb_strlen($index[$index_loc][0])) {
321
					$index[$index_loc] = $ligne;
322
				} else if(mb_strlen($ligne[0]) == mb_strlen($index[$index_loc][0])) {
323
					if(preg_match("/(.)* ([a-zA-Z0-9]{2})/", $ligne[0])) {
324
						$index[$index_loc] = $ligne;
325
					}
326
				}
327
			}
328
		}
329
 
330
		return array_values($index);
331
	}
332
 
333
	private function trierListeLieux($a, $b) {
334
		return strcmp($a[0], $b[0]);
335
	}
336
 
2538 aurelien 337
	private function testerCoordonneesWgsFrance($latitude, $longitude) {
338
		$coord_france = false;
339
		if ($latitude != '' && $longitude != '') {
340
			if ($latitude < 51.071667 && $latitude > 41.316667) {
341
				if ($longitude < 9.513333 && $longitude > -5.140278) {
342
					$coord_france = true;
343
				}
344
			}
345
		}
346
		return $coord_france;
347
	}
2557 aurelien 348
 
349
	private function traiterLigneResultatRequeteGeocodingMondiale($res) {
350
 
351
		$lat = "";
352
		$lon = "";
353
		$nom = "";
354
		$code ="";
355
		$code_pays="";
356
		$type="";
357
 
358
		if(!empty($res)) {
359
			$code_pays = $this->extraireCodePays($res);
360
 
361
			$nom = $res['intitule'];
362
			$lat = $res['centre_lat'];
363
			$lon = $res['centre_lng'];
364
			$type = $code_pays == 'FR' || $code_pays == 'FX' ? 'france' : 'monde';
365
 
366
			$retour = array(
367
					'lat' => (float)$lat,
368
					'lng' => (float)$lon,
369
					'nom' => $nom,
370
					'code_zone' => $code,
371
					'code_pays' => $code_pays,
372
					'type' => $type
373
			);
374
		}
375
 
376
		return $retour;
377
	}
378
 
379
	private function convertirUtmVersLatLong($x, $y, $sector) {
380
		$convertisseur = new gPoint();
381
		$convertisseur->setUTM($x, $y, $sector);
382
		$convertisseur->convertTMtoLL();
383
 
384
		$lat_long = array();
385
		$lat_long['lat'] = str_replace(',','.',$convertisseur->Lat());
386
		$lat_long['lng'] = str_replace(',','.',$convertisseur->Long());
387
		return $lat_long;
388
	}
389
 
390
	private function extraireCodePays($infos_pays) {
391
		$code = "";
392
		if(!empty($infos_pays['codeIso31661'])) {
393
			$code = $infos_pays['codeIso31661'];
394
		} elseif(!empty($infos_pays['codeIso31662'])) {
395
			// Si pas de code pays, le sous code peut exister et content le code pays
396
			// sous la forme codepays-souscode
397
			$code = substr($infos_pays['codeIso31662'], 0, 2);
398
		}
399
		return $code;
400
	}
401
 
402
	private function estUneChaineRequeteValide($lieu) {
403
		return (strlen($lieu) > 0) && ($lieu != '%');
404
	}
405
 
406
	private function formaterLigneResultat($ligne) {
407
		if(empty($ligne['code']) || $ligne['code'] == null) {
408
			$res = array($ligne['nom'], '');
409
		} else {
410
			$res = array($ligne['nom'].' ('.substr(sprintf('%02s', $ligne['code']),0,2).')', $ligne['code']);
411
		}
412
		return $res;
413
	}
414
 
415
	/**
416
	 * Remplace les * par % pour faire des recherches floues
417
	 * Remplace les + par _ (nginx envoie des "+" dans l'URL à la place des espaces)
418
	 * Remplace les espaces et les - par _ car les noms de communes peuvent avoir des espaces ou des tirets
419
	 * @param string $lieu
420
	 * @return string le lieu formaté pour la recherche
421
	 */
422
	private function formaterChaineLieuPourRequete($lieu) {
423
		$lieu = ltrim($lieu);
424
		$lieu = preg_replace('/\*+/', '%', $lieu);
425
		$lieu = str_replace(['+', ' ', '-'], '_', $lieu);
426
		return $lieu;
427
	}
2538 aurelien 428
}