reading a feature class
The reason is that Joseph Bechara's post included a jpeg image at 16KB, and an HTML component breaking the list requirement to post in plain text only. Multiple copies of his logo then clog up the system. The real answer to is problem is that ogrDrivers() will show that he does not have the driver (certainly not in the CRAN binaries). GDB is a proprietary format that can be used if you builds GDAL/OGR against ESRI binaries, not otherwise (certainly for Personal GDB, and effectively for File GDB). The data must be saved from ArcGIS in an open format. Hope this clarifies, Roger
On Wed, 9 Jul 2014, Frede Aakmann T?gersen wrote:
From: Frede Aakmann T??gersen Sent: 9. juli 2014 15:12 To: 'JOSEPH BECHARA'; 'Michael Sumner' Cc: 'RsigGeo' Subject: FW: [R-sig-Geo] reading a feature class Hi I am getting some information like this from the list: Your mail to 'R-sig-Geo' with the subject RE: [R-sig-Geo] reading a feature class Is being held until the list moderator can review it for approval. The reason it is being held: Message body is too big: 76484 bytes with a limit of 50 KB Either the message will get posted to the list, or you will receive notification of the moderator's decision. If you would like to cancel this posting, please visit the following URL:
https://stat.ethz.ch/mailman/confirm/r-sig-geo/5859051c379864eae2b43f7c1fcf9c94159715a0 Cutting and resending: Hi Well both Michael and I are just guessing because of lack of information. I cannot really help you before I can try to read your file on my system (both windows and linux) before I can be of any use to help you. I suppose that Michael expressed the same. One thing that comes to my mind is that reading vector formats into R requires lots of physical memory. If your computer do not have that it is probably swapping to hard disk (hence working along for ever). Can you check that? If so you might as well kill the process and think on other approaches. Yes, other software may handle this more efficient than R. 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<mailto:frtog at vestas.com> http://www.vestas.com<http://www.vestas.com/> Company reg. name: Vestas Wind Systems A/Sas well This e-mail is subject to our e-mail disclaimer statement. Please refer to www.vestas.com/legal/notice<http://www.vestas.com/legal/notice> If you have received this e-mail in error please contact the sender. From: JOSEPH BECHARA [mailto:jbechara at lri-lb.org] Sent: 9. juli 2014 14:44 To: 'Michael Sumner' Cc: Frede Aakmann T??gersen; 'RsigGeo' Subject: RE: [R-sig-Geo] reading a feature class I tried your suggestion and it is working. I got the below message but it have been running for almost 1 h and didn???t close the ring yet. bb<-readOGR(dsn="All_poly.gdb",layer="All_lays1") OGR data source with driver: OpenFileGDB Source: "All_poly.gdb", layer: "All_lays1" with 9067252 features and 5 fields Feature type: wkbMultiPolygon with 2 dimensions JOSEPH BECHARA Wildfire Management Component Coordinator Ashrafieh 1020 Building, 8th floor, Mathaf, Beirut, Lebanon Tel.: +961 (1) 426630-1-2 Fax: +961 (1) 426631 Mail: jbechara at lri-lb.org<mailto:jbejjani at lri-lb.org> www.lri-lb.org<http://www.lri-lb.org/> From: Michael Sumner [mailto:mdsumner at gmail.com] Sent: Wednesday, July 09, 2014 2:00 PM To: Joseph Bechara Cc: Frede Aakmann T??gersen; RsigGeo Subject: RE: [R-sig-Geo] reading a feature class dd is the layer name/s, so literally I meant to use it to find what "All_lays1" really should be. (I have no experience with this format thankfully). So see for yourself what ogrListLayers reports, then the dsn for readOGR is the same, but it also wants a layer name as the second arg. (I probably need to find one of these "geodatabases" things and try for myself before I can help.) dsn is the data source (the mothership) layer is the one or more layer-thingys in that ship (shapefiles have just one layer, hence the most available example breaks the more general mould) Clear? See ?readOGR for examples On 9 Jul 2014 19:59, "JOSEPH BECHARA" <jbechara at lri-lb.org<mailto:jbechara at lri-lb.org>> wrote: Hi Michael, It worked partially. I still have an error message dd<-ogrListLayers(dsn="C:/sophie/Allpoints/All_poly.gdb") bb<-readOGR(dsn=dd,layer="All_lays1") Error in ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv) : Cannot open file In addition: Warning messages: 1: In if (nchar(dsn) == 0) stop("empty name") : the condition has length > 1 and only the first element will be used 2: In if (nchar(dsn) == 0) stop("empty name") : the condition has length > 1 and only the first element will be used JOSEPH BECHARA Wildfire Management Component Coordinator Ashrafieh 1020 Building, 8th floor, Mathaf, Beirut, Lebanon Tel.: +961 (1) 426630-1-2 Fax: +961 (1) 426631<tel:%2B961%20%281%29%20426631> Mail: jbechara at lri-lb.org<mailto:jbejjani at lri-lb.org> www.lri-lb.org<http://www.lri-lb.org/> [[alternative HTML version deleted]]
Roger Bivand Department of Economics, Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 91 00 e-mail: Roger.Bivand at nhh.no