Excluindo recursos selecionados usando o PyQGIS?

Eu tenho um shapefile com apenas um DN de campo . Gostaria de excluir todos os recursos que não possuem valor específico do campo. Digamos que o DN possa ser 1,2 ou 3.

Até agora, pude descobrir como selecionar Recursos e como excluir recursos, mas não consegui combiná-los. Digamos que eu queira que apenas os recursos com DN = 3 permaneçam:

Seleção dos recursos encontrados aqui , no entanto, não sei se ele faz o que deveria fazer, pois não encontrei a maneira de verificar o que está dentro.

select = layer.getFeatures(QgsFeatureRequest().setFilterExpression(u'"DN"!=3'))

excluir recursos encontrados aqui

res = Ilayer.dataProvider().deleteFeatures()

Tanto quanto eu sei, o deleteFeatures()requer IDs dos recursos, mas não tenho idéia de como obtê-los a partir do selectque é QgsFeatureIterator.

Você pode percorrer o iterador e obter .id()todos os recursos nele:

with edit(layer):
    # build a request to filter the features based on an attribute
    request = QgsFeatureRequest().setFilterExpression('"DN" != 3')

    # we don't need attributes or geometry, skip them to minimize overhead.
    # these lines are not strictly required but improve performance
    request.setSubsetOfAttributes([])
    request.setFlags(QgsFeatureRequest.NoGeometry)

    # loop over the features and delete
    for f in layer.getFeatures(request):
        layer.deleteFeature(f.id())

Ou com QGIS <2.12

request = QgsFeatureRequest().setFilterExpression('"DN" != 3')
request.setSubsetOfAttributes([])
request.setFlags(QgsFeatureRequest.NoGeometry)
ids = [f.id() for f in layer.getFeatures(request)]
layer.startEditing()
for fid in ids:
    layer.deleteFeature(fid)
layer.commitChanges()