Subversion Repositories eFlore/Projets.eflore-projets

Rev

Rev 743 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 743 Rev 1091
Line 1... Line 1...
1
<?php
1
<?php
2
 
-
 
-
 
2
// declare(encoding='UTF-8');
3
/**
3
/**
-
 
4
 * @api {get} /nasa-srtm/altitude?lat=:lat&lon=:lon Donne l'altitude.
-
 
5
 * @apiExample Exemple d'usage :
-
 
6
 * curl -i "http://api.tela-botanica.org/service:eflore:0.1/nasa-srtm/altitude?lon=3.86589&lat=44.73702"
-
 
7
 * @apiName GetNasaSrtm
-
 
8
 * @apiGroup nasa-srtm
4
 * Classe unique du web service nasa-srtm qui détermine l'altitude d'un point à partir de ses coordonnees
9
 * @apiGroupDescription  Le web service nasa-srtm/altitude 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
10
 * 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).
11
 * (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
12
 * 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.
13
 * 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
14
 * 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
15
 * (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
16
 * 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
17
 * 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.
18
 * 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
 *
19
 *
-
 
20
 * @apiParam {Number} lat Latitude au format décimal (séparateur ".").
-
 
21
 * @apiParam {Number} lon Longitude au format décimal (séparateur ".").
-
 
22
 *
-
 
23
 * @apiSuccess {Number} altitude Altitude en mètre du point correspondant aux coordonées passées en paramètre.
-
 
24
 * @apiSuccess {Number} latitude Latitude au format décimal (séparateur ".").
-
 
25
 * @apiSuccess {Number} longitude Longitude au format décimal (séparateur ".").
-
 
26
 *
-
 
27
 * @category   eFlore
-
 
28
 * @package    Services
-
 
29
 * @subpackage Nasa-srtm
-
 
30
 * @version    0.1
-
 
31
 * @author     Mathias CHOUET <mathias@tela-botanica.org>
-
 
32
 * @author     Jean-Pascal MILCENT <jpm@tela-botanica.org>
-
 
33
 * @author     Aurelien PERONNET <aurelien@tela-botanica.org>
-
 
34
 * @license    GPL v3 <http://www.gnu.org/licenses/gpl.txt>
-
 
35
 * @license    CECILL v2 <http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt>
-
 
36
 * @copyright  1999-2014 Tela Botanica <accueil@tela-botanica.org>
26
 */
37
 */
27
 
-
 
28
 
-
 
