Subversion Repositories eFlore/Applications.cel

Rev

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

Rev Author Line No. Line
2538 aurelien 1
<?php
2914 mathias 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
 */
2538 aurelien 18
class RechercheInfosZoneGeo extends Cel {
2557 aurelien 19
 
2565 aurelien 20
	public $maxTimeoutRequeteMondiale = 0;
21
 
2557 aurelien 22
	/**
23
	 * Cas de la demande par lat et lon
24
	 * array(
25
	 * 		'code_pays' => 'FR',
26
	 * 		'code_zone' => 'INSEE-D:35800',
27
	 * 		'type' => 'france',
28
	 * 		'nom' => 'Taulignan-les-oies'
29
	 * 	);
30
	 *
31
	 * 	array(
32
	 * 		'code_pays' => 'RU',
33
	 * 		'code_zone' => 'OSM-ID:4511574',
34
	 * 		'type' => 'monde',
35
	 * 		'nom' => 'Moscou'
36
	 * 	);
37
	 *
38
	 * Cas de la demande par nom
39
	 * array(
40
     * 		'lat': 40.372934,
41
 	 *		'lng': -3.690601,
42
 	 *		'nom': 'Madrid',
43
 	 *		'code_zone': 'OSM-ID:558474',
44
     *		'code_pays': 'ES',
45
     *		'type': 'monde'
46
	 * 	);
47
	 *
48
	 * 	array(
49
     * 		'lat': 48.860799,
50
 	 *		'lng': 2.3457986,
51
 	 *		'nom': 'Paris',
52
 	 *		'code_zone': '75056',
53
     *		'code_pays': 'FR',
54
     *		'type': 'france'
55
	 * 	);
56
	 *
57
	 */
2538 aurelien 58
 
59
	//TODO: migrer tout ce qui concerne la zone géo dans cette classe
2914 mathias 60
	public function obtenirInfosPourCoordonnees($coordonnees) {
61
		$infos_coord = null;
62
		// Test facile qui permet d'abord de tenter une localisation bien plus rapide si les coordonnées
63
		// sont dans la bounding box approximative de la France
64
		if($this->testerCoordonneesWgsFrance($coordonnees['latitude'], $coordonnees['longitude'])) {
65
			$infos_coord = $this->chercherInfosCommune('osm', $coordonnees['latitude'], $coordonnees['longitude']);
66
			if ($infos_coord == null) {
67
				// Sinon recherche par pays
68
				$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);
69
			}
70
		} else {
71
			// Recherche par pays immédiate si en dehors de la bouding box française
72
			$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);
73
		}
74
		return $infos_coord;
2557 aurelien 75
	}
2914 mathias 76
 
77
	// gère aussi les id_zone_geo @TODO renommer la méthode et internationaliser
78
	// la gestion des codes
79
	public function obtenirInfosPourNom($nom, $pays, $code_departement=null, $id_zone=null) {
2564 aurelien 80
		$nom = trim($nom);
81
		$pays = trim($pays);
2914 mathias 82
		$code_departement = trim($code_departement);
83
 
2557 aurelien 84
		$retour = false;
85
		// Dans le cas de la france on possède directement une table
86
		// avec les centroïdes des communes
2914 mathias 87
		if($pays == 'FR' || $code_departement != "") {
88
			$retour = $this->chercherCentroideCommuneBdd($nom, $code_departement, $id_zone);
2557 aurelien 89
		} else {
90
			// Sinon méthode habituelle, un test simple suivi d'un test plus complexe
2564 aurelien 91
			// On cherche donc d'abord en France
2914 mathias 92
			// @TODO c'est quoi ce if() mal ordonné ?
2564 aurelien 93
			if($pays == "") {
2914 mathias 94
				$retour = $this->chercherCentroideCommuneBdd($nom, $code_departement, $id_zone);
2557 aurelien 95
			}
2564 aurelien 96
			// Si on ne trouve rien on teste dans le monde entier
2914 mathias 97
			if($retour == false) {
98
				$retour = $this->effectuerRequeteGeocodingMondiale($nom, '2,3,4,5,6,7,8', 1, $pays, $id_zone);
2564 aurelien 99
			}
2557 aurelien 100
		}
101
 
102
		return $retour;
2914 mathias 103
	}
104
 
105
	/**
106
	 * Renvoie tous les codes de zones pour le groupe $groupe
107
	 */
108
	protected function obtenirZonesDuGroupe($groupe) {
109
		$requete = 'SELECT valeur as id_zone FROM cel_groupes_zones '
110
			. 'WHERE id_groupe = ' . Cel::db()->proteger($groupe)
111
			. ' -- '.__FILE__.':'.__LINE__;
112
 
113
		$liste_zones = Cel::db()->requeter($requete);
114
		return $liste_zones;
115
	}
