¿Eliminar características seleccionadas usando PyQGIS?

Tengo un shapefile con solo un campo DN . Me gustaría eliminar todas las funciones que no tienen un valor específico del campo. Digamos que el DN puede ser 1,2 o 3.

Hasta ahora pude averiguar cómo seleccionar Funciones y cómo eliminar funciones, pero no logré combinarlas. Digamos que quiero que solo queden las características con DN = 3:

Selección de las características que se encuentran aquí , sin embargo, no sé si hace lo que se supone que debe hacer, ya que no encontré la forma de verificar lo que hay dentro.

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

eliminar funciones encontradas aquí

res = Ilayer.dataProvider().deleteFeatures()

Por lo que yo sé, el deleteFeatures()requiere ID de las características, pero no tengo ni idea de cómo obtenerlos a partir de selectla cual es QgsFeatureIterator.

Puede recorrer el iterador y obtener el .id()para cada característica en él:

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())

O con 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()