<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi All,<br>
<br>
I'd like to write the results of a kriging call to an ascii grid using
GSTAT -- but I would like to write BOTH the predictions and the
variances to grids. Is there a more elegant (and less dangerous) way to
do it than what I have below? <br>
<br>
<br>
depth_uk <- krige(DEPTH~slope.asc, depth, slope, vgm_depth_r)<br>
<br>
# write the predictions<br>
<br>
write.asciigrid(depth_uk "c:/junk/rk/predictions.asc")<br>
<br>
# replace predictions with variances and then write the variances<br>
<br>
depth_uk$var1.pred<-depth_uk$var1.var<br>
write.asciigrid(depth_uk, "c:/junk/rk/.asc")<br>
<br>
PS this sample code comes from Tomislav Hengl's page
(<a class="moz-txt-link-freetext" href="http://spatial-analyst.net/regkriging.php">http://spatial-analyst.net/regkriging.php</a>)<br>
<br>
Zev<br>
<div class="moz-signature">-- <br>
<title>Untitled Document</title>
<meta http-equiv="Content-Type" content="text/html; ">
<font face="Arial, Helvetica, sans-serif, Trebuchet MS" size="-1">Zev
Ross<br>
<strong><font color="#660000">ZevRoss Spatial Analysis</font></strong><br>
303 Fairmount Ave<br>
Ithaca, NY 14850<br>
(607) 277-0004 (phone)<br>
(866) 877-3690 (fax toll-free)<br>
<a class="moz-txt-link-abbreviated" href="mailto:zev@zevross.com">zev@zevross.com</a><br>
<a class="moz-txt-link-abbreviated" href="http://www.zevross.com">www.zevross.com</a></font><br>
</div>
</body>
</html>