116
 
117
	/**
118
	 * Renvoie des infos pour un groupes de zones, déclaré dans la table
119
	 * cel_groupe_zones : toutes les infos des zones concernées, plus le carré
120
	 * englobant (bounding box) et le centroïde de l'union des zones
121
	 */
122
	public function obtenirInfosPourGroupeZonesFrance($groupe) {
123
		$groupe = trim($groupe);
124
		$retour = false;
125
 
126
		$zones = $this->obtenirZonesDuGroupe($groupe);
127
		$retour = array(
128
			'groupe' => array(),
129
			'zones' => array()
130
		);
131
		foreach ($zones as $zone) {
132
			$idZone = $zone['id_zone'];
133
			$retour['zones'][] = $this->obtenirInfosPourNom(null, null, null, $idZone);
134
		}
135
 
136
		// recherche bbox (de Bouygues Telecom) @TODO faire ça au dessus, en une passe ?
137
		$lngmin = 99999;
138
		$latmin = 99999;
139
		$lngmax = -99999;
140
		$latmax = -99999;
141
		foreach ($retour['zones'] as $zone) {
142
			$lat = $zone['lat'];
143
			$lng = $zone['lng'];
144
			if ($lat < $latmin) { $latmin = $lat; }
145
			if ($lat > $latmax) { $latmax = $lat; }
146
			if ($lng < $lngmin) { $lngmin = $lng; }
147
			if ($lng > $lngmax) { $lngmax = $lng; }
148
		}
149
		$retour['groupe']['bbox'] = array(
150
			'lat_min' => $latmin,
151
			'lat_max' => $latmax,
2915 mathias 152
			'lng_min' => $lngmin,
153
			'lng_max' => $lngmax
2914 mathias 154
		);
2915 mathias 155
		// @TODO suggérer un niveau de zoom en fonction de la plus grande dimension
156
		// de la bounding-box (ne marchera pas dans tous les cas mais ce serait
157
		// déjà pas mal)
2914 mathias 158
		$retour['groupe']['centroide'] = array(
159
			'lat' => ($latmin + $latmax) / 2,
2915 mathias 160
			'lng' => ($lngmin + $lngmax) / 2
2914 mathias 161
		);
162
 
163
		return $retour;
164
	}
165
 
166
	public function chercherInfosCommune($projet, $latitude, $longitude) {
167
 
168
		$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;
169
		$url_service = $this->config['cel']['url_service_geo_local'].$valeurs_params;
170
		$url_service = str_replace(',', '.', $url_service);
171
 
172
		$ch = curl_init($url_service);
173
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
174
		$reponse = curl_exec($ch);
175
 
176
		$infos_coord = $this->traiterReponseServiceCommune($reponse);
177
		curl_close($ch);
178
 
179
		return $infos_coord;
180
	}
181
 
2538 aurelien 182
	public function chercherInfosPays($latitude, $longitude) {
2914 mathias 183
 
2538 aurelien 184
		$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;
185
		$url_service = $this->config['cel']['url_service_geo_mondial'].$valeurs_params;
2914 mathias 186
		$url_service = str_replace(',', '.', $url_service);
187
 
188
		$ch = curl_init($url_service);
189
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
190
		$reponse = curl_exec($ch);
191
 
192
		$res = json_decode($reponse, true);
2557 aurelien 193
		$retour = null;
194
 
2914 mathias 195
		$retour = array(
196
					'code_zone' => '',
197
					'code_pays' => '',
198
					'nom' => '',
2557 aurelien 199
					'type' => ''
200
				);
2914 mathias 201
 
202
		if(!empty($res)) {
203
			// Cas où l'on a recherché des coordonnées en France
204
			// On en profite pour remplir plus de données
205
			// que le simple pays si elles sont présentes
206
			if(isset($res[8]) && isset($res[8]['codeInsee'])) {
207
				$retour = array(
2557 aurelien 208
						'code_zone' => $res[8]['codeInsee'],
2914 mathias 209
						'code_pays' => 'FR',
210
						'nom' => $res[8]['nom'],
211
						'type' => 'france'
212
				);
213
			} else {
214
				// Cas de la recherche en dehors de France
215
				// La zone de plus haut niveau est toujours un pays
216
				// (car le niveau de zone est limité à 2)
217
				$infos_pays = $res[min(array_keys($res))];
218
				// La zone de niveau le plus bas est la "localité"
219
				// la plus précise qu'on a pu trouver
220
				$infos_localite = $res[max(array_keys($res))];
221
 
222
				// Cas où l'on a trouvé un code pays
223
				if(!empty($infos_pays['codeIso31661'])) {
224
					$retour = array(
225
							'code_pays' => $infos_pays['codeIso31661'],
226
							'nom' => ''
227
					);
228
				} elseif(!empty($infos_pays['codeIso31662'])) {
229
					// Quelquefois la zone de plus haut niveau est une région ou une province
230
					// heureusement son code est de forme XX-YY ou XX est le code iso du pays !
231
					$retour = array(
232
							'code_pays' => substr($infos_pays['codeIso31662'], 0, 2),
233
							'nom' => ''
234
					);
235
				}
236
 
237
				// Pas de sens de prendre un nom de localité si on a pas de pays
238
				if(!empty($retour['code_pays'])) {
239
					// Type sert à savoir quelle fonction de mise à jour on appellera
240
					$retour['type'] = 'monde';
2557 aurelien 241
					$retour['nom'] = $infos_localite['intitule'];
2914 mathias 242
					$retour['code_zone'] = $infos_localite['codeZoneGeo'];
243
				}
244
			}
2557 aurelien 245
		}
2914 mathias 246
 
247
		return $retour;
248
	}
