search for: krig

hello I try to build DEM using Krig function of fields package. And I get this error message. here is the procedure I followed: > dt<-read.table("/usr/local/bartin/stat/topostat1",header=F,sep="|") > names(dt) <-c("x","y","z") > coord<-cbind(dt$x,dt$y) > elev...

Hi, I have a 4 column file (long/lat/elev/variable) and I tried to fit the values of my variable to the XYZ space using Tps and I keep getting the following message: Warning messages: 1: GCV search gives a minumum at the endpoints of the grid search in: Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, 2: GCV search gives a minumum at the endpoints of the grid search in: Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV.one, Krig.fgcv.one, 3: Failed to find the upcrossing in: find.upcross(Krig.fs2hat, info, upcross.level = shat.pure.error^2...

I'm new at this, so please bear with me. I am trying to get Krig (actually tps) results into a data frame at a set of points. Frankly, all I really need is to query the Krig on a specific XY pair and return a result. I can carry the rest. I am using tps from the fields package. Here's the code I use to create the Krig object: =========== Snip! ==========...

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Mike J Smith > Sent: 09 June 2005 09:58 > To: r-help at stat.math.ethz.ch > Subject: [R] krig.image help > > > Hi > > I have recently been experimenting with the use of kriging, primarily > through Goldensoftware's Surfer. Having had a look around at the > packages available for R I note a couple that are available so have > initially added "fields"....

...files installing data files installing man source files installing indices not zipping data installing help >>> Building/Updating help pages for package 'fields' Formats: text html latex example BD text html latex example Krig text html latex example Tps text html latex example US text html latex example US.dat text html latex Wtransform.image text...

...imply cannot get the black points ("+") to display the actual value. I'm guessing the problem is somewhere in the second to last line of code (starts with "pts<-..."). I have attached the two text files needed to run the code. http://www.nabble.com/file/p21734824/R_Hk_Krig_File_Log.txt R_Hk_Krig_File_Log.txt http://www.nabble.com/file/p21734824/Predict_Location_XY.txt Predict_Location_XY.txt Eric library(gstat) K.dat<-read.table("C:/temp/R_Hk_Krig_File_Log.txt", header = TRUE) my.dat<-data.frame(K.dat) attach(my.dat) coordinates(my.dat)=~X+Y pred...

....2/../../.. -lfrtbegin -lg2c -lm -lgcc_s ** R ** data ** help >>> Building/Updating help pages for package 'fields' Formats: text html latex example BD text html latex example Exponential text html latex Krig.Amatrix text html latex example Krig text html latex example Krig.engine.default text html latex example Krig.null.function text html latex RMprecip text...

Dear list, I get the following error message when installing the package GRASS. missing link(s): interp.new missing link(s): vgram Krig predict.Krig predict.se.Krig See also below for the full installation message. What does this mean exactly? Regards, Ulrich ------------------------------------------------------------------ > R CMD INSTALL --configure-args=--grass5=/usr/local/grass5/grass5bin/grass5 -l ~/.lib/R/library/...

I have a question to everybody. After log10 transfprmation, I have done ordinary kriging in gstat in R? I need to back trnasform the prediction values to orgiginal scale. How do I do this in gstat in R? Thanks Zia -- Zia Uddin Ahmed 915 Brad Field Hall Department of Crop and Soil Cornell University Ithaca NY 14850 USA

...nal login failure: rjhjr [127.0.0.1] Jan 13 18:51:02 kongemord dovecot: child 1320 (auth) killed with signal 6 messages Jan 13 18:51:02 kongemord kernel: pid 1320 (dovecot-auth), uid 1006: exited on signal 6 debug.log Jan 13 18:51:02 kongemord slapd[1103]: conn=14 op=1 SRCH base="ou=imap,dc=krig,dc=net" scope=2 deref=0 filter="(&(objectClass=posixAccount)(uid=rjhjr))" Jan 13 18:51:02 kongemord slapd[1103]: conn=14 op=1 SRCH attr=uid userPassword Jan 13 18:51:02 kongemord slapd[1103]: conn=14 op=1 SEARCH RESULT tag=101 err=0 nentries=1 text= Jan 13 18:51:02 kongemord slap...

Hi, I am trying to use the Tps procedure to visualize 3 variables in XYZ space, but keep getting the error messages below: Warning messages: 1: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, : GCV search gives a minimum at the endpoints of the grid search 2: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV.one, Krig.fgcv.one, : GCV search gives a minimum at the endpoints of the grid search This is the code I am using: tp...

A majority of my data makes a kriged map perfectly using an exponential model for the semivariogram to fit my data and then going through the commands variofit() to define the model and then krige.conv() to use the model to predict values in a grid. But?one set of my data appears to be linearly correlated for the first 5000 met...

Hi ! Sorry for my ignorance. I have two questions: 1) which is in R the function to make a covariogram of spatial data ? 2) does anyone know if there exist the possibility in R of performing kriging with an arbitrary covariance function ? If I have well understood the Krig procedure offers only three possibilities (expcov, gauscov and sphercov) that are the common models for stationary and isotropic processes. Thanks, Paola. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

