> Hi all,
>
> I am new to R and would like to get some advice. I performed kernel
> density estimation using spatstat. The point data (i.e. drawing) is
> in ESRI shapefile format. It was read in using "readShapePoints" of
> maptools. After the analysis, the result was exported using
> "writeA
AFAIK, 4326 is the PROJ.4 code for "+proj=longlat". You should find that
proj4string(data)
shows you something like that string above, and
summary(data)
will report that it is not projected.
In general, please provide at least the results of summary on your data
object/s.
BTW, just because
You ought to be able to read that using readGDAL (rgdal) which will
return you a SpatialGridDataFrame. You then have methods for image and
spplot that will work:
## NOT TESTED:
library(rgdal)
asc <- "file.asc" ## or whatever it is
g <- readGDAL(asc)
summary(g)
image(g)
## OR
spplot(g)
Just want to follow up from my earlier post
Since the voronoi polygon layer and the point layer share the same study area,
is there a way to do a 'spatial joint' such that I can find out the points fall
into which polygon?
And to have a dataset like the following:
Caseid xc yc voronoi.pol
Dear R-users,
I have a dataset with over 5000 pts and I know their x and y coordinates.
Based on the study area, I created 10 random points and used these points to
create voronoi to dissect the space.
genxy.vm<-voronoi.mosaic(genxy$x, genxy$y, duplicate="remove")
genxy.vp<-voronoi.polygons(gen
sorry for spamming all of you with my problems...but here's another one:
I have ascii grids from ESRI ArcGIS in the following format:
ncols 42
nrows 36
xllcorner -10.311630235357
yllcorner 35.110501339805
cellsize 1
NODATA_value -
Is there an already existing com
Hi there
I have an ordinary kriging result (gstat), which is in a projected
coordinate system of course (I used
proj4string(data)=CRS("+init=epsg:4326"). But I will need the map in
geographical coordinates (lon, lat...) for visualisation. How can I
back-project the map?
Thanks
Thomas
__
Thomas