Displaying 20 results from an estimated 2000 matches similar to: "Bayesian Block Kriging?"
2004 Feb 23
1
border of a polygon in contour.kriging - geoR
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",type="n",asp=1)
>
2003 Jun 26
1
krige in gstat() package
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
2008 Dec 04
2
Plotting a kriging on a map
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
2010 Apr 07
1
kriging problem - very urgent
Hi everybody,
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
2011 Jan 05
1
Prediction error for Ordinary Kriging
Hi ALL,
Can you please help me on how to determine the prediction error for ordinary
kriging?Below are all the commands i used to generate the OK plot:
rsa2 <- readShapeSpatial("residentialsa", CRS("+proj=tmerc
+lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0
+ellps=intl +units=m +no_defs"))
x2 <- readShapeSpatial("ptna2",
2012 Oct 30
2
issues with krige function
Greetings all,
Ran into a strange problem with the krige function from geoR. The
problem that I am having is that while the krige function seems to
work well, the resulting predicted values are all NAs. Given the size
of the datasets I am working with can't attach it, but I can provide
snippets of the datasets.
> casedata
station year month day obs mpe bias type
2008 Jul 20
2
Indicator Kriging?
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
2011 Oct 05
2
kriging shapefiles
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 deleted]]
2011 Jan 06
1
Cross validation for Ordinary Kriging
ear ALL,
The last part of my thesis analysis is the cross validation. Right now I am
having difficulty using the cross validation of gstat. Below are my commands
with the tsport_ace as the variable:
nfold <- 3
part <- sample(1:nfold, 69, replace = TRUE)
sel <- (part != 1)
m.model <- x2[sel, ]
m.valid <- x2[-sel, ]
t<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0))
cv69
2004 Mar 24
2
geoR - help for bayesian modelling
Hi,
I am trying to do a bayesian prediction for soil pollution data above
a certain threshold, using geoR.
Everything is working fine until i am doing the krig.bayes. I tried to
do the prediction on a grid 67 by 113 cells and my computer is
freezing to death. At larger numbers of cells it tells me after a while
that it reaches the max. memory of 511 Mb. My computer has only
512 Mb of RAM.
2013 Apr 04
5
help with kriging interpolation
All,
I am new to using R and know some basics. I wish to use kriging in R to
do the following:
given data Y =f(X1,X2,X3,.....,Xn) --1000+ irregular measured data set.
I would like to be able to get a single value y given sinle input set
(x1,x2,x3,...xn)
A google search on this takes me lierally to the same example on involving
analysis with soil sampling and I cannot figure out how to
2002 Oct 16
0
Ordinary and simple 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)
krige.conv: model with constant mean
krige.conv: Kriging
2002 Oct 16
0
[Fwd: Ordinary and simple kriging]
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)
>
>
2010 Aug 19
0
2d kriging with anisotropy on an irregular network (RandomFields Package)
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 network
Here is
2010 Jan 07
1
kriging with geoR package
Hi all,
I have to draw a map of microorganisms repartition in a sample. I use the
"image" function of geoR package to build the map but there are three
problems:
1) the legend is Horizontal (with or without vert.leg=TRUE)
2) I want to plot the position of organisms and level curves of density on
the map but I don't find the script to insert these data on the map.
3)the grey range
2004 Jan 08
1
Using geoR krige.conv
Hi,
I'm using the function krige.conv from package geoR but I'm getting a
warning that I do not understand:
Warning messages:
1: NaNs produced in: sqrt(variance)
2: NaNs produced in: sqrt(variance)
I'm performimg an ordinary kriging with a log transform (lambda=0).
Does anyone knows where this comes from ?
Thanks and regards
EJ
2008 Dec 10
2
exporting rast from R to GRASS
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 module r.mapcalc. I can't do it. I tried use
2003 Jul 24
1
geoR size limit problem
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
2003 Jan 17
3
kriging in R
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
2006 Jan 05
1
Memory limitation in GeoR - Windows or R?
Dear Aaron,
I am really a tool user and not a tool maker (actually an ecologist
doing some biostatistics)... so, I take the liberty of sending a copy of
this e-mail to the r-help list where capable computer persons and true
statisticians may provide more relevant information and also to Paulo
Ribeiro and Peter Diggle, the authors of geoR..
I really feel that your huge matrix cannot be