Multidimensional Scaling (MDS) is certainly an option to be considered.
Its suitability depends, of course, on the specific application.
In my case, which is conceptually similar to Martin's, that option was
discarded.
And I think I can explain why whith this single image:
Tugores Ferra pilar.tugo...@ba.ieo.es wrote:
From: Pilar Tugores Ferra pilar.tugo...@ba.ieo.es
Subject: Re: [R-sig-Geo] kriging as fish swim, not as crows fly
To: Martin Renner martin.ren...@stonebow.otago.ac.nz
Cc: r-sig-geo@stat.math.ethz.ch
Date: Wednesday, January 27, 2010, 5:35 PM
Hi, Martin
Ooops, you're right, sorry.
This shows that nobody tried to use it before, haha :)
Now I have uploaded the source files
(http://www.geeitema.org/guenmap/index.jsp?opcion=resultadosidioma=en).
Please let me know if you have any trouble. I have just tested it in
Linux. It compiled and seem to
Renner
Sent: Wednesday, January 27, 2010 1:07 AM
To: r-sig-geo@stat.math.ethz.ch
Subject: [R-sig-Geo] kriging as fish swim, not as crows fly
Hi All,
I want to kirg fish and seabird densities within an estuary which has
several arms. Since neither organisms cross land, the appropriate
, 2010 1:07 AM
To: r-sig-geo@stat.math.ethz.ch
Subject: [R-sig-Geo] kriging as fish swim, not as crows fly
Hi All,
I want to kirg fish and seabird densities within an estuary which has
several arms. Since neither organisms cross land, the appropriate
distances would not be euclidian
@stat.math.ethz.ch
Asunto: Re: [R-sig-Geo] kriging as fish swim, not as crows fly
Not kriging as such, but check out the soap-film smoothing in package mgcv:
http://www.maths.bath.ac.uk/~sw283/simon/papers/soap.pdf
FWIW, there are binning methods with MCMC in the package
tripEstimation that have similar
Hi Facu,
great, this is just what I have been looking for. You say that your
modifications to geoR are open source. On your website I found the geoR.dll
file, but no source code. The .dll is of little use platforms other than
windows. Is there any way I could get the source?
Thank you for
Hi All,
I want to kirg fish and seabird densities within an estuary which has several
arms. Since neither organisms cross land, the appropriate distances would not
be euclidian but over-water (as fish swim). There are several papers,
describing this problem and how to deal with it (see below),
Tobin Cara wrote:
Hello,
I have recently read an interesting article about integrating Limited Area
Models (LAMs) into kriging with external drift for temperature (Libert� et al.
link below).
www.wmo.int/pages/prog/www/IMOP/.../IOM.../P2(05)_Perini_Italy.doc
This URL is incomplete. Please
Hello,
I have recently read an interesting article about integrating Limited Area
Models (LAMs) into kriging with external drift for temperature (Libertà et al.
link below).
www.wmo.int/pages/prog/www/IMOP/.../IOM.../P2(05)_Perini_Italy.doc
As I understand, it seems that the authors generated
Hi Greg,
Variogram modelling is slower with large data sets, but 8-10.000
observations should not be a problem, unless you need the results
extremely fast. On my computer (3 years old) it takes about 4 seconds
with 8.000 random observations, using the variogram function in gstat.
Time
Hi, I have a very non-specific question about the number of sample points that
can be used
for developing experimental variograms and kriging in R/gstat/etc. Does anyone
have experience
using a very large number of data points in the kriging processes with R? It
has been a few
years since my
Hello,
I am taking the log of precipitation values and therefore many are now NA
values. I want to continue to krig my precipitation matrix.
Is there a way to ignore these values with kriging.
My attempt with is.nan still gives:
Erreur : dimensions do not match: locations 105 and data 12
If you are interpolating precipitation you probably should not ignore
the zeros. If you want to log transform your values, perhaps you can
use log(x+1) instead of log(x). Robert
On Tue, Nov 17, 2009 at 10:43 AM, Tobin Cara cara.to...@epfl.ch wrote:
Hello,
I am taking the log of precipitation
Hello,
Thank you all for your previous help. I began using R 2 weeks ago, and I am
getting somewhere finally. I have been able to run universal kriging with a
Digital Elevation Model trend. Now, does anyone have experience with kriging
more than one trend variable?
I assume you have to have
Tobin Cara wrote:
Hello,
Thank you all for your previous help. I began using R 2 weeks ago, and I am
getting somewhere finally. I have been able to run universal kriging with a
Digital Elevation Model trend. Now, does anyone have experience with kriging
more than one trend variable?
I
Dear all,
I have borehole data from which I want to interpolate the basis of a
faulted sandstone formation using kriging. The structural analysis
clearly shows the effect of faults with a directional variogram that is
nonstationary perpendicular to them, but I am a bit at a loss as to how
Hi all
Please help and correct me, to predict Z(S0) at (0.5,0.5) given Z(S1) = 3 at
(0,0), Z(S2) = 5 at(0,1), Z(S3) = 6(1,0) and Z(S4) = 4(1,1). Let
γ (h)=h^2 h1
= 1 h=1 using geo.
I did this: coords-matrix(c(0,0,1,1,0,1,0,1), nrow=4, ncol=2)
data-c(3,5,6,4)
Hi Kabeli,
I never saw Brian Vinyard slides 40! :-)
It is accessible on a internet site?
bests
milton
On Tue, Aug 25, 2009 at 6:56 PM, KABELI MEFANE kabelimef...@yahoo.co.ukwrote:
Hi all
Please help and correct me, to predict Z(S0) at (0.5,0.5) given Z(S1) = 3
at (0,0), Z(S2) = 5
Maybe something like this:
http://www.ars.usda.gov/sp2UserFiles/ad_hoc/1200SpatialWorkshop/01VinyardOverview.pdf
- Original Message -
From: milton ruser milton.ru...@gmail.com
Date: Tuesday, August 25, 2009 7:50 pm
Subject: Re: [R-sig-Geo] Kriging
Hi Kabeli,
I never saw Brian
On Sunday 21 June 2009, Ebrahim Jahanshiri wrote:
It has been long that I wanted to suggets this for automatic trend
detection based on our previous conversations with Edzer and Anne. I
found two ways that seem to be reasonable and have potential for
automizing the trend detection ( I got
It has been long that I wanted to suggets this for automatic trend
detection based on our previous conversations with Edzer and Anne. I
found two ways that seem to be reasonable and have potential for
automizing the trend detection ( I got these from my conversations
with Margaret Oliver and Dick
Hello all,
I have been using gstat to krige temperature data with elevation as
external drift.
feb01.meantemp.krig-krige(MEANTEMP~elevation, locations=
feb01.meantemp, newdata=elevation, model=feb01.meantemp.r)
I keep encountering the warning message:
Error: cannot allocate vector of size
Hi Paulo,
you suggest to add conversion functions to geoR later on, but I (as the
average lazy user) would ask the following: if in the following geoR call
kc.s - krige.conv(s100, loc=gr.s, krige=krige.control(obj=ml.s))
gr.s is of class SpatialPixels[DataFrame] or SpatialGrid[DataFrame],
Hello,
After kriging with the use of the krige.conv function, I would like to
export my result under the grid format. In order to view it into a GIS.
It was suggested to me to use the writeGDAL function into the rgdal
package. However It looks like I have a problem of supported formats
Oceanográfico de Baleares
Muelle de Poniente s/n
07015 Palma de Mallorca (España)
Tel.: (34) 971 401561
-Mensaje original-
De: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] En nombre de Ashton Shortridge
Enviado el: 27 August 2008 22:22
Para: r-sig-geo@stat.math.ethz.ch
Asunto: Re: [R-sig-Geo
Thanks Edzer,
I've requested Cressie's book from our library (just waiting on it).
My main concern was the many 0 counts. I also was not enthusiastic about
odd transformations which then require appropriate back-transforms (I
imagine the back transform of the kriging variance gets messy)
Dave,
Transformation to a continuous distribution when the data follow a
discrete distribution is always messy, and the back-transform may get worse.
While you're at the library, try to pick up Diggle Ribeiro's
Model-based geostatistics; they describe a model-based approach that
extends
a
copy of your lecture notes?
Thanks,
Jin
-Original Message-
From: Hengl, T. [mailto:[EMAIL PROTECTED]
Sent: Monday, 7 July 2008 6:23
To: Li Jin
Cc: r-sig-geo@stat.math.ethz.ch
Subject: RE: [R-sig-Geo] kriging [SEC=UNCLASSIFIED]
Jin,
Do not get me wrong. I support your effort
everyone interested in RK or the like.
Best regards,
Jin
-Original Message-
From: Hengl, T. [mailto:[EMAIL PROTECTED]
Sent: Friday, 4 July 2008 5:55
To: Li Jin
Cc: r-sig-geo@stat.math.ethz.ch
Subject: RE: [R-sig-Geo] kriging [SEC=UNCLASSIFIED]
Dear Jin,
I really think that this list
of
Agricultural Sciences immediately and delete this email.
Fra: [EMAIL PROTECTED] på vegne af Edzer Pebesma
Sendt: to 03-07-2008 13:33
Til: Hengl, T.
Cc: r-sig-geo@stat.math.ethz.ch; Dave Depew
Emne: Re: [R-sig-Geo] kriging
Hengl, T. wrote:
I agree with Paulo - gstat
: Re: [R-sig-Geo] kriging
Hengl, T. wrote:
I agree with Paulo - gstat can work with any linear model including the
transforms of the original predictors e.g.:
Z ~ X + X^2 + Y + Y^2etc.
The problem is that gstat implements the so-called
Kriging-with-external-trend algorithm
Dear geo experts,
has anyone looked into kriging of spatial processes that do not live in a 2D
continuum but instead are constrained to a network/graph (e.g. a street
grid) ?
Clearly, distances need to be redefined but more than that, the covariance
matrix is a very different animal.
Thanks!
Dear Dave,
I separate fitting of the deterministic (trend) and residual part of the
universal kriging model all
the time. Adding OK of residuals to the trend is fine, as long as the
regression model is estimated
using GLS (but many do it even if they use only OLS; the difference is often
hi Tom,
That was my impression from reading some introductory texts... I'll have
to see how the mgcv package fits the polynomial function to the
data...it isn't clear to me at first glance how it is accomplished.
Many thanks for your advice.
Dave
Tomislav Hengl wrote:
Dear Dave,
I separate
=EUB:NOTICE:LBNA22904:EN:HTML
-Original Message-
From: [EMAIL PROTECTED] on behalf of Dave Depew
Sent: Mon 6/16/2008 10:54 PM
To: Paulo Justiniano Ribeiro Jr
Cc: r-sig-geo@stat.math.ethz.ch
Subject: Re: [R-sig-Geo] kriging
Ok,
What about higher order polynomials? I have fitted one using a gam
Thanks Tom,
I've been able to fit a polynomial function to the data quite well. The
residuals are behaving (i.e normal distribution and no skewness of
variance). I'm assuming this means that I could krige the residuals
(Ordinary K?) and then add the trend back to the predicted residual
grid?
Hi all,
I have a data set that I would like to krige to interpolate between
transects. There is a non-linear trend between two of the variables...my
impression from reading the gstat help file is that there must be a
linear relationship between the data to use universal kriging?
Second, would
Dave,
what is necessary for UK is a relation expressed by a linear model, not
necessaraly a linear relation between the variables.
e.g. you could have a second degree polinomial and still work within the
scope of universal kriging.
On Mon, 16 Jun 2008, Dave Depew wrote:
Hi all,
I have a data
Ok,
What about higher order polynomials? I have fitted one using a gam to
the data which which helps to normalize the residuals, and reduce the
variance of the residuals.
Is it simply a matter of plugging in the function into the gstat command
line? Or is it simpler to krig the residuals and
On Mon, 5 May 2008, Edzer Pebesma wrote:
Dave Depew wrote:
Thanks,
This worked.
I'm still confused why the if else statement didn't work...
If one wanted to do conditional arithmetic would a for statement bee
needed?
e.g.
meuse.grid[[class]] = for(i in 1:length(meuse.grid[[dist]])){
Dave Depew wrote:
Hi all,
I've got a question regarding kriging outputs. I have an interpolated
dataset which due to the nugget effect contains some negative values
as the predictions. I would like to truncate these @ 0, rather than
having them as a negative prediction.
I've tried something
Thanks,
This worked.
I'm still confused why the if else statement didn't work...
If one wanted to do conditional arithmetic would a for statement bee needed?
e.g.
meuse.grid[[class]] = for(i in 1:length(meuse.grid[[dist]])){
if (meuse.grid[[dist]]0.5)
Dear members,
I have tried to export a kriging map to arcgis as asciigrid or image.
I have used the functions write.asciigrid and writeRast6sp(grass), in
both cases any success; In the former when exporting it, I got the
following message Asciigrid does not support grids with non-square
cells. I
analysis, properly conducted, is a delicate dissection of
uncertainties, a surgery of suppositions. ~M.J.Moroney
-Oorspronkelijk bericht-
Van: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Namens Robert Helber
Verzonden: donderdag 26 juli 2007 22:25
Aan: R geo
Onderwerp: [R-sig-Geo
Hi,
i've tryed to delete the \ backslash,
now i've a different error, at the same line :
grd - GridTopology(cellcentre.offset=c(G$west+(G$ewres/2), G$south
+(G$nsres/2)), cellsize=c(G$ewres, G$nsres), cells.dim=c(G$cols, G
$rows));
Errore in validObject(.Object) : invalid class GridTopology
On Sat, 20 Jan 2007, epifanio wrote:
Hi,
i've tryed to delete the \ backslash,
now i've a different error, at the same line :
grd - GridTopology(cellcentre.offset=c(G$west+(G$ewres/2), G$south
+(G$nsres/2)), cellsize=c(G$ewres, G$nsres), cells.dim=c(G$cols, G
$rows));
Errore in
hi i've some problem to do a tutorial on the kriging
interpolation
i found instruction on how to interpolate the srtm data to
increase the resolution :
http://grass.itc.it/newsletter/GRASS_OSGeo_News_vol4.pdf
at page 20 ... at the line :
grd -
On Fri, 19 Jan 2007, epifanio wrote:
hi i've some problem to do a tutorial on the kriging
interpolation
i found instruction on how to interpolate the srtm data to
increase the resolution :
http://grass.itc.it/newsletter/GRASS_OSGeo_News_vol4.pdf
at page
49 matches
Mail list logo