Thank you all for your kind help.
Following examples from the link provided by Edzer, I tested the difference
between:
1. projection specified
coordinates(data.file.dev) = ~LON+LAT
proj4string(data.file.dev)=CRS("+proj=longlat")
vgm1 <- variogram(sqrt(mud)~1, data.file.dev)
model.1 <- fit.variogr
On Tue, 28 Apr 2009, Alina Sheyman wrote:
I'm trying to combine two shapefiles using
zipmaps5 <- spRbind(zipmaps1,zipmaps2)
and getting the following error message
Error in spRbind(as(obj, "SpatialPolygons"), as(x, "SpatialPolygons")) :
non-unique polygon IDs
Does anyone know how I can get ri
I'm trying to combine two shapefiles using
zipmaps5 <- spRbind(zipmaps1,zipmaps2)
and getting the following error message
Error in spRbind(as(obj, "SpatialPolygons"), as(x, "SpatialPolygons")) :
non-unique polygon IDs
Does anyone know how I can get rid of duplicate polygons in a shapefile?
th
While it should be straightforward to calculate great circle
distances, designing a rectangular grid on a sphere is not. There's a
solution for a global grid using hexagons (Teanby, 2006), but I
haven't seen it implemented in R (yet). My gut feeling is that for
Australia there's a standard
Paul, as mentioned here:
https://stat.ethz.ch/pipermail/r-sig-geo/2008-October/004456.html
gstat does use spherical distances when data are in geographical
coordinates, however the usual variogram models are typically not
non-negative definite on the sphere, and no appropriate models are
availabl
In my case, I do not see the reason why the current region and the
raster region should be
different. Everything is within a for() loop, in which I
set the region,
create a MASK,
create a raster with r.mapcalc (hence the same region than the current
region)
import the raster to R
It's not crit
Agustin Lobo wrote:
>
> Great!
> The URL is
> http://n2.nabble.com/R-sig-geo-f2731867.html
>
> but it's readonly (I think), we cannot post from there, Posting would be
> the real advantage
> of using Nabble. Note that anything posted from Nabble goes through the
> list, Nabble
> is just a fron
Great!
The URL is
http://n2.nabble.com/R-sig-geo-f2731867.html
but it's readonly (I think), we cannot post from there, Posting would be
the real advantage
of using Nabble. Note that anything posted from Nabble goes through the
list, Nabble
is just a front end, users have to be subscribed to the l
On Tue, 28 Apr 2009, Agustin Lobo wrote:
Hi!
I'm getting this message at reading a grass raster file:
rastdum <- readRAST6("nomdum",cat=F)
raster map/current region mismatch detected in components:
colsrows origin.northing origin.easting
TRUETRUE
Hi!
I'm getting this message at reading a grass raster file:
rastdum <- readRAST6("nomdum",cat=F)
raster map/current region mismatch detected in components:
colsrows origin.northing origin.easting
TRUETRUE FALSETRUE
set plugin=TR
Hi Jin Li,
I seem to remember that gstat can deal with Great circle distances, but
I could be wrong as I have never used them before. You say that a
certain projection does not produce satisfactory results, how have you
defined satisfactory?
cheers,
Paul
jin...@ga.gov.au wrote:
Dear all,
11 matches
Mail list logo
