977 |
jpm |
1 |
<?php
|
|
|
2 |
class CartoGroupage {
|
978 |
jpm |
3 |
|
|
|
4 |
private static $seuilClusterisation = 100;
|
|
|
5 |
private static $zoomDefaut = 3;
|
|
|
6 |
private static $zoomMaxClustering = 12;
|
|
|
7 |
private static $pasZoomDefaut = 1;
|
|
|
8 |
private static $pasZoomMaxClustering = 0.05;
|
|
|
9 |
private static $profondeurMin = 0;
|
|
|
10 |
private static $profondeurMax = 8;
|
1003 |
jpm |
11 |
|
978 |
jpm |
12 |
private static $pasCorrectionCentre = null;
|
1003 |
jpm |
13 |
private static $coefficientReductionPas = null;
|
978 |
jpm |
14 |
private static $coefficientProfondeurMax = null;
|
1003 |
jpm |
15 |
|
979 |
jpm |
16 |
private static $nbElements = array('stations' => 0,'communes' => 0, 'observations' => 0);
|
978 |
jpm |
17 |
|
977 |
jpm |
18 |
private static $listeNoeudsSelectionnes = array();
|
1003 |
jpm |
19 |
|
979 |
jpm |
20 |
private static $pointsDejaTraites = array();
|
978 |
jpm |
21 |
|
|
|
22 |
/*
|
|
|
23 |
+---------+---------+
|
|
|
24 |
| | |
|
|
|
25 |
| A | B |
|
|
|
26 |
| | |
|
|
|
27 |
+---------*---------+
|
|
|
28 |
| | |
|
|
|
29 |
| D | C |
|
|
|
30 |
| | |
|
|
|
31 |
+---------+---------+
|
1003 |
jpm |
32 |
|
978 |
jpm |
33 |
Quatres cadrans sont considérés par le quad tree
|
1003 |
jpm |
34 |
* = centre de la fenetre
|
978 |
jpm |
35 |
*/
|
979 |
jpm |
36 |
public static function creerGroupesQuadtree(&$markers, $neLat, $neLng, $swLat, $swLng, $zoom = 3) {
|
1003 |
jpm |
37 |
|
978 |
jpm |
38 |
if(count($markers) > self::$seuilClusterisation) {
|
1003 |
jpm |
39 |
|
978 |
jpm |
40 |
self::calculerProfondeurMax($zoom);
|
|
|
41 |
self::calculerPasCorrectionCentre($zoom);
|
|
|
42 |
|
|
|
43 |
$noeudRacine = array('nbrePoints' => count($markers), 'points' => $markers);
|
1003 |
jpm |
44 |
self::attribuerAuCadran($noeudRacine, $neLat, $neLng, $swLat, $swLng);
|
|
|
45 |
|
978 |
jpm |
46 |
} else {
|
|
|
47 |
foreach($markers as $marker) {
|
|
|
48 |
$points = array($marker);
|
|
|
49 |
$noeudSimple = array('points' => $points, 'nbrePoints' => 1);
|
|
|
50 |
self::$listeNoeudsSelectionnes[] = self::ajouterGroupeOuPoint($noeudSimple);
|
977 |
jpm |
51 |
}
|
978 |
jpm |
52 |
}
|
1003 |
jpm |
53 |
|
978 |
jpm |
54 |
return self::$listeNoeudsSelectionnes;
|
|
|
55 |
}
|
1003 |
jpm |
56 |
|
978 |
jpm |
57 |
private function calculerCoefficientReductionPas() {
|
|
|
58 |
if(self::$coefficientReductionPas == null) {
|
|
|
59 |
self::$coefficientReductionPas = (self::$pasZoomMaxClustering - self::$pasZoomDefaut)/(self::$zoomMaxClustering - self::$zoomDefaut);
|
977 |
jpm |
60 |
}
|
1003 |
jpm |
61 |
|
978 |
jpm |
62 |
return self::$coefficientReductionPas;
|
977 |
jpm |
63 |
}
|
1003 |
jpm |
64 |
|
978 |
jpm |
65 |
private function calculerPasCorrectionCentre($zoom) {
|
|
|
66 |
self::$pasCorrectionCentre = ($zoom - self::$zoomDefaut) * self::calculerCoefficientReductionPas() + self::$pasZoomDefaut;
|
977 |
jpm |
67 |
}
|
1003 |
jpm |
68 |
|
978 |
jpm |
69 |
private function calculerCoefficientProfondeurMax() {
|
|
|
70 |
if(self::$coefficientProfondeurMax == null) {
|
|
|
71 |
self::$coefficientProfondeurMax = (self::$profondeurMax - self::$profondeurMin)/(self::$zoomMaxClustering - self::$zoomDefaut);
|
977 |
jpm |
72 |
}
|
1003 |
jpm |
73 |
|
978 |
jpm |
74 |
return self::$coefficientProfondeurMax;
|
977 |
jpm |
75 |
}
|
1003 |
jpm |
76 |
|
|
|
77 |
private function calculerProfondeurMax($zoom) {
|
|
|
78 |
if($zoom > self::$zoomDefaut) {
|
|
|
79 |
//self::$profondeurMax = round(($zoom - self::$zoomDefaut) * self::calculerCoefficientProfondeurMax() + self::$profondeurMin,0);
|
|
|
80 |
self::$profondeurMax = $zoom-1;
|
977 |
jpm |
81 |
} else {
|
978 |
jpm |
82 |
self::$profondeurMax = 1;
|
977 |
jpm |
83 |
}
|
|
|
84 |
}
|
1003 |
jpm |
85 |
|
979 |
jpm |
86 |
public static function getNbElements() {
|
|
|
87 |
return self::$nbElements;
|
|
|
88 |
}
|
1003 |
jpm |
89 |
|
978 |
jpm |
90 |
/**
|
1003 |
jpm |
91 |
*
|
978 |
jpm |
92 |
* @param mixed $noeud Le noeud à traiter par le quadtree
|
1003 |
jpm |
93 |
* @param float $neLat Latitude du coin nord est de la fenetre
|
|
|
94 |
* @param float $neLng Longitude du coin nord est de la fenetre
|
|
|
95 |
* @param float $swLat Latitude du coin sud ouest de la fenetre
|
|
|
96 |
* @param float $swLng Longitude du coin sud ouest de la fenetre
|
978 |
jpm |
97 |
* @param int $profondeur profondeur courante de l'arbre
|
|
|
98 |
*/
|
|
|
99 |
private static function attribuerAuCadran(&$noeud, $neLat, $neLng, $swLat, $swLng, $profondeur = 0) {
|
|
|
100 |
$latCentre = round((($neLat+$swLat)/2)/self::$pasCorrectionCentre,0)*self::$pasCorrectionCentre;
|
|
|
101 |
$lngCentre = round((($neLng+$swLng)/2)/self::$pasCorrectionCentre,0)*self::$pasCorrectionCentre;
|
1003 |
jpm |
102 |
|
977 |
jpm |
103 |
foreach ($noeud['points'] as &$point) {
|
979 |
jpm |
104 |
self::$nbElements['observations']++;
|
1003 |
jpm |
105 |
if (isset($point['type_emplacement'])) {
|
|
|
106 |
self::$nbElements[$point['type_emplacement']]++;
|
|
|
107 |
}
|
979 |
jpm |
108 |
unset($point['type_emplacement']);
|
|
|
109 |
$cadran = self::obtenirCadranPourPoint($latCentre, $lngCentre, $point);
|
|
|
110 |
self::ajouterFils($noeud,$cadran,$point);
|
977 |
jpm |
111 |
}
|
1003 |
jpm |
112 |
|
977 |
jpm |
113 |
$profondeur++;
|
1003 |
jpm |
114 |
|
|
|
115 |
if ($profondeur <= self::$profondeurMax) {
|
|
|
116 |
(isset($noeud['A'])) ? self::attribuerAuCadran($noeud['A'], $neLat, $lngCentre , $latCentre, $swLng, $profondeur) : '';
|
|
|
117 |
(isset($noeud['B'])) ? self::attribuerAuCadran($noeud['B'], $neLat, $neLng, $latCentre, $lngCentre, $profondeur) : '';
|
|
|
118 |
(isset($noeud['C'])) ? self::attribuerAuCadran($noeud['C'], $latCentre, $neLng, $swLat, $lngCentre, $profondeur) : '';
|
|
|
119 |
(isset($noeud['D'])) ? self::attribuerAuCadran($noeud['D'], $latCentre, $lngCentre, $swLat, $swLng, $profondeur) : '';
|
977 |
jpm |
120 |
}
|
1003 |
jpm |
121 |
|
978 |
jpm |
122 |
if(self::estUnParentFeuilles($noeud)) {
|
|
|
123 |
self::$listeNoeudsSelectionnes[] = self::ajouterGroupeOuPoint($noeud);
|
977 |
jpm |
124 |
}
|
|
|
125 |
}
|
1003 |
jpm |
126 |
|
978 |
jpm |
127 |
private function obtenirCadranPourPoint($latCentre,$lngCentre, &$point) {
|
|
|
128 |
if ($point['lng'] < $lngCentre) {
|
|
|
129 |
if ($point['lat'] > $latCentre) {
|
|
|
130 |
$cadran = 'A';
|
|
|
131 |
} else {
|
|
|
132 |
$cadran = 'D';
|
|
|
133 |
}
|
|
|
134 |
} else {
|
|
|
135 |
if ($point['lat'] > $latCentre) {
|
|
|
136 |
$cadran = 'B';
|
|
|
137 |
} else {
|
|
|
138 |
$cadran = 'C';
|
|
|
139 |
}
|
1003 |
jpm |
140 |
}
|
978 |
jpm |
141 |
return $cadran;
|
|
|
142 |
}
|
1003 |
jpm |
143 |
|
|
|
144 |
private static function ajouterFils(&$noeud, $cadran, &$point) {
|
|
|
145 |
$fils = array('points' => array(), 'nbrePoints' => 1, 'latMoyenne' => 0, 'lngMoyenne' => 0);
|
|
|
146 |
$fils['points'][] = $point;
|
|
|
147 |
$fils['nbrePoints'] = isset($noeud[$cadran]['nbrePoints']) ? $noeud[$cadran]['nbrePoints']++ : 1;
|
|
|
148 |
$fils['latMoyenne'] += $point['lat'];
|
|
|
149 |
$fils['lngMoyenne'] += $point['lng'];
|
|
|
150 |
$noeud[$cadran] = $fils;
|
978 |
jpm |
151 |
}
|
1003 |
jpm |
152 |
|
979 |
jpm |
153 |
private static function ajouterGroupeOuPoint(&$noeud) {
|
977 |
jpm |
154 |
$groupe = array();
|
|
|
155 |
if ($noeud['nbrePoints'] > 1) {
|
|
|
156 |
$groupe['lat'] = $noeud['latMoyenne']/$noeud['nbrePoints'];
|
|
|
157 |
$groupe['lng'] = $noeud['lngMoyenne']/$noeud['nbrePoints'];
|
978 |
jpm |
158 |
$groupe['id'] = 'GROUPE:'.$groupe['lat'].';'.$groupe['lng'];
|
977 |
jpm |
159 |
$groupe['nbreMarqueur'] = $noeud['nbrePoints'];
|
|
|
160 |
} else {
|
|
|
161 |
$groupe = $noeud['points'][0];
|
|
|
162 |
}
|
|
|
163 |
return $groupe;
|
|
|
164 |
}
|
1003 |
jpm |
165 |
|
979 |
jpm |
166 |
private static function estUnParentFeuilles(&$noeud) {
|
1003 |
jpm |
167 |
return self::estUneFeuille($noeud['A']) &&
|
|
|
168 |
self::estUneFeuille($noeud['B']) &&
|
|
|
169 |
self::estUneFeuille($noeud['C']) &&
|
978 |
jpm |
170 |
self::estUneFeuille($noeud['D']);
|
|
|
171 |
}
|
1003 |
jpm |
172 |
|
|
|
173 |
private static function estUneFeuille(&$noeud) {
|
|
|
174 |
return $noeud == null || (!isset($noeud['A']) && !isset($noeud['B']) && !isset($noeud['C']) && !isset($noeud['D']));
|
977 |
jpm |
175 |
}
|
|
|
176 |
}
|
|
|
177 |
?>
|