Subversion Repositories eFlore/Applications.cel

Rev

Rev 3857 | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
1500 delphine 1
<?php
2414 jpm 2
// declare(encoding='UTF-8');
1500 delphine 3
/**
2533 aurelien 4
 * Traitement des codes INSEE et des codes pays de la base de données CEL
5
 * Permet d'affecter des codes INSEE et des codes pays aux obs possédant des coordonnées mais pas de commune ou de pays.
1500 delphine 6
 *
2414 jpm 7
 * Utilisation :
2533 aurelien 8
 * - Pour ? : <code>/opt/lamp/bin/php cli.php code_zone_geo -a tout</code>
9
 * - Pour ? : <code>/opt/lamp/bin/php cli.php code_zone_geo -a colonne</code>
10
 * - Pour ? : <code>/opt/lamp/bin/php cli.php code_zone_geo -a payssanscorrespondance</code>
11
 * - Pour ? : <code>/opt/lamp/bin/php cli.php code_zone_geo -a codeinseesanscorrespondance</code>
1500 delphine 12
 *
2414 jpm 13
 * @category   CEL
14
 * @package    Scripts
2533 aurelien 15
 * @subpackage Zones Geo
2414 jpm 16
 * @author     Mathias CHOUET <mathias@tela-botanica.org>
17
 * @author     Aurelien PERONNET <aurelien@tela-botanica.org>
18
 * @license    GPL v3 <http://www.gnu.org/licenses/gpl.txt>
19
 * @license    CECILL v2 <http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt>
20
 * @copyright  1999-2014 Tela Botanica <accueil@tela-botanica.org>
1500 delphine 21
 */
2414 jpm 22
include_once dirname(__FILE__).'/bibliotheque/Dao.php';
23
 
