Subversion Repositories eFlore/Projets.eflore-projets

Rev

Rev 1091 | Go to most recent revision | Blame | Compare with Previous | Last modification | View Log | RSS feed

<?php
/**
 * @api {get} /nasa-srtm/altitude?lat=:lat&lon=:lon Donne l'altitude.
 * @apiExample Exemple d'usage :
 * curl -i "http://api.tela-botanica.org/service:eflore:0.1/nasa-srtm/altitude?lon=3.86589&lat=44.73702"
 * @apiName GetNasaSrtm
 * @apiGroup nasa-srtm
 * @apiGroupDescription  Le web service nasa-srtm/altitude détermine l'altitude d'un point à partir de ses coordonnees.
 * Les coordonnees sont definies dans des fichiers au format HGT dans un dossier specifique
 * (dont le chemin est defini dans le fichier de configuration propre au web service).
 * Les ressources utilisees sont les donnees issues du programme SRTM-3 de la NASA qui couvre
 * l'ensemble terrestre du monde. La precision des points dont on renvoie l'altitude est de 90 metres.
 * Chaque fichier couvre un secteur de 1 degre sur 1 degre et contient un tableau de 1201 lignes
 * (axe des latitudes) sur 1201 colonnes (axe des longitudes) contenant l'altitude en metres
 * correspondant a des point precis. L'ecart entre chaque entree des tableaux est constant, ce qui
 * permet de calculer la latitude et la longitude de chaque point. L'altitude du point le plus proche
 * de celui passe en parametres sera renvoyee au client.
 *
 * @apiParam {Number} lat Latitude au format décimal (séparateur ".").
 * @apiParam {Number} lon Longitude au format décimal (séparateur ".").
 *
 * @apiSuccess {Number} altitude Altitude en mètre du point correspondant aux coordonées passées en paramètre.
 * @apiSuccess {Number} latitude Latitude au format décimal (séparateur ".").
 * @apiSuccess {Number} longitude Longitude au format décimal (séparateur ".").
 *
 * @category   eFlore
 * @package    Services
 * @subpackage Nasa-srtm
 * @version    0.1
 * @author     Mathias CHOUET <mathias@tela-botanica.org>
 * @author     Jean-Pascal MILCENT <jpm@tela-botanica.org>
 * @author     Aurelien PERONNET <aurelien@tela-botanica.org>
 * @license    GPL v3 <http://www.gnu.org/licenses/gpl.txt>
 * @license    CECILL v2 <http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt>
 * @copyright  1999-2014 Tela Botanica <accueil@tela-botanica.org>
 */
class Altitude {

        private $parametres = array();
        private $ressources = array();
        private $coordonnees = null;
        private $fichierSrtm = '';

        const LONGUEUR_COTE = 1201;
        const OFFSET = 2;


        public function consulter($ressources, $parametres) {
                $this->ressources = $ressources;
                $this->parametres = $parametres;
                $retour = null;
                try {
                        $this->traiterCoordonnees();
                        $this->rechercherFichierSrtm();
                        $this->recupererAltitude();
                        $retour = $this->coordonnees;
                } catch (Exception $erreur) {
                        $retour = $erreur->getMessage();
                }
                return $retour;
        }

        private function traiterCoordonnees() {
                if ($this->estParametreExistant('lat') && $this->estParametreExistant('lon')) {
                        $longitude = $this->parametres['lon'];
                        $latitude  = $this->parametres['lat'];
                        if ($this->estUnFloat($longitude) && $this->estUnFloat($latitude)) {
                                $this->verifierValiditeCoordonnees($longitude, $latitude);
                        } else {
                                $message = "La valeur des coordonnées longitude ou latitude n'est pas correcte. ".
                                " Elle doit être pour les deux paramètres une valeur décimale.";
                                throw new Exception($message, RestServeur::HTTP_CODE_MAUVAISE_REQUETE);
                        }
                } else {
                        $message = "Tous les paramètres passés dans l'URL ne correspondent pas à ceux attendus. ".
                        "Le web service nécessite de lui fournir une longitude et une latitude pour son bon fonctionnement.";
                        throw new Exception($message, RestServeur::HTTP_CODE_CONTENU_REQUIS);
                }
        }

