[R] border of a polygon in contour.kriging - geoR

Patrick Giraudoux patrick.giraudoux at univ-fcomte.fr
Mon Feb 23 13:43:53 CET 2004

Dear all,

When a conventional kriging and then a contour plot is limited with a polygon (as possible with krige.conv and contour.kriging), the
polygon border is displayed in black by default.

> kc<-krige.conv(CZdata,loc=pred.grid,borders=czpoly,krige=krige.control(obj.m=ls))
> plot(CZcoord,xlab="x",ylab="y",type="n",asp=1)
> contour(kc,loc=pred.grid,add=T)

I would like to be capable to get the contour lines only (limited by the polygon, but without the polygon border displayed). I tried
to pass the argument border=NULL, unsuccessfully:

> contour(kc,loc=pred.grid,add=T,border=NULL)
Error in contour(as.double(x), as.double(y), z, as.double(levels), labels,  :
        dimension mismatch
In addition: Warning message:
Replacement length not a multiple of the elements to replace in matrix(...)

Any hint?

Thanks in advance,

Patrick Giraudoux

More information about the R-help mailing list