2533 aurelien 24
class CodeZoneGeo extends Script {
2414 jpm 25
 
1500 delphine 26
	protected $tableauTaxon;
27
	protected $dao;
28
	protected $observations;
2533 aurelien 29
 
30
	protected $url_service_geo_tpl;
31
	protected $url_service_pays_tpl;
2414 jpm 32
 
1523 aurelien 33
	protected $mode_verbeux = false;
2414 jpm 34
 
1500 delphine 35
	// Paramêtres autorisées lors de l'appel au script en ligne de commande
36
	protected $parametres_autorises = array(
37
		'-n' => array(true, true, 'Nom du fichier ou du dossier à traiter'));
38
 
39
	public function executer() {
2533 aurelien 40
 
41
		$this->url_service_geo_tpl = Config::get('urlZoneGeoTpl');
42
		$this->url_service_pays_tpl = Config::get('urlPaysTpl');
2414 jpm 43
 
1500 delphine 44
		$this->dao = new Dao();
45
		// Récupération de paramétres
46
		// Lancement de l'action demandée
47
		$cmd = $this->getParametre('a');
1523 aurelien 48
		$this->mode_verbeux = $this->getParametre('v');
2414 jpm 49
 
1500 delphine 50
		switch ($cmd) {
51
			case 'tout' :
1517 aurelien 52
				// à effectuer manuellement au premier lancement du script
1500 delphine 53
				$this->creerColonne();
2533 aurelien 54
				$this->calculerToutSansCorrespondanceDepuisLeDebut();
1500 delphine 55
				break;
56
			case 'colonne' :
57
				$this->creerColonne();
58
				break;
2533 aurelien 59
			case 'codeinseesanscorrespondance' :
1523 aurelien 60
				$this->calculerCodesInseeSansCorrespondance();
61
				break;
2533 aurelien 62
			case 'codeinseesanscorrespondancedepuisledebut' :
63
				$this->calculerCodesInseeSansCorrespondanceDepuisLeDebut();
1500 delphine 64
				break;
2533 aurelien 65
			case 'payssanscorrespondance' :
66
				$this->calculerPaysSansCorrespondance();
67
				break;
68
			case 'payssanscorrespondancedepuisledebut' :
69
				$this->calculerPaysSansCorrespondanceDepuisLeDebut();
70
				break;
71
			case 'toutsanscorrespondance' :
72
				// devrait être être l'option utilisée dans le cron quotidien (avec l'option -v 0)
73
				$this->calculerToutSansCorrespondance();
74
				break;
2536 aurelien 75
			case 'toutsanscorrespondancedepuisledebut' :
2533 aurelien 76
				$this->calculerToutSansCorrespondanceDepuisLeDebut();
77
				break;
1500 delphine 78
			default :
79
				$this->traiterErreur('Erreur : la commande "%s" n\'existe pas!', array($cmd));
80
		}
81
	}
82
 
2533 aurelien 83
	protected function creerColonne() {
1523 aurelien 84
		$this->dao->creerColonneCodeInseeCalcule();
1500 delphine 85
	}
2414 jpm 86
 
2533 aurelien 87
	protected function calculerCodesInseeSansCorrespondance() {
88
		$liste_coordonnees = $this->dao->rechercherCoordonneesFrancaisesSansCorrespondances();
1517 aurelien 89
		$this->traiterCoordonnees($liste_coordonnees);
90
	}
2533 aurelien 91
 
92
	protected function calculerCodesInseeSansCorrespondanceDepuisLeDebut() {
93
		$liste_coordonnees = $this->dao->rechercherCoordonneesFrancaisesSansCorrespondanceDepuisLeDebut();
94
		$this->traiterCoordonnees($liste_coordonnees);
95
	}
96
 
97
	protected function calculerPaysSansCorrespondance() {
98
		$liste_coordonnees = $this->dao->rechercherPaysSansCorrespondance();
99
		$this->traiterCoordonnees($liste_coordonnees);
100
	}
101
 
102
	protected function calculerPaysSansCorrespondanceDepuisLeDebut() {
103
		$liste_coordonnees = $this->dao->rechercherPaysSansCorrespondanceDepuisLeDebut();
104
		$this->traiterCoordonnees($liste_coordonnees);
105
	}
106
 
107
	protected function calculerToutSansCorrespondance() {
2559 aurelien 108
		$this->affecterPaysFranceAuxCommunes();
2533 aurelien 109
		$liste_coordonnees = $this->dao->rechercherToutSansCorrespondance();
110
		$this->traiterCoordonnees($liste_coordonnees);
111
	}
2414 jpm 112
 
2533 aurelien 113
	protected function calculerToutSansCorrespondanceDepuisLeDebut() {
2559 aurelien 114
		$this->affecterPaysFranceAuxCommunes();
2533 aurelien 115
		$liste_coordonnees = $this->dao->rechercherSansCorrespondanceDepuisLeDebut();
1517 aurelien 116
		$this->traiterCoordonnees($liste_coordonnees);
117
	}
2559 aurelien 118
 
119
	protected function affecterPaysFranceAuxCommunes() {
120
		if ($this->mode_verbeux) {
121
			echo "-------------------------------------------------------------------\n";
122
			echo "  Affectation du pays France aux communes ayant un code INSEE mais pas de pays ... ";
123
		}
124
 
125
		$maj = $this->dao->affecterPaysFranceAuxCommunesRenseignees();
126
 
127
		if ($this->mode_verbeux) {
128
			echo "OK \n";
129
			echo "  ".$maj." commune".($maj > 1 ? 's' : '')." affectée".($maj > 1 ? 's' : '')."\n";
130
			echo "-------------------------------------------------------------------\n";
131
		}
132
	}
2414 jpm 133
 
1517 aurelien 134
	private function traiterCoordonnees($liste_coordonnees) {
135
		$debut = microtime(true);
1523 aurelien 136
		$nb_coordonnees_modifiees = 0;
137
		$nb_coordonnees_ignorees = 0;
1517 aurelien 138
		$total = count($liste_coordonnees);
2414 jpm 139
 
140
		if ($this->mode_verbeux) {
1523 aurelien 141
			echo "-------------------------------------------------------------------\n";
2533 aurelien 142
			echo "  Calcul des codes INSEE, du code pays et modification (".count($liste_coordonnees)." coordonnées en tout) ... \n";
1523 aurelien 143
			echo "  Enrichissement des champs zone geo et ce_zone_geo vides ... \n";
144
			echo "-------------------------------------------------------------------\n";
145
		}
2414 jpm 146
 
1500 delphine 147
		foreach ($liste_coordonnees as $coordonnees) {
1523 aurelien 148
			$infos_coord = $this->obtenirInfosPourCoordonnees($coordonnees);
149
			if ($infos_coord != null) {
2533 aurelien 150
				// Le type renvoyé permet de connaitre ce qu'on peut modifier
151
				if($infos_coord['type'] == 'france') {
152
					// Cas d'une commune française
153
					$coordonnee_a_traiter = array(
154
						'latitude' => $coordonnees['latitude'],
155
						'longitude' => $coordonnees['longitude'],
156
						'code_insee' => $infos_coord['code_insee'],
157
						'nom' => $infos_coord['nom']
158
					);
159
					$this->dao->modifierCodeInseeEtZoneGeo($coordonnee_a_traiter);
160
				} else if($infos_coord['type'] == 'pays') {
161
					// Cas d'un lieu à l'étranger
162
					$coordonnee_a_traiter = array(
163
							'latitude' => $coordonnees['latitude'],
164
							'longitude' => $coordonnees['longitude'],
165
							'code_pays' => $infos_coord['code_pays'],
166
							'nom' => $infos_coord['nom']
167
					);
168
					$this->dao->modifierPays($coordonnee_a_traiter);
169
				}
1523 aurelien 170
				$nb_coordonnees_modifiees++;
171
			} else {
172
				$nb_coordonnees_ignorees++;
1500 delphine 173
			}
2414 jpm 174
			if ($this->mode_verbeux) {
1523 aurelien 175
				$this->afficherAvancement('  Coordonnées traitées ', $nb_coordonnees_modifiees);
176
			}
1500 delphine 177
		}
2414 jpm 178
 
1517 aurelien 179
		$fin = microtime(true);
2414 jpm 180
		if ($this->mode_verbeux) {
1523 aurelien 181
			echo "\n";
182
			echo "-------------------------------------------------------------------\n";
2533 aurelien 183
			echo "  Fin de la mise à jour des codes INSEE calculés, du pays et de l'enrichissement, \n";
1523 aurelien 184
			echo "  ".($fin - $debut)." secondes écoulées \n";
185
			echo "  ".$nb_coordonnees_modifiees." code insee calculés et modifiés \n";
186
			echo "  ".$nb_coordonnees_ignorees." coordonnées ignorées \n";
187
			echo "-------------------------------------------------------------------\n";
188
			echo "\n";
189
		}
1500 delphine 190
	}
2414 jpm 191
 
1523 aurelien 192
	private function obtenirInfosPourCoordonnees($coordonnees) {
193
		$infos_coord = null;
2533 aurelien 194
		// Test facile qui permet d'abord de tenter une localisation bien plus rapide si les coordonnées
195
		// sont dans la bounding box approximative de la France
196
		if($this->testerCoordonneesWgsFrance($coordonnees['latitude'], $coordonnees['longitude'])) {
1523 aurelien 197
			$infos_coord = $this->chercherInfosCommune('osm', $coordonnees['latitude'], $coordonnees['longitude']);
198
			if ($infos_coord == null) {
2533 aurelien 199
				// Sinon recherche par pays
200
				$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);
1523 aurelien 201
			}
2533 aurelien 202
		} else {
203
			// Recherche par pays immédiate si en dehors de la bouding box française
204
			$infos_coord = $this->chercherInfosPays($coordonnees['latitude'], $coordonnees['longitude']);
1523 aurelien 205
		}
206
		return $infos_coord;
207
	}
