Subversion Repositories eFlore/Applications.cel

Rev

Rev 2565 | Rev 2915 | 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,
152
			'lon_min' => $lngmin,
153
			'lon_max' => $lngmax
154
		);
155
		$retour['groupe']['centroide'] = array(
156
			'lat' => ($latmin + $latmax) / 2,
157
			'lon' => ($lngmin + $lngmax) / 2
158
		);
159
 
160
		return $retour;
161
	}
162
 
163
	public function chercherInfosCommune($projet, $latitude, $longitude) {
164
 
165
		$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;
166
		$url_service = $this->config['cel']['url_service_geo_local'].$valeurs_params;
167
		$url_service = str_replace(',', '.', $url_service);
168
 
169
		$ch = curl_init($url_service);
170
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
171
		$reponse = curl_exec($ch);
172
 
173
		$infos_coord = $this->traiterReponseServiceCommune($reponse);
174
		curl_close($ch);
175
 
176
		return $infos_coord;
177
	}
178
 
2538 aurelien 179
	public function chercherInfosPays($latitude, $longitude) {
2914 mathias 180
 
2538 aurelien 181
		$valeurs_params = '?lat='.$latitude.'&lon='.$longitude;
182
		$url_service = $this->config['cel']['url_service_geo_mondial'].$valeurs_params;
2914 mathias 183
		$url_service = str_replace(',', '.', $url_service);
184
 
185
		$ch = curl_init($url_service);
186
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
187
		$reponse = curl_exec($ch);
188
 
189
		$res = json_decode($reponse, true);
2557 aurelien 190
		$retour = null;
191
 
2914 mathias 192
		$retour = array(
193
					'code_zone' => '',
194
					'code_pays' => '',
195
					'nom' => '',
2557 aurelien 196
					'type' => ''
197
				);
2914 mathias 198
 
199
		if(!empty($res)) {
200
			// Cas où l'on a recherché des coordonnées en France
201
			// On en profite pour remplir plus de données
202
			// que le simple pays si elles sont présentes
203
			if(isset($res[8]) && isset($res[8]['codeInsee'])) {
204
				$retour = array(
2557 aurelien 205
						'code_zone' => $res[8]['codeInsee'],
2914 mathias 206
						'code_pays' => 'FR',
207
						'nom' => $res[8]['nom'],
208
						'type' => 'france'
209
				);
210
			} else {
211
				// Cas de la recherche en dehors de France
212
				// La zone de plus haut niveau est toujours un pays
213
				// (car le niveau de zone est limité à 2)
214
				$infos_pays = $res[min(array_keys($res))];
215
				// La zone de niveau le plus bas est la "localité"
216
				// la plus précise qu'on a pu trouver
217
				$infos_localite = $res[max(array_keys($res))];
218
 
219
				// Cas où l'on a trouvé un code pays
220
				if(!empty($infos_pays['codeIso31661'])) {
221
					$retour = array(
222
							'code_pays' => $infos_pays['codeIso31661'],
223
							'nom' => ''
224
					);
225
				} elseif(!empty($infos_pays['codeIso31662'])) {
226
					// Quelquefois la zone de plus haut niveau est une région ou une province
227
					// heureusement son code est de forme XX-YY ou XX est le code iso du pays !
228
					$retour = array(
229
							'code_pays' => substr($infos_pays['codeIso31662'], 0, 2),
230
							'nom' => ''
231
					);
232
				}
233
 
234
				// Pas de sens de prendre un nom de localité si on a pas de pays
235
				if(!empty($retour['code_pays'])) {
236
					// Type sert à savoir quelle fonction de mise à jour on appellera
237
					$retour['type'] = 'monde';
2557 aurelien 238
					$retour['nom'] = $infos_localite['intitule'];
2914 mathias 239
					$retour['code_zone'] = $infos_localite['codeZoneGeo'];
240
				}
241
			}
2557 aurelien 242
		}
2914 mathias 243
 
244
		return $retour;
245
	}
246
 
