Subversion Repositories eFlore/Projets.eflore-projets

Rev

Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
736 alex 1
<?php
2
 
743 alex 3
/**
4
 * Classe unique du web service nasa-srtm qui détermine l'altitude d'un point à partir de ses coordonnees
5
 * Les coordonnees sont definies dans des fichiers au format HGT dans un dossier specifique
6
 * (dont le chemin est defini dans le fichier de configuration propre au web service).
7
 * Les ressources utilisees sont les donnees issues du programme SRTM-3 de la NASA qui couvre
8
 * l'ensemble terrestre du monde. La precision des points dont on renvoie l'altitude est de 90 metres.
9
 * Chaque fichier couvre un secteur de 1 degre sur 1 degre et contient un tableau de 1201 lignes
10
 * (axe des latitudes) sur 1201 colonnes (axe des longitudes) contenant l'altitude en metres
11
 * correspondant a des point precis. L'ecart entre chaque entree des tableaux est constant, ce qui
12
 * permet de calculer la latitude et la longitude de chaque point. L'altitude du point le plus proche
13
 * de celui passe en parametres sera renvoyee au client.
14
 *
15
 * Parametres du web service :
16
 *   - latitude : latitude du point dont on recherche les coordonnees
17
 *   - longitude : longitude du point dont on recherche les coordonnees
18
 *
19
 * @package framework-0.4
20
 * @author Alexandre GALIBERT <alexandre.galibert@tela-botanica.org>
21
 * @license GPL v3 <http://www.gnu.org/licenses/gpl.txt>
22
 * @license CECILL v2 <http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt>
23
 * @version $Id$
24
 * @copyright 2013 Tela Botanica (accueil@tela-botanica.org)
25
 *
26
 */
27
 
28
 
