Re: [R-sig-Geo] azymuth error, polygons and Lisbon projection

2007-07-26 Thread Roger Bivand
On Thu, 26 Jul 2007, Marta Rufino wrote: > Hello, > > > I have four doubts: > > 1) > why gzAzymuth is giving me Northern direction in this case, for example? gzAzymuth() in the maptools package is for geographical coordinates, as its help page says, and yours seem to be projected. > > require(s

[R-sig-Geo] azymuth error, polygons and Lisbon projection

2007-07-26 Thread Marta Rufino
Hello, I have four doubts: 1) why gzAzymuth is giving me Northern direction in this case, for example? require(sp); require(rgdal) kk=data.frame(x=c(575552.8,576295.9), y=c(4103900,4103611)) coordinates(kk)=~x+y plot(kk, typ="o", col=4) gzAzimuth(coordinates(kk)[1,], coordinates(kk)[2,]) 2) h