249
 
250
	public function traiterReponseServiceCommune($reponse) {
251
		$retour = null;
252
		$reponse = json_decode($reponse);
253
		// cas du service lion 1906 qui peut renvoyer plusieurs communes (on prend la première)
254
		if (is_array($reponse)) {
255
			$reponse = $reponse[0];
256
		}
257
		if (isset($reponse->codeINSEE)) {
258
			// Type sert à savoir quelle fonction de mise à jour on appellera
2557 aurelien 259
			$retour = array(
260
						'code_zone' => $reponse->codeINSEE,
2914 mathias 261
						'code_pays' => 'FR',
262
						'nom' => $reponse->nom,
2557 aurelien 263
						'type' => 'france'
2914 mathias 264
					);
265
		}
266
		return $retour;
2538 aurelien 267
	}
2914 mathias 268
 
269
	// @TODO gérer $id_zone (inactif pour l'instant)
270
	protected function effectuerRequeteGeocodingMondiale($nom, $niveau = '2,3,4,5,6,7,8', $limite = 1, $pays = null, $id_zone=null) {
271
		$url_sans_pays = $this->config['cel']['url_service_geo_mondial'].'?masque='.urlencode($nom).'&niveau='.$niveau.'&limite='.$limite;
2564 aurelien 272
		$url = $url_sans_pays.(!empty($pays) ? '&pays='.urlencode($pays) : '');
273
 
2914 mathias 274
		$ch = curl_init($url);
2565 aurelien 275
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
2914 mathias 276
		curl_setopt($ch, CURLOPT_TIMEOUT, $this->maxTimeoutRequeteMondiale);
2565 aurelien 277
		$reponse = curl_exec($ch);
278
 
279
		$res = @json_decode($reponse, true);
2557 aurelien 280
 
2914 mathias 281
		$retour = array();
2557 aurelien 282
 
283
		if(!empty($res)) {
2914 mathias 284
			// cas de la localisation directe d'une commune
2557 aurelien 285
			if($limite == 1) {
286
				$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);
287
			} else {
288
				// cas de l'autocompletion
289
				$retour = array();
290
				foreach($res as $ligne) {
291
					$ligne_decodee = $this->traiterLigneResultatRequeteGeocodingMondiale($ligne);
292
					$param = array('nom' => $ligne_decodee['nom'], 'code' => $ligne_decodee['code_pays']);
293
					$retour[] = $this->formaterLigneResultat($param);
294
				}
295
			}
2564 aurelien 296
		} else if($limite == 1) {
297
			// Dans le cas où un pays a été demandé mais pas de résultat, on rentente la recherche sans
298
			// préciser de pays
299
			$res = @json_decode(file_get_contents($url_sans_pays), true);
300
 
2914 mathias 301
			if(!empty($res)) {
2564 aurelien 302
				$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);
303
			}
2914 mathias 304
		}
305
 
306
		return $retour;
2557 aurelien 307
	}
2914 mathias 308
 
309
	/**
310
	 * Renvoie les infos d'une commune, incluant les coordonnées du centroïde,
311
	 * en fonction du nom et/ou du département, et/ou de l'id_zone (plus fiable)
312
	 */
