[STATSGRASS] Problems when overlaying plots

Raffaele Morelli raffaele.morelli at gmail.com
Fri Oct 20 10:56:13 EDT 2006


Hello.

I am running grass & R and want to overlay maps and graphs created using
vector attributes stored in a MySQL database.

This is the essential script I am working with

...
province<-readVECT6("M_Province",remove.duplicates=FALSE,ignore.stderr=TRUE)
sites<-readVECT6("M_Province",type="centroid",remove.duplicates=FALSE,
ignore.stderr=TRUE)

par(fig=c(0, 1, 0, 1))
plot(c(0,10), c(0,10), type="n", xlab="v", ylab="c",axes=T )
box()

par(new=T)
plot(province)

par(new=T,oma=c(0, 2, 0.5, 1.5))
plot(sites)


After execution the sites vector is wrongly displayed, meaning that
centroids are sligthly translated from the original position.
Neverthless they are right when using points(sites) instead of plot(sites).

Am I missing something?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-stats/attachments/20061020/bb9ce745/attachment.html


More information about the grass-stats mailing list