247
	public function traiterReponseServiceCommune($reponse) {
248
		$retour = null;
249
		$reponse = json_decode($reponse);
250
		// cas du service lion 1906 qui peut renvoyer plusieurs communes (on prend la première)
251
		if (is_array($reponse)) {
252
			$reponse = $reponse[0];
253
		}
254
		if (isset($reponse->codeINSEE)) {
255
			// Type sert à savoir quelle fonction de mise à jour on appellera
2557 aurelien 256
			$retour = array(
257
						'code_zone' => $reponse->codeINSEE,
2914 mathias 258
						'code_pays' => 'FR',
259
						'nom' => $reponse->nom,
2557 aurelien 260
						'type' => 'france'
2914 mathias 261
					);
262
		}
263
		return $retour;
2538 aurelien 264
	}
2914 mathias 265
 
266
	// @TODO gérer $id_zone (inactif pour l'instant)
267
	protected function effectuerRequeteGeocodingMondiale($nom, $niveau = '2,3,4,5,6,7,8', $limite = 1, $pays = null, $id_zone=null) {
268
		$url_sans_pays = $this->config['cel']['url_service_geo_mondial'].'?masque='.urlencode($nom).'&niveau='.$niveau.'&limite='.$limite;
2564 aurelien 269
		$url = $url_sans_pays.(!empty($pays) ? '&pays='.urlencode($pays) : '');
270
 
2914 mathias 271
		$ch = curl_init($url);
2565 aurelien 272
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
2914 mathias 273
		curl_setopt($ch, CURLOPT_TIMEOUT, $this->maxTimeoutRequeteMondiale);
2565 aurelien 274
		$reponse = curl_exec($ch);
275
 
276
		$res = @json_decode($reponse, true);
2557 aurelien 277
 
2914 mathias 278
		$retour = array();
2557 aurelien 279
 
280
		if(!empty($res)) {
2914 mathias 281
			// cas de la localisation directe d'une commune
2557 aurelien 282
			if($limite == 1) {
283
				$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);
284
			} else {
285
				// cas de l'autocompletion
286
				$retour = array();
287
				foreach($res as $ligne) {
288
					$ligne_decodee = $this->traiterLigneResultatRequeteGeocodingMondiale($ligne);
289
					$param = array('nom' => $ligne_decodee['nom'], 'code' => $ligne_decodee['code_pays']);
290
					$retour[] = $this->formaterLigneResultat($param);
291
				}
292
			}
2564 aurelien 293
		} else if($limite == 1) {
294
			// Dans le cas où un pays a été demandé mais pas de résultat, on rentente la recherche sans
295
			// préciser de pays
296
			$res = @json_decode(file_get_contents($url_sans_pays), true);
297
 
2914 mathias 298
			if(!empty($res)) {
2564 aurelien 299
				$retour = $this->traiterLigneResultatRequeteGeocodingMondiale($res);
300
			}
2914 mathias 301
		}
302
 
303
		return $retour;
2557 aurelien 304
	}
2914 mathias 305
 
306
	/**
307
	 * Renvoie les infos d'une commune, incluant les coordonnées du centroïde,
308
	 * en fonction du nom et/ou du département, et/ou de l'id_zone (plus fiable)
309
	 */
