[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