The OGR Simple Features Library is a C++ open source library (and command line tools) providing access to vector file formats.

learn more… | top users | synonyms

7
votes
2answers
155 views
+50

Conflict between ArcPy and GDAL/OGR?

TL;DR I can't get an installation free version of GDAL/OGR to work from within ArcMap 10.0. I get an ImportError: DLL load failed on _mod = imp.load_module('_gdal', fp, pathname, description). ...
0
votes
1answer
15 views

Adding a Vertex To Exisiting Feature(linear) Using OGR Python

I am working on a project that involves finding the intersection of linears in a shapefile, if that vertex(the found intersection) doesn't already exist for each(or one) of the features, then that ...
2
votes
0answers
77 views

Can OGR handle Esri's Geodatabase extensible markup language (XML) files?

The USDA Forest Service distributes GIS data as Esri Geodatabase extensible markup language (XML) files. I downloaded one in the hopes that I would be able to convert it to something less esoteric ...
1
vote
0answers
27 views

Why does this python OGR code for creating a layer and then adding a field not work?

Spot the difference! So this works... fp = "C:\\LiDAR\\MyLayer.shp" layer = "MyLayer" driver = ogr.GetDriverByName("ESRI Shapefile") data_source = driver.CreateDataSource(fp) srs = ...
0
votes
0answers
15 views

Exporting ArcSDE feature class to shapefile using OGR?

I'd like to export data from ArcSDE for SQL Server to SHP. I've tried first to run this python script in QGIS: http://trac.osgeo.org/gdal/browser/trunk/autotest/ogr/ogr_sde.py But finally I've got ...
0
votes
1answer
31 views

Using OGR Intersection Method Python

I'm currently working to find the intersections between linears in a shapefile. I'm a newbie to the GDAL/OGR library. I've done some googling and came up with what seems to be a useful example(see ...
3
votes
1answer
307 views

Latest version of ogr2ogr fails with Layer Creation

I am able to successfully run ogr2ogr to add S57 data to PostgresDB. The version I use is GDAL 1.7.0b2, FWTools 2.4.7 and release date 2010/01/19 I upgraded my version to GDAL 1.9.1 and release date ...
2
votes
1answer
17 views

OGR DeleteDataSource not working properly in Python

I have a small Python script which is supposed to delete a shapefile in case the file is empty. import sys from osgeo import ogr drv = ogr.GetDriverByName("ESRI Shapefile") ds = drv.Open(inData) if ...
0
votes
2answers
86 views

MultiLineString to separate individual lines using Python with GDAL/OGR, Fiona, Shapely

There are lots of MultiLineString features in a polyline shapefile. The problem is to split these lines into separate individual lines. There is a similar question with PostGIS, but I need to code in ...
1
vote
0answers
11 views

QGIS view bounding box on Add vector layer GeoGSON protocol?

Our server can serve GeoJSON with an optional bounding box HTTP parameter. e.g. http://our_gis_server.org.au/layername?bbox=xmin,ymin,xmax,ymax In QGIS I can add this vector layer as a Protocol and ...
1
vote
1answer
22 views

Extracting z values from shapefile using OGR Api

Is there an example of how we can extract the z value from a shp file given the x,y coordinate value (namely utm) using gdal/ogr c++ lib? alternatively, would it be easier if we convert the shp file ...
6
votes
1answer
442 views

Geometry errors when importing data from ESRI file geodatabase to PostGIS

When I run the following command to import by data from an ESRI fgdb to postgis... ogr2ogr -f "PostgreSQL" -a_srs "EPSG:28355" -t_srs "EPSG:3857" PG:"dbname=tp_test user=postgres port=54321 ...
5
votes
1answer
328 views

GDAL.RasterizeLayer not ported to the C# bindings?

