search for: kriging

I'm performing ordinary and simple kriging from a set of non-negative values: >krige.control.sk<-krige.control(type.krige="sk",obj.model=my.variogram.model,beta=my.variogram.model$beta) >sk<-krige.conv(cadiz.geo,data=amostragem.cadiz$dia1.std,locations=as.matrix(cadiz.polygrid[,c(1,2)]),krige=krige.control.sk) kri...

Juan Zwolinski wrote: > > I'm performing ordinary and simple kriging from a set of non-negative > values: > > > >krige.control.sk<-krige.control(type.krige="sk",obj.model=my.variogram.model,beta=my.variogram.model$beta) > > >sk<-krige.conv(cadiz.geo,data=amostragem.cadiz$dia1.std,locations=as.matrix(cadiz.polygrid[,c(1,2)])...

Hi all, I tried to produce some kriged surfaces with geoR (latest version). The size of the grid should be around 900 x 650 cells (what I find is not a very big grid), and the number of points is around 2500. The command krige.conv stopped after arround 5 min saying it can not allocate a vector with around 1.5 billion units. Sounds reasonable. Is there a workaround? How would I partition the

Dear List I am using the RandomFields package, and I have a problem when 2d-kriging, with an anisotropy, some values from an irregular network. It works well when : - 2d-kriging, without any anisotropy, some data from an irregular network - 2d-kriging, with and without anisotropy, some data from a regular network - 3d-kriging, with and without anisotropy, some data from a regular...

Hello, I'm interested in doing Bayesian kriging using R. I see that the package geoR has a function that will allow one to do this (krige.bayes). However, my data are not in the form of points, but rather they are blocks that represent spatial averages (i.e., the number of fishing hooks per month in a given lat x long square). I am there...

...36083 149 319675 2011 8 28 0.00 0.000000000 0.00000000 COOP 36.13056 -80.54806 > v.fit = vgm(phi=1.2914,model="Mat",range=14.1611,nugget=0.23,kappa=83.1047) > krige.pred.out = krige(formula=bias~1,locations=~lon+lat,data=casedata,newdata=gridmeta,model=v.fit) [using ordinary kriging] > krige.pred.out lon lat var1.pred var1.var 1 -86.28771 33.91335 NaN NaN 2 -86.27388 33.94722 NaN NaN 3 -86.26003 33.98111 NaN NaN 4 -86.24615 34.01500 NaN NaN 5 -86.23225 34.04890 NaN NaN 6 -86.21835 34.08280 NaN...

Dear all, When a conventional kriging and then a contour plot is limited with a polygon (as possible with krige.conv and contour.kriging), the polygon border is displayed in black by default. > kc<-krige.conv(CZdata,loc=pred.grid,borders=czpoly,krige=krige.control(obj.m=ls)) > plot(CZcoord,xlab="x",ylab="y&quot...

HI, I wonder does anyone have experience with doing sequential gaussian simulation with krige() function in gstat? I find it VERY slow compared to use krige() to achieve kriging function itself.. I wonder why, is that because it has to model the variogram, and do the kriging separately for each point to be simulated? so it would be N times slower to achieve the simulation than the kriging if the number of points to be estimated is N?? Any hints is welcome and appreciate...

Hi, everybody! i created a imagem by kriging using geoR package. I imported points from GRASS("zn", after converted to geodata "zn_geo"), the border "zn_border" and a raster mask. Then i interpolated the points by kriging and created a raster image. Now, i need export this image back to GRASS to use it in the mo...

Dear all, I would like to estimate rainfall with the krige function. First, I produce a polygon of my region of interest (where poly_test.txt is a x,y suite of points defining a polygon, obtained with the grass v.out.ascii command) : user>border_limite<-read.table("/home/lionel/rwork/poly_test.txt",header=FALSE) user>polygone<-list(x=border_limite[,1],y=border_limite[,2])

...out the same as the linear model except now the output parameters are tausq=10.02, sigmasq=318909, and phi = 3714567. And calling krige.conv() again, the predicted values are in the range that I would expect, and the kriged map looks fine. I?m not sure if this has something to do with how the kriging is using the model beyond the 5 km. It seems like I need to be able to set the function to only apply the linear weighting to the linearly correlated portion (data less than 5km away), and the rest to 0, but I?m not sure how to do that (or maybe I?m completely going down the wrong track). I?m...

Hi, all, Have anyone used kringing included in R? How is it? Does it handle anisotropy data well? How does it compare with Kriging in Arc/Info? or other geostatistics software customized to do kriging or other geostatistics functions? I tried Easykriging, a geostatistics tool developed for Matlab. It has very nice GUI, but it does not provide library which i can call in my programs. so it is good for doing it mannually for sm...

Hi list, Well, this time I’ve a doubt with mapping generation. I was already able to read and plot shapefiles, plot point on this map. All this without any problems, but now I want to do something that I think, only Golden Software Surfer is capable of. I would like to plot a kriging result on the existing map (map script below). Well I looked for on the net, but I was not able to be sure of which packages I’ll need to perform the kriging (maybe geoR?) and then plot it in the area map. Is R able to do that? What packages I need to know to do so? I just need a starting point t...

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080730/00e14657/attachment.pl>

Hi! Im new to R and I need to interpolate a shapefile using kriging. I've been able to plot/read the shapefile using the package maptools or rgdal. I've searched the internet for sample codes but most of the kriging codes that I've found done in R is done using txtfiles or CSVs.  An example could be of great help. Thanks. [[alternative HTML version del...

Hello All! I like to do some indicator kriging with R. So far I used geoR for simple and ordinary kriging. Does anybody know which package I should use? Thanks for your help! Sascha! -- Jetzt dabei sein: http://www.shortview.de/wasistshortview.php?mc=sv_ext_mf at gmx

Full_Name: Albrecht Gebhardt Version: 0.90.0 OS: Linux + Tru64 Submission from: (NULL) (143.205.61.72) This is where example(krige) from library(sgeostat) stops with R 0.90.0: krige> image(grid$x, grid$y, grid$krige$zhat, add = T) Error in length(x) == nrow(z) : comparison (1) is possible only for vector types in previous versions of R the above image() call worked, because image() did the

R Users: ...been working with the sp and gstat packages for the past couple of days in an effort to analyze a set of ~ 200 soil samples collected from various eastings, northings, and depths and containing a wide range of measured hydrocarbon concentrations. Thus far, I've managed to import the data, log-transform the concentrations, assign coordinates, generate and fit a variogram model and

...I have a longitude vector and a latitude one. Associated to these coordinates, i have a matrix with some data at some coordinates but not all. Lon <- seq(136.025,144.975,0.05) Lat <- rev(seq(-66.975,-65.525,0.05)) dim(z) <- c(Lon,Lat) And i have tried to apply to these data a kriging function. But first i need to reshape these 3 variables to have a dataframe with these and with the corresponding z values with the good coordinates. Do you know how i can do that? And i will have the same problem with the object that return the krige function (gstat package). It gives me a dat...

Hi, I'm looking for software that can perform kriging on systems with dimensionality higher than 3, say d=5. Are anyone aware of packages in R that can do this? Thanks, Eivind Sm??rgrav ------------------------------------------------------------------- The information contained in this message may be CONFIDENTIAL and is intended for the addresse...