310
	protected function chercherCentroideCommuneBdd($commune, $departement=null, $id_zone=null) {
2557 aurelien 311
		$commune_formatee = str_replace(array(' ', '-'), '_', trim($commune));
2914 mathias 312
		// clauses
313
		$clauses = array();
314
		if ($commune != null) {
315
			$clauses[] = 'nom LIKE '.Cel::db()->proteger($commune_formatee);
316
		}
317
		if ($departement != null) {
318
			$departement = trim($departement, "*");
319
			if (strlen($departement) > 2) {
320
				// @TODO piège : si deux communes ont le même nom dans deux
321
				// DROM-COM, la troncature à 2 chiffres sera insuffisante
322
				$departement = substr($departement, 0, 2);
323
			}
324
			$clauses[] = 'code LIKE '.Cel::db()->proteger($departement.'%');
325
		}
326
		if ($id_zone != null) {
327
			$clauses[] = 'id_zone_geo = '.Cel::db()->proteger($id_zone);
328
		}
329
		$clauses = implode(' AND ', $clauses);
330
 
331
		$requete = 'SELECT utm_x, utm_y, utm_secteur, code FROM cel_zones_geo '
332
				. 'WHERE ' . $clauses
333
				. ' -- '.__FILE__.':'.__LINE__;
334
 
335
		$commune_coordonnees = Cel::db()->requeter($requete);
336
 
337
		$retour = false;
338
		if ($commune_coordonnees && is_array($commune_coordonnees) && count($commune_coordonnees) > 0) {
339
			$lat_lon = $this->convertirUtmVersLatLong($commune_coordonnees[0]['utm_x'],$commune_coordonnees[0]['utm_y'],$commune_coordonnees[0]['utm_secteur']);
340
 
341
			$retour = array(
342
					'lat' => (float) $lat_lon['lat'],
343
					'lng' => (float) $lat_lon['lng'],
344
					'nom' => $commune,
345
					'code_zone' => $commune_coordonnees[0]['code'],
2557 aurelien 346
					'code_pays' => 'FR',
2914 mathias 347
					'type' => 'france'
348
			);
349
		}
350
		return $retour;
2557 aurelien 351
	}
352
 
353
	public function obtenirListeInfosPourNom($lieu) {
354
		$lieux_retour = $this->obtenirListeCommunes($lieu);
355
 
356
		// S'il y a moins de 50 résultats on tente des les fusionner avec une recherche mondiale
357
		// avec 50 - n éléments (où n est la taille du tableau de communes)
358
		if(count($lieux_retour) < 50) {
359
			$lieux_mondiaux = $this->effectuerRequeteGeocodingMondiale($lieu, '2,3,4,5,6,7,8', 50 - count($lieux_retour));
360
			usort($lieux_mondiaux, array($this, 'trierListeLieux'));
361
			$lieux_retour += $lieux_mondiaux;
362
		}
363
 
364
		// Suppression de possibles doublons
365
		$lieux_retour = $this->reduireListeLieux($lieux_retour);
366
 
367
		return $lieux_retour;
368
	}
369
 
370
	public function obtenirListeCommunes($lieu) {
2914 mathias 371
		$lieu_formate = $this->formaterChaineLieuPourRequete($lieu);
372
		$retour = array();
373
 
374
		if ($this->estUneChaineRequeteValide($lieu_formate)) {
375
			$requete = 'SELECT DISTINCT nom, code '.
376
					'FROM cel_zones_geo '.
377
					'WHERE nom LIKE '.Cel::db()->proteger($lieu_formate.'%').' '.
378
					'ORDER BY nom '.
379
					'LIMIT 50 '.
380
					' -- '.__FILE__.':'.__LINE__;
381
 
382
			$liste_lieux = Cel::db()->requeter($requete);
383
			if ($liste_lieux) {
384
				foreach ($liste_lieux as $lieu_trouve) {
385
					$retour[] = $this->formaterLigneResultat($lieu_trouve);
386
				}
387
			}
388
		}
2557 aurelien 389
		return $retour;
390
	}
391
 
392
	/**
393
	 * Fonctions utilitaires
394
	 *
395
	*/
396
 
2914 mathias 397
	protected function reduireListeLieux($tableau) {
2557 aurelien 398
		$index = array();
399
		foreach($tableau as $ligne) {
400
			$index_loc = trim(preg_replace("/\([a-zA-Z0-9]{2}\)$/", "", $ligne[0]));
401
			if(!isset($index[$index_loc])) {
402
				$index[$index_loc] = $ligne;
403
			} else {
404
				if(mb_strlen($ligne[0]) > mb_strlen($index[$index_loc][0])) {
405
					$index[$index_loc] = $ligne;
406
				} else if(mb_strlen($ligne[0]) == mb_strlen($index[$index_loc][0])) {
407
					if(preg_match("/(.)* ([a-zA-Z0-9]{2})/", $ligne[0])) {
408
						$index[$index_loc] = $ligne;
409
					}
410
				}
411
			}
412
		}
413
 
414
		return array_values($index);
415
	}
416
 
