Return to Revisions

You can use this R script

You will need these packages:

• library(rgdal)
• library (raster)

For example for a map called "Prodes1":

Prodes1<-readOGR(dsn="E:/PRODES/PDigital2014_22768_shp", layer="PDigital2014_22768__pol", dropNULLGeometries=TRUE)

Prodes1$Grd_ranks<-rank(Prodes1$mainclass) #Creating a numeric column for rasterize

i<-1
Prodes1$value<-rank(Prodes1$mainclass)
for (i in 1:(length(Prodes1$mainclass))) {
  if (Prodes1$mainclass[i]=="DESFLORESTAMENTO") Prodes1$value[i]<-0 else 
    if (Prodes1$mainclass[i]=="FLORESTA") Prodes1$value[i]<-1 else
      if (Prodes1$mainclass[i]=="HIDROGRAFIA") Prodes1$value[i]<-2 else
        if (Prodes1$mainclass[i]=="NAO_FLORESTA") Prodes1$value[i]<-3 else
          if (Prodes1$mainclass[i]=="NUVEM") Prodes1$value[i]<-4 else
            if (Prodes1$mainclass[i]=="RESIDUO") Prodes1$value[i]<-5 else Prodes1$mainclass[i]<-Prodes1$classe[i])
}                            

#Creating a ROI to receive the raster data
ROI2=raster() 
extent(ROI2)<-extent(Prodes1)
res(ROI2)<-0.002232143
proj4string(ROI2)<-CRS(proj4string(Prodes1))
r<-rasterize(Prodes1, ROI2, "value", fun="first")

#Saving the raster TIFF

writeRaster(r, filename="Prodes.tiff")