Hallo everybody.
Using the krige() function of the gstat package for ordinary and universal
kriging and the predict() function of the same package for ordinary and
universal cokriging I obtain the prediction of the variable and also the
variance of prediction. Is possible to do the same with the ra
Hallo everybody.
Using the krige() function of the gstat package for ordinary and universal
kriging and the predict() function of the same package for ordinary and
universal cokriging I obtain the prediction of the variable and also the
variance of prediction. Is possible to do the same with the ra
Good point, thanks!
zerodist computes Euclidian distances, meaning you should use have
projected coordinates to do something useful at the poles.
I'll put it on the TODO list to compute great circle distances when data
are in long lat, hopefully for the next sp, which might come out pretty
soon.
Hi everyone,
I would like to know how the distance for "zero" parameter in zerodist/remove.duplicates function (sp package) should be entered when coordinates for SpatialPoints are in lat/long. Function documentation has no mention about
you have to use either map units (decimal degrees for
Sorry ... my mistake ... the rectangle is slightly rotated ... here is the
difference.
Arnaud
2013/6/26 Arnaud Mosnier
> Hi,
>
> Maybe I made an error somewhere, but I can't find where !
> There is differences when I calculated the area of a rectangular polygon
> and the value calculated by th
Hi,
Maybe I made an error somewhere, but I can't find where !
There is differences when I calculated the area of a rectangular polygon
and the value calculated by the function Polygon.
Sorry for the large precision numbers, but I tested if it was due to a
rounding error.
Hope you can enlighten m
On Wed, 26 Jun 2013, Renger van Nieuwkoop wrote:
Hi
I tried running the following code using the shapefiles for the city of
Zurich
require("rgdal") # requires sp, will use proj.4 if installed
require("maptools")
zh = readOGR(".",layer="Zurich")
and got following error:
OGR data source with
In general, no you cannot say that this would be correct. It depends
on where and how extensive the region is, so if you can tell us that I
could put together something to show how valid it would be. You should
reproject to a local equal area projection and then it's going to be
valid, but even the
How can I do boxcox transformation for explanatory variables and use it in a
multiple linear regression?
Kind regards,
ISRAEL O. IKOYI
MSc. Thesis Research Student
Wageningen University
Department of Soil Quality
P.O. Box 47
NL-6700 AA Wageningen
The Netherlands
___
Dear All,
I try to plot some variogram models in a panel by show.vgms in gstat.
my program is:
show.vgms(min =0, max =20, n = 50, sill =c(v.ols$psill[2],v.wls$psill[2],
v.gls$psill[2],v.rml$sigmasq,v.reml$sigmasq,v.opvme$psill[2],v.opv$psill[2],
v.opvab$psill[2],v.opvlin$psill[2]/3),range =c(
Hi,
Good news.
For me, what was missing for the moment with spatial tasks was something
that could coerce linear spatial objects to (i)graph objects so as to pull
out all the potentiel of igraph package.
I can remember some attempts by Barry Rownlingson with a buildTopo function
http://rpubs.com
Hi
I tried running the following code using the shapefiles for the city of Zurich
require("rgdal") # requires sp, will use proj.4 if installed
require("maptools")
zh = readOGR(".",layer="Zurich")
and got following error:
OGR data source with driver: ESRI Shapefile
Source: ".", layer: "Zurich"
wi
12 matches
Mail list logo