313
	protected function chercherCentroideCommuneBdd($commune, $departement=null, $id_zone=null) {
2557 aurelien 314
		$commune_formatee = str_replace(array(' ', '-'), '_', trim($commune));
2914 mathias 315
		// clauses
316
		$clauses = array();
317
		if ($commune != null) {
318
			$clauses[] = 'nom LIKE '.Cel::db()->proteger($commune_formatee);
319
		}
320
		if ($departement != null) {
321
			$departement = trim($departement, "*");
322
			if (strlen($departement) > 2) {
323
				// @TODO piège : si deux communes ont le même nom dans deux
324
				// DROM-COM, la troncature à 2 chiffres sera insuffisante
325
				$departement = substr($departement, 0, 2);
326
			}
327
			$clauses[] = 'code LIKE '.Cel::db()->proteger($departement.'%');
328
		}
329
		if ($id_zone != null) {
330
			$clauses[] = 'id_zone_geo = '.Cel::db()->proteger($id_zone);
331
		}
332
		$clauses = implode(' AND ', $clauses);
333
 
334
		$requete = 'SELECT utm_x, utm_y, utm_secteur, code FROM cel_zones_geo '
335
				. 'WHERE ' . $clauses
336
				. ' -- '.__FILE__.':'.__LINE__;
337
 
338
		$commune_coordonnees = Cel::db()->requeter($requete);
339
 
340
		$retour = false;
341
		if ($commune_coordonnees && is_array($commune_coordonnees) && count($commune_coordonnees) > 0) {
342
			$lat_lon = $this->convertirUtmVersLatLong($commune_coordonnees[0]['utm_x'],$commune_coordonnees[0]['utm_y'],$commune_coordonnees[0]['utm_secteur']);
343
 
344
			$retour = array(
345
					'lat' => (float) $lat_lon['lat'],
346
					'lng' => (float) $lat_lon['lng'],
347
					'nom' => $commune,
348
					'code_zone' => $commune_coordonnees[0]['code'],
2557 aurelien 349
					'code_pays' => 'FR',
2914 mathias 350
					'type' => 'france'
351
			);
352
		}
353
		return $retour;
2557 aurelien 354
	}
355
 
356
	public function obtenirListeInfosPourNom($lieu) {
357
		$lieux_retour = $this->obtenirListeCommunes($lieu);
358
 
359
		// S'il y a moins de 50 résultats on tente des les fusionner avec une recherche mondiale
360
		// avec 50 - n éléments (où n est la taille du tableau de communes)
361
		if(count($lieux_retour) < 50) {
362
			$lieux_mondiaux = $this->effectuerRequeteGeocodingMondiale($lieu, '2,3,4,5,6,7,8', 50 - count($lieux_retour));
363
			usort($lieux_mondiaux, array($this, 'trierListeLieux'));
364
			$lieux_retour += $lieux_mondiaux;
365
		}
366
 
2925 mathias 367
		// Suppression de possibles doublons @WARNING marche pas !!!
368
		//$lieux_retour = $this->reduireListeLieux($lieux_retour);
2557 aurelien 369
 
370
		return $lieux_retour;
371
	}
372
 
373
	public function obtenirListeCommunes($lieu) {
2914 mathias 374
		$lieu_formate = $this->formaterChaineLieuPourRequete($lieu);
375
		$retour = array();
376
 
377
		if ($this->estUneChaineRequeteValide($lieu_formate)) {
378
			$requete = 'SELECT DISTINCT nom, code '.
379
					'FROM cel_zones_geo '.
380
					'WHERE nom LIKE '.Cel::db()->proteger($lieu_formate.'%').' '.
381
					'ORDER BY nom '.
382
					'LIMIT 50 '.
383
					' -- '.__FILE__.':'.__LINE__;
2925 mathias 384
 
2914 mathias 385
			$liste_lieux = Cel::db()->requeter($requete);
386
			if ($liste_lieux) {
387
				foreach ($liste_lieux as $lieu_trouve) {
388
					$retour[] = $this->formaterLigneResultat($lieu_trouve);
389
				}
390
			}
391
		}
2557 aurelien 392
		return $retour;
393
	}
394
 
395
	/**
396
	 * Fonctions utilitaires
2925 mathias 397
	 */
398
 
399
	// qui a écrit cette m*rde ?
2914 mathias 400
	protected function reduireListeLieux($tableau) {
2557 aurelien 401
		$index = array();
402
		foreach($tableau as $ligne) {
2925 mathias 403
			// merci de mettre des commentaires !!!
404
			// ça enlève les 2 premières lettres ou chiffres du nom de lieu => wtf ?
2557 aurelien 405
			$index_loc = trim(preg_replace("/\([a-zA-Z0-9]{2}\)$/", "", $ligne[0]));
406
			if(!isset($index[$index_loc])) {
407
				$index[$index_loc] = $ligne;
408
			} else {
409
				if(mb_strlen($ligne[0]) > mb_strlen($index[$index_loc][0])) {
410
					$index[$index_loc] = $ligne;
411
				} else if(mb_strlen($ligne[0]) == mb_strlen($index[$index_loc][0])) {
412
					if(preg_match("/(.)* ([a-zA-Z0-9]{2})/", $ligne[0])) {
413
						$index[$index_loc] = $ligne;
414
					}
415
				}
416
			}
417
		}
418
 
419
		return array_values($index);
420
	}
