Hi!
If you have a Spatial Points Data Frame, you can use the following function
to create a grid:
create1 <- function(obj) {
# Function that creates a new_data object if one is missing
convex_hull = chull(coordinates(obj)[,1],coordinates(obj)[,2])
convex_hull = c(convex_hull,
how to create prediction grid like meuse grid
*Thanks & Regards*
*C Dineshkumar*
*Bsc Agriculture*
*M.Tech Remote Sensing and GIS*
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
Hello,
I have a problem with the kriging function. I get this error:
Error in solve.default(matrix(A, n + 1, n + 1)) :
system is computationally singular: reciprocal condition number =
5.53457e-17
But it has worked with other data and now I don't understand why this
error.
Here is a data
Try initial fitting with fewer lags; I get a fit to your data with lags=3.
S Ellison
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Diego Ubuntu
Sent: 06 July 2013 16:09
To: r-help@r-project.org
Subject: [R] Kriging
I'm trying to get a kriging surface using these data, but I only get this
error message. Any ideas on why or how to solve it.
d
x y e
1 551595.2 18040622.0
2 599591.7 18208232.0
3 615604.7 18209032.0
4 612337.6 17709892.1
5 614203.6 17618163.0
6
On 10/05/2011 09:07 AM, Leynnard Rey Matillano wrote:
Hi! Im new to R and I need to interpolate a shapefile using kriging. I've
been able to plot/read the shapefile using the package maptools or rgdal.
I've searched the internet for sample codes but most of the kriging codes
that I've found
Hi! Im new to R and I need to interpolate a shapefile using kriging. I've been
able to plot/read the shapefile using the package maptools or rgdal. I've
searched the internet for sample codes but most of the kriging codes that I've
found done in R is done using txtfiles or CSVs. An example
It is perhaps easier than you think, see the example in gstat:
library(gstat)
?krige
After this is run meuse is a SpatialPointsDataFrame:
coordinates(meuse) = ~x+y
Any point shapefile read with readOGR from rgdal (or the alternative
functions in maptools) will also be SpatialPointsDataFrames,
Please consider using the R-sig-geo list for kriging questions.
In fact, using str() would show that your z is badly impacted by your misuse
of the c() function - this was almost certainly not what you wanted to do.
Constructing your coordinates by using seq() is always a bad idea, as you
need
Hi everybody,
I have a longitude vector and a latitude one. Associated to these coordinates,
i have a matrix with some data at some coordinates but not all.
Lon - seq(136.025,144.975,0.05)
Lat - rev(seq(-66.975,-65.525,0.05))
dim(z) - c(Lon,Lat)
And i have tried to apply to these
hi all.
If someone have the same problems this is the answer:
-vertical legend :
legend.krige(x.leg=c(X,X),y.leg=c(X,X),kr$pred,vert=TRUE,col=gray(seq(.7,0,l=10)))
-sample's positions on the map:
###coords.dat=table$coords### like in
Hi all,
I have to draw a map of microorganisms repartition in a sample. I use the
image function of geoR package to build the map but there are three
problems:
1) the legend is Horizontal (with or without vert.leg=TRUE)
2) I want to plot the position of organisms and level curves of density on
Hei,
I have two spatial datasets Sa and Sb, both with lat-lon coordinates and
from same geographic area, but from different localities within the area
(independent samples). Sa is biotoc data, Sb is some environmental
parameter (fertility). I 'know' that Sb affects Sa, but wonder on which
scale. I
Hi Robert,
A package I use for spatial statistics is gstat. It supports a method
called co-kriging that seems to fit your needs. To get a good overview
of the spatial capabilities of R a good place to start is the Spatial
Task View (http://cran.r-project.org/web/views/Spatial.html). If you
14 matches
Mail list logo