[R-sig-Geo] problem with disc ( ) in spatstat

Barry Rowlingson b.rowlingson at lancaster.ac.uk
Mon Dec 27 16:40:25 CET 2010

On Mon, Dec 27, 2010 at 3:07 PM, gianni lavaredo
<gianni.lavaredo at gmail.com> wrote:
> Dear Researchers,
> is There some problem with DISC() in spatstat?
> disc =  disc(radius=10, centre=c(774000.8,7249592.3), npoly=126)
> I am try also to run the example but I have always a problem

 Works for me - again you've not posted any error messages! Note its a
bad idea to call something by the same name as a function. Anyway,
here's how it goes for me - this creates and plots a nice disc called

 > library(spatstat)
Loading required package: mgcv
This is mgcv 1.6-1. For overview type `help("mgcv-package")'.
Loading required package: deldir
deldir 0.0-12

     Please note: The process for determining duplicated points
     has changed from that used in version 0.0-9 (and previously).

spatstat 1.20-3
Type ‘help(spatstat)’ for an overview of spatstat
     ‘latest.news()’ for news on latest version
     ‘licence.polygons()’ for licence information on polygon calculations
 > d =  disc(radius=10, centre=c(774000.8,7249592.3), npoly=126)
 > d
window: polygonal boundary
enclosing rectangle: [773990.8, 774010.8] x [7249582, 7249602] units
 > plot(d)

 - a disc appears!

More information about the R-sig-Geo mailing list