2.2 Collecte des données

2.2.1 Réseau d’échantillonnage

Les parcelles d’échantillonage couvrent une surface de 30 x 30 mètres et correspondent en taille, forme et position aux pixels des images Landsat. Les parcelles d’échantillonnage sont régulièrement réparties sur le territoire du Togo, sur une grille d’une taille de maille de 480 mètres (tous les 16 pixels Landsat) ou 4,3 parcelles d’échantillonnage par km². Sur l’ensemble du territoire togolais, cela donne un réseau de 250 000 parcelles.

Si nécessaire pour autres utilisation, la taille de maille de 480 mètres permet d’enlargir la maille à 679 ou 960 mètres. D’autre part, il est également possible de comprimer la maille à 340, 240, 170, 120, 85, … mètres.

La grille d’échantillonnage est basée sur l’étendue du Togo (alignée avec la grille Landsat) et la taille du maillage. Les points d’échantillonnage résultants sont attribués avec les coordonnées x et y et l’ID correspondant. La grille de points résultante est enregistrée sous forme de shapefile (télécharger archive ZIP).

Script R: 01_SSTS/02_BdD/_src/create-grid.R

###############################################################################
# create-grid.R: créer une grille de points d'observation SSTS
# -----------------------------------------------------------------------------
# Bern University of Applied Sciences
# Oliver Gardi, <oliver.gardi@bfh.ch>
# 13 Mai 2020

# Définitions des variables ===================================================

OUT.DIR <- paste0(DIR.SST.BDD, "/01_reseau-SSTS")
if(!dir.exists(OUT.DIR)) dir.create(OUT.DIR)

RES <- 480   # Résolution de la grille (mètres)

# Grille d'observation ========================================================

# Coordonnées min/max en ligne avec images Landsat
x.min <- RES * extent(TGO)@xmin %/% RES
x.max <- RES * extent(TGO)@xmax %/% RES
y.min <- RES * extent(TGO)@ymin %/% RES
y.max <- RES * extent(TGO)@ymax %/% RES

# Creation de la grille
frame.points <- SpatialPoints(expand.grid(seq(x.min, x.max, by=RES), 
                                          seq(y.min, y.max, by=RES)), 
                              proj4string=UTM.31)[TGO]

# Ajouter attribues PLOTID, xcoords et ycoords
frame.points$xcoords <- frame.points@coords[,1]
frame.points$ycoords <- frame.points@coords[,2]
frame.points$PLOTID  <- paste0(str_pad(frame.points@xcoords, 7, "left", "0"), "_", 
                               str_pad(frame.points@ycoords, 7, "left", "0"))

# Sauveguarder comme fichier Shapefile
writeOGR(frame.points, dsn=OUT.DIR, layer="TGO_frame_480m", 
         driver="ESRI Shapefile", overwrite=TRUE)