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 |
?>
|