readShapePoints issue
On Mon, 26 Sep 2011, Raffaele Morelli wrote:
Hi, I am working with R and using pgsql2shp to export shape files from pg. The shapes are created by pgsql2shp but when I try to import them in a R enviromnent I get this error Error in res[i, ] <- Map$Shapes[[i]]$verts : replacement has length zero So, I can load the shape files in Qgis but not in R. Moreover, if I load one of those files in Qgis and then export it again... this file is correctly loaded by R. I really can't guess what's wrong. Any idea?
Use readOGR() in rgdal, not readShapeSpatial() in maptools. You are also generating objects with no coordinates, so the real issue is upstream, in PostGIS or the intermediate utility. Why not build your own PostGIS driver for OGR, and read directly with readOGR(): http://wiki.intamap.org/index.php/PostGIS Hope this clarifies, Roger
Regards Raffaele -- L'unica speranza di catarsi, ammesso che ne esista una, resta affidata all'istinto di ribellione, alla rivolta non isterilita in progetti, alla protesta violenta e viscerale.
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Roger Bivand Department of Economics, NHH Norwegian School of Economics, 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