2011/9/26 Roger Bivand <roger.biv...@nhh.no>
> 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<http://wiki.intamap.org/index.php/PostGIS>
>
> Hope this clarifies,
>
yes, readOGR did the job.
Thank you
-r
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
