[R-sig-Geo] Error importing ESRI Shapefile

Alec Friedhoff hrbuilder at hotmail.com
Tue Mar 10 15:46:40 CET 2015


Hi Frede,

My apologies for the incomplete question. I did post on stackoverflow and stackexchange, but have deleted the latter post (I've never visited onlinedevs.com). For what it's worth, I do update questions I've asked on sites like stackoverflow if I find answers elsewhere, but I do understand your general annoyance.

Here's some more information:

Output of traceback():
4: .Call("make_Polygonlist", iG, gComments[[i]], PACKAGE = "rgdal")
3: stopifnot(is.list(srl))
2: Polygons(.Call("make_Polygonlist", iG, gComments[[i]], PACKAGE = "rgdal"), 
       ID = as.character(fids[i]))
1: readOGR("/geo/layers", "tracts")

Message I receive upon loading of rgdal library:
Loading required package: sp
rgdal: version: 0.9-1, (SVN revision 518)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.10.1, released 2013/08/26
Path to GDAL shared files: /usr/share/gdal/1.10
Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
Path to PROJ.4 shared files: (autodetected)

ogrDrivers() reveals I do have Shapefile support:
    name write
1      AeronavFAA FALSE
2          ARCGEN FALSE
3          AVCBin FALSE
4          AVCE00 FALSE
5             BNA  TRUE
6         CouchDB  TRUE
7             CSV  TRUE
8             DGN  TRUE
9            DODS FALSE
10            DXF  TRUE
11         EDIGEO FALSE
12  ElasticSearch  TRUE
13 ESRI Shapefile  TRUE
14     Geoconcept  TRUE
...

Output of sessionInfo():
R version 3.1.2 (2014-10-31)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8     LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                  LC_ADDRESS=C               LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] rgdal_0.9-1 sp_1.0-17  

loaded via a namespace (and not attached):
[1] grid_3.1.2      lattice_0.20-30 tools_3.1.2  


Output of ogrListLayers():
[1] "tracts"
attr(,"driver")
[1] "ESRI Shapefile"
attr(,"nlayers")
[1] 1


Output of ogrInfo():
Source: "/geo/layers", layer: "tracts"
Driver: ESRI Shapefile number of rows 73682 
Feature type: wkbPolygon with 2 dimensions
Extent: (-179.1473 17.88124) - (179.7785 71.39048)
CRS: +proj=longlat +datum=WGS84 +no_defs  
LDID: 0 
Number of fields: 34 
         name type length typeName
1        FIPS    4     11   String
2        SQMI    2     16     Real
...

Any help is greatly appreciated.

Best,
Alec

> From: frtog at vestas.com
> To: hrbuilder at hotmail.com; r-sig-geo at r-project.org
> Date: Tue, 10 Mar 2015 09:49:55 +0100
> Subject: RE: [R-sig-Geo] Error importing ESRI Shapefile
> 
> Hi Alec
> 
> It is very impolite to cross post you questions. Besides the r-sig-geo you have probably posted on
> 
> 1) http://onlinedevs.com/problem-reading-shapefile-using-r/ 
> 2) http://gis.stackexchange.com/questions/138217/problem-reading-shapefile-using-r 
> 3) http://stackoverflow.com/questions/28957014/error-reading-in-esri-shapefile-using-r 
> 
> Are those sites sharing information? Do people on those sites see my specific answer to your question on r-sig-geo?
> 
> Also please do not post in html but in plain text. Set your email client to post in text when posting to r-sig-geo.
> 
> Nobody can really help you. You give no information from which we can see why you get this error message.
> 
> To make you help us to help you please do something along this route. If you get an error please issue an
> 
>> traceback()
> 
> at command prompt in the R shell. Now do this at R prompt:
> 
>> library(rgdal)
> 
> and post the output from that call probably something like this:
> 
> " Loading required package: sp
> rgdal: version: 0.9-1, (SVN revision 518)
> Geospatial Data Abstraction Library extensions to R successfully loaded
> Loaded GDAL runtime: GDAL 1.11.1, released 2014/09/24
> Path to GDAL shared files: c:/Programmer/R/R-3.1.2/library/rgdal/gdal
> GDAL does not use iconv for recoding strings.
> Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
> Path to PROJ.4 shared files: c:/Programmer/R/R-3.1.2/library/rgdal/proj"
> 
>> sessionInfo()
> 
> Post the output!
> 
>> ogrDrivers()
> 
>> dsn <- "/path/to/layer"
> 
>> file.choose(dsn)
> 
> Can you see any shapefiles?
> 
>> ogrListLayers(dsn)
> 
> Can you see any layer names?
> 
>> ogrInfo(dsn=dsn, layer="layer")
> 
> Any information from this call?
> 
> 
> Yours sincerely / Med venlig hilsen
> 
> 
> Frede Aakmann Tøgersen
> Specialist, M.Sc., Ph.D.
> Plant Performance & Modeling
> 
> Technology & Service Solutions
> T +45 9730 5135
> M +45 2547 6050
> frtog at vestas.com
> http://www.vestas.com
> 
> Company reg. name: Vestas Wind Systems A/S
> This e-mail is subject to our e-mail disclaimer statement.
> Please refer to www.vestas.com/legal/notice
> If you have received this e-mail in error please contact the sender. 
> 
>> -----Original Message-----
>> From: R-sig-Geo [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of Alec
>> Friedhoff
>> Sent: 10. marts 2015 07:21
>> To: r-sig-geo at r-project.org
>> Subject: [R-sig-Geo] Error importing ESRI Shapefile
>> 
>> Hi all,
>> I'm trying to read in a Shapefile using the rgdal library, and am unable to do
>> so.When using the following code, geo <-
>> readOGR("/path/to/layer","layer"),I'm met with the error "Error in
>> stopifnot(is.list(srl)) : infinite label point"
>> This one has me stumped. My searches haven't yielded any useful
>> information about this error. I suspect there are some bad geometries in the
>> underlying file. I have access to ArcMap and tried running the
>> RepairGeometry_management method of arcpy, but this has no effect.
>> How can I go about understanding this error and ultimately fixing it?
>> Many thanks!
>> -HR
>> [[alternative HTML version deleted]]
>> 
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
 		 	   		  
	[[alternative HTML version deleted]]



More information about the R-sig-Geo mailing list