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]]