An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20071213/b0a04b2c/attachment.pl>
maybe a new bug in the elide() funcftion for maptools’ version 0.7-2
2 messages · zhijie zhang, Roger Bivand
On Thu, 13 Dec 2007, zhijie zhang wrote:
Dear Roger, There maybe a new bug in the new elide() funcftion for maptools' version 0.7-2, which was correct in the former version. See below.
Again correct - please simply replace the elide()ed object data slot with the original. Roger
#Example in the help page
data(meuse)
coordinates(meuse) <- c("x", "y")
proj4string(meuse) <- CRS("+init=epsg:28992")
data(meuse.riv)
river_polygon <- Polygons(list(Polygon(meuse.riv)), ID="meuse")
rivers <- SpatialPolygons(list(river_polygon))
proj4string(rivers) <- CRS("+init=epsg:28992")
##maptools' version 0.7-2, the transformed dataset dropped two variables in
the rightmost: landuse and dist.m.
meuse1 <- elide(meuse, bb=bbox(rivers), scale=TRUE)
meuse[1,]
coordinates cadmium copper lead zinc elev dist om ffreq soil lime landuse dist.m 1 (181072, 333611) 11.7 85 299 1022 7.909 0.00135803 13.6 1 1 1 Ah 50
meuse1[1,]
coordinates x y cadmium copper lead zinc elev dist om ffreq soil lime 1 (0.23093, 0.660129) 181072 333611 11.7 85 299 1022 7.909 0.00135803 13.6 1 1 1 ##maptools' version 0.6-19 + original elide() function: i can't found the former version for maptools package with elide() function, but i use the elide() function that you sent me for the 1st time to show its results. You can try the former version. meuse1 <- elide(meuse, bb=bbox(rivers), scale=TRUE)
meuse[1,]
coordinates cadmium copper lead zinc elev dist om ffreq soil lime landuse dist.m 1 (181072, 333611) 11.7 85 299 1022 7.909 0.00135803 13.6 1 1 1 Ah 50
meuse1[1,] #it didn't drop any variables
coordinates cadmium copper lead zinc elev dist om ffreq soil lime landuse dist.m 1 (0.23093, 0.660129) 11.7 85 299 1022 7.909 0.00135803 13.6 1 1 1 Ah 50 The reason maybe the latest version remains the original x/y coordinates that resulted in overlaying the rightmost two variables, while the former version dropped the original x/y coordinates . So i think two additional colunms maybe needed there. Thanks.
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: Roger.Bivand at nhh.no