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
Paulo Justiniano Ribeiro Jr
2004-Feb-23 13:51 UTC
[R] border of a polygon in contour.kriging - geoR
Dear Patrick border=NULL should work but there is a bug. I'm fixing this right now and a new version will be available later today P.J. On Mon, 23 Feb 2004, Patrick Giraudoux wrote:> 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 > > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://www.stat.math.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html > >Paulo Justiniano Ribeiro Jr Departamento de Estat?stica Universidade Federal do Paran? Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 361 3471 Fax: (+55) 41 361 3141 e-mail: pj at est.ufpr.br http://www.est.ufpr.br/~paulojus