#read ascii file header <- scan("Typologie des climats francais.asc", nlines=6, what = character()) num1<- seq(1,length(header),2) num2<- seq(2, length(header),2) header.df<-data.frame(labels=header[num1], values=header[num2]) View(header.df) header.df<- data.frame(labels=header[num1], values=as.numeric(header[num2])) #read the cells' information from the ascii file read.table("Typologie des climats francais.asc", header=F, skip=6)-> mat as.matrix(mat)-> mat dim(mat) raster(mat) -> rasterClimat rasterClimat xmin <- header.df[3,2] ymin <- header.df[4,2] xmax<- header.df[5,2]*header.df[1,2]+xmin ymax<-header.df[5,2]*header.df[2,2]+ymin extent(rasterClimat)<- extent(xmin, xmax, ymin, ymax) plot(rasterClimat)