On Tue, 26 Jun 2007, Agustin Lobo wrote:
Still some problems:
sessionInfo()
R version 2.5.0 (2007-04-23)
i386-pc-mingw32
locale:
LC_COLLATE=Spanish_Spain.1252;LC_CTYPE=Spanish_Spain.1252;LC_MONETARY=Spanish_Spain.1252;LC_NUMERIC=C;LC_TIME=Spanish_Spain.1252
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "datasets" "methods"
"base"
other attached packages:
maptools foreign spatstat mgcv sp
"0.6-13" "0.8-20" "1.11-7" "1.3-25" "0.9-14"
limit1 <- read.shape("limit_PNSN2000.shp")
limit1owin <- as(limit2,"owin")
Error in inherits(x, "factor") : cannot get a slot ("contains") from an
object of type "NULL"
limit1owin <- as(as(limit2, "SpatialPolygons"), "owin")
Checking 3 polygons...1, 2, [Checking polygon with 2749 edges...]
3.
done.
Checking for cross-intersection between 3 polygons...1, Polygons 1 and 3
cross over
2.
done.
Error in owin(poly = opls) : Polygon data contain overlaps between polygons
plot(limit2)
I can send you limit2.rda, is that ok?
Yes, but as you can see, the problem seems to be that your polygons
overlap. That probably needs resolving first, or in the limit2 object
inside R.
Roger
Agus
Roger Bivand escribió:
On Mon, 25 Jun 2007, Edzer J. Pebesma wrote:
> Roger Bivand wrote:
> > On Mon, 25 Jun 2007, Agustin Lobo wrote:
> >
> > > I've tried (among other things):
> > > > limit1 <- read.shape("limit_PNSN2000.shp")
> > > Shapefile type: Polygon, (5), # of Shapes: 2
> > > > limit1poly <- Map2poly(limit1)
> > > > limit1owin <- as(limit1poly,"owin")
> > > Error in .classEnv(thisClass) : unable to find an environment
> > containing > class "polylist"
> > > > > limit2 <- readOGR("limit_PNSN2000.shp",layer="limit_PNSN2000")
> > > OGR data source with driver: ESRI Shapefile
> > > Source: "limit_PNSN2000.shp", layer: "limit_PNSN2000"
> > > with 2 rows and 1 columns
> >
> > Then just:
> >
> > limit1owin <- as(limit2,"owin")
> >
> > ought to work with current maptools, or possibly coerce:
> >
> > limit1owin <- as(as(limit2, "SpatialPolygons"), "owin")
> >
> the latter form shouldn't be needed thanks to S4 inheritance.
In principle yes, especially with fresh sp and maptools, but I'm not sure
that my code in earlier versions allowed inheritance properly, and we
don't know which package versions are being used. Hence the general
request to state sessionInfo() output in all questions to the list.
Roger
> --
> Edzer
>
--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: [EMAIL PROTECTED]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo