layer.GetFeature() failing

[jivan]

Works correctly:

from osgeo import ogr

geopackage = 'rosarito-osm-20170131.gpkg'
layer_name = 'amenities_all_points'

driver = ogr.GetDriverByName("GPKG")
dataSource = driver.Open(geopackage)
print dataSource
layer = dataSource.GetLayer(layer_name)

for feature in layer:
    print feature
layer.ResetReading()

Returns None for every feature:

from osgeo import ogr

geopackage = 'rosarito-osm-20170131.gpkg'
layer_name = 'amenities_all_points'

driver = ogr.GetDriverByName("GPKG")
dataSource = driver.Open(geopackage)
print dataSource
layer = dataSource.GetLayer(layer_name)

for i in range(layer.GetFeatureCount()):
    feature = layer.GetFeature(i)
    print feature

[jivan]

Geopackage file with data from OSM.

[bitner]

Just to add a bit of context, this appears to be the case with geopackages. If you convert the gpkg to a shapefile and rerun the two examples that Jivan pointed out both approaches work as expected.

[Even Rouault]

In 37266:

GPKG: make GetFeature() works on non conformant tables that have no integer primary key field (fixes #6799), and be robust to non standard column types

[Even Rouault]

In 37267:

GPKG: make GetFeature() works on non conformant tables that have no integer primary key field (fixes #6799), and be robust to non standard column types

[Even Rouault]

The tables in the geopackage are non compliant with the specification, as they lack a INTEGER PRIMARY KEY field

​http://www.geopackage.org/spec/#feature_user_tables: """ Requirement 29

A GeoPackage MAY contain tables or updateable views containing vector features. Every such feature table or view in a GeoPackage SHALL have a column with column type INTEGER and PRIMARY KEY AUTOINCREMENT column constraints

"""

I've nevertheless fixed GetFeature() to work in the case such a column is missing