        private function estParametreExistant($nomParametre) {
                return in_array($nomParametre, array_keys($this->parametres));
        }

        private function estUnFloat($variable) {
                return (preg_match("/^(-)?\d+(\.\d+)?$/", $variable) == 1);
        }

        private function verifierValiditeCoordonnees($longitude, $latitude) {
                $longitude = floatval($longitude);
                $latitude  = floatval($latitude);
                $longitudeMax = Config::get("limite_longitude");
                $latitudeMax  = Config::get("limite_latitude");
                if (abs($longitude) > $longitudeMax || abs($latitude) > $latitudeMax) {
                        $message = "Les coordonnées passées en paramètres désignent un point qui se trouve ".
                        "en dehors des limites du monde. Elles doivent être comprises entre -{$longitudeMax} ".
                        "et $longitudeMax sur l'axe des longitudes, et entre -{$latitudeMax} et {$latitudeMax} ".
                        "sur l'axe des latitudes.";
                        throw new Exception($message, RestServeur::HTTP_CODE_MAUVAISE_REQUETE);
                } else {
                        $this->coordonnees = new StdClass();
                        $this->coordonnees->longitude = $longitude;
                        $this->coordonnees->latitude  = $latitude;
                }
        }

        private function rechercherFichierSrtm() {
                $nomFichierSrtm = $this->construireNomFichierSrtm();
                if (!file_exists($nomFichierSrtm)) {
                        $message = "Erreur interne : la ressource « $nomFichierSrtm » demandée n'a pas pu être trouvée sur le serveur.";
                        throw new Exception($message, restServeur::HTTP_CODE_ERREUR);
                } else {
                        $this->fichierSrtm = $nomFichierSrtm;
                }
        }

        private function construireNomFichierSrtm() {
                $latitudeEntier = abs(floor($this->coordonnees->latitude));
                if ($latitudeEntier < 10) {
                        $latitudeEntier = "0".$latitudeEntier;
                }
                $suffixeLatitude = $this->coordonnees->latitude < 0 ? "S" : "N";
                $longitudeEntier = abs(floor($this->coordonnees->longitude));
                if ($longitudeEntier < 10) {
                        $longitudeEntier = "00".$longitudeEntier;
                } elseif ($longitudeEntier < 100) {
                        $longitudeEntier = "0".$longitudeEntier;
                }
                $suffixeLongitude = $this->coordonnees->longitude < 0 ? "W" : "E";
                $dossierSrtm = Config::get('dossier_srtm').DS;
                $nomFichierSrtm = $dossierSrtm.$suffixeLatitude.$latitudeEntier.$suffixeLongitude.$longitudeEntier.".hgt.zip";
                return $nomFichierSrtm;
        }

        private function recupererAltitude() {
                $zip = zip_open($this->fichierSrtm);
                $fichier = zip_read($zip);
                $donnees = zip_entry_read($fichier, zip_entry_filesize($fichier));
                zip_close($zip);

                $xDepart = floor($this->coordonnees->longitude);
                $yDepart = floor($this->coordonnees->latitude);
                $longitude = $this->coordonnees->longitude;
                $latitude = $this->coordonnees->latitude;
                $positionX = (self::LONGUEUR_COTE-1) * ($longitude - $xDepart);
                $positionY = (self::LONGUEUR_COTE-1) * (1 - $latitude + $yDepart);
                $positionX = ($positionX + 0.5 > ceil($positionX)) ? ceil($positionX) : floor($positionX);
                $positionY = ($positionY + 0.5 > ceil($positionY)) ? ceil($positionY) : floor($positionY);

                $binaire = substr($donnees, ($positionY * self::LONGUEUR_COTE + $positionX) * self::OFFSET, self::OFFSET);
                $this->coordonnees->altitude = current(unpack("n*", $binaire));
                if (!$this->coordonnees->altitude) {
                        $this->coordonnees->altitude = 0;
                }
        }
}