2914 mathias 417
	protected function trierListeLieux($a, $b) {
2557 aurelien 418
		return strcmp($a[0], $b[0]);
419
	}
2914 mathias 420
 
421
	protected function testerCoordonneesWgsFrance($latitude, $longitude) {
422
		$coord_france = false;
423
		if ($latitude != '' && $longitude != '') {
424
			if ($latitude < 51.071667 && $latitude > 41.316667) {
425
				if ($longitude < 9.513333 && $longitude > -5.140278) {
426
					$coord_france = true;
427
				}
428
			}
429
		}
430
		return $coord_france;
2538 aurelien 431
	}
2557 aurelien 432
 
2914 mathias 433
	protected function traiterLigneResultatRequeteGeocodingMondiale($res) {
434
 
435
		$lat = "";
436
		$lon = "";
437
		$nom = "";
438
		$code ="";
439
		$code_pays="";
440
		$type="";
441
 
442
		if(!empty($res)) {
443
			$code_pays = $this->extraireCodePays($res);
444
 
445
			$nom = $res['intitule'];
446
			$lat = $res['centre_lat'];
447
			$lon = $res['centre_lng'];
448
			$type = $code_pays == 'FR' || $code_pays == 'FX' ? 'france' : 'monde';
449
 
450
			$retour = array(
451
					'lat' => (float)$lat,
452
					'lng' => (float)$lon,
453
					'nom' => $nom,
454
					'code_zone' => $code,
455
					'code_pays' => $code_pays,
456
					'type' => $type
457
			);
458
		}
459
 
460
		return $retour;
2557 aurelien 461
	}
462
 
2914 mathias 463
	protected function convertirUtmVersLatLong($x, $y, $sector) {
464
		$convertisseur = new gPoint();
465
		$convertisseur->setUTM($x, $y, $sector);
466
		$convertisseur->convertTMtoLL();
467
 
468
		$lat_long = array();
469
		$lat_long['lat'] = str_replace(',','.',$convertisseur->Lat());
470
		$lat_long['lng'] = str_replace(',','.',$convertisseur->Long());
471
		return $lat_long;
2557 aurelien 472
	}
473
 
2914 mathias 474
	protected function extraireCodePays($infos_pays) {
475
		$code = "";
476
		if(!empty($infos_pays['codeIso31661'])) {
477
			$code = $infos_pays['codeIso31661'];
478
		} elseif(!empty($infos_pays['codeIso31662'])) {
479
			// Si pas de code pays, le sous code peut exister et content le code pays
480
			// sous la forme codepays-souscode
481
			$code = substr($infos_pays['codeIso31662'], 0, 2);
482
		}
483
		return $code;
2557 aurelien 484
	}
2914 mathias 485
 
486
	protected function estUneChaineRequeteValide($lieu) {
487
		return (strlen($lieu) > 0) && ($lieu != '%');
2557 aurelien 488
	}
2914 mathias 489
 
490
	protected function formaterLigneResultat($ligne) {
2557 aurelien 491
		if(empty($ligne['code']) || $ligne['code'] == null) {
492
			$res = array($ligne['nom'], '');
493
		} else {
494
			$res = array($ligne['nom'].' ('.substr(sprintf('%02s', $ligne['code']),0,2).')', $ligne['code']);
2914 mathias 495
		}
496
		return $res;
2557 aurelien 497
	}
498
 
2914 mathias 499
	/**
500
	 * Remplace les * par % pour faire des recherches floues
501
	 * Remplace les + par _ (nginx envoie des "+" dans l'URL à la place des espaces)
502
	 * Remplace les espaces et les - par _ car les noms de communes peuvent avoir des espaces ou des tirets
503
	 * @param string $lieu
504
	 * @return string le lieu formaté pour la recherche
505
	 */
506
	protected function formaterChaineLieuPourRequete($lieu) {
507
		$lieu = ltrim($lieu);
508
		$lieu = preg_replace('/\*+/', '%', $lieu);
509
		$lieu = str_replace(['+', ' ', '-'], '_', $lieu);
510
		return $lieu;
2557 aurelien 511
	}
2538 aurelien 512
}