Subversion Repositories eFlore/Applications.cel

Rev

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

Rev Author Line No. Line
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 
?>