29
class Altitude {
38
class Altitude {
30
	
39
 
31
	private $parametres = array();
40
	private $parametres = array();
32
	private $ressources = array();
41
	private $ressources = array();
33
	private $coordonnees = null;
42
	private $coordonnees = null;
34
	private $fichierSrtm = '';
43
	private $fichierSrtm = '';
35
	
44
 
36
	const LONGUEUR_COTE = 1201;
45
	const LONGUEUR_COTE = 1201;
37
	const OFFSET = 2;
46
	const OFFSET = 2;
38
	
47
 
39
	
48
 
40
	public function consulter($ressources, $parametres) {
49
	public function consulter($ressources, $parametres) {
41
		$this->ressources = $ressources;
50
		$this->ressources = $ressources;
42
		$this->parametres = $parametres;
51
		$this->parametres = $parametres;
43
		$retour = null;
52
		$retour = null;
44
		try {
53
		try {
Line 49... Line 58...
49
		} catch (Exception $erreur) {
58
		} catch (Exception $erreur) {
50
			$retour = $erreur->getMessage();
59
			$retour = $erreur->getMessage();
51
		}
60
		}
52
		return $retour;
61
		return $retour;
53
	}
62
	}
54
	
63
 
55
	private function traiterCoordonnees() {
64
	private function traiterCoordonnees() {
56
		if ($this->estParametreExistant('lat') && $this->estParametreExistant('lon')) {
65
		if ($this->estParametreExistant('lat') && $this->estParametreExistant('lon')) {
57
			$longitude = $this->parametres['lon'];
66
			$longitude = $this->parametres['lon'];
58
			$latitude  = $this->parametres['lat'];
67
			$latitude  = $this->parametres['lat'];
59
			if ($this->estUnFloat($longitude) && $this->estUnFloat($latitude)) {
68
			if ($this->estUnFloat($longitude) && $this->estUnFloat($latitude)) {
Line 67... Line 76...
67
			$message = "Tous les paramètres passés dans l'URL ne correspondent pas à ceux attendus. ".
76
			$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.";
77
			"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);
78
			throw new Exception($message, RestServeur::HTTP_CODE_CONTENU_REQUIS);
70
		}
79
		}
71
	}
80
	}
72
	
81
 
73
	private function estParametreExistant($nomParametre) {
82
	private function estParametreExistant($nomParametre) {
74
		return in_array($nomParametre, array_keys($this->parametres));
83
		return in_array($nomParametre, array_keys($this->parametres));
75
	}
84
	}
76
	
85
 
77
	private function estUnFloat($variable) {
86
	private function estUnFloat($variable) {
78
		return (preg_match("/^(-)?\d+(\.\d+)?$/", $variable) == 1);
87
		return (preg_match("/^(-)?\d+(\.\d+)?$/", $variable) == 1);
79
	}
88
	}
80
	
89
 
81
	private function verifierValiditeCoordonnees($longitude, $latitude) {
90
	private function verifierValiditeCoordonnees($longitude, $latitude) {
82
		$longitude = floatval($longitude);
91
		$longitude = floatval($longitude);
83
		$latitude  = floatval($latitude);
92
		$latitude  = floatval($latitude);
84
		$longitudeMax = Config::get("limite_longitude");
93
		$longitudeMax = Config::get("limite_longitude");
85
		$latitudeMax  = Config::get("limite_latitude");
94
		$latitudeMax  = Config::get("limite_latitude");
Line 93... Line 102...
93
			$this->coordonnees = new StdClass();
102
			$this->coordonnees = new StdClass();
94
			$this->coordonnees->longitude = $longitude;
103
			$this->coordonnees->longitude = $longitude;
95
			$this->coordonnees->latitude  = $latitude;
104
			$this->coordonnees->latitude  = $latitude;
96
		}
105
		}
97
	}
106
	}
98
	
107
 
99
	private function rechercherFichierSrtm() {
108
	private function rechercherFichierSrtm() {
100
		$nomFichierSrtm = $this->construireNomFichierSrtm();
109
		$nomFichierSrtm = $this->construireNomFichierSrtm();
101
		if (!file_exists($nomFichierSrtm)) {
110
		if (!file_exists($nomFichierSrtm)) {
102
			$message = "Erreur interne : certaines ressources demandées n'ont pas pu être trouvées sur le serveur.";
111
			$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);
112
			throw new Exception($message, restServeur::HTTP_CODE_ERREUR);
104
		} else {
113
		} else {
105
			$this->fichierSrtm = $nomFichierSrtm;
114
			$this->fichierSrtm = $nomFichierSrtm;
106
		}
115
		}
107
	}
116
	}
108
	
117
 
109
	private function construireNomFichierSrtm() {
118
	private function construireNomFichierSrtm() {
110
		$latitudeEntier = abs(floor($this->coordonnees->latitude));
119
		$latitudeEntier = abs(floor($this->coordonnees->latitude));
111
		if ($latitudeEntier < 10) {
120
		if ($latitudeEntier < 10) {
112
			$latitudeEntier = "0".$latitudeEntier;
121
			$latitudeEntier = "0".$latitudeEntier;
113
		}
122
		}
Line 121... Line 130...
121
		$suffixeLongitude = $this->coordonnees->longitude < 0 ? "W" : "E";
130
		$suffixeLongitude = $this->coordonnees->longitude < 0 ? "W" : "E";
122
		$dossierSrtm = Config::get('dossier_srtm').DS;
131
		$dossierSrtm = Config::get('dossier_srtm').DS;
123
		$nomFichierSrtm = $dossierSrtm.$suffixeLatitude.$latitudeEntier.$suffixeLongitude.$longitudeEntier.".hgt.zip";
132
		$nomFichierSrtm = $dossierSrtm.$suffixeLatitude.$latitudeEntier.$suffixeLongitude.$longitudeEntier.".hgt.zip";
124
		return $nomFichierSrtm;
133
		return $nomFichierSrtm;
125
	}
134
	}
126
	
135
 
127
	private function recupererAltitude() {
136
	private function recupererAltitude() {
128
		$zip = zip_open($this->fichierSrtm);
137
		$zip = zip_open($this->fichierSrtm);
129
		$fichier = zip_read($zip);
138
		$fichier = zip_read($zip);
130
		$donnees = zip_entry_read($fichier, zip_entry_filesize($fichier));
139
		$donnees = zip_entry_read($fichier, zip_entry_filesize($fichier));
131
		zip_close($zip);
140
		zip_close($zip);
132
				
141
 
133
		$xDepart = floor($this->coordonnees->longitude);
142
		$xDepart = floor($this->coordonnees->longitude);
134
		$yDepart = floor($this->coordonnees->latitude);
143
		$yDepart = floor($this->coordonnees->latitude);
135
		$longitude = $this->coordonnees->longitude;
144
		$longitude = $this->coordonnees->longitude;
136
		$latitude = $this->coordonnees->latitude;
145
		$latitude = $this->coordonnees->latitude;
137
		$positionX = (self::LONGUEUR_COTE-1) * ($longitude - $xDepart);
146
		$positionX = (self::LONGUEUR_COTE-1) * ($longitude - $xDepart);
138
		$positionY = (self::LONGUEUR_COTE-1) * (1 - $latitude + $yDepart);
147
		$positionY = (self::LONGUEUR_COTE-1) * (1 - $latitude + $yDepart);
139
		$positionX = ($positionX + 0.5 > ceil($positionX)) ? ceil($positionX) : floor($positionX);
148
		$positionX = ($positionX + 0.5 > ceil($positionX)) ? ceil($positionX) : floor($positionX);
140
		$positionY = ($positionY + 0.5 > ceil($positionY)) ? ceil($positionY) : floor($positionY);
149
		$positionY = ($positionY + 0.5 > ceil($positionY)) ? ceil($positionY) : floor($positionY);
141
		
150
 
142
		$binaire = substr($donnees, ($positionY * self::LONGUEUR_COTE + $positionX) * self::OFFSET, self::OFFSET);
151
		$binaire = substr($donnees, ($positionY * self::LONGUEUR_COTE + $positionX) * self::OFFSET, self::OFFSET);
143
		$this->coordonnees->altitude = current(unpack("n*", $binaire));
152
		$this->coordonnees->altitude = current(unpack("n*", $binaire));
144
		if (!$this->coordonnees->altitude) {
153
		if (!$this->coordonnees->altitude) {
145
			$this->coordonnees->altitude = 0;
154
			$this->coordonnees->altitude = 0;
146
		}
155
		}
147
	}
156
	}
148
	
157
 
149
}
158
}
Line 150... Line 159...
150
 
159
 
151
?>
160
?>