Dear useRs,

I am trying to convert the predictions of a kriging model into a  
shapefile but I am getting some errors and I am getting nowere with  
my solutions...

This is the code I am running:

a <- data.frame(Id=seq(1,length(pred.grid[,1]),1),X=pred.grid[, 
1],Y=pred.grid[,2])
a_dbf <- data.frame(Id=seq(1,length(pred.grid[,1]),1), data=kc2$predict)
shp_1 <- convert.to.shapefile(a, a_dbf, field="Id", type=1)

were pre.grid is determined by this code:

min_x <- 142794
max_x <- 152121
min_y <- 485080
max_y <- 508887
pred.grid <- expand.grid(seq(min_x, max_x, 50), seq(min_y, max_y, 50))

and kc2 is an object returned by the aplication of the function  
krige.conv:

kc2 <- krige.conv(dados_g, loc=pred.grid, krige=krige.control 
(obj.m=vario_fit2))




Another question:

as you can see I am predicting my values to a square set of points...  
my question is how can I can I generate a set of point that match a  
specific area (because I have a specific limite(area) in an esri's  
shapefile).



Thanks in advance.
Carlos



Carlos GUERRA

Gabinete de Sistemas de Informacao Geografica
Escola Superior Agraria de Ponte de Lima
Mosteiro de Refoios do Lima
4990-706 Ponte de Lima

Tlm: +351 91 2407109
Tlf: +351 258 909779

Be a Mac user... update your self!


        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to