<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=iso-8859-1"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:70.85pt 70.85pt 56.7pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Dear All,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am facing problem to obtain the nested variogram model for the 3D lithological data. I have 3D data, in which the coordinate are x,y,z (z corresponds to the depth) and the lithofacie information for different materials as a binary values (Ind01, Ind02, Ind03, Ind04, Ind05) . Now I would like to develop sequential indicator simulation to obtain the facies distribution for each of the category. I have used gstat to develop the vertical and horizontal varigram but unable to fit as a nested variogram for both vertical and horizontal direction. Is there any one help me to combine the vertical and horizontal variogram in one nested structure. <o:p></o:p></p><p class=MsoNormal>The data I have used from the greengoblin data. A dropbox link has been placed at the end of this post to obtain the data.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The code I have used as follows to obtain the vertical and horizontal variogram for the Ind03 facies (medium sand)<o:p></o:p></p><p class=MsoNormal>green_gstat <- as.data.frame(green_goblin)<o:p></o:p></p><p class=MsoNormal>coordinates(green_gstat) = ~X+Y+Z<o:p></o:p></p><p class=MsoNormal>names (green_gstat)<o:p></o:p></p><p class=MsoNormal>class(green_gstat) <o:p></o:p></p><p class=MsoNormal># vertical variogram <o:p></o:p></p><p class=MsoNormal>ind03_vert <- variogram((Ind03)~X+Y+Z, green_gstat, cutoff = 250, width= 2, alpha= 0, beta = 90)<o:p></o:p></p><p class=MsoNormal># fitting the variogram model<o:p></o:p></p><p class=MsoNormal>ind03_fit <- fit.variogram (ind03_vert, vgm(1, "Exp", 500,1))<o:p></o:p></p><p class=MsoNormal># plotting<o:p></o:p></p><p class=MsoNormal>plot(ind03_vert, ind03_fit)<o:p></o:p></p><p class=MsoNormal><img width=544 height=248 style='width:5.6666in;height:2.5833in' id="Picture_x0020_2" src="cid:image001.jpg@01D2CF27.47841B80" alt=Rplot01><o:p></o:p></p><p class=MsoNormal># estimating the horizontal variogram in omni direction<o:p></o:p></p><p class=MsoNormal>Ind03_hori_omni <- variogram(ind03, cutoff = 15000, width = 1000, alpha = 0)<o:p></o:p></p><p class=MsoNormal># fiting variogram mdoel<o:p></o:p></p><p class=MsoNormal>ind03_fit_dir <- vgm(.188, "Exp", 1500, .05, anis = c (45, .25))<o:p></o:p></p><p class=MsoNormal># obtain the fit model informaiotn <o:p></o:p></p><p class=MsoNormal>Ind03_fit_vrio_hor <- fit.variogram(Ind03_hori_omni,vgm(.188, "Exp", 1500, .05, anis = c (45, .25)))<o:p></o:p></p><p class=MsoNormal># plotting the data<o:p></o:p></p><p class=MsoNormal>plot(Ind03_hori_omni,ind03_fit_dir)<o:p></o:p></p><p class=MsoNormal><img width=521 height=320 style='width:5.425in;height:3.3333in' id="Picture_x0020_1" src="cid:image002.jpg@01D2CF27.47841B80" alt=hori><o:p></o:p></p><p class=MsoNormal># now I have tried to go for the nested variogram approach, typically what I have found there, that one can create the nested model like this from the same variogram, but using diffent fitting options. <o:p></o:p></p><p class=MsoNormal>ind03_vertnested <- fit.variogram (<b>ind03_vert</b>, vgm (1, "Exp", 500, add.to = vgm (.188, "Exp", .05, anis = c(45, 0.25))))<o:p></o:p></p><p class=MsoNormal>here I have only used the ind03_vert (variogram), but I also have a horizontal variogram model (Ind03_hori_omni). So, how I can incorporate two model into a single variogram model. A feedback in this regard would be highly appreciate. Or, is there any other option to create a 3D variogram modeling. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>To have the data, please follow this dropbox link for the easy access. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Link to the data- <a href="https://www.dropbox.com/s/6xl8zh7ofntlna6/GreenGoblin.csv?dl=0">https://www.dropbox.com/s/6xl8zh7ofntlna6/GreenGoblin.csv?dl=0</a><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><b><span style='font-size:12.0pt;color:#1F497D;mso-fareast-language:DE'>--<o:p></o:p></span></b></p><p class=MsoNormal><b><span style='font-size:12.0pt;color:#1F497D;mso-fareast-language:DE'>A.B.M Firoz</span></b><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'><br><b>Researcher- GIS & Hydrological Modeling<o:p></o:p></b></span></p><p class=MsoNormal><b><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'><o:p> </o:p></span></b></p><p class=MsoNormal><b><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'>ITT</span></b><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'>- Institute for Technology and Resources Management <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'>in the Tropics and Subtropics<o:p></o:p></span></p><p class=MsoNormal><b><span style='color:#B00004;mso-fareast-language:DE'>Technology</span></b><span style='color:#E76B0C;mso-fareast-language:DE'> </span><b><span style='color:#FC4F08;mso-fareast-language:DE'>Arts</span></b><b><span style='color:#F59236;mso-fareast-language:DE'> </span></b><b><span style='color:#BE1687;mso-fareast-language:DE'>Sciences</span></b><span style='color:black;mso-fareast-language:DE'><br></span><b><span style='font-size:10.0pt;color:black;mso-fareast-language:DE'>TH Köln</span></b><span style='font-size:10.0pt;color:black;mso-fareast-language:DE'> - University of Applied Sciences <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;color:black;mso-fareast-language:DE'><o:p> </o:p></span></p><p class=MsoNormal><span lang=DE style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'>T: +49 221 8275-2059<o:p></o:p></span></p><p class=MsoNormal><span lang=DE style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'>F: +49 221 8275-2736<br>E : </span><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'><a href="mailto:abm.firoz@fh-koeln.de"><span lang=DE style='color:blue'>abm.firoz@fh-koeln.de</span></a></span><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'> <span lang=DE><o:p></o:p></span></span></p><p class=MsoNormal><span lang=DE style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;color:#1F497D;mso-fareast-language:DE'>Kalk- Campus<br>RobertStrasse 2<br>51105 Köln Germany</span><b><span style='font-size:13.5pt;color:#1F497D;mso-fareast-language:DE'><o:p></o:p></span></b></p><p class=MsoNormal style='margin-bottom:.25in'><span lang=DE style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333;mso-fareast-language:DE'><a href="http://www.tt.fh-koeln.de/"><span lang=EN-US style='color:blue'>www.tt.fh-koeln.de</span></a></span><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333;mso-fareast-language:DE'><br></span><span lang=DE style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333;mso-fareast-language:DE'><a href="http://www.th-koeln.de/"><span lang=EN-US style='color:blue'>www.th-koeln.de</span></a></span><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>