Prova 2

library(sp)

library(gstat)

library(lattice)

library(latticeExtra)

caminho ="C:/Users/kamil/Desktop/GEOLOGIA/P2-Geoestat/"

#Planilha com coordenadas utm e valores da variável em tela;

data=read.table(header=TRUE,paste(caminho,"duckhead.txt",sep = ""))

edit(data)

coluna=data["V"]

coord=subset(data,select=c("UTMX","UTMY"))

edit(coord)

amostra=cbind(coord,coluna)

colnames(amostra)[3]="V"

edit(amostra)

#1

caminho="C:/Users/kamil/Desktop/GEOLOGIA/P2-Geoestat/"

vgm1=variogram(V~1,locations=~UTMX+UTMY, data=amostra)

f<-fit.variogram(vgm1,vgm(0.01,"Sph",190,0.1))

ff<-variogramLine(f,maxdist=650,n=700,min=1.0)

plot(ff,col="orange")

points(vgm1[,2],vgm1[,3], col="red")

print(f)

#3

s.grid=GridTopology(c(406900,87900),c(60,40),c(50,50))

s.grid=SpatialPoints(s.grid)

m=vgm( 0.033081720,"Sph", 143.3975, 0.003302009)

pred=krige(V~1, ~UTMX+UTMY, model = m,data = amostra, newd = s.grid)

dfpred=as.data.frame(pred)

mz=matrix(dfpred[,3], nrow=50, ncol=50,byrow=FALSE)

nmz=matrix(nrow=50, ncol=50)

for(i in 1:50)

  for(j in 1:50)

  {nmz[i,j]=mz[i,51-j]}

persp(x = seq(407000,409940, by = 60), y=seq(88600, 90599,by=40),nmz,xlab="UTMX",ylab="UTMY", main="Diagrama de Bloco do Vadium (V)", theta=30, phi=30, r=100, d=20, scale=TRUE, col="Orange")