2011/9/26 Roger Bivand <Roger.Bivand@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]]

