Hello Didier,
Many thanks for your replay
Now is everything working on, when I did it manually, by using your suggested command
I have last question in this topic. Namely, what settings are given as default (about Anamorphosis) in
simtub function
When I set
NA in parameter
anam and
flag.back will be equal
FALSE, what could be done with raw data during simulation? Then I also obtain results with different scale (more symmetric/gaussian). If I set
anam=model.anam and
flag.back=TRUE, I have to also use
anam.y2z function (it seems works the same like
flag.back=TRUE) after calculated
simtub? Below I pasted my code to make my questions clear.
Thank you a lot for your time and given answers.
Best regards,
Ania
- Code: Select all
library(RGeostats)
name=file.choose()
B4=read.csv(name,header=T,na="MISS")
class(B4)
val=db.create(B4)
val=db.locate(val,4:6,"x")
val=db.locate(val,7,"z")
model.anam <- anam.fit(val,type="gaus",nbpoly=100,draw=T)
val <- anam.z2y(val,anam=model.anam)
grid3D=db.create(flag.grid=TRUE,nx=c(48,55,6),x0=c(20,10,3),dx=c(1,1,3),ndim=3,nvar=1,autoname=TRUE)
grid3D=migrate(val,grid3D)
vario3D_4dir=vario.calc(val,lag=5,nlag=15,dir=c(0,45,90,135))
plot(vario3D_4dir,npairdw=T,npairpt=T)
x11()
m4dir=model.auto(vario3D_4dir,struct=melem.name(c(1,2,3)))
plot(m4dir,add=T,vario=vario3D_4dir)
x11()
moving.neigh3D=neigh.create(ndim=3,type=2,flag.sector=T,flag.aniso=T,flag.rotation=F,nmini=1,nmaxi=3,radius=c(120,120,100))
result3D=simtub(val,grid3D,model=m4dir,neigh=moving.neigh3D,anam=model.anam,flag.back=TRUE,seed=232131,nbtuba=100,nbsimu=100)
result3D <- anam.y2z(result3D,anam=model.anam)
plot(result3D,col=topo.colors(100),pos.legend=5,xlab="x [m]",ylab="y [m]")