Hi list,

Well, this time I’ve a doubt with mapping generation.

I was already able to read and plot shapefiles, plot point on this map. All
this without any problems, but now I want to do something that I think, only
Golden Software Surfer is capable of.

I would like to plot a kriging result on the existing map (map script
below). Well I looked for on the net, but I was not able to be sure of which
packages I’ll need to perform the kriging (maybe geoR?) and then plot it in
the area map.

Is R able to do that? What packages I need to know to do so?

I just need a starting point to go forward.

 

I need to generate something like this:

http://i288.photobucket.com/albums/ll162/Godrigos/Kriging.jpg

(but without the top right area, just the different grayscale kriging area).


 

library(xlsReadWrite)

library(maps)

library(mapdata)

library(maptools)

Pontos<-read.xls('Pontos.xls',sheet=1,rowNames=T)

Prof<-read.xls('Pontos.xls',sheet=2,rowNames=T)

Batimetria<-readShapeSpatial(‘Batimetria.shp')

Municipios<-readShapeSpatial('Municipios.shp')

Recifes<-readShapeSpatial('Recifes.shp')

par(mar=c(4,0,0,0))

map('worldHires','brazil',ylim=c(,),xlim=c(,),type='n')

plot(Batimetria,ylab='',xlab='',border='lightgray',add=T)

plot(Recifes,ylab='',xlab='',lty='dashed',add=T)

plot(Municipios,ylab='',xlab='',col=rgb(229,229,229,max=255),add=T)

axis(1,xaxp=c(,,))

axis(2,yaxp=c(,,))

points(Pontos$long,Pontos$lat,pch=20)

text(Pontos$long,Pontos$lat,rownames(Pontos),pos=1,cex=0.5)

text(Prof$Long,Prof$Lat,rownames(Prof),col='darkgray',cex=0.8,srt=-24)

 

Ps.: Sorry for hiding the maps coordinates. I had to, a demand of my
coworkers.

 

Thanks in advance.

___________________________________
MSc.  <mailto:r.aluizio@gmail.com> Rodrigo Aluizio
Centro de Estudos do Mar/UFPR
Laboratório de Micropaleontologia


	[[alternative HTML version deleted]]

