Subversion Repositories eFlore/Applications.cel

Rev

Rev 1106 | Go to most recent revision | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 1106 Rev 1142
1
<?php
1
<?php
2
class CartoGroupage {
2
class CartoGroupage {
-
 
3
 
-
 
4
	const MARQUEUR_GROUPE = 'GROUPE';
-
 
5
	const MARQUEUR_COMMUNE = 'COMMUNE';
-
 
6
	const MARQUEUR_STATION = 'STATION';
3
 
7
	
4
	private static $seuilClusterisation = 200;
8
	private static $seuilClusterisation = 200;
5
	private static $zoomDefaut = 3;
9
	private static $zoomDefaut = 3;
6
	private static $zoomMaxClustering = 12;
10
	private static $zoomMaxClustering = 12;
7
	private static $pasZoomDefaut = 1;
11
	private static $pasZoomDefaut = 1;
8
	private static $pasZoomMaxClustering = 0.05;
12
	private static $pasZoomMaxClustering = 0.05;
9
	private static $profondeurMin = 0;
13
	private static $profondeurMin = 0;
10
	private static $profondeurMax = 8;
14
	private static $profondeurMax = 8;
11
	
15
	
12
	private static $pasCorrectionCentre = null;
16
	private static $pasCorrectionCentre = null;
13
	private static $coefficientReductionPas = null;	
17
	private static $coefficientReductionPas = null;	
14
	private static $coefficientProfondeurMax = null;
18
	private static $coefficientProfondeurMax = null;
15
	
19
	
16
	private static $nbElements = array('stations' => 0,'communes' => 0, 'observations' => 0);
20
	private static $nbElements = array('stations' => 0,'communes' => 0, 'observations' => 0);
17
 
-
 
18
	private static $listeNoeudsSelectionnes = array();
21
	private static $listeNoeudsSelectionnes = array();	
19
	
-
 
20
	private static $pointsDejaTraites = array();
-
 
21
	
-
 
22
	private static $bornesMax = array('latMin' => null, 'lngMin' => null, 'latMax' => null, 'lngMax' => null);
22
	private static $bornesMax = array('latMin' => null, 'lngMin' => null, 'latMax' => null, 'lngMax' => null);
23
 
-
 
-
 
23
	private static $id_traites = array();
24
	/*
24
	/*
25
	  +---------+---------+
25
	  +---------+---------+
26
	  |         |         |
26
	  |         |         |
27
	  |    A    |    B    |
27
	  |    A    |    B    |
28
	  |         |         |
28
	  |         |         |
29
	  +---------*---------+
29
	  +---------*---------+
30
	  |         |         |
30
	  |         |         |
31
	  |    D    |    C    |
31
	  |    D    |    C    |
32
	  |         |         |
32
	  |         |         |
33
	  +---------+---------+
33
	  +---------+---------+
34
	  
34
	  
35
	  Quatres cadrans sont considérés par le quad tree
35
	  Quatres cadrans sont considérés par le quad tree
36
	  * = centre de la fenetre 
36
	  * = centre de la fenetre 
37
	 */
37
	 */
38
	public static function creerGroupesQuadtree(&$markers, $neLat, $neLng, $swLat, $swLng, $zoom = 3) {
38
	public static function creerGroupesQuadtree(&$markers, $neLat, $neLng, $swLat, $swLng, $zoom = 3) {
39
		if(count($markers) > self::$seuilClusterisation) {
39
		if(count($markers) > self::$seuilClusterisation) {
40
						
40
						
41
			self::calculerProfondeurMax($zoom);
41
			self::calculerProfondeurMax($zoom);
42
			self::calculerPasCorrectionCentre($zoom);
42
			self::calculerPasCorrectionCentre($zoom);
43
 
43
 
44
			$noeudRacine = array('nbrePoints' => count($markers), 'points' => $markers);
44
			$noeudRacine = array('nbrePoints' => count($markers), 'points' => $markers);
45
			self::attribuerAuCadran($noeudRacine, $neLat, $neLng, $swLat, $swLng);	
45
			self::attribuerAuCadran($noeudRacine, $neLat, $neLng, $swLat, $swLng);	
46
			
46
			
47
		} else {
47
		} else {
48
			foreach($markers as $marker) {				
48
			foreach($markers as &$marker) {	
-
 
49
				if(!self::estUnPointAExclure($marker)) {
-
 
50
					$emplacement = self::formaterPointPourAjout(&$marker);
49
				self::mettreAJourBornes(&$marker);
51
					self::mettreAJourBornes(&$marker);
50
				$points = array($marker);
52
					$points = array($marker);
51
				$noeudSimple = array('points' => $points, 'nbrePoints' => 1);
53
					$noeudSimple = array('points' => $points, 'nbrePoints' => 1);
-
 
54
					self::$nbElements[$emplacement]++;
-
 
55
					self::$listeNoeudsSelectionnes[] = self::ajouterGroupeOuPoint($noeudSimple);
-
 
56
				}
52
				self::$nbElements['observations']++;
57
				self::$nbElements['observations']++;
53
				$emplacement = isset($marker['type_emplacement']) ? $marker['type_emplacement'] : self::obtenirTypeEmplacementParId(&$marker);
-
 
54
				self::$nbElements[$emplacement]++;
-
 
55
				unset($marker['type_emplacement']);
-
 
56
				self::$listeNoeudsSelectionnes[] = self::ajouterGroupeOuPoint($noeudSimple);
-
 
57
			}
58
			}
58
		}
59
		}
59
			
-
 
60
		return self::$listeNoeudsSelectionnes;
60
		return self::$listeNoeudsSelectionnes;
61
	}
61
	}
62
	
62
	
63
	private function calculerCoefficientReductionPas() {
63
	private function calculerCoefficientReductionPas() {
64
		if(self::$coefficientReductionPas == null) {
64
		if(self::$coefficientReductionPas == null) {
65
			self::$coefficientReductionPas = (self::$pasZoomMaxClustering - self::$pasZoomDefaut)/(self::$zoomMaxClustering - self::$zoomDefaut);
65
			self::$coefficientReductionPas = (self::$pasZoomMaxClustering - self::$pasZoomDefaut)/(self::$zoomMaxClustering - self::$zoomDefaut);
66
		}
66
		}
67
		
67
		
68
		return self::$coefficientReductionPas;
68
		return self::$coefficientReductionPas;
69
	}
69
	}
70
	
70
	
71
	private function calculerPasCorrectionCentre($zoom) {
71
	private function calculerPasCorrectionCentre($zoom) {
72
		self::$pasCorrectionCentre = ($zoom - self::$zoomDefaut) * self::calculerCoefficientReductionPas() + self::$pasZoomDefaut;
72
		self::$pasCorrectionCentre = ($zoom - self::$zoomDefaut) * self::calculerCoefficientReductionPas() + self::$pasZoomDefaut;
73
	}
73
	}
74
	
74
	
75
	private function calculerCoefficientProfondeurMax() {
75
	private function calculerCoefficientProfondeurMax() {
76
		if(self::$coefficientProfondeurMax == null) {
76
		if(self::$coefficientProfondeurMax == null) {
77
			self::$coefficientProfondeurMax = (self::$profondeurMax - self::$profondeurMin)/(self::$zoomMaxClustering - self::$zoomDefaut);
77
			self::$coefficientProfondeurMax = (self::$profondeurMax - self::$profondeurMin)/(self::$zoomMaxClustering - self::$zoomDefaut);
78
		}
78
		}
79
		
79
		
80
		return self::$coefficientProfondeurMax;
80
		return self::$coefficientProfondeurMax;
81
	}
81
	}
82
	
82
	
83
	private function calculerProfondeurMax($zoom) {	
83
	private function calculerProfondeurMax($zoom) {	
84
		if($zoom > self::$zoomDefaut) {	
84
		if($zoom > self::$zoomDefaut) {	
85
			self::$profondeurMax = round(($zoom - self::$zoomDefaut) * self::calculerCoefficientProfondeurMax() + self::$profondeurMin,0);
85
			self::$profondeurMax = round(($zoom - self::$zoomDefaut) * self::calculerCoefficientProfondeurMax() + self::$profondeurMin,0);
86
		} else {
86
		} else {
87
			self::$profondeurMax = 1;
87
			self::$profondeurMax = 1;
88
		}
88
		}
89
	}
89
	}
90
	
90
	
91
	public static function getNbElements() {
91
	public static function getNbElements() {
92
		return self::$nbElements;
92
		return self::$nbElements;
93
	}
93
	}
94
	
94
	
95
	public static function mettreAJourBornes(&$point) {
95
	public static function mettreAJourBornes(&$point) {
-
 
96
		
-
 
97
		if(!is_numeric($point['lat']) || abs($point['lat']) > 90) return;
-
 
98
		if(!is_numeric($point['lng']) || abs($point['lng']) > 180) return;
-
 
99
		
96
		self::$bornesMax['latMin'] = (is_numeric($point['lat']) && $point['lat'] < self::$bornesMax['latMin'] || self::$bornesMax['latMin'] == null) ? $point['lat'] : self::$bornesMax['latMin'] ;
100
		self::$bornesMax['latMin'] = ($point['lat'] < self::$bornesMax['latMin'] || self::$bornesMax['latMin'] == null) ? $point['lat'] : self::$bornesMax['latMin'] ;
97
		self::$bornesMax['lngMin'] = (is_numeric($point['lng']) && $point['lng'] < self::$bornesMax['lngMin'] || self::$bornesMax['lngMin'] == null) ? $point['lng'] : self::$bornesMax['lngMin'] ;
101
		self::$bornesMax['lngMin'] = ($point['lng'] < self::$bornesMax['lngMin'] || self::$bornesMax['lngMin'] == null) ? $point['lng'] : self::$bornesMax['lngMin'] ;
98
		self::$bornesMax['latMax'] = (is_numeric($point['lat']) && $point['lat'] > self::$bornesMax['latMax'] || self::$bornesMax['latMax'] == null) ? $point['lat'] : self::$bornesMax['latMax'] ;
102
		self::$bornesMax['latMax'] = ($point['lat'] > self::$bornesMax['latMax'] || self::$bornesMax['latMax'] == null) ? $point['lat'] : self::$bornesMax['latMax'] ;
99
		self::$bornesMax['lngMax'] = (is_numeric($point['lng']) && $point['lng'] > self::$bornesMax['lngMax'] || self::$bornesMax['lngMax'] == null) ? $point['lng'] : self::$bornesMax['lngMax'] ;
103
		self::$bornesMax['lngMax'] = ($point['lng'] > self::$bornesMax['lngMax'] || self::$bornesMax['lngMax'] == null) ? $point['lng'] : self::$bornesMax['lngMax'] ;
100
	}
104
	}
101
	
105
	
102
	public static function getBornes() {
106
	public static function getBornes() {
103
		return self::$bornesMax;
107
		return self::$bornesMax;
104
	}
108
	}
105
	
109
	
106
	/**
110
	/**
107
	 * 
111
	 * 
108
	 * @param mixed $noeud Le noeud à traiter par le quadtree
112
	 * @param mixed $noeud Le noeud à traiter par le quadtree
109
	 * @param float $neLat Latitude du coin nord est de la fenetre 
113
	 * @param float $neLat Latitude du coin nord est de la fenetre 
110
	 * @param float $neLng Longitude du coin nord est de la fenetre 
114
	 * @param float $neLng Longitude du coin nord est de la fenetre 
111
	 * @param float $swLat Latitude du coin sud ouest de la fenetre 
115
	 * @param float $swLat Latitude du coin sud ouest de la fenetre 
112
	 * @param float $swLng Longitude du coin sud ouest de la fenetre 
116
	 * @param float $swLng Longitude du coin sud ouest de la fenetre 
113
	 * @param int $profondeur profondeur courante de l'arbre
117
	 * @param int $profondeur profondeur courante de l'arbre
114
	 */
118
	 */
115
	private static function attribuerAuCadran(&$noeud, $neLat, $neLng, $swLat, $swLng, $profondeur = 0) {
119
	private static function attribuerAuCadran(&$noeud, $neLat, $neLng, $swLat, $swLng, $profondeur = 0) {
116
			
120
			
117
		$latCentre = round((($neLat+$swLat)/2)/self::$pasCorrectionCentre,0)*self::$pasCorrectionCentre;
121
		$latCentre = round((($neLat+$swLat)/2)/self::$pasCorrectionCentre,0)*self::$pasCorrectionCentre;
118
		$lngCentre = round((($neLng+$swLng)/2)/self::$pasCorrectionCentre,0)*self::$pasCorrectionCentre;
122
		$lngCentre = round((($neLng+$swLng)/2)/self::$pasCorrectionCentre,0)*self::$pasCorrectionCentre;
119
						
123
						
-
 
124
		foreach ($noeud['points'] as &$point) {			
-
 
125
			if(!self::estUnPointAExclure($point)) {
120
		foreach ($noeud['points'] as &$point) {
126
				$emplacement = self::formaterPointPourAjout(&$point);
121
				self::mettreAJourBornes(&$point);
-
 
122
				self::$nbElements['observations']++;
-
 
123
				$emplacement = isset($point['type_emplacement']) ? $point['type_emplacement'] : self::obtenirTypeEmplacementParId(&$point);
-
 
124
				self::$nbElements[$emplacement]++;
-
 
125
				unset($point['type_emplacement']);
127
				self::mettreAJourBornes(&$point);
126
				$cadran = self::obtenirCadranPourPoint($latCentre, $lngCentre, $point);
128
				$cadran = self::obtenirCadranPourPoint($latCentre, $lngCentre, $point);
-
 
129
				self::ajouterFils($noeud,$cadran,$point);
-
 
130
				self::$nbElements[$emplacement]++;
-
 
131
			}
127
				self::ajouterFils($noeud,$cadran,$point);
132
			self::$nbElements['observations']++;
128
		}
133
		}
129
		
134
		
130
		$profondeur++;
135
		$profondeur++;
131
		
136
		
132
		if($profondeur <= self::$profondeurMax) {
137
		if($profondeur <= self::$profondeurMax) {
133
			(isset($noeud['A']) && $noeud['A'] != null) ? self::attribuerAuCadran($noeud['A'], $neLat, $lngCentre , $latCentre, $swLng, $profondeur) : '';
138
			(isset($noeud['A']) && $noeud['A'] != null) ? self::attribuerAuCadran($noeud['A'], $neLat, $lngCentre , $latCentre, $swLng, $profondeur) : '';
134
			(isset($noeud['B']) && $noeud['B'] != null) ? self::attribuerAuCadran($noeud['B'], $neLat, $neLng, $latCentre, $lngCentre, $profondeur) : '';
139
			(isset($noeud['B']) && $noeud['B'] != null) ? self::attribuerAuCadran($noeud['B'], $neLat, $neLng, $latCentre, $lngCentre, $profondeur) : '';
135
			(isset($noeud['C']) && $noeud['C'] != null) ? self::attribuerAuCadran($noeud['C'], $latCentre, $neLng, $swLat, $lngCentre, $profondeur) : '';
140
			(isset($noeud['C']) && $noeud['C'] != null) ? self::attribuerAuCadran($noeud['C'], $latCentre, $neLng, $swLat, $lngCentre, $profondeur) : '';
136
			(isset($noeud['D']) && $noeud['D'] != null) ? self::attribuerAuCadran($noeud['D'], $latCentre, $lngCentre, $swLat, $swLng, $profondeur) : '';
141
			(isset($noeud['D']) && $noeud['D'] != null) ? self::attribuerAuCadran($noeud['D'], $latCentre, $lngCentre, $swLat, $swLng, $profondeur) : '';
137
		}
142
		}
138
		
143
		
139
		if(self::estUnParentFeuilles($noeud)) {
144
		if(self::estUnParentFeuilles($noeud)) {
140
			self::$listeNoeudsSelectionnes[] = self::ajouterGroupeOuPoint($noeud);
145
			self::$listeNoeudsSelectionnes[] = self::ajouterGroupeOuPoint($noeud);
141
		}
146
		}
142
	}
147
	}
143
	
148
	
-
 
149
	private static function estUnPointAExclure(&$point) {
144
	private static function obtenirTypeEmplacementParId($point) {
150
		return self::estSensible($point) && 
-
 
151
		self::coordonneesCommuneSontNulles($point);
-
 
152
	}
-
 
153
	
-
 
154
	private static function coordonneesCommuneSontNulles(&$point) {
145
		$tableau_point_id = explode(':',$point['id'],2);
155
		$coord_nulles = ($point['wgs84_latitude'] == null ||
-
 
156
		$point['wgs84_latitude'] == '' ||
-
 
157
		$point['wgs84_longitude'] == null ||
146
		switch($tableau_point_id[0]) {
158
		$point['wgs84_longitude'] == '');
-
 
159
		return $coord_nulles;
-
 
160
	}
-
 
161
	
-
 
162
	private static function coordonneesSontNulles(&$point) {
-
 
163
		
-
 
164
		$coord_nulles = ($point['coord_x'] == '000null' || 
147
			case 'STATION':
165
				$point['coord_x'] == '' ||
-
 
166
				$point['coord_y'] == '000null' ||
-
 
167
				$point['coord_y'] == '');
-
 
168
		return $coord_nulles;
-
 
169
	}
-
 
170
	
-
 
171
	private static function estSensible(&$point) {
148
				$type_emplacement = 'stations';
172
		$sensible = isset($point['mots_cles_texte']) && substr_count($point['mots_cles_texte'], 'sensible') != 0;
-
 
173
		return $sensible;
-
 
174
	}
-
 
175
	
149
			break;
176
	private static function formaterNomStation(&$point) {
-
 
177
		$station = '';
150
			case 'COMMUNE':
178
		if($point['station'] != '' && $point['station'] != '000null') {
151
				$type_emplacement = 'communes';
179
			$station = $point['station'];
-
 
180
		} else {
-
 
181
			$station = $point['location'].(($point['id_location'] != '' && $point['id_location'] != '000null') ? '('.$point['id_location'].')' : '');
-
 
182
		}
-
 
183
		
-
 
184
		return $station;
-
 
185
	}
-
 
186
	
-
 
187
	private static function formaterPointPourAjout(&$point) {
-
 
188
		
-
 
189
		if(isset($point['type_emplacement'])) {
152
			break;
190
			return $point['type_emplacement'];
-
 
191
		}
-
 
192
		
-
 
193
		if(self::coordonneesSontNulles($point) || self::estSensible($point)) {
-
 
194
				$point['id'] = self::MARQUEUR_COMMUNE.':'.$point['wgs84_latitude'].'|'.$point['wgs84_longitude'];
-
 
195
				$point['type_emplacement'] = 'communes';
-
 
196
				$point['lat'] = (float)$point['wgs84_latitude'];
-
 
197
				$point['lng'] = (float)$point['wgs84_longitude'];
-
 
198
		} else {
-
 
199
			$point_allege = array();
-
 
200
			$point_allege['id'] = self::MARQUEUR_STATION.':'.$point['coord_x'].'|'.$point['coord_y'];
-
 
201
			$point_allege['type_emplacement'] = 'stations';
-
 
202
			$point_allege['nom'] = self::formaterNomStation($point);
-
 
203
			$point_allege['lat'] = (float)$point['coord_x'];
-
 
204
			$point_allege['lng'] = (float)$point['coord_y'];
-
 
205
			
-
 
206
			$point = $point_allege;
-
 
207
		}
153
		}
208
		
154
		return $type_emplacement;
209
		return $point['type_emplacement'];
155
	}
210
	}
156
	
211
	
157
	private function obtenirCadranPourPoint($latCentre,$lngCentre, &$point) {
212
	private function obtenirCadranPourPoint($latCentre,$lngCentre, &$point) {
158
		if ($point['lng'] < $lngCentre) {
213
		if ($point['lng'] < $lngCentre) {
159
			if ($point['lat'] > $latCentre) {
214
			if ($point['lat'] > $latCentre) {
160
					$cadran = 'A';
215
					$cadran = 'A';
161
				} else {
216
				} else {
162
					$cadran = 'D';
217
					$cadran = 'D';
163
				}
218
				}
164
		} else {
219
		} else {
165
			if ($point['lat'] > $latCentre) {
220
			if ($point['lat'] > $latCentre) {
166
				$cadran = 'B';
221
				$cadran = 'B';
167
			} else {
222
			} else {
168
				$cadran = 'C';
223
				$cadran = 'C';
169
			}
224
			}
170
		}		
225
		}		
171
		return $cadran;
226
		return $cadran;
172
	}
227
	}
173
	
228
	
174
	private static function ajouterFils(&$noeud, $cadran, &$point) {
229
	private static function ajouterFils(&$noeud, $cadran, &$point) {
175
		if(!isset($noeud[$cadran])) {
230
		if(!isset($noeud[$cadran])) {
176
			$noeud[$cadran] = array('points' => array(),'nbrePoints' => 0, 'latMoyenne' => 0, 'lngMoyenne' => 0);
231
			$noeud[$cadran] = array('points' => array(),'nbrePoints' => 0, 'latMoyenne' => 0, 'lngMoyenne' => 0);
177
		}	
232
		}	
178
		$noeud[$cadran]['points'][] = $point;
233
		$noeud[$cadran]['points'][] = $point;
179
		$noeud[$cadran]['nbrePoints']++;
234
		$noeud[$cadran]['nbrePoints']++;
180
		$noeud[$cadran]['latMoyenne'] += $point['lat'];
235
		$noeud[$cadran]['latMoyenne'] += $point['lat'];
181
		$noeud[$cadran]['lngMoyenne'] += $point['lng'];
236
		$noeud[$cadran]['lngMoyenne'] += $point['lng'];
182
	}
237
	}
183
	
238
	
184
	private static function ajouterGroupeOuPoint(&$noeud) {
239
	private static function ajouterGroupeOuPoint(&$noeud) {
185
		$groupe = array();
240
		$groupe = array();
186
		if ($noeud['nbrePoints'] > 1) {
241
		if ($noeud['nbrePoints'] > 1 && isset($noeud['latMoyenne']) && isset($noeud['lngMoyenne'])) {
187
			$groupe['lat'] = $noeud['latMoyenne']/$noeud['nbrePoints'];
242
			$groupe['lat'] = $noeud['latMoyenne']/$noeud['nbrePoints'];
188
			$groupe['lng'] = $noeud['lngMoyenne']/$noeud['nbrePoints'];
243
			$groupe['lng'] = $noeud['lngMoyenne']/$noeud['nbrePoints'];
189
			$groupe['id'] = 'GROUPE:'.$groupe['lat'].';'.$groupe['lng'];
244
			$groupe['id'] = 'GROUPE:'.$groupe['lat'].';'.$groupe['lng'];
190
			$groupe['nbreMarqueur'] = $noeud['nbrePoints'];
245
			$groupe['nbreMarqueur'] = $noeud['nbrePoints'];
191
		} else {
246
		} else {
192
			$groupe = $noeud['points'][0];
247
			$groupe = $noeud['points'][0];
193
		}
248
		}
194
		return $groupe;
249
		return $groupe;
195
	}
250
	}
196
	
251
	
197
	private static function estUnParentFeuilles(&$noeud) {
252
	private static function estUnParentFeuilles(&$noeud) {
198
		return  self::estUneFeuille($noeud['A']) && 
253
		return  self::estUneFeuille($noeud['A']) && 
199
				self::estUneFeuille($noeud['B']) && 
254
				self::estUneFeuille($noeud['B']) && 
200
				self::estUneFeuille($noeud['C']) && 
255
				self::estUneFeuille($noeud['C']) && 
201
				self::estUneFeuille($noeud['D']);
256
				self::estUneFeuille($noeud['D']);
202
	}
257
	}
203
	
258
	
204
	private static function estUneFeuille(&$noeud) {		
259
	private static function estUneFeuille(&$noeud) {		
205
		return $noeud == null || 
260
		return $noeud == null || 
206
		(!isset($noeud['A']) || $noeud['A'] == null) && 
261
		(!isset($noeud['A']) || $noeud['A'] == null) && 
207
		(!isset($noeud['B']) || $noeud['B'] == null) && 
262
		(!isset($noeud['B']) || $noeud['B'] == null) && 
208
		(!isset($noeud['C']) || $noeud['C'] == null) && 
263
		(!isset($noeud['C']) || $noeud['C'] == null) && 
209
		(!isset($noeud['D']) || $noeud['D'] == null);
264
		(!isset($noeud['D']) || $noeud['D'] == null);
210
	}
265
	}
211
}
266
}
212
?>
267
?>