421
 
2914 mathias 422
	protected function trierListeLieux($a, $b) {
2557 aurelien 423
		return strcmp($a[0], $b[0]);
424
	}
2914 mathias 425
 
426
	protected function testerCoordonneesWgsFrance($latitude, $longitude) {
427
		$coord_france = false;
428
		if ($latitude != '' && $longitude != '') {
429
			if ($latitude < 51.071667 && $latitude > 41.316667) {
430
				if ($longitude < 9.513333 && $longitude > -5.140278) {
431
					$coord_france = true;
432
				}
433
			}
434
		}
435
		return $coord_france;
2538 aurelien 436
	}
2557 aurelien 437
 
2914 mathias 438
	protected function traiterLigneResultatRequeteGeocodingMondiale($res) {
439
 
440
		$lat = "";
441
		$lon = "";
442
		$nom = "";
443
		$code ="";
444
		$code_pays="";
445
		$type="";
446
 
447
		if(!empty($res)) {
448
			$code_pays = $this->extraireCodePays($res);
449
 
450
			$nom = $res['intitule'];
451
			$lat = $res['centre_lat'];
452
			$lon = $res['centre_lng'];
453
			$type = $code_pays == 'FR' || $code_pays == 'FX' ? 'france' : 'monde';
454
 
455
			$retour = array(
456
					'lat' => (float)$lat,
457
					'lng' => (float)$lon,
458
					'nom' => $nom,
459
					'code_zone' => $code,
460
					'code_pays' => $code_pays,
461
					'type' => $type
462
			);
463
		}
464
 
465
		return $retour;
2557 aurelien 466
	}
467
 
2914 mathias 468
	protected function convertirUtmVersLatLong($x, $y, $sector) {
469
		$convertisseur = new gPoint();
470
		$convertisseur->setUTM($x, $y, $sector);
471
		$convertisseur->convertTMtoLL();
472
 
473
		$lat_long = array();
474
		$lat_long['lat'] = str_replace(',','.',$convertisseur->Lat());
475
		$lat_long['lng'] = str_replace(',','.',$convertisseur->Long());
476
		return $lat_long;
2557 aurelien 477
	}
478
 
2914 mathias 479
	protected function extraireCodePays($infos_pays) {
480
		$code = "";
481
		if(!empty($infos_pays['codeIso31661'])) {
482
			$code = $infos_pays['codeIso31661'];
483
		} elseif(!empty($infos_pays['codeIso31662'])) {
484
			// Si pas de code pays, le sous code peut exister et content le code pays
485
			// sous la forme codepays-souscode
486
			$code = substr($infos_pays['codeIso31662'], 0, 2);
487
		}
488
		return $code;
2557 aurelien 489
	}
2914 mathias 490
 
491
	protected function estUneChaineRequeteValide($lieu) {
492
		return (strlen($lieu) > 0) && ($lieu != '%');
2557 aurelien 493
	}
2914 mathias 494
 
495
	protected function formaterLigneResultat($ligne) {
2557 aurelien 496
		if(empty($ligne['code']) || $ligne['code'] == null) {
497
			$res = array($ligne['nom'], '');
498
		} else {
499
			$res = array($ligne['nom'].' ('.substr(sprintf('%02s', $ligne['code']),0,2).')', $ligne['code']);
2914 mathias 500
		}
501
		return $res;
2557 aurelien 502
	}
503
 
2914 mathias 504
	/**
505
	 * Remplace les * par % pour faire des recherches floues
506
	 * Remplace les + par _ (nginx envoie des "+" dans l'URL à la place des espaces)
507
	 * Remplace les espaces et les - par _ car les noms de communes peuvent avoir des espaces ou des tirets
508
	 * @param string $lieu
509
	 * @return string le lieu formaté pour la recherche
510
	 */
511
	protected function formaterChaineLieuPourRequete($lieu) {
512
		$lieu = ltrim($lieu);
513
		$lieu = preg_replace('/\*+/', '%', $lieu);
514
		$lieu = str_replace(['+', ' ', '-'], '_', $lieu);
515
		return $lieu;
2557 aurelien 516
	}
2538 aurelien 517
}