736 alex 29
class Altitude {
30
 
31
	private $parametres = array();
32
	private $ressources = array();
33
	private $coordonnees = null;
34
	private $fichierSrtm = '';
35
 
36
	const LONGUEUR_COTE = 1201;
37
	const OFFSET = 2;
38
 
39
 
40
	public function consulter($ressources, $parametres) {
41
		$this->ressources = $ressources;
42
		$this->parametres = $parametres;
43
		$retour = null;
44
		try {
45
			$this->traiterCoordonnees();
46
			$this->rechercherFichierSrtm();
47
			$this->recupererAltitude();
48
			$retour = $this->coordonnees;
49
		} catch (Exception $erreur) {
50
			$retour = $erreur->getMessage();
51
		}
52
		return $retour;
53
	}
54
 
55
	private function traiterCoordonnees() {
743 alex 56
		if ($this->estParametreExistant('lat') && $this->estParametreExistant('lon')) {
57
			$longitude = $this->parametres['lon'];
58
			$latitude  = $this->parametres['lat'];
736 alex 59
			if ($this->estUnFloat($longitude) && $this->estUnFloat($latitude)) {
60
				$this->verifierValiditeCoordonnees($longitude, $latitude);
61
			} else {
62
				$message = "La valeur des coordonnées longitude ou latitude n'est pas correcte. ".
63
				" Elle doit être pour les deux paramètres une valeur décimale.";
64
				throw new Exception($message, RestServeur::HTTP_CODE_MAUVAISE_REQUETE);
65
			}
66
		} else {
67
			$message = "Tous les paramètres passés dans l'URL ne correspondent pas à ceux attendus. ".
68
			"Le web service nécessite de lui fournir une longitude et une latitude pour son bon fonctionnement.";
69
			throw new Exception($message, RestServeur::HTTP_CODE_CONTENU_REQUIS);
70
		}
71
	}
72
 
73
	private function estParametreExistant($nomParametre) {
74
		return in_array($nomParametre, array_keys($this->parametres));
75
	}
76
 
77
	private function estUnFloat($variable) {
78
		return (preg_match("/^(-)?\d+(\.\d+)?$/", $variable) == 1);
79
	}
80
 
81
	private function verifierValiditeCoordonnees($longitude, $latitude) {
82
		$longitude = floatval($longitude);
83
		$latitude  = floatval($latitude);
84
		$longitudeMax = Config::get("limite_longitude");
85
		$latitudeMax  = Config::get("limite_latitude");
86
		if (abs($longitude) > $longitudeMax || abs($latitude) > $latitudeMax) {
87
			$message = "Les coordonnées passées en paramètres désignent un point qui se trouve ".
88
			"en dehors des limites du monde. Elles doivent être comprises entre -{$longitudeMax} ".
89
			"et $longitudeMax sur l'axe des longitudes, et entre -{$latitudeMax} et {$latitudeMax} ".
90
			"sur l'axe des latitudes.";
91
			throw new Exception($message, RestServeur::HTTP_CODE_MAUVAISE_REQUETE);
92
		} else {
93
			$this->coordonnees = new StdClass();
94
			$this->coordonnees->longitude = $longitude;
95
			$this->coordonnees->latitude  = $latitude;
96
		}
97
	}
98
 
99
	private function rechercherFichierSrtm() {
100
		$nomFichierSrtm = $this->construireNomFichierSrtm();
101
		if (!file_exists($nomFichierSrtm)) {
102
			$message = "Erreur interne : certaines ressources demandées n'ont pas pu être trouvées sur le serveur.";
103
			throw new Exception($message, restServeur::HTTP_CODE_ERREUR);
104
		} else {
105
			$this->fichierSrtm = $nomFichierSrtm;
106
		}
107
	}
108
 
109
	private function construireNomFichierSrtm() {
110
		$latitudeEntier = abs(floor($this->coordonnees->latitude));
111
		if ($latitudeEntier < 10) {
112
			$latitudeEntier = "0".$latitudeEntier;
113
		}
114
		$suffixeLatitude = $this->coordonnees->latitude < 0 ? "S" : "N";
115
		$longitudeEntier = abs(floor($this->coordonnees->longitude));
116
		if ($longitudeEntier < 10) {
117
			$longitudeEntier = "00".$longitudeEntier;
118
		} elseif ($longitudeEntier < 100) {
119
			$longitudeEntier = "0".$longitudeEntier;
120
		}
121
		$suffixeLongitude = $this->coordonnees->longitude < 0 ? "W" : "E";
122
		$dossierSrtm = Config::get('dossier_srtm').DS;
123
		$nomFichierSrtm = $dossierSrtm.$suffixeLatitude.$latitudeEntier.$suffixeLongitude.$longitudeEntier.".hgt.zip";
124
		return $nomFichierSrtm;
125
	}
126
 
127
	private function recupererAltitude() {
128
		$zip = zip_open($this->fichierSrtm);
129
		$fichier = zip_read($zip);
130
		$donnees = zip_entry_read($fichier, zip_entry_filesize($fichier));
131
		zip_close($zip);
132
 
133
		$xDepart = floor($this->coordonnees->longitude);
134
		$yDepart = floor($this->coordonnees->latitude);
135
		$longitude = $this->coordonnees->longitude;
136
		$latitude = $this->coordonnees->latitude;
137
		$positionX = (self::LONGUEUR_COTE-1) * ($longitude - $xDepart);
138
		$positionY = (self::LONGUEUR_COTE-1) * (1 - $latitude + $yDepart);
139
		$positionX = ($positionX + 0.5 > ceil($positionX)) ? ceil($positionX) : floor($positionX);
140
		$positionY = ($positionY + 0.5 > ceil($positionY)) ? ceil($positionY) : floor($positionY);
141
 
142
		$binaire = substr($donnees, ($positionY * self::LONGUEUR_COTE + $positionX) * self::OFFSET, self::OFFSET);
143
		$this->coordonnees->altitude = current(unpack("n*", $binaire));
144
		if (!$this->coordonnees->altitude) {
145
			$this->coordonnees->altitude = 0;
146
		}
147
	}
148
 
149
}
150
 
151
?>