I am trying to use the OSGeo.GDAL.Gdal.RasterizeLayer(…) from GDAL 1.9.2 and it requires some SWIGTYPE parameters: OSGeo.GDAL.Gdal.RasterizeLayer( OSGeo.GDAL.Dataset dataset, int bands, ...
0
votes
1answer
41 views

Creating an OVF file for connecting Oracle and QGIS

I'm trying to create an OGR file that can be opened in QGIS that pulls in vector data from an Oracle Database. I've tried a few combinations of things but can't seem to get anything to work. I've few ...
2
votes
2answers
638 views

OSR Python testSR.ImportFromEPSG(4326) returns empty spatial reference

This is empty using Python .ImportFromEPSG(4326): import osr #this fails: testSR = osr.SpatialReference() testSR.ImportFromEPSG(4326) print testSR.ExportToPrettyWkt() #why did the import from EPSG ...
13
votes
7answers
3k views

Free tool to simplify parcel shapefiles?

Is there any tool available which let's me run millions of parcel shapefiles and simplify them?
1
vote
2answers
172 views

I need to extract a specific layer from a Geospatial PDF using GDAL/OGR

I have GDAL/OGR 1.10 installed on a 64bit Windows 7 machine. I have some Geospatial PDFs that I want to crop to the outline of the map and convert to Rasters. However I cannot clip to the neatline. ...
4
votes
1answer
187 views

Loading a multipoint shapefile in R

This issue has only been raised in the community once, two years ago, to the best of my knowledge. I think it's high time to raise the issue again! readOGR, from the rgdal package, cannot handle ...
3
votes
2answers
2k views

ogr2ogr srs options: where is the well known definition determined?

In the ogr2ogr documentation it says that when using -a_srs srs_def: Srs_def can be a full WKT definition (hard to escape properly), or a well known definition (ie. EPSG:4326) or a file with a WKT ...
0
votes
1answer
61 views

Creating records in output shapefile using Python

I have a list (output from a previous loop) containing attributes from an input shapefile. What I am stuck with is inputting those attributes to a new shapefile. I have plotted the point data using ...
1
vote
1answer
30 views

ogr SimplifyPreserveTopology does not keep the Topology

I use ogr2ogr to simplify a layer of polyline, from the document: -simplify tolerance: (starting with GDAL 1.9.0) distance tolerance for simplification. This method will preserve topology, in ...
2
votes
0answers
45 views

Append ArcGIS feature layer data using GDAL/OGR API

ArcGIS servers are configured with a maximum records count that they're allowed to return following a REST query. This maximum is usually 1000 records, which means that GDAL/OGR clients should ...
0
votes
1answer
83 views

Writing shapefile with projection defined crashes fiona

I've just started trying to use Fiona to write out shapefiles. I found the example code below, which seems to work fine from shapely.geometry import mapping, Polygon import fiona # Here's an example ...
3
votes
2answers
137 views

Join a CSV file to shapefile using gdal/ogr?

I have a shapefile with several attributes, for example YEAR, COUNTY, and AREA. I also have a CSV file with more fields that I want in the shapefile, such as POPULATION. Both the shapefile and the CSV ...
0
votes
0answers
26 views

SetSpatialFilter worked unexpected when using ogr by python

I am using ogr python binding to read the shapfile, this is the code: ..... bbox="POLYGON ((20.18164062 30.63867188,20.63867188 30.63867188,20.63867188 30.18164062,20.18164062 30.18164062,20.18164062 ...
0
votes
1answer
37 views

How to filter date values using OGR 1.11.0?

With OGR 1.10.1 this works: ogrinfo test2.shp -sql 'SELECT COUNT(*) FROM test2 WHERE DATE < "1970-01-01"' Note that DATE is a column of type date (not string). (Here's test data if you want to ...
1
vote
0answers
11 views

Asking about OFTInteger data type in OGR/GDAL?

I created the shape file by using the feature field type: OFTInteger. But when I open the file in Arcgis (ArcCatalog) the data type is: Double (not integer). Do you know how to fix this (I want it ...
2
votes
2answers
141 views

OSM linestring shown as point

I have loaded the following file http://openstreetmapdata.com/data/coastlines to postgreSQL. In visual studio with gdal/ogr and geos enabled: SELECT geom FROM waterpoly WHERE geom IS NOT NULL LIMIT ...
2
votes
2answers
33 views

How to replicate ArcGIS coordinate transformation with ogr2ogr?

I need to be able to replicate with ogr2ogr a coordinate transformation that was executed in ArcGIS for Desktop. Both ArcGIS and ogr2ogr are based on EPSG standard definition. There's not only one ...
0
votes
0answers
23 views

Comparing datums with GDAL/OGR

I have a bit of code in Java that takes in a SpatialReference object and checks that it is one of the acceptable datums (WGS84, NAD83, NAD83 HARN, NAD27), but because of the lack of standardization in ...
0
votes
1answer
31 views

TileStache can't find spatial ref for vector layer

I have a query based vector layer in my TileStache configuration, and it's hitting this block of code if layer.GetSpatialRef() is None and driver_name != 'SQLite': raise KnownUnknown('Couldn\'t ...
8
votes
2answers
1k views

Processing SQLite files using ArcGIS Desktop?

I'm using qgis and grass for my work. I've produced a shapefile which is over 100mbs (s a result of some statistical analysis). If i convert it to sqlite format (spatially enabled) the data shrinks ...
0
votes
1answer
52 views

no viewable feature when creating line with ogr in python

the python code below produces the expected line shapefile viewable in QGIS import os from osgeo import ogr drv=ogr.GetDriverByName('ESRI Shapefile') if os.path.exists('test.shp'): ...
2
votes
1answer
456 views

ogr2ogr specify lat/long order on conversion to GeoJSON

When using ogr2ogr to convert a polyon from a .shp file to GeoJSON, I noticed that it writes the coordinates in order lat, long. I noticed this because whenever I tried to view the GeoJSON on a map, ...
0
votes
1answer
67 views

VSI File API support in ogr2ogr

It seems that support for streaming input is wishy-washy in ogr2ogr. For example, why can I do this: curl "http://data.usgin.org/arizona/wfs?service=WFS&version=1.1.0& ...
2
votes
1answer
189 views

Ogr2Ogr: Failed to create GPX file output.gpx

When i try to export an shapefile to gpx using ogr2ogr. ogr2ogr -f GPX output.gpx D:\mypaste\input.shp -sql "SELECT Id AS name FROM layername" It gives me the following error "Failed to create GPX ...
1
vote
1answer
36 views

List the tab/shp files properties within folders/sub folders w.r.t relative folder path [closed]

Currently i am using this command line For /R "U:\01 My Proj" %f in (*.tab) do ogrinfo -al -so "%f" >>U:\New.txt to extract all the tab file related information as a text.Instead of giving the ...
0
votes
1answer
66 views

Get the outer ring inner ring for a given polygon using ogr python binding

I am trying to read the geometry of some shapefiles, now I meet some problem when I try to identify the polygon with holes, since a polygon may(or not) have one hole(or multiple holes). This is the ...
1
vote
1answer
91 views

Union multiple shapefiles with overlapping shapes

I have several shapefiles that overlap and I would like to union their geometries to get rid of those overlaps. I know how to do this with OGR and Python for single geometries. from osgeo import ogr ...
0
votes
2answers
1k views

No ogr module despite having gdal installed?

I am trying to use the mapping component of the wonderful PySAL package, but for some reason I do not have the ogr module. I was under the impression that it came packaged with gdal. Am I missing ...
0
votes
0answers
26 views

GeoJson attributes encoding limitation

I have problem when trying to encode attribute values in latin1 using ogr GeoJson driver. Accent char are not displayed correctly. Is the geoJson encoding is limited to UTF-8? I cant find any info ...
1
vote
1answer
24 views

ogr: what is the advantage of using `geom.Clone()` over `geom`?

What is the advantage of using geom.Clone() over geom? Simple example: geom = link_feat.GetGeometryRef() newfeature.SetGeometry(geom.Clone()) # or newfeature.SetGeometry(geom) ...
0
votes
1answer
37 views

How to get feature for a given bbox from shapefile by ogr

I have some features point,line,polygon with shapefile format, now I want to get features for a given boundbox, is it possible to remove the features out of the bbox? And for a polygon, I think it is ...
1
vote
1answer
78 views

GDAL - How to get gemetry type of a MapInfo TAB file?

I am using the OGR MapInfo File driver to read a MapInfo TAB file. However, the function GetGeomType() returns 0 which means 'wkbUnknown` even if the features in the TAB are only multi-polygons. Is ...
0
votes
0answers
61 views

aggregate counts on a grid and determine grid adjacency

So here is the problem I am working on. I have a bunch of lat, long coordinates for event occurrences. I want to create a grid and then count how many points fall into each grid pixel. Once I have ...
0
votes
1answer
85 views

create polygon adjacency matrix using python

I am trying to create an adjacency matrix from a set of polygons. That is, I have a bunch of polygons and I want to identify which polygons have a common edge or "touch" each other. Once I find this ...
1
vote
1answer
65 views

My class using ogr to calculate distance is returning incorrect results away from the equator

In an answer to this question ustroetz kindly provided some sample code for use of ogr to do lat/lon distance and point in polygon calculations. I've created a Python class based on his example to do ...
3
votes
1answer
110 views

Are the functions in OGRPolygon accessible through the python bindings?

I'm currently trying to use the OGR python bindings to process some shapefiles. I want to obtain the number of interior rings in a polygon feature. Unfortunately I'm unable to use the documented ...
3
votes
1answer
56 views

writeOGR column limit error: “Creating Name field failed”

I have a SpatialPointsDataFrame with 1057 elements and 4760 columns and want to export it to a shapefile. When I try to do this with writeOGR() I get an error: > ...
0
votes
0answers
15 views

ogr: confusing results from `SetSpatialFilter()`

My layer is derived from a data source like this: l_131 = bounds.ExecuteSQL("select * from 'nam_l' where 'ND_5'= 131 ") feat = l_131.GetNextFeature() feat.GetFeatureCount() >>> 200 It is a ...