...s. I include the data file and the syntax file for reference. Again, if I use the fields function, as is, I get the message: Error in svd(tempM) : error 159 in dsvdc using traceback, I get: > traceback() 4: stop(paste("error ", z$info, " in dsvdc")) 3: svd(tempM) 2: Krig(x, Y, cov.function = rad.cov, m = m, decomp = decomp, scale.type = scale.type, outputcall = Tpscall, p = p, ...) 1: Tps(bvolcap, bdsm) > if I change the occurrence of svd in the fields package to La.svd, I get the error message > bout <- Tps( bvolcap, bdsm) Error in "[&...

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. What RAM capacity should i look for to do a 150 x 250 cell gri...

Hello to everybody ! I don't have a question directly connected with R, but it takes starting point from the use of the R function Krig for kriging interpolation and from the idea of making a different implemetation of the function for my purposes. Doas anyone know where can I find free C or C++ codes of kriging ? As always, thanks a lot in advance, Paola. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

R Experts, I'm looking for some help with the geoR package. I'm trying to krig some data without using a global neighborhood. I would like to set my moving neighborhood to a distance, say 100 meters, where I know my data is spatially correlated. I have tried the ksline function, but that only allows my moving neighborhood to be set to a number of data points. But, si...

...0,3), ylim=c(0,3), col="blue", pch=19) I would like to interpolate the quantities on grid1 (red) and grid2 (blue) on the same grid (black). This scenario is very common in geophysical data and models. I only found: - functions in package akima which are designed for irregular grids - krigging in package fields, which also requires irregular spaced data - approx or spline which works in 1D and which I could apply line by line and column by column and use a mean of both estimates I am sure there are plenty of functions already available to do this but searching R-help and the pack...

...5,2 15,6 43,4 64,5 15 0,2 0 575091 4089921 183 1,2 13,2 64,6 70,6 16,2 0,2 3,4 The area where this data must be interpolated in the natural park with an irregular shape inside of a grid of 2190 lines and 2282 columns with utm coordinates. I have tried "Krig" from "fields" package. A created a matrix with independant data (a) and a vector with dependant data (b). a<-matrix(c(clima$X,clima$Y,clima$Altitud),,3) b<-clima$PrepJul Interpolation<-Krig(a,b) mapa<-predict.surface(interpolacion, nx=2190, ny=2282) str(mapa) num [...

...p-OP2003Seppos.AICsp #[1] -4 (OP2003Seppos.lf0.p<-exp(OP2003Seppos.lf0$beta)/(1+exp(OP2003Seppos.lf0$beta))) #P non spatial #[1] 0.9717596 (OP2003Seppos.lf.p<-exp(OP2003Seppos.lf$beta)/(1+exp(OP2003Seppos.lf$beta))) #P spatial #[1] 0.9717596 It must what have an important influence at kriging, because it shows as following: OP2003Sepposbin.krig<-glsm.krige(OP2003Seppos.tune,loc=OP2003Seppospro.pred.grid,bor=OP2003Sepposbor) #glsm.krige: Prediction for a generalised linear spatial model #There are 50 or mode advices (use warnings() to see them) #> warnings() #Warning messages...