Questions tagged [ogr]
The OGR Simple Features Library is a C++ open source library (and command line tools) providing access to vector file formats.
2
votes
2answers
33 views
OGR: Load layer with sql statement and write to new data set
I am having difficulties after loading a shapefile with a sql statement with ogr.
The layer and its features seem to be correctly loaded. Still I am not able to write them into a new shapefile or pass ...
1
vote
0answers
22 views
Select features of a shape which are not in another shape
I have a shape with some features (let's say 100) and a second shape with some of this features (50) buffered (5 meters inside buffer). I want to select and write in a new shape every feature of the ...
0
votes
0answers
26 views
Export OGR feature style data from QGIS for use in GeoPDF
I have a fully formated map in QGIS from which I wish to make a GeoPDF.
I followed this and this tutorial and generated a simple version of my map without much problems.
I am however looking for an ...
2
votes
1answer
27 views
GDAL SQLite geometry error
I am having the following issue writing points to a SQLite file using ogr with SPATIALITE enabed:
ERROR 1: Cannot insert feature with geometry of type POINTZ in column GEOMETRY. Type POINT expected
...
0
votes
0answers
13 views
gdal FlattenTo2D makes geom null
I have a df that has geoms as ogr geom in 3D. I was trying to convert to 2D but it empties out the geom.
ipdb> df['default_geom'][0]
<osgeo.ogr.Geometry; proxy of <Swig Object of ...
0
votes
0answers
21 views
What is the geom_type at CreateLayer() in ogr for xyz-GeoJSON?
My GeoJSON-File is a Multipolygon which creates several triangles with a X, Y and Z coordinate. Z determines the height of the specified point for X/Y.
If i try to render my GeoJSON into a shapefile ...
1
vote
1answer
32 views
OGR SetAttributeFilter with multiple fields?
I am trying to open a layer using OGR in Python and filter the features using SetAttributeFilter() where two fields are greater than zero.
from osgeo import ogr
driver = ogr.GetDriverByName('ESRI ...
0
votes
1answer
45 views
How to extract FileGDB's domains without ArcGIS?
I'm trying to get an ESRI File Geodatabase domains by using the GDAL\OGR library. By doing so, I wish to edit this datasource using QGIS' Attributes Form with designated value map widgets that point ...
0
votes
0answers
31 views
GDAL OGR Polygon.Contains(Point) returns false but point is contained in QGIS map
I'm sure I must be missing something here but for the life of me I can't see what it might be. I have been struggling for the last hour with an issue with the Contains() method in OGR.
I am iterating ...
0
votes
1answer
43 views
Cascading union OGR Multipolygons?
I have 2 OGR Multipolygons that intersect, hence I want to cascade Union them. I tried this but it returns me empty and doesn't union it. Seems this works only on Polygons and not Mul;tipolygons. What'...
1
vote
1answer
95 views
Drawing graticule with GDAL/OGR?
I'm trying to generate a graticule with user-defined range of latitude, longitude, and interval; however, my code only displays one line of latitude and one line of longitude when loading it in QGIS. ...
2
votes
0answers
54 views
Writing selected features to shapefile using OGR and Python?
I have two shapefiles that intersect: one containing a line, and the other containing polygons. I'm trying to use OGR and Python to create a buffer around the line and select polygons within it, then ...
0
votes
1answer
48 views
Clip shapefile based on another shapefile [closed]
I have two shapefiles with me, one big and one small. I want to clip part of big shapefile where the small shapefile is located. I tried QGIS vector intersect by geometry, but that is returning me ...
0
votes
0answers
17 views
QGIS with MSSQLSpatial : no geometry. OgrInfo Corrupt data
I am having this issue while I am trying to connect to my MSSQLSpatial database.
So, I've created this features with just one record with geometry. I can see this feature on Microsoft Management ...
0
votes
1answer
38 views
How to convert an ASCII file (with lat,lon coordinates) to a shapefile using OGR?
I have an ASCII file with the following columns : ID, longitude, latitude, value
Is it possible to use OGR (ogr2ogr or ogrinfo) in command line in order to open the ASCII file and convert it to a ...
1
vote
1answer
30 views
Find layer name of sp object
After making some changes to an sp object I want to save it as a shapefile using writeOGR()so that I don't have to repeat the cleaning every time.
I haven't found a way to identify the layer name of ...
0
votes
1answer
66 views
“OGR: not enough data error” when using st_transform in R
I'm trying to use st_transform (from the R package sf) on a sf object and am getting the following error when projecting from WGS 1984 to Albers Equal Area:
> sf_obj_prj <- sf_obj %>%
...
0
votes
3answers
57 views
ogr2ogr CSV conversion with two points into column names
I have a CSV with ":" in all column names. I am trying to convert it with ogr2ogr to another format. I don't need all fields
example for formating:
#0:imgID,1:camSN,2:lat,3:lon,...,57:fields
...
0
votes
0answers
11 views
QgsVectorLayer, configuring ogr uri for a shapefile (to set field datatypes)
I'd like to load a .shp and set the data type of one of the fields. I thought something like this would work:
uri = 'C:/myfile.shp?field=myfield:integer'
vlay = QgsVectorLayer(uri ,'some_name','ogr')...
2
votes
1answer
48 views
When is it worthwhile to reproject in QGIS? [closed]
My QGIS project has various raster and (complex) vector layers from different sources, in different CRSes. Using reprojection on the fly, I don't have to reproject them. But should I, especially for ...
0
votes
2answers
47 views
How to efficiently extract specific metadata for many shapefiles?
I have a directory with 50 subdirectories containing a total of 300 shapefiles which I would like to query for their feature count and output to a text file.
Thus far I am able to get the feature ...
2
votes
0answers
37 views
QGIS parsing layer source
I periodically need to refresh my QGIS map based on changed underlying data sources (currently locally-saved KML, GPKG, but trying to be general) using PyQGIS. Following Make a Qgis layer update from ...
0
votes
0answers
65 views
Alter and update a shapefile with a field from an ASCII file using ogr
I have a vector file "map.shp" with a field ID.
I have an attribute table (data.txt) with the following fields: ID, value1, value2. The field ID is the same in both files.
I want to update map.shp by ...
2
votes
1answer
71 views
Populating new field in shapefile with values from list?
I'm looking to populate a newly created field in a shapefile that from a list using basic python and osgeo tools. The values are integers and in a python list that is the length of the number of ...
1
vote
1answer
45 views
Getting Field by name from feature using GDAL?
I have recently started using GDAL in C#, I want to get a value from a field of a feature using the name of that field. In the docs here it says that the operator [char*] is supported, however it's ...
0
votes
1answer
57 views
How to specify GPKG Driver in GDAL/OGR under C++?
Working under CentOS (Linux), we have run up against the 4GB shapefile (dbf) limitation in our Monte Carlo model. We've decided to switch to GeoPackage to output our point attributes. I cannot find a ...
1
vote
1answer
30 views
Add FID to attribute table using ogrinfo
I want to add the fid to the attribute table using ogrinfo in OSGeo4W Shell.
Here is my code:
ogrinfo layer.shp -sql "ALTER TABLE layer ADD COLUMN ID integer"
ogrinfo layer.shp -dialect SQLite -...
0
votes
1answer
19 views
How could a GeoJSON look that OGR recognizes more than one layer?
I load a GeoJSON with OGR and apply the method GetLayerCount().
How must the GeoJSON look like that GetLayerCount() > 1?
driver = ogr.GetDriverByName("GeoJSON")
dataset = driver.Open(my.geojson, 0)
# ...
2
votes
1answer
94 views
Can't delete features from TAB layer table in QGIS if no geometry
I have some layers containing DAT,MAP and TAB files, when I loaded these layers into QGIS some layers have check boxes before layer name and some layers won't like showing in below image.
In the ...
1
vote
1answer
58 views
Convert an OGR File to a Raster
i tried out this recipe from the GDAL/OGR cookbook:
from osgeo import gdal, ogr
# Define pixel_size and NoData value of new raster
pixel_size = 25
NoData_value = -9999
# Filename of input OGR file
...
0
votes
1answer
29 views
How to free the memory of return value of OGRGeometry::Buffer()?
I have in C++ (gdal 2.3.2)
OGRPolygon poly;
OGRGeometry* geom = poly.Buffer(1e-12, 1);
now how do I free the memory allocated for geom after done using it?
0
votes
0answers
116 views
Merge two shapefiles into one new shapefile with OGR/GDAL library in C#?
I am currently coding a C#/WPF application that lets you choose two shapefiles and create a new merged one, with the content of both shapefiles.
My problem is that I am quite new (1 month with gis) ...
3
votes
1answer
51 views
How to export a unique ID field of a geopackage
When I am converting my geopackage file in shapefile, the unique ID field is not exported. I could duplicate it before exportation, but I wanted to know if it could be done in one step, in order not ...
0
votes
1answer
179 views
Converting self-intersecting polygons into distinct polygons
How to convert self-intersecting polygons into distinctive, well behaving polygons?
My polygon inputs consist of self-intersecting polygons, which fails the st_isvalid test, ie: the below statement ...
0
votes
0answers
129 views
OGR Python: Copying features from one datasource to another works with shapefile but not with geopackage
I have a shp/dbf EPSG4674 input datasource, and I want to copy features from one layer of this input datasource, create a new output datasource and write these features to that output datasource.
As ...
0
votes
1answer
209 views
Getting OGR to work on a Mac
I downloaded GDAL on my Mac (running Mojave 10.14.1) through Homebrew and I'm having trouble getting OGR to work with Postgresql. When I try the following to create a geojson from gpx files I have, I ...
0
votes
1answer
92 views
How can I write a zipped shapefile with ogr2ogr and /vsizip/?
After reading a bit of GDAL Virtual File Systems docs, I am trying to write a zipped Shapefile using this gdal command involving the use of /vsizip/:
ogr2ogr -f "ESRI Shapefile" /vsizip/test.zip\test....
0
votes
2answers
182 views
QGIS SERVER - convert WFS GetFeature GML as Shapefile with ogr2ogr
I'd like to use ogr2ogr with the WFS driver to convert the GML response of my QGIS Server 3.4.2 GetFeature layer into Shapefile. I have to do this from a WebGIS JavaScript application I am developing, ...
1
vote
1answer
43 views
Keep only one geometry type when clipping polygons to sqlite with ogr2ogr
I have a regular grid of polygons and want to clip them using ogr2ogr. While using sqlite files as destination I stumbled over some wierd behavior:
When clipping polygons exactly at their borders, a ...
2
votes
1answer
64 views
Can't open specific layers with Qgis nor OGR
I'm working with mutliple sets of layers in order to create a new one.
As I need to load them directly from the internet, I use both /vsizip/ and /vsicurl/
But I'm stuck at loading one of them.
Here ...
0
votes
1answer
110 views
ogr2ogr polygon to point with point inside the polygon
Can ogr2ogr be used to convert a polygon to a point but the point must remain inside the polygon? It does not matter where the point is as long as it is inside the object polygon.
0
votes
0answers
195 views
Combine GeoJSON files into bigger shapes with external JSON to merge on
I have a shapefile of New York ZCTAs here
I also have these neighborhood definitions that map zip codes to New York City Boroughs.
I'd like to take my "all of New York zip code" data and:
Clump up ...
1
vote
0answers
137 views
OGR SetAttributeFilter does not let datasource to update
I am trying to update an attribute value for a certain feature in a shapefile. To filter out that feature, I use SetAttributeFilter. This returns the feature I want, and SetFeature also works as ...
1
vote
0answers
31 views
Which polygons intersect rectangle in PHP/Linux GDAL?
I'm developing a website where I want to test which polygons from a shapefile intersect a given rectangle. Using PHP, I thought about using exec() function, calling a GDAL program, like ogrinfo. ...
1
vote
1answer
93 views
QGIS cannot save to SpatiaLite after vector difference?
In QGIS I have a spatialite DB full of roads. I'd like to do a Difference operation (Vector > Geometry Tools > Difference) with another layer in order to keep roads only available in a given ...
0
votes
1answer
189 views
Get a bounding box of a geometry that crosses the antimeridian using OGR
I'm trying to get a bounding box of the geometry that crosses the antimeridian using OGR:
from osgeo import ogr
geom = ogr.CreateGeometryFromJson('{"type":"MultiPolygon","coordinates":[[[[-180.0,-40....
1
vote
1answer
98 views
C++ OGR GDAL Set spatial reference of new layer
When creating a new layer using OGR / GDAL in C++, how do i set the spatial reference?
OGRLayer *poLayer2;
poLayer2 = poDS2->CreateLayer( "layerName", NULL, wkbPoint, NULL );
if( poLayer2 == NULL ...
4
votes
1answer
93 views
Get a list of all supported EPSG code from OGR
I need get a a list of all supported EPSG code from OGR/GDAL. I need this becaouse i'm developing a C++ application map, and I need this information to put it in a user form.
I read the OGR ...
1
vote
1answer
86 views
Get path of feature classes in feature datasets without arcpy
I'm looking for a way to get the full path of a feature class within feature datasets (either in personal or file geodatabases) without using arcpy.
My first attempt was to use ogr which works fine ...
0
votes
1answer
94 views
GDAL/BUILD can't get result files Geojson to mbtiles
i try to 2.55gb of GeoJson (maybe 200 over layers in one file) to vector Mbtiles, but fail to create .mbtiles.
so i try convert each layers(.shp) to geojson(over 200 files),
then create each ....
