Subversion Repositories Sites.obs-saisons.fr

<?


if (!isset($_POST['valid_station']) || sizeof($erreurs)>0) {

?>

Définissez vos stations d'observation :

<br>
<form id="addr" action="#" method="get" onsubmit="return showAddress(this);" style="float:left;">
                <ul style="list-style-type:none;margin:0;padding:0;">
                        <li>
                                <label for="lieu">Selectionnez la commune la plus proche de votre station :</label>
                                <input id="lieu" type="text" size="30" name="lieu" value=""/>
                                <input type="submit" value="Localiser" />
                        </li>
                        <li>
                                Ensuite, positionnez la pastille rouge sur votre station et donnez lui un nom.
                        </li
                </ul>
        </fieldset>
</form>

<div id="map" style="width:600px;height:400px;border:2px solid grey;margin:10px 0;clear:both;"></div>

<?php

  //Récupération des stations existantes

  //FONCTIONNEMENT_TESTE_ET_OK


  // les stations de l'observateur :

  $requete_stations = mysql_query("select  STATION_ID, STATION_NOM, STATION_LATITUDE,  STATION_LONGITUDE, STATION_ALTITUDE, STATION_ENVIRONNEMENT_ID, ENVIRONNEMENT_NOM from STATION, ENVIRONNEMENT  where STATION.PARTICIPANT_ID=".$_SESSION['participant']." AND ENVIRONNEMENT_ID=STATION_ENVIRONNEMENT_ID ORDER BY STATION_NOM");



  $nb_requete_stations = @mysql_num_rows($requete_stations);






  $stations=array();


  echo "<br>";
  echo "</br>\n";
  echo "<b>Vos stations enregistrées (Cliquez sur le nom de la station pour ajouter de nouvelles espèces) :</b> <br>";


  ?>
  <form name="form_station" method="post" action="<?= $url_page ?>">
  <input type="hidden" name="etape" value="<?= $etape ?>">
  <input type="hidden" name="<?= $provenance ?>" value="1">
  <input type="hidden" name="<?= $provenance ?>etape" value="<?= $etape ?>">
  <?php

  echo "<table id='tabstations'>";



  if ($nb_requete_stations > 0) {


  // Liste des stations

  $MOIS_EN_COURS = '09';
  $MOIS_RENTREE = '09'; // FIXME : quelle difference avec la variable precedente ?
  $mois = date('m');
  $annee = date('Y');

  //$annee=$annee+1;
 /*
  if ($mois > $MOIS_RENTREE) {
                $annee = $annee+1;
  }
  */

  while ($station = mysql_fetch_assoc($requete_stations)) {


          $stations [] = $station;
          echo "<tr>";
          echo "<td><strong>Station : <a title=\"Ajouter de nouvelles espèces\" href=\"".$url_page."&station_id=".$station['STATION_ID']."\">".$station['STATION_NOM']."</a><strong></td>";



         // echo "<td>".$station['STATION_LATITUDE']."</td>";
          //echo "<td>".$station['STATION_LONGITUDE']."</td>";
          //echo "<td>".$station['STATION_ALTITUDE']."</td>";
          //echo "<td>"."<input type=\"submit\" class=\"submit\" name=\"".$station['STATION_ID']."\" value=\"Choisir cette station\"></td>";

          // A t-on des mesures associe a cette station ? Si oui on ne peut pas la modifier ni supprimer

          $requete_mesure =  mysql_query("SELECT MESURE.SEQUENCE_ID FROM   MESURE,SEQUENCE,STATION where MESURE.SEQUENCE_ID=SEQUENCE.SEQUENCE_ID AND STATION.STATION_ID=SEQUENCE.STATION_ID AND  STATION.PARTICIPANT_ID=".$_SESSION['participant']." AND STATION.STATION_ID=".$station['STATION_ID']);
          $nb_requete_mesure = mysql_num_rows($requete_mesure);

      $obs_existe=false;

      if ($nb_requete_mesure  == 0) {  // Pas de mesure on peut modifier
          echo "<td>"."<input type=\"submit\" class=\"submit\"  name=\"".$station['STATION_ID']."\" value=\"Modifier\"></td>";
      }
      else {
           $obs_existe=true;
          echo "<td></td>";
      }

          echo "</tr>";




         $requete_obs_annee = mysql_query("select distinct SEQUENCE.SEQUENCE_ID, ESPECE_NOM_VERNACULAIRE, STATION_NOM from SEQUENCE,ESPECE,STATION,MESURE ".
        "where ESPECE.ESPECE_ID=SEQUENCE.ESPECE_ID".
        " and MESURE.SEQUENCE_ID=SEQUENCE.SEQUENCE_ID".
        " and SEQUENCE.PARTICIPANT_ID=".$_SESSION['participant'].
        " and STATION.STATION_ID=". $station['STATION_ID'] .
        " and STATION.STATION_ID=SEQUENCE.STATION_ID".
        " AND MESURE_DATE!='0000-00-00' and MESURE_DATE>'".($annee-2).'-'.($mois)."-01'".
        " and MESURE_DATE<'".($annee).'-'.($mois)."-01'".
        " order by STATION_NOM, ESPECE_NOM_VERNACULAIRE");

        $nb_requete_obs_annee = mysql_num_rows($requete_obs_annee);

        if (($nb_requete_obs_annee  == 0) && ($obs_existe)) {
                echo "<tr>";
                echo "<td>Observations archivées </td>";
                echo "</tr>";

        }

        while ($obs = mysql_fetch_assoc($requete_obs_annee)) {
                echo "<tr>";
                echo "<td>Espèce : <a title=\"Ajouter/modifier des observations\" href=\"".$url_page."&sequence_id=".$obs['SEQUENCE_ID']."\">".$obs['ESPECE_NOM_VERNACULAIRE']."</a></td>";
                echo "</tr>";
        }

         echo "<tr>";
         echo "</tr>";
         echo "<tr>";
         echo "</tr>";

  }



        // On retient la premiere station pour centrage

        $init_lat=$stations[0]['STATION_LATITUDE'];
        $init_lon=$stations[0]['STATION_LONGITUDE'];


  }
  else {

        // Pas de stations : affichage centre carte france
        // FONCTIONNEMENT_TESTE_ET_OK


        $init_lat=45.545;
        $init_lon=3.249722;

        echo "<td><strong>Pas de station enregistr&eacute;e</strong></td>";
        echo "<td></td>";
        echo "<td></td>";
        echo "<td></td>";
        echo "<td></td>";
        echo "<td></td>";
        echo "</tr>";



  }

  echo "</table>";
  echo "</form>";





//$espece_encours=$obs[1];
$commune_encours=$obs[2];



/***** Selection des stations *****/

//Affichage des éventuelles erreurs
afficherErreurs($erreurs);

?>




<form id="form_coordonnee" action="<?= $url_page ?>" method="post" style="float:left;">
                                <input id="ll_latitude" name="ll_latitude"  type="hidden" value="" />
                                <input id="ll_latitude_dms" name="ll_latitude_dms"  type="hidden" value="" />
                                <input id="ll_longitude" name="ll_longitude"  type="hidden" value="" />
                                <input id="ll_longitude_dms" name="ll_longitude_dms"  type="hidden" value="" />
                                <input id="ll_altitude" name="ll_altitude"  type="hidden" value="" />
                                <input id="form_coordonnee_valider" name="valid_station" class="submit" type="hidden" value="Poursuivre" />
                                <input type="hidden" name="etape" value="3">
                                <input type="hidden" name="provenance" value="<?= $provenance ?>">
                        </li>
                </ul>

</form>



<?php

include_once 'modules/configuration/car_config.inc.php' ;

if (PARTICIPANT_EST_ADULTE) {
?>
<script type="text/javascript" src="http://maps.google.com/maps?file=api&amp;v=2&amp;key=ABQIAAAAnm2MSMmwsDuoh7THcVDnlBQH9AL26AfEfaraxwo1G6_adyVjPRRQE4Gx2u7MomK9DJWf9kMLJ9y-ZA" type="text/javascript"></script>
<?php
}
else {
?>
<script type="text/javascript" src="http://maps.google.com/maps?file=api&amp;v=2&amp;key=ABQIAAAAFMANWib9GSm7T4W6AN9qzxRmHYqXwt_BUwczny2vFmnHz3sr3BTy_SMRJ-NzWXDUynHHSWQ7_aEZ6Q" type="text/javascript"></script>
<?php
}


// Centrage initial : premiere station ou ville d'inscription

   echo "<script type=\"text/javascript\">";
   echo "\n//<![CDATA[\n";
   echo "var init_lat = ".$init_lat.";\n";
   echo "var init_long = ".$init_lon.";\n";

?>


        var largeur_carte = 690;
        var haureur_retrait = 330;

        var init_zoom = 10;

// Contient la carte GoogleMap
var map;
// Contient l'utilitaire permettant de trouver une adresse.
var EfGeocoder;
// Point de l'adresse saisie et à passer au GeoCoder
var addrpnt;
// L'icône à afficher sur la carte
var icon;

function load() {
        // Vérification de la compatibilité du navigateur
        if (GBrowserIsCompatible()) {
                // Création de la carte
                map = new GMap2(document.getElementById("map"));
                // Positionnement de la carte
                map.setCenter(new GLatLng(init_lat, init_long), init_zoom);
                // Types de carte disponibles :
                map.addMapType(G_NORMAL_MAP);
                map.addMapType(G_SATELLITE_MAP);
                map.addMapType(G_PHYSICAL_MAP);
                // Bouton pour changer de type de carte (mixte, satellite, plan)
                map.addControl(new GMapTypeControl(1));
                // Définition du type de la carte par défaut
                map.setMapType(G_HYBRID_MAP);
                // Bouton large de déplancement et zoom apparaissant dans le coin haut gauche
                map.addControl(new GLargeMapControl());
                // Pour hierarchiser les types de carte --> pas encore dans la version stable...
                //var mapControl = new GHierarchicalMapTypeControl();
                // Définition des relations des types de menus
                //mapControl.clearRelationships();
                //mapControl.addRelationship(G_SATELLITE_MAP, G_HYBRID_MAP, 'Labels', false);
                // Ajout du controle après avoir spécifié les relations
                //map.addControl(mapControl);
                // Vue d'ensemble apparaissant dans le coin bas droit
                var EfOverViewMap = new GOverviewMapControl();
                map.addControl(EfOverViewMap);
                var mini = EfOverViewMap.getOverviewMap();
                // ?
                map.enableContinuousZoom();
                // Zoom avec la molette de la souris
                map.enableScrollWheelZoom();
                // Zoom avec le double clic de la souris
                map.enableDoubleClickZoom();

                // Création du GéoCoder
                EfGeocoder = new GClientGeocoder() ;

                // Icone station enregistrees (vert)
                icon = new GIcon();
                icon.image = "http://www.google.com/mapfiles/ms/micons/green-dot.png";
                icon.shadow = "http://www.google.com/mapfiles/shadow50.png";
        icon.iconSize = new GSize(34, 34);
        icon.shadowSize = new GSize(37, 34);
        icon.iconAnchor = new GPoint(16, 34);
        // Pour affichage infobulle
        icon.infoWindowAnchor = new GPoint(5,1);


                // Icone nouvelle  station (rouge) ou station en modification
                Nicon = new GIcon();
                Nicon.image = "http://www.google.com/mapfiles/ms/micons/red-dot.png";
                Nicon.shadow = "http://www.google.com/mapfiles/shadow50.png";
        Nicon.iconSize = new GSize(34, 34);
        Nicon.shadowSize = new GSize(37, 34);
        Nicon.iconAnchor = new GPoint(16, 34);
        // Pour affichage infobulle
        Nicon.infoWindowAnchor = new GPoint(5,1);




        // Creates a marker at the given point with the given number label

                function createMarker(point, text, altitude, envir) {
                  var marker = new GMarker(point,icon);
                  var html = text;
                  GEvent.addListener(marker, "click", function() {
                        map.setCenter(marker.getPoint());
                        afficherInfo(marker.getPoint());
                    afficherMasque(marker,text,altitude,envir);
                        });
                  return marker;
                }

        // Affichage des stations presentes TODO : et si plusieurs  stations  au meme endroit ? Utiliser cluster ?

        <? foreach($stations as $station) { ?>
          var point = new GLatLng(<?= $station['STATION_LATITUDE'] ?>, <?= $station['STATION_LONGITUDE'] ?>);
          map.addOverlay(createMarker(point, "<?= $station['STATION_NOM'] ?>","<?= $station['STATION_ALTITUDE'] ?>", "<?= $station['ENVIRONNEMENT_NOM'] ?>"));
        <? } ?>


        }
}



function showAddress(form) {
        if (form.lieu.value != '') {
                var lieu = form.lieu.value;
                EfGeocoder.getLatLng(lieu, function(position) {
                        if (!position) {
                                alert('Le lieu "'+lieu+'" est introuvable! Veuillez essayer un autre nom.')
                        } else {
                                placerMarkeur(position);
                                map.setCenter(position, 14) ;
                        }
                        // Vidage du champ du formulaire
                        form.lieu.value = '';
                });
        } else {
                var position = new GLatLng(form.lati.value, form.longi.value);
                placerMarkeur(position);
                map.setCenter(position, 14);
                // Vidage des champs du formulaire
                form.lati.value = '';
                form.longi.value = '';
        }
        return false;
}

// Ajout d'un nouveau marqueur dragable , depuis le geocodage uniquement (blocage du clic depuis la carte)

function placerMarkeur(position) {
        if (addrpnt) {
                map.removeOverlay(addrpnt);
        }
        addrpnt = new GMarker(position, {icon: Nicon, draggable: true,  title: 'Click sur la carte!'}) ;
        addrpnt.enableDragging() ;
        map.addOverlay(addrpnt) ;
        GEvent.addListener(addrpnt,'dragend',function() {
                        map.setCenter(addrpnt.getPoint());
                        afficherInfo(addrpnt.getPoint());
                    afficherMasque(addrpnt,"");

                }) ;

        GEvent.addListener(addrpnt, "click", function() {
                        map.setCenter(addrpnt.getPoint());
                        afficherInfo(addrpnt.getPoint());
                    afficherMasque(addrpnt,"");
                        });

        afficherInfo(position);
}

// Masque de mise a jour station :

function afficherMasque(position,nom,altitude,envir) {

        // On recupere les données du formulaire (caché) , notamment pour traiter simplement l'appel asynchrone de recup altitude

                                         var html = "<table>";
                                         if (nom.length>0) {
                                                  html = html + "<tr><td>Station :</td> <td>" + nom + "</td> </tr>";




                                                        html = html + "<tr><td>Environnement :</td> <td>" + envir + "</td> </tr>";

                                                        html = html + "</td></tr>";



                                         }
                                         else {
                                                 html = html + "<tr><td>Station</td> <td><input type='text' id='nom_station' value=''/> </td> </tr>";
                                                 <?
                                                                $requete_environnement = mysql_query("select ENVIRONNEMENT_ID, ENVIRONNEMENT_NOM from ENVIRONNEMENT");
                                                        ?>

                                                        html = html + "<tr><td>Environnement :</td>";
                                                        html = html + "<td><select id='environnement_id' name='environnement'>";

                                                        <? while ($env = mysql_fetch_row($requete_environnement)) { ?>
                                                        html = html + "<option value='<?= $env[0] ?>'><?= $env[1] ?>";
                                                        <? } ?>

                                                        html = html + "</select></td></tr>";

                                         }




                                         html = html + "<tr><td>Latitude :</td> <td>" +document.getElementById('ll_latitude').value + "</td> </tr>" ;
                         html = html + "<tr><td>Latitude :</td> <td>" +document.getElementById('ll_latitude').value + "</td> </tr>" ;
                         html = html + "<tr><td>Longitude :</td> <td>" + document.getElementById('ll_longitude').value + "</td> </tr>";
                         if (altitude) {
                                html = html + "<tr><td>Altitude :</td> <td>" + altitude + "</td> </tr>";
                         }
                         else {
                                html = html + "<tr><td>Altitude :</td> <td>" + document.getElementById('ll_altitude').value + "</td> </tr>";
                         }
                         if (nom.length==0) {
                                                html = html + "<tr><td></td><td><input type='button' value='Enregistrer' onclick='saveData()'/></td></tr>";
                                         }


              position.openInfoWindow(html);

}


// Sauvegarde station

// TODO : Retrouver marker depuis lat/lon

 function saveData() {

      var nom_station = escape(document.getElementById("nom_station").value);
      var lat_station = document.getElementById("ll_latitude").value;
      var lon_station = document.getElementById("ll_longitude").value;
      var alt_station = document.getElementById("ll_altitude").value;
          var environnement_id = document.getElementById("environnement_id").value;

      var url = "modules/observations/sauvegarde_station.php?nom_station=" + nom_station + "&action=C" + "&lat_station=" + lat_station +
               "&lon_station=" + lon_station +  "&environnement_id=" + environnement_id + "&alt_station=" + alt_station + "&participant_id=<?= $_SESSION['participant'] ?>";

      GDownloadUrl(url, function(data, responseCode) {
        if (responseCode == 200 && data.length <= 1) {
                        window.location = "<?= $url_page ?>";
        }
      });
    }

function afficherInfo(position) {
        var form = document.getElementById('form_coordonnee');
        afficherCoordonnees(position.lat().toString(), position.lng().toString());
}

function afficherCoordonnees(lati, longi) {
        var ellipsoid = 21;// 21 = WGS-84 (voir geo_constants)
        var xtm = 0;// 0 = UTM et 1= MTM
        var Coordonnee = calculer_coordonnee(lati, longi, ellipsoid, xtm);

        // Nous affichons des informations.
        document.getElementById('ll_latitude').value =  Coordonnee.latitude.toFixed(4);
        document.getElementById('ll_longitude').value = Coordonnee.longitude.toFixed(4);
        document.getElementById('ll_latitude_dms').value = Coordonnee.latitude_dms;
        document.getElementById('ll_longitude_dms').value = Coordonnee.longitude_dms;

        // Calcul altitude
        var script = document.createElement('script');
    document.body.appendChild(script);
    script.src = 'http://ws.geonames.org/' + "srtm3" + 'JSON?lat=' + lati + '&lng=' + longi + '&callback=GetAltitude';


}

function GetAltitude(result ) {
                        document.getElementById('ll_altitude').value = result.srtm3;
}

window.onload = load
window.onunload = GUnload

/* A supprimer ? */

/* Source : Code provenant du site http://pages.globetrotter.net/roule/utmgoogle.htm
*/

var deg2rad = Math.PI / 180;

/*
* Fonction de transformation de coordonnées géographiques entre ellipsoide.
* Source : http://home.hiwaay.net/~taylorc/bookshelf/math-science/geodesy/datum/transform/molodensky/
* Licence : "The source code in the listing below may be copied and reused without restriction, but it is offered AS-IS with NO WARRANTY."
*
* Parameters:
*     from:     The geodetic position to be translated. (objet : GeodeticPosition)
*     from_a:   The semi-major axis of the "from" ellipsoid. (double)
*     from_f:   Flattening of the "from" ellipsoid. (double)
*     from_esq: Eccentricity-squared of the "from" ellipsoid. (double)
*     da:       Change in semi-major axis length (meters); "to" minus "from"
*     df:       Change in flattening; "to" minus "from"  (double)
*     dx:       Change in x between "from" and "to" datum.  (double)
*     dy:       Change in y between "from" and "to" datum.  (double)
*     dz:       Change in z between "from" and "to" datum.  (double)
*
* Paramêtres :
*       from :          l'objet GeodesiquePosition à transformer (contient latitute, longitude et hauteur)
*       from_a :        le demi grand axe (=a) de l'ellipsoïde utilisé pour la position géodésique contenu dans l'objet "from".
*       from_f :        l'applatissement (=f) = 1/valeur... de l'ellipsoïde de l'objet "from"
*       from_esq :      l'excentricité au carré de l'ellipsoïde de l'objet "from"
*       da :            différence de demi grand axe (en metres) entre l'ellipsoïde "to" et "from".
*       df :            différence d'applatissement entre l'ellipsoïde "to" et "from".
*       dx :            différence pour l'axe x entre l'ellipsoïde "to" et "from".
*       dy :            différence pour l'axe y entre l'ellipsoïde "to" et "from".
*       dz :            différence pour l'axe z entre l'ellipsoïde "to" et "from".
*/
function transform(from, from_a, from_f, from_esq, da, df, dx, dy, dz) {
        var slat = Math.sin (from.lat);
        var clat = Math.cos (from.lat);
        var slon = Math.sin (from.lon);
        var clon = Math.cos (from.lon);
        var ssqlat = slat * slat;
        // "a divided by b"
        var adb = 1.0 / (1.0 - from_f);
        var dlat, dlon, dh;

        var rn = from_a / Math.sqrt (1.0 - from_esq * ssqlat);
        var rm = from_a * (1. - from_esq) / Math.pow ((1.0 - from_esq * ssqlat), 1.5);

        dlat = (((((-dx * slat * clon - dy * slat * slon) + dz * clat) + (da * ((rn * from_esq * slat * clat) / from_a))) + (df * (rm * adb + rn / adb) * slat * clat))) / (rm + from.h);

        dlon = (-dx * slon + dy * clon) / ((rn + from.h) * clat);

        dh = (dx * clat * clon) + (dy * clat * slon) + (dz * slat) - (da * (from_a / rn)) + ((df * rn * ssqlat) / adb);

        // Retour des données sous forme d'objet
        var GeodeticPosition =  {       lon: from.lon + dlon,
                                                                lat: from.lat + dlat,
                                                                h: from.h + dh
                                                        };

        return GeodeticPosition;
}

function test_transform() {
        var PositionNtf = { lon: 7.7372, lat: 48.6, h: 0 };
        var PositionWGS84 =     transform(PositionNtf, 6378249.2, 0.00340755, 0.0068034881, -112.2, 0.00005474, -168, -60, 320);
        alert('Longitude : ' + PositionWGS84.lon + 'Latitude : ' + PositionWGS84.lat + 'Hauteur : ' + PositionWGS84.h);
}

// Indiquer l'ellipsoid par son index dans le table et pour le xtm mettre 0 pour l'UTM  et 1 pour le MTM (?)
function geo_constants(ellipsoid, xtm) {
        // returns ellipsoid values
        ellipsoid_axis = new Array(); // valeur du demi grand axe (a) de l'ellipsoïde
        ellipsoid_eccen = new Array(); // valeur de l'excentricité au carré (e²) de l'ellipsoïde
        ellipsoid_axis[0] = 6377563.396; ellipsoid_eccen[0] = 0.00667054;  //airy
        ellipsoid_axis[1] = 6377340.189; ellipsoid_eccen[1] = 0.00667054;  // mod airy
        ellipsoid_axis[2] = 6378160; ellipsoid_eccen[2] = 0.006694542;  //aust national
        ellipsoid_axis[3] = 6377397.155; ellipsoid_eccen[3] = 0.006674372;  //bessel 1841
        ellipsoid_axis[4] = 6378206.4; ellipsoid_eccen[4] = 0.006768658;  //clarke 1866 == NAD 27 (TBC)
        ellipsoid_axis[5] = 6378249.145; ellipsoid_eccen[5] = 0.006803511;  //clarke 1880
        ellipsoid_axis[6] = 6377276.345; ellipsoid_eccen[6] = 0.00637847;  //everest
        ellipsoid_axis[7] = 6377304.063; ellipsoid_eccen[7] = 0.006637847;  // mod everest
        ellipsoid_axis[8] = 6378166; ellipsoid_eccen[8] = 0.006693422;  //fischer 1960
        ellipsoid_axis[9] = 6378150; ellipsoid_eccen[9] = 0.006693422;  //fischer 1968
        ellipsoid_axis[10] = 6378155; ellipsoid_eccen[10] = 0.006693422;  // mod fischer
        ellipsoid_axis[11] = 6378160; ellipsoid_eccen[11] = 0.006694605;  //grs 1967
        ellipsoid_axis[12] = 6378137; ellipsoid_eccen[12] = 0.00669438;  //  grs 1980
        ellipsoid_axis[13] = 6378200; ellipsoid_eccen[13] = 0.006693422;  // helmert 1906
        ellipsoid_axis[14] = 6378270; ellipsoid_eccen[14] = 0.006693422;  // hough
        ellipsoid_axis[15] = 6378388; ellipsoid_eccen[15] = 0.00672267;  // int24
        ellipsoid_axis[16] = 6378245; ellipsoid_eccen[16] = 0.006693422;  // krassovsky
        ellipsoid_axis[17] = 6378160; ellipsoid_eccen[17] = 0.006694542;  // s america
        ellipsoid_axis[18] = 6378165; ellipsoid_eccen[18] = 0.006693422;  // wgs-60
        ellipsoid_axis[19] = 6378145; ellipsoid_eccen[19] = 0.006694542;  // wgs-66
        ellipsoid_axis[20] = 6378135; ellipsoid_eccen[20] = 0.006694318;  // wgs-72
        ellipsoid_axis[21] = 6378137; ellipsoid_eccen[21] = 0.00669438;  //wgs-84 (et NAD83 : pourquoi mis avec WGS84 ?)

        if (ellipsoid == 0) {
                ellipsoid = 22;
        }

        --ellipsoid; // table indexed differently

        if (xtm == 1) {
                // WAS: if (ellipsoid > 22)
                scaleTm = 0.9999;
                eastingOrg = 304800.;
        } else {
                scaleTm = 0.9996;
                eastingOrg = 500000.;
        }

        // Retour des données sous forme d'objet
        var ellipsoid = {       axis: ellipsoid_axis[ellipsoid],
                                                eccentricity: ellipsoid_eccen[ellipsoid],
                                                eastingOrg: eastingOrg,
                                                scaleTm: scaleTm
                                        };
        return ellipsoid;
}

/**This routine determines the correct UTM letter designator for the given latitude
* returns 'Z' if latitude is outside the UTM limits of 84N to 80S
*/
function get_zoneletter(lat) {
        var zoneletter;
        if ((84 >= lat) && (lat >= 72)) zoneletter = 'X';
        else if ((72 > lat) && (lat >= 64)) zoneletter = 'W';
        else if ((64 > lat) && (lat >= 56)) zoneletter = 'V';
        else if ((56 > lat) && (lat >= 48)) zoneletter = 'U';
        else if ((48 > lat) && (lat >= 40)) zoneletter = 'T';
        else if ((40 > lat) && (lat >= 32)) zoneletter = 'S';
        else if ((32 > lat) && (lat >= 24)) zoneletter = 'R';
        else if ((24 > lat) && (lat >= 16)) zoneletter = 'Q';
        else if ((16 > lat) && (lat >= 8)) zoneletter = 'P';
        else if (( 8 > lat) && (lat >= 0)) zoneletter = 'N';
        else if (( 0 > lat) && (lat >= -8)) zoneletter = 'M';
        else if ((-8> lat) && (lat >= -16)) zoneletter = 'L';
        else if ((-16 > lat) && (lat >= -24)) zoneletter = 'K';
        else if ((-24 > lat) && (lat >= -32)) zoneletter = 'J';
        else if ((-32 > lat) && (lat >= -40)) zoneletter = 'H';
        else if ((-40 > lat) && (lat >= -48)) zoneletter = 'G';
        else if ((-48 > lat) && (lat >= -56)) zoneletter = 'F';
        else if ((-56 > lat) && (lat >= -64)) zoneletter = 'E';
        else if ((-64 > lat) && (lat >= -72)) zoneletter = 'D';
        else if ((-72 > lat) && (lat >= -80)) zoneletter = 'C';
        else zoneletter = chr(32 + 66); //This is here as an error flag to show that the Latitude is outside the UTM limits
        return zoneletter;
}

function parseCoordinate(coordinate, type, format, spaced) {
        var coordd;
        var coordm;
        var coords;
        var coordh = 0;

        if (coordinate.search(/(^ *-|[WOS])/i) >= 0) {
                coordh = -1;
        }
        if (coordinate.search(/(^ *\+|[NE])/i) >= 0) {
                coordh = 1;
        }

        // if (coordinate.search(/[EW]/i) >= 0 && !type) { type = 'lon'; }
        // if (coordinate.search(/[NS]/i) >= 0 && !type) { type = 'lat'; }

        coordinate = coordinate.replace(/,/g, '.');  // french commas
        // not sure really needed.
        coordinate = coordinate.replace(/[NESWO+\-]/gi,' ');  // add also: °, ', "  --  or replace everything that is not a number or a dot

        // alert("coordinate = " + coordinate);
        if (coordinate.search(/[0-9]/i) < 0) {
                alert("Can't parse input field");
                coordd = "";
                return;
        }

        // http://www.gpsvisualizer.com/calculators
        // http://www.javascriptkit.com/javatutors/redev2.shtml
        var parts = coordinate.match(/([0-9\.\-]+)[^0-9\.]*([0-9\.]+)?[^0-9\.]*([0-9\.]+)?/);
        //                           /                          (one or more)             /
        //                            ([0-9\.\-]+)              digits, dot or -, one set
        //                                        [^0-9\.]*            separator
        //                                                 ([0-9\.]+)?       digits, dot, zero or one set
        //                                                            [^0-9\.]*     separator
        //                                                                     ([0-9\.]+)?  digits, dot, zero or one set

        // *: 0|more   +: 1|more   ?: 0|1        http://www.javascriptkit.com/javatutors/redev2.shtml
        if (!parts[1]) {
                alert("Can't parse input field");
                coordd = '';
                return;
        } else {
                coordd = parts[1];
                if (parts[2]) {
                        coordm = parts[2];
                } else {
                        coordm = 0;
                }
                if (parts[3]) {
                        coords = parts[3];
                } else {
                        coords = 0;
                }
                // n = parseFloat(parts[1]);
                // if (parts[2]) { n = n + parseFloat(parts[2])/60; }
                // if (parts[3]) { n = n + parseFloat(parts[3])/3600; }
        }

        // Retour des données sous forme d'objet
        var Coordonnee =        {       coordd: coordd,
                                                        coordm: coordm,
                                                        coords: coords,
                                                        coordh: coordh
                                                };
        return Coordonnee;
}

function validate_dms(latd, latm, lats, latb, lonm, lond, lons, lonb) {
        var valid = 1;
        if (Math.abs(Number(latd)) >= 90)
                valid = 0;
        if (Number(latm) >= 60)
                valid = 0;
        if (Number(lats) >= 60)
                valid = 0;
        if (Math.abs(Number(lond)) >= 180)
                valid = 0;
        if (Number(lonm) >= 60)
                valid = 0;
        if (Number(lons) >= 60)
                valid = 0;

        return(valid);
}

// convert decimal degrees to dms
function convertir_en_dms(lat, lon) {
        var latbrg = 1;
        var lonbrg = 2;
        if (lat < 0)
                latbrg = 2
        if (lon < 0)
                lonbrg = 1;
        // LEW: have to round here, else could end up with 60 seconds :-)
        var tlat = Math.abs(lat) + 0.5 / 360000;  // round up 0.005 seconds (1/100th)
        var tlon = Math.abs(lon) + 0.5 / 360000;

        var tlatdm = Math.abs(lat) + 0.5 / 60000;  // round up 0.0005 minutes (1/1000th)
        var tlondm = Math.abs(lon) + 0.5 / 60000;

        var deglat = Math.floor(tlat);
        var t = (tlat - deglat) * 60;
        var minlat = Math.floor(t);

        var minlatdm = Math.floor((tlatdm - Math.floor(tlatdm)) * 60 * 1000) / 1000;

        var seclat = (t - minlat) * 60;
        seclat = Math.floor(seclat * 100) / 100;  //  works in js 1.4
        // seclat = seclat.toFixed(2);  // 2 decimal places js 1.5 and later

        var deglon = Math.floor(tlon);

        t = (tlon - deglon) * 60;
        var minlon = Math.floor(t);

        var minlondm = Math.floor((tlondm - Math.floor(tlondm)) * 60 * 1000) / 1000;

        var seclon = (t - minlon) * 60;
        seclon = Math.floor(seclon * 100) / 100;  // js 1.4
        // seclon = seclon.toFixed(2);  // js 1.5 and later

        var latb = '';
        if (latbrg > 0) {
                latb = 'N';// 1 = N (nord)
        } else {
                latb = 'S';// 2 = S (sud)
        }
        var latdms = deglat + "° " + minlat + "' " + seclat + "\"";

        var lonb = '';
        if (lonbrg > 0) {
                lonb = 'E'; // 2 = E (est)
        } else {
                lonb = 'W';// 1 = W (west)
        }
        var londms = deglon + "° " + minlon + "' " + seclon + "\"";

        var chaine_latlon_dms = latdms + '  ' + latb + '  ' + londms + '  ' + lonb;

        // Retour des données sous forme d'objet
        var LatLongDms =        {       lat: lat,
                                                        lat_degre: Math.floor(tlatdm),
                                                        lat_min: minlatdm,
                                                        lat_s: seclat,
                                                        lat_direction:latb,
                                                        lat_dms:latdms,
                                                        lon: lon,
                                                        lon_degre: Math.floor(tlondm),
                                                        lon_min: minlondm,
                                                        lon_s: seclon,
                                                        lon_direction:lonb,
                                                        lon_dms:londms,
                                                        ll_dms:chaine_latlon_dms
                                                };
        return LatLongDms;
}

function calculer_coordonnee(lati, longi, ellips, xtm) {
        // Récupération de l'ellipsoide
        var ellipsoid = geo_constants(ellips, xtm);
        var axis = ellipsoid.axis;
        var eccent = ellipsoid.eccentricity;
        var scaleTm = ellipsoid.scaleTm;
        var eastingOrg = ellipsoid.eastingOrg;

        // Nous parsons la latitude saisie
        var CoordonneeLat = parseCoordinate(lati);
        var latb = '';// vide
        if (CoordonneeLat.coordh != 0) {
                if (CoordonneeLat.coordh > 0) {
                        latb = 'N';// 1 = N (nord)
                } else {
                        latb = 'S';// 2 = S (sud)
                }
        }
        var latd = CoordonneeLat.coordd;
        var latm = CoordonneeLat.coordm;
        var lats = CoordonneeLat.coords;

        // Nous parsons la longitude saisie
        var CoordonneeLon = parseCoordinate(longi);
        var lonb = '';// vide
        if (CoordonneeLon.coordh != 0) {
                if (CoordonneeLon.coordh > 0) {
                        lonb = 'E'; // 2 = E (est)
                } else {
                        lonb = 'W';// 1 = W (west)
                }
        }
        var lond = CoordonneeLon.coordd;
        var lonm = CoordonneeLon.coordm;
        var lons = CoordonneeLon.coords;

        // cope with blank fields
        if (latd == '' || lond == '')   {
                alert("Latitude and longitude degrees must be entered");
                return;
        }

        // Indication de la direction par défaut pour la latitude et la longitude en France
        if (latb == '') {
                latb = 'N';// 1 = N (Nord)
        }
        if (lonb == '') {
                lonb = 'E';// 2 = E (Est)
        }

        // Validation
        var valid = validate_dms(latd, latm, lats, latb, lonm, lond, lons, lonb);
        if (valid == 0) {
                alert("Invalid degrees, minutes or seconds");
                return;
        }

        var lat = Number(latd);
        lat = lat + Number(latm) / 60;
        lat = lat + Number(lats) / 3600;
        if (latb == 'S') {  // 2 = S
                lat = lat * -1;
        }
        var lon = Number(lond);
        lon = lon + Number(lonm) / 60;
        lon = lon + Number(lons) / 3600;
        if (lonb == 'W') { // 1 = W
                lon = lon * -1;
        }

        if (lat >= 84 || lat <= -80) {
                alert("UTM latitudes should be between 84N and 80S\nCalculation will proceed anyway as locator will be valid.");
        }

        var k0 = scaleTm;
        var latrad = lat * deg2rad;
        var longrad = lon * deg2rad;
        var zonenum = Math.floor((lon + 180) / 6) + 1;
        // @dc (180 - (-70.5))/3 - 76
        if (eastingOrg == 304800.) {
                zonenum = Math.floor((180 - lon) / 3) - 76;  // MTM, only in Quebec
                if (zonenum < 3 || zonenum > 10) {
                        alert("MTM zone numbers only confirmed for 3-10, province of Quebec\nContinuing anyway");
                }
        }

        if (lat >= 56.0 && lat < 64.0 && lon >= 3.0 && lon < 12.0 ) {
                zonenum = 32;
        }
        // Zones speciales pour Svalbard
        if( lat >= 72.0 && lat < 84.0 ) {
                if (lon >= 0.0  && lon <  9.0 ) zonenum = 31;
                else if ( lon >= 9.0  && lon < 21.0 ) zonenum = 33;
                else if ( lon >= 21.0 && lon < 33.0 ) zonenum = 35;
                else if ( lon >= 33.0 && lon < 42.0 ) zonenum = 37;
        }

        var lonorig = (zonenum - 1) * 6 - 180 + 3;  //+3 puts origin in middle of zone
        // @dc 180 - (7+76) * 3 - 1.5
        if (eastingOrg == 304800.) {
                lonorig = 180 - (zonenum + 76) * 3 - 1.5;
        }
        var lonorigrad = lonorig * deg2rad;

        //Récupération de la lettre du fuseau xTM
        var letter = get_zoneletter(lat);

        // Calcul xTM
        var eccPrimeSquared = (eccent) / (1 - eccent);
        var N = axis / Math.sqrt(1 - eccent * Math.sin(latrad) * Math.sin(latrad));
        var T = Math.tan(latrad) * Math.tan(latrad);
        var C = eccPrimeSquared * Math.cos(latrad) * Math.cos(latrad);
        var A = Math.cos(latrad) * (longrad - lonorigrad);
        var M = axis * ((1 - eccent / 4 - 3 * eccent * eccent / 64 - 5 * eccent * eccent * eccent / 256) * latrad - (3 * eccent / 8 + 3 * eccent * eccent / 32 + 45 * eccent * eccent *eccent / 1024) * Math.sin(2 * latrad) + (15 * eccent * eccent / 256 + 45 * eccent * eccent * eccent / 1024) * Math.sin(4 * latrad) - (35 * eccent * eccent * eccent / 3072) * Math.sin(6 * latrad));

        var easting = (k0 * N * (A + (1 - T + C) * A * A * A / 6 + (5 - 18 * T + T * T + 72 * C - 58 * eccPrimeSquared) * A * A * A * A * A / 120) + eastingOrg);
        var northing = (k0 * (M + N * Math.tan(latrad) * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * A * A * A * A / 24 + (61 - 58 * T + T * T + 600 * C - 330 * eccPrimeSquared) * A * A * A * A * A * A / 720)));
        if (lat < 0) {
                northing += 10000000.0; //10000000 meter offset for southern hemisphere
        }

        // xTM
        // Arrondi au nombre supérieur
        // alert("easting = " + easting);
        preciseEasting = easting;
        easting = Math.floor(easting + .5);
        // alert("easting = " + easting);
        // alert("northing = " + northing);
        preciseNorthing = northing;
        northing = Math.floor(northing + .5);
        // alert("northing = " + northing);
        var chaine_xtm = '';
        var x_tm = '';
        if (eastingOrg == 304800.) {
                chaine_xtm = zonenum + '  ' + easting + 'm E  ' + northing + 'm N';
                x_tm = 'MTM';
        } else {
                chaine_xtm = zonenum + letter + '  ' + easting + 'm E  ' + northing + 'm N';
                x_tm = 'UTM';
        }

        // Latitude/Longitude
        var chaine_latlon = lat + ' ' + latb + ' ' + lon + ' ' + lonb;
        var LatLongDms = convertir_en_dms(lat,lon);

        // Retour des données sous forme d'objet
        var Coordonnee = {      xtm: x_tm,
                                                nord: northing,
                                                est: easting,
                                                zone_numero: zonenum,
                                                zone_lettre: letter,
                                                xtm_chaine: chaine_xtm,
                                                latitude:lat,
                                                longitude:lon,
                                                ll_chaine: chaine_latlon,
                                                latitude_dms:LatLongDms.lat_dms,
                                                longitude_dms:LatLongDms.lon_dms,
                                                ll_chaine_dms: LatLongDms.ll_dms
                                        };
        return Coordonnee;
}
//]]>

</script>

<?
}
?>