2414 jpm 208
 
1523 aurelien 209
	private function chercherInfosCommune($projet, $latitude, $longitude) {
2533 aurelien 210
 
211
		$noms_params = array('{projet}', '{latitude}', '{longitude}');
212
		$valeurs_params = array($projet, $latitude, $longitude);
213
 
214
		$url_service = str_replace($noms_params, $valeurs_params, $this->url_service_geo_tpl);
215
		$url_service = str_replace(',', '.', $url_service);
2414 jpm 216
 
1500 delphine 217
		$ch = curl_init($url_service);
218
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
219
		$reponse = curl_exec($ch);
2414 jpm 220
 
1523 aurelien 221
		$infos_coord = $this->traiterReponseServiceCommune($reponse);
1518 aurelien 222
		curl_close($ch);
2414 jpm 223
 
1523 aurelien 224
		return $infos_coord;
1518 aurelien 225
	}
2533 aurelien 226
 
227
	private function chercherInfosPays($latitude, $longitude) {
228
		$noms_params = array('{latitude}', '{longitude}');
229
		$valeurs_params = array($latitude, $longitude);
230
		$url_service = str_replace($noms_params, $valeurs_params, $this->url_service_pays_tpl);
231
		$url_service = str_replace(',', '.', $url_service);
1518 aurelien 232
 
2533 aurelien 233
		$ch = curl_init($url_service);
234
		curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
235
		$reponse = curl_exec($ch);
236
 
237
		$res = json_decode($reponse, true);
238
		$retour = null;
239
 
240
		if(!empty($res)) {
241
			// Cas où l'on a recherché des coordonnées en France
242
			// On en profite pour remplir plus de données
243
			// que le simple pays si elles sont présentes
244
			if(isset($res[8]) && isset($res[8]['codeInsee'])) {
245
				$retour = array(
246
					'code_insee' => $res[8]['codeInsee'],
247
					'nom' => $res[8]['intitule'],
248
					'type' => 'france'
249
				);
250
			} else {
251
				// Cas de la recherche en dehors de France
252
				// La zone de plus haut niveau est toujours un pays
253
				// (car le niveau de zone est limité à 2)
254
				$infos_pays = $res[min(array_keys($res))];
255
				// La zone de niveau le plus bas est la "localité"
256
				// la plus précise qu'on a pu trouver
257
				$infos_localite = $res[max(array_keys($res))];
258
 
259
				// Cas où l'on a trouvé un code pays
260
				if(!empty($infos_pays['codeIso31661'])) {
261
					$retour = array(
262
						'code_pays' => $infos_pays['codeIso31661'],
263
						'nom' => ''
264
					);
265
				} elseif(!empty($infos_pays['codeIso31662'])) {
266
					// Quelquefois la zone de plus haut niveau est une région ou une province
267
					// heureusement son code est de forme XX-YY ou XX est le code iso du pays !
268
					$retour = array(
269
						'code_pays' => substr($infos_pays['codeIso31662'], 0, 2),
270
						'nom' => ''
271
					);
272
				}
273
 
274
				// Pas de sens de prendre un nom de localité si on a pas de pays
275
				if(!empty($retour['code_pays'])) {
276
					// Type sert à savoir quelle fonction de mise à jour on appellera
277
					$retour['type'] = 'pays';
278
					$retour['nom'] = $infos_localite['intitule'];
279
				}
280
			}
281
		}
282
		return $retour;
283
	}
