Dear All, I want to do regression kriging and cokriging with the following datasets: a point data set, DEM, EVI, landcover map, TWI, TRI and soil map. All these maps are spatialgriddataframes and have the same extend. The landcover map has 9 classes.
Question1. How can i remove the ninth class in R because it does not have to be included in geostatical analysis. Question 2 how can i use the map created in Q1 to clip the other 5 aboventioned maps (eg DEM etc) or how can i create a mask from the map in Q1. Please find the summary of landcover: summary(landcover) Object of class SpatialGridDataFrame Coordinates: min max x 464469.4 813469.4 y 8105042.0 8965042.0 Is projected: NA proj4string : [NA] Number of points: 2 Grid attributes: cellcentre.offset cellsize cells.dim x 464969.4 1000 349 y 8105542.0 1000 860 Data attributes: Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 1.00 4.00 5.00 5.31 6.00 9.00 20735.00 May you please help. Regards, Rise.M [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo