Catégorie : Non classé

Ce qu'il faut savoir avant de démarrer :

1. Avant de commencer à charger les données sur le GPS, il faut savoir que les GPS Garmin utilisent le fomat .img pour stocker les cartes. 
   Pour en savoir plus sur ce format.
2. Le Garmin GPSMAP 64/64s permet d'ajouter plusieurs cartes au format .img contrairement à certains modèles plus anciens de GPS Garmin.

Introduction

Revenons en à nos moutons, grâce aux contributions de différentes personnes ou organisations, on peut facilement récupérer une carte OpenStreetMap compatible avec les GPS Garmin depuis un simple navigateur. En effet, des sites comme https://extract.bbbike.org que nous utilisons régulièrement pour préparer nos sorties terrain permet la récupération très rapide d'un fichier .img OSM de n'importe quelle zone de la planète 🌍.

Toutefois, ce n'est pas le seul à proposer ce service, il existe une multitude de solutions plus ou moins similaires. Vous pouvez consulter une liste non exhaustive par ici. Prenez le temps de bien identifier votre besoin pour choisir la carte qui répondra la mieux à vos attentes en terme de représentation de l'information.

L'exemple ci-dessous présente le site https://extract.bbbike.org

`

Le site minimaliste est assez bien conçu. Voici la démarche à suivre : 

1. Choisir le format désiré, dans notre cas Garmin!
   A noter qu'il existe différentes cartes Garmin, elles correspondent à des rendus différents de la donnée OSM (couleurs, pictos,…). En en selectionnant une vous aurez un aperçu du rendu sur le côté.
2. Renseigner son e-mail afin de recevoir le lien vers le fichier .img une fois qu'il est généré par le serveur 
3. Donner un nom à la zone choisie
4. Choisir l'emprise de la zone qui vous intéresse
5. Lancer la demande et attendre le mail avec le lien vers le fichier gmapsupp.img

​Une fois le fichier téléchargé 

1. Brancher le GPS à son ordinateur
2. Ouvrir l'espace disque correspondant au GPS
3. Copier le fichier gmapsupp.img dans le répertoire ./Garmin
4. Déconnecter le GPS de l'ordinateur
5. Lancer le GPS et admirer le résultat!
6. Pour ajouter plusieurs cartes, vous prouvez renommer le fichier gmapsupp.img. Par exemple : senegal.img, mali.img,…

Passer d'une carte à l'autre depuis le GPS :

1. Accéder à la carte en utilisant le bouton "PAGE"
2. Ensuite, appuyer sur "MENU"
3. Choisir "PARAMETRAGE CARTE" et appuyer sur "ENTER"
4. En bas du menu, aller sur "SELECTIONNER CARTE"
5. En utilisant la touche "ENTER" vous pouvez choisir d'activer ou non les cartes
6. Pour optimiser les performance d'affichage de la carte, il est préférable de désactiver les cartes non utilisées

Nota : Autre possibilité qui n'est pas développée dans cette article, vous pouvez créer vous même votre carte à partir des données OSM.

En complément : 

• Installer les cartes Open Street Map (OSM) dans le logiciel Basecamp de Garmin

Toujours dans le cadre de mes travaux de recherches, je devais comparer deux façons de calculer un NDVI cumulé annuel, 1) en intégrant toutes mes valeurs bi-mensuelles sur toute l'année, 2) en intégrant mes valeurs correspondant uniquement à la saison des pluies, c'est à dire environ de Mai à Octobre (contexte africain, au Niger).

Mes données d'entrée sont donc deux rasters calculés:

• raster1 = la manière 1 et comptant autant de bandes que d'années (ici 11)
• raster2 = la manière 2 et comptant autant de bandes que d'années (ici 11)

​Le résultat que je souhaitais obtenir en sortie était quelque chose comme ceci :

mais avec un graphique par année (ou par bande).

J'ai donc décidé de faire cela sur …. R (^^).

Voici donc le code que j'ai mis en place. Dans mon cas j'avais besoin de le faire tourner sur plusieurs zones donc sur plusieurs images, il est donc fait pour tourner sur plusieurs paires de fichiers, mais en supprimer la première boucle vous pouver le faire tourner sur uniquement deux fichiers.

######### Comparaison de rasterS multibandes #######

# input : 2 Rasters multibandes format ENVI (dans ce cas si, mais fonctionne avec autres formats raster)
# output : graphique au format pdf de la "corrélation" de la bande 1…n d'un raster1 avec la bande 1…n d'un raster2
# Leroux Louise, CIRAD, UMR TETIS 2014

######## Initialisation #########
## Chargement de packages nécéssaires

library(raster)
library(sp)
library(Hmisc)

## Initialisation de l'espace de travail
setwd("D:/MyFolder")

##Importation des données raster##
# Pour le premier raster
liste.file1<-list.files(path=".", pattern="motif1")
liste.file1.hdr<-liste.file1[-grep(« .hdr », liste.file1, fixed=T)] #Récupération du nom sans l'extension
 

# Pour le second raster
liste.file2<-list.files(path=".", pattern="motif2")
liste.file2.hdr<-liste.file2[-grep(« .hdr », liste.file2, fixed=T)]  #Récupération du nom sans l'extension

## Boucle sur les couples de fichiers ##
for (j in 1:length(liste.file1.hdr)){
  name1<-liste.file1.hdr[j]
  name2<-liste.file2.hdr[j]
  nband<-nbands(raster(name1)) #variable stockant le nombre de bandes dans l'image
  cell<-ncell(raster(name1)) #variable stockant le nombre de pixel d'une bande
  mat1<-matrix(0, ncol=nband, nrow=cell)
  mat2<-matrix(0, ncol=nband, nrow=cell)
  
  # Stockage des bandes du raster 1 dans un matrice
  for (i in 1:nband){
    mat1[,i]<-getValues(raster(name1, band=i))
  }
  
  # Stockage des bandes du raster 2 dans un matrice
  for (l in 1:nband){
    mat2[,l]<-getValues(raster(name2, band=l))
  }
  
  ## Création du graphique ##
  out_name<-paste(out_name, ".pdf",sep="")
  pdf(out_name, height=10,width=10)
  ye<-seq(2000,2010,1) # Séquence pour nommer les graphiques, ici, 1 bande= 1 année
  par(mfrow=c(4,3)) # Préparation de la fenêtre graphique : 4 lignes, 3 colonnes pour mettre un graphique par bande dans chaque case

  # Boucle sur les bandes
  for (k in 1:nband){
    df <- data.frame(mat1[,i],mat2[,k])
    x <- densCols(mat1[,k],mat2[,k], colramp=colorRampPalette(c("black", "white")))
    # Création d'une palette pour représenter la concentration des individus
    df$dens <- col2rgb(x)[1,] + 1L
    cols <-  colorRampPalette(c("#000099", "#00FEFF", "#45FE4F","#FCFF00", "#FF9400", "#FF3100"))(256)
    df$col <- cols[df$dens]
    plot(mat2[,k]~mat1[,k], data=df, pch=20, col=col, xlab=name_xlab, ylab=name_ylab, main=ye[k], cex.lab=0.75)
    reg<-lm(mat2[,k]~mat1[,k])
    abline(coef=coef(reg), lwd=1)
  }
  dev.off()
}

Et voici le résultat en sortie pour l'un de mes fichiers :