####################PREMIERES ETAPES################################################
# ---Clear the workspace 
rm(list=ls())
 
#Load libraries
library(raster)
library(rgdal)
library(stringr)
 
###################DEFINIR DES ENTREES###############################################
#set working directory. Il faut definir ici le directoire. Il faut mettre / au lieu de \
wd <- "C:/Users/olofsson/Desktop/STR_example/"
setwd(wd)
 
# Read a raster, GeoTiff or something. Il faut mettre le nom du raster de stratification. 
# Il faut que le raster soit dans le working directory
Raster <- raster('Col_strat_buffer_subset_UTM.tif')
 
##Resolution de la carte utilise
Res <- 30
 
##Fichier avec les superficies par estrate
Output <- "Col_strat_buffer_wh.csv"

###NE TOUCHER CE QUI EST ECRIT CI DESSUS##################################################
############################CALCUL DU NOMBRE DE PIXELS PER CLASSE DE LA CARTE#############
t <- Raster@data@attributes
map_area <- as.data.frame(t)
map_area$Ha
map_area$Ha <- map_area$Count*Res*Res/10000
write.csv(map_area, file= Output, row.names= FALSE)