Error importing ESRI Shapefile

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]]

	[[alternative HTML version deleted]]