Mapa de concentração

O mapa foi feito através do programa de krigagem ordinária.

Programa:

data(jura)

jura.grid <- juragrid.dat

coordinates(jura.grid)=~Xloc+Yloc

gridded(jura.grid) = TRUE

class(jura.grid)

coordinates(jura.pred)=~Xloc+Yloc

class(jura.pred)

m <- vgm(12.524976,"Sph",1.183541,1.304963)

x <- krige(Co~1, jura.pred, jura.grid, model=m)

class(x)

image(x)

spplot(x[1],main='Mapa de concentração de cobalto (ppm)',col.regions=bpy.colors(20),xlab='Xloc (km)',ylab='Yloc (km)')

Resultado: