Simulações Gaussianas

k) Nove Simulações Gaussianas do Atributo em Mapas de Concentração.

s.grid<-GridTopology(c(178260,329460),c(160,160),c(22.5,30)) s.grid<-SpatialPoints(s.grid)gridded(s.grid)<-TRUEdata(meuse.alt)m<- vgm(0.0045326051,"Exp",241.1351, 0.0007736096) xx <- krige(log(alt)~1, ~x+y, model = m, data = meuse.alt, newd = s.grid )X11()xx <- krige(log(alt)~1, ~x+y, model = m, data = meuse.alt, newd = s.grid, nsim=9)spplot(xx["sim1"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim1")X11()spplot(xx["sim2"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim2")X11()spplot(xx["sim3"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim3")X11()spplot(xx["sim4"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim4")X11()spplot(xx["sim5"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim5")X11()spplot(xx["sim6"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim6")X11()spplot(xx["sim7"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim7")X11()spplot(xx["sim8"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim8")X11()spplot(xx["sim9"],xlab="Xloc(0-4)",ylab="Yloc(0-4)",main="Valores Simulados-sim9")X11()















l) Mapa de Concentração da Mediana do Atributo Simulado.

s.grid<-GridTopology(c(178260,329460),c(160,160),c(30,30)) s.grid<-SpatialPoints(s.grid)gridded(s.grid)<-TRUEdata(meuse.alt)#set.seed(9999)
m<- vgm(0.0045326051,"Exp",241.1351, 0.0007736096)xx <- krige(log(alt)~1, ~x+y, model = m, data = meuse.alt, newd = s.grid, nsim=30)
X11()
aux=matrix(1:27000,900,30)
for (i in 1:900)
for (j in 1:30)
{ aux[i,j]=xx[[j]][i]}
#obtendo as medianas
for (i in 1:900)
{xx[[1]][i]=median(aux[i,])}
#plotando o mapa
spplot(xx[,1], xlab="xloc(0-4)", ylab="yloc(0-4)", main="Mediana")