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