284
 
1523 aurelien 285
	private function traiterReponseServiceCommune($reponse) {
286
		$retour = null;
1500 delphine 287
		$reponse = json_decode($reponse);
1518 aurelien 288
		// cas du service lion 1906 qui peut renvoyer plusieurs communes (on prend la première)
2414 jpm 289
		if (is_array($reponse)) {
1518 aurelien 290
			$reponse = $reponse[0];
291
		}
1500 delphine 292
		if (isset($reponse->codeINSEE)) {
2533 aurelien 293
			// Type sert à savoir quelle fonction de mise à jour on appellera
1523 aurelien 294
			$retour = array('code_insee' => $reponse->codeINSEE,
2533 aurelien 295
							'nom' => $reponse->nom,
296
							'type' => 'france');
2414 jpm 297
		}
1523 aurelien 298
		return $retour;
1500 delphine 299
	}
2533 aurelien 300
 
301
	private function testerCoordonneesWgsFrance($latitude, $longitude) {
302
		$coord_france = false;
303
		if ($latitude != '' && $longitude != '') {
304
			if ($latitude < 51.071667 && $latitude > 41.316667) {
305
				if ($longitude < 9.513333 && $longitude > -5.140278) {
306
					$coord_france = true;
307
				}
308
			}
309
		}
310
		return $coord_france;
1500 delphine 311
	}
2414 jpm 312
}