Hi,I am sorry for not having control checked the code before posting (in my earlier mail). I have done this now. I just have a simple question now. Why is the overlay plot (the final plot ) not falling into the intended position (as defined by the viewport vp2)?#Setting up the initial example raster image library(grid) x <- y <- c(-15,seq(-4*pi, 4*pi, len=27),15) r <- sqrt(outer(x^2, y^2, "+")) z <- cos(r^2)*exp(-r/6) image <- (z - min(z))/diff(range(z)) xrange <- range(x) yrange <- range(y)
vp1<-viewport(x = unit(0.5, "npc"), y = unit(0.5, "npc"), width = unit(1, "npc"), height = unit(1, "npc"), default.units = "npc", just = c("left","bottom"),name="vp1") pushViewport(vp1) plot(x, y, ann=FALSE,xlim=xrange, ylim=yrange,xaxs="i", yaxs="i") rasterImage(image,xrange[1], yrange[1],xrange[2], yrange[2],interpolate=FALSE) lines(c(-10,4),c(-12,-12),col="red",lwd=2) lim<-par("usr") lim limplt<-par("plt") limplt #par(new=TRUE,plt=c(limplt[1],limplt[2],limplt[3])) #convertX(unit(0:1,"npc"),"native") x1<-convertX(unit(-10,"native"),"npc") x2<-convertX(unit(5,"native"),"npc") y1<-convertY(unit(-12,"native"),"npc") y2<-convertY(unit(3,"native"),"npc") wx<-(5-(-10))/(15-(-15)) wy<-(3-(-12))/(15-(-15)) xw<-convertWidth(unit(wx,"native"),"npc") yw<-convertHeight(unit(wy,"native"),"npc") vp2<-viewport(x = unit(x1, "npc"), y = unit(y1, "npc"), width = unit(xw, "npc"), height = unit(yw, "npc"), default.units = "npc", name="vp2") pushViewport(vp2) #data with a different scale xs<-seq(0.01,0.1,0.01)ys<-xs^2 #points(xs,ys,type='b',col='red',newpage=FALSE) plot(xs,ys,type='b',col='red') #fills the whole screen instead of just the intended viewport Thanks,Ravi [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.