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, c
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
https:/
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 ex
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
> Su
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 51424
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 fo
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,
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 coul
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 t
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 the
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
>image(kr,col="gray"(seq(.7,0,l=10)),xlim=c(-1
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
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
ha